Construction of BSseq Object from Preprocessed Methylation Data

This R script aggregates multiple preprocessed methylation count files into a single BSseq object using the DSS and bsseq 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 formattingv2.r), structured in DSS-compatible format with columns: chr, pos, N, and X.Process Overview:Loads 40 paired .txt files representing different samples and conditions.Constructs a BSseq object using makeBSseqData(), associating each dataset with a sample ID.Saves the full R environment to Loaded_Dataframe.RData for downstream analysis (e.g., in experimentaldesign.r).Inputs:40 DSS-formatted methylation files (W1__DSS_Format.txt, ..., W40__DSS_Format.txt), each with:chr: chromosomepos: positionN: total read coverageX: number of methylated readsOutputs:A BSseq object containing all samples.Saved R workspace as Loaded_Dataframe.RData (contains the BSobj).Software Requirements:R (≥ 4.0)Packages: DSS, bsseqUsage Notes:Ensure all 40 DSS format files are present in the working directory before running.The sample labels used in makeBSseqData() (e.g., "D6C1", "D6D1", etc.) should align with your experimental design matrix for consistency in downstream modeling.The saved .RData file is expected by scripts performing differential methylation analysis, such as experimentaldesign.r.