Views
No views yet
.RData File Structure.RData file contains exactly three objects:CpG_matCpG_mat is a data.frame with one row per CpG model in the file.CpG: CpG identifier.chr: chromosome for the CpG.p.0, p.1: CpG genomic coordinates stored in the file.Corr: correlation/performance value for the CpG model.R2: R-squared value for the CpG model.pval: p-value for the CpG model.SNPsSNPs is a named list. The list names are CpG IDs and match CpG_mat$CpG.SNPs[[<CpG>]] entry is a data.frame containing the SNP model for that CpG, with columns:CHR: chromosome of the SNP.SNP: SNP identifier, usually an rsID.P0: genomic position of the SNP in GRCh38.A0: effect alleleA1: other allele.weight: model weights.weight column is stored as a one-column numeric matrix inside the data frame.elastic_net.MatrixLDMatrixLD is a named list. The list names are CpG IDs and match CpG_mat$CpG.MatrixLD[[<CpG>]] entry is a square numeric matrix giving pairwise LD/correlation values for the SNPs used in the corresponding CpG model.nrow(CpG_mat) == length(SNPs) == length(MatrixLD).CpG_mat$CpG, names(SNPs), and names(MatrixLD) align.cg, SNPs[[cg]] contains the SNP model and MatrixLD[[cg]] contains the LD matrix for the same CpG model.nrow(SNPs[[cg]]) matches nrow(MatrixLD[[cg]]), and the SNP row names match the LD matrix row/column names.1load("chr1_87.RData")
2
3str(CpG_mat)
4names(SNPs)[1:3]
5
6cg <- CpG_mat$CpG[1]
7str(SNPs[[cg]])
8dim(MatrixLD[[cg]])