Question: "MethylKit" package for WGBS data
0
gravatar for shrinka.genetics
6 weeks ago by
shrinka.genetics0 wrote:

Hello I am working with WGBS data and using MethyKit package. I have 12 samples in 3 groups. I used codes like

> file.list = list ("JChr21-modd.txt", "BChr21-modd.txt", "SChr21-modd.txt", "NChr21-modd.txt","KChr21-modd.txt","EChr21-modd.txt","CChr21-modd.txt","PChr21-modd.txt", "IChr21-modd.txt","RChr21-modd.txt","LChr21-modd.txt","MChr21-modd.txt")

> myobj = methRead (file.list, sample.id = list ("JChr21-modd", "BChr21-modd", "SChr21-modd", "NChr21-modd","KChr21-modd","EChr21-modd", "CChr21-modd", "PChr21-modd","IChr21-modd","RChr21-modd", "LChr21-modd", "MChr21-modd"), assembly = "hg19", treatment = c(2,2,2,2,1,1,1,1,0,0,0,0), context = "CpG")

> meth=unite(myobj, destrand=FALSE)  uniting...

> PCASamples (meth)

Here they are giving me PC1 Vs PC2 plot. My question is how can I get PC1 Vs PC3, PC2 Vs PC3 etc others PCA plots, by changing arguments or what?

Thanks Shrinka

sequencing next-gen • 113 views
ADD COMMENTlink modified 6 weeks ago by Friederike6.2k • written 6 weeks ago by shrinka.genetics0
1
gravatar for Friederike
6 weeks ago by
Friederike6.2k
United States
Friederike6.2k wrote:

My question is how can I get PC1 Vs PC3, PC2 Vs PC3 etc others PCA plots, by changing arguments or what?

Yes.

If you type ?PCASamples you should see the documentation of that function. This will explain that there is a comp parameter, which can be used to define the components to be plotted. Without having used the function I'll guess that that's what you're looking for.

ADD COMMENTlink written 6 weeks ago by Friederike6.2k
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: 1955 users visited in the last hour