isglobal-brge / meal Goto Github PK

> GRset_EPICs
class: GenomicRatioSet 
dim: 705657 25 
metadata(0):
assays(1): Beta
rownames(705657): cg14817997 cg26928153 ... cg07660283 cg09226288
rowData names(0):
colnames(25): AT102_VAT AT180_VAT ... CAT84_CF CAT86_CF
colData names(1): X1
Annotation
  array: IlluminaHumanMethylationEPIC
  annotation: ilm10b4.hg19
Preprocessing
  Method: Matrix converted with makeGenomicRatioSetFromMatrix
  minfi version: NA
  Manifest version: NA

> expr_set
ExpressionSet (storageMode: lockedEnvironment)
assayData: 16255 features, 28 samples 
  element names: exprs 
protocolData: none
phenoData
  sampleNames: AT102_VAT AT180_VAT ... CAT86_CF (28 total)
  varLabels: Individual Gender ... agegroup (14 total)
  varMetadata: labelDescription
featureData
  featureNames: ENSG00000000003 ENSG00000000005 ... ENSG00000278845
    (16255 total)
  fvarLabels: chromosome start end strand
  fvarMetadata: labelDescription
experimentData: use 'experimentData(object)'
Annotation:

multi <- createMultiDataSet()
multi <- add_genexp(multi, expr_set)
multi <- add_methy(multi, GRset_EPICs)

> commonSamples(multi)
Object of class 'MultiDataSet'
 . assayData: 2 elements
    . expression: 16255 features, 25 samples 
    . methylation: 705657 features, 25 samples 
 . featureData:
    . expression: 16255 rows, 4 cols (chromosome, ..., end)
    . methylation: 705657 rows, 5 cols (seqnames, ..., width)
 . rowRanges:
    . expression: YES
    . methylation: YES
 . phenoData:
    . expression: 25 samples, 15 cols (Individual, ..., agegroup)
    . methylation: 25 samples, 2 cols (X1, id)

> methExprs <- correlationMethExprs(multi)
Warning messages:
1: In .Seqinfo.mergexy(x, y) :
  The 2 combined objects have no sequence levels in common. (Use
  suppressWarnings() to suppress this warning.)
2: In correlationMethExprs(multi.filt) :
  There are no expression probes in the range of the cpgs. An empty data.frame will be returned.