"MethylKit" package for WGBS data
0
0
Entering edit mode
4.1 years ago

Hello I am working on WGBS data using MethylKIT package. I am trying to analyse few samples whose dimension are not same (see below). What kind of measure can be taken to handle this type of problem

library(methylKit) file.list = list ("JChr21-modd.txt","SChr21-modd.txt","Delta21-mod.txt","BChr21-modd.txt","NChr21-modd.txt") myobj = methRead (file.list, sample.id = list ("JChr21-modd","SChr21-modd","Delta21-mod", "BChr21-modd", "NChr21-modd"), assembly = "hg19", treatment = c(1,1,0,0,0), context = "CpG") meth=unite(myobj, destrand=FALSE) uniting... Error in unite(myobj, destrand = FALSE) : no base were united. try adjusting 'min.per.group’.

dim(Jchr21)

[1] 857692 6

dim(Delta21)

[1] 760888 6

next-gen sequencing • 1.1k views
ADD COMMENT
0
Entering edit mode

You are posting these kinds of question for weeks (maybe months now), don't you think it is time to talk to someone who can guide you at your institution locally, given that you have trouble all the way? Online communities are great for specific questions but hands-on guidance which are apparently need at this point is not in the scope of them. I really say this in your best interest.

ADD REPLY
0
Entering edit mode

It seems methylkit was not able to find anything common between 5 samples. Is your input file tab delimited? If yes try providing seperator as sep = "\t" .

ADD REPLY

Login before adding your answer.

Traffic: 1823 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6