Question: Generation of Dendogram on WGBS data
0
gravatar for shrinka.genetics
11 weeks ago by
shrinka.genetics0 wrote:

Hi I am working on WGBS data and using MethylKit package. Methykit has its own function to make a dendogram on CpG methylation, but I need bootstraping on it. I used these codes.

library(methylKit) file.list = list (........) myobj = methRead (file.list, sample.id = list (..........), assembly = "hg19", treatment = c(), context = "CpG") meth=unite(myobj, destrand=FALSE) data<-clusterSamples(meth, dist= "correlation", method= "ward", plot=TRUE) library(pvclust) result <- pvclust(data, nboot=10000)

These are the usage of these arguments

1) clusterSamples(.Object, dist="correlation", method="ward", sd.filter=TRUE,sd.threshold=0.5, filterByQuantile=TRUE, plot=TRUE,chunk.size)

2) pvclust(data, method.hclust="average", method.dist="correlation", use.cor="pairwise.complete.obs", nboot=1000, parallel=FALSE, r=seq(.5,1.4,by=.1), store=FALSE, weight=FALSE, iseed=NULL, quiet=FALSE)

I have tried various ways, sometimes by putting meth directly or sometime that data object and tried it, but every time getting error something or other.

Anyone could you please help me regarding this? Any suggestion, that how can I modify my codes to have my result.

sequencing • 138 views
ADD COMMENTlink modified 20 days ago by iwwilson0 • written 11 weeks ago by shrinka.genetics0
0
gravatar for iwwilson
20 days ago by
iwwilson0
iwwilson0 wrote:

This article may be of use (I'm about to work through it myself). Check out Additional File 7 - they included the code used to pass MethylKit's dendrogram to pvclust:

https://bmcgenomics.biomedcentral.com/articles/10.1186/s12864-019-5578-4

ADD COMMENTlink modified 20 days ago • written 20 days ago by iwwilson0
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 2008 users visited in the last hour