[BioC] methyAnalysis - changes to GenoSet
Tim Triche, Jr.
tim.triche at gmail.com
Tue Mar 12 16:32:19 CET 2013
Right, my original intent was simply to coerce the SummarizedExperiment-based containers from methylumi and minfi into a genoset for Pan's package.
That reminds me, I need to update the documentation for 450kProbe to note that it is essentially obsolete, since both methylumi and minfi can emit SummarizedExperiments. The only mop-up now is for me to send in a few minfi patches for processing. The innards of minfi are tidier than... Others.
--t
On Mar 12, 2013, at 6:57 AM, Kasper Daniel Hansen <kasperdanielhansen at gmail.com> wrote:
> Lavina,
>
> I don't really have time to really look into all of this, but note the
> function mapToGenome which associates a MethylSet with genomic
> coordinates. So you can do
>
> gmSet = mapToGenome(mSet)
>
> and now you have
>
> granges(gmSet) # locations
> getMeth(gmSet) # methylation channel
>
> etc. Note that the transformation is non-invertible because a few
> CpGs have no associated genomic locations and they are dropped (well,
> you can use drop=FALSE).
>
> Kasper
>
> On Mon, Mar 11, 2013 at 5:13 PM, Lavinia Gordon
> <lavinia.gordon at mcri.edu.au> wrote:
>> Argh!, thanks Martin, sorry about that.
>>
>> So in that case I think I've worked out a (rather ugly) way to force minfi data into a MethyGenoSet object for use within the package methyAnalysis:
>>
>> Library(minfi)
>> # mset is minfi object, e.g. used in dmpFinder code
>> myexprs <- getM(mSet)
>> mymeth <- getMeth(mSet)
>> myunmeth <- getUnmeth(mSet)
>>
>> # using Tim Triches helpful code, thanks!
>> library(IlluminaHumanMethylation450kprobe)
>> data(IlluminaHumanMethylation450kprobe)
>> library(GenomicRanges)
>> chs = levels(IlluminaHumanMethylation450kprobe$chr)
>> names(chs) = paste('chr',chs,sep='')
>> CpGs.450k = with(IlluminaHumanMethylation450kprobe,
>> GRanges(paste('chr',chr,sep=''),
>> IRanges(start=site, width=2, names=Probe_ID),
>> strand=strand))
>>
>> CpGs.450k.df <- as.data.frame(CpGs.450k)
>>
>> mylocData <- RangedData(ranges=IRanges(start=CpGs.450k.df$start, end=CpGs.450k.df$end, names=row.names(CpGs.450k.df)), space=CpGs.450k.df$seqnames)
>>
>> mymethy.obj <- new("MethyGenoSet", locData=mylocData, assayData=assayDataNew(exprs=new("matrix"), methylated=new("matrix"), unmethylated=new("matrix"), detection=new("matrix")))
>>
>> # make sure the CpG ids in mylocData agree with what is in CpGs.450k.df
>>
>> exprs(mymethy.obj) <- myexprs
>> methylated(mymethy.obj) <- mymeth
>> unmethylated(mymethy.obj) <- myunmeth
>> # had to add this afterwards
>>
>> methyGenoSet.sm <- smoothMethyData(mymethy.obj, winSize = 250)
>>
>> I haven't run through all of the MethyAnalysis functions so the MethyGenoSet object may need more information added to it, in other slots, for it to work properly.
>>
>>
>> Lavinia Gordon
>> Senior Research Officer
>> Quantitative Sciences Core, Bioinformatics
>>
>> Murdoch Childrens Research Institute
>> The Royal Children's Hospital
>> Flemington Road Parkville Victoria 3052 Australia
>> T 03 8341 6221
>> www.mcri.edu.au
>>
>>
>> -----Original Message-----
>> From: Martin Morgan [mailto:mtmorgan at fhcrc.org]
>> Sent: Saturday, 9 March 2013 1:49 AM
>> To: Lavinia Gordon
>> Cc: bioconductor at r-project.org; dupan at gmal.com
>> Subject: Re: [BioC] methyAnalysis - changes to GenoSet
>>
>> On 03/07/2013 08:23 PM, Lavinia Gordon wrote:
>>> Dear Bioc users,
>>> I have just tried the example code supplied with the package methyAnalysis and am getting an error:
>>>
>>>> ###################################################
>>>> ### code chunk number 3: Slide-window smoothing of the DNA
>>>> methylation data ###################################################
>>>> methyGenoSet.sm <- smoothMethyData(exampleMethyGenoSet, winSize =
>>>> 250)
>>> Smoothing Chromosome 21 ...
>>>
>>> Note: Method with signature âGenoSet#character#ANY#ANYâ chosen for
>>> function â[â, target signature âMethyGenoSet#character#character#missingâ.
>>> "MethyGenoSet#ANY#ANY#ANY" would also be valid Warning message:
>>> The ranges method on a GenoSet is depricated. Please use space(locData(x)) or seqnames(locData(x)) as appropriate for RangedData or GRanges.
>>>> methyGenoset.sm
>>> Error: object 'methyGenoset.sm' not found
>>
>> typo -- methyGenoSet.sm
>>
>> Martin
>>
>>>> sessionInfo()
>>> R version 2.15.2 (2012-10-26)
>>> Platform: x86_64-redhat-linux-gnu (64-bit)
>>>
>>> locale:
>>> [1] LC_CTYPE=en_US.UTF-8
>>> [2] LC_NUMERIC=C
>>> [3] LC_TIME=en_US.UTF-8
>>> [4] LC_COLLATE=en_US.UTF-8
>>> [5] LC_MONETARY=en_US.UTF-8
>>> [6] LC_MESSAGES=en_US.UTF-8
>>> [7] LC_PAPER=C
>>> [8] LC_NAME=C
>>> [9] LC_ADDRESS=C
>>> [10] LC_TELEPHONE=C
>>> [11] LC_MEASUREMENT=en_US.UTF-8
>>> [12] LC_IDENTIFICATION=C
>>>
>>> attached base packages:
>>> [1] grid stats graphics grDevices
>>> [5] utils datasets methods base
>>>
>>> other attached packages:
>>> [1] methyAnalysis_1.0.0 org.Hs.eg.db_2.8.0
>>> [3] RSQLite_0.11.2 DBI_0.2-5
>>> [5] AnnotationDbi_1.20.5 Biobase_2.18.0
>>> [7] IRanges_1.16.6 BiocGenerics_0.4.0
>>>
>>> loaded via a namespace (and not attached):
>>> [1] affy_1.36.1 affyio_1.26.0
>>> [3] annotate_1.36.0 BiocInstaller_1.8.3
>>> [5] biomaRt_2.14.0 Biostrings_2.26.3
>>> [7] biovizBase_1.6.2 bitops_1.0-5
>>> [9] BSgenome_1.26.1 cluster_1.14.3
>>> [11] colorspace_1.2-1 dichromat_2.0-0
>>> [13] genefilter_1.40.0 GenomicFeatures_1.10.2
>>> [15] GenomicRanges_1.10.7 genoset_1.10.1
>>> [17] Gviz_1.2.1 Hmisc_3.10-1
>>> [19] KernSmooth_2.23-9 labeling_0.1
>>> [21] lattice_0.20-13 lumi_2.10.0
>>> [23] MASS_7.3-23 Matrix_1.0-11
>>> [25] methylumi_2.4.0 mgcv_1.7-22
>>> [27] munsell_0.4 nleqslv_2.0
>>> [29] nlme_3.1-108 parallel_2.15.2
>>> [31] plyr_1.8 preprocessCore_1.20.0
>>> [33] RColorBrewer_1.0-5 RCurl_1.95-4.1
>>> [35] Rsamtools_1.10.2 rtracklayer_1.18.2
>>> [37] scales_0.2.3 splines_2.15.2
>>> [39] stats4_2.15.2 stringr_0.6.2
>>> [41] survival_2.37-4 tools_2.15.2
>>> [43] XML_3.95-0.2 xtable_1.7-1
>>> [45] zlibbioc_1.4.0
>>>
>>> With regards,
>>>
>>> Lavinia Gordon
>>> Senior Research Officer
>>> Quantitative Sciences Core, Bioinformatics
>>>
>>> Murdoch Childrens Research Institute
>>> The Royal Children's Hospital
>>> Flemington Road Parkville Victoria 3052 Australia T 03 8341 6221
>>> www.mcri.edu.au<http://www.mcri.edu.au/>
>>>
>>>
>>> ______________________________________________________________________
>>> This email has been scanned by the Symantec Email Security.cloud service.
>>> For more information please visit http://www.symanteccloud.com
>>> ______________________________________________________________________
>>> [[alternative HTML version deleted]]
>>>
>>>
>>>
>>> _______________________________________________
>>> Bioconductor mailing list
>>> Bioconductor at r-project.org
>>> https://stat.ethz.ch/mailman/listinfo/bioconductor
>>> Search the archives:
>>> http://news.gmane.org/gmane.science.biology.informatics.conductor
>>
>>
>> --
>> Computational Biology / Fred Hutchinson Cancer Research Center
>> 1100 Fairview Ave. N.
>> PO Box 19024 Seattle, WA 98109
>>
>> Location: Arnold Building M1 B861
>> Phone: (206) 667-2793
>>
>> ______________________________________________________________________
>> This email has been scanned by the Symantec Email Security.cloud service.
>> For more information please visit http://www.symanteccloud.com
>>
>> If you have any question, please contact MCRI IT Helpdesk for further assistance.
>> ______________________________________________________________________
>>
>> ______________________________________________________________________
>> This email has been scanned by the Symantec Email Security.cloud service.
>> For more information please visit http://www.symanteccloud.com
>>
>> _______________________________________________
>> Bioconductor mailing list
>> Bioconductor at r-project.org
>> https://stat.ethz.ch/mailman/listinfo/bioconductor
>> Search the archives: http://news.gmane.org/gmane.science.biology.informatics.conductor
>
> _______________________________________________
> Bioconductor mailing list
> Bioconductor at r-project.org
> https://stat.ethz.ch/mailman/listinfo/bioconductor
> Search the archives: http://news.gmane.org/gmane.science.biology.informatics.conductor
More information about the Bioconductor
mailing list