Construction of BSseq Object from Preprocessed Methylation Data

This R script aggregates multiple preprocessed methylation count files into a single <code>BSseq</code> object using the <b>DSS</b> and <b>bsseq</b> packages. The generated object is suitable for differential methylation analysis in bisulfite sequencing (BS-seq) experiments.The script expects input files produced from a binomial filtering step (e.g., from a script like <code>formattingv2.r</code>), structured in DSS-compatible format with columns: <code>chr</code>, <code>pos</code>, <code>N</code>, and <code>X</code>.<b>Process Overview:</b>Loads 40 paired <code>.txt</code> files representing different samples and conditions.Constructs a <code>BSseq</code> object using <code>makeBSseqData()</code>, associating each dataset with a sample ID.Saves the full R environment to <code>Loaded_Dataframe.RData</code> for downstream analysis (e.g., in <code>experimentaldesign.r</code>).<b>Inputs:</b>40 DSS-formatted methylation files (<code>W1__DSS_Format.txt</code>, ..., <code>W40__DSS_Format.txt</code>), each with:<code>chr</code>: chromosome<code>pos</code>: position<code>N</code>: total read coverage<code>X</code>: number of methylated reads<b>Outputs:</b>A <code>BSseq</code> object containing all samples.Saved R workspace as <code>Loaded_Dataframe.RData</code> (contains the <code>BSobj</code>).<b>Software Requirements:</b>R (≥ 4.0)Packages: <code>DSS</code>, <code>bsseq</code><b>Usage Notes:</b>Ensure all 40 DSS format files are present in the working directory before running.The sample labels used in <code>makeBSseqData()</code> (e.g., <code>"D6C1"</code>, <code>"D6D1"</code>, etc.) should align with your experimental design matrix for consistency in downstream modeling.The saved <code>.RData</code> file is expected by scripts performing differential methylation analysis, such as <code>experimentaldesign.r</code>.