WGBS Data Analysis Related
1
0
Entering edit mode
2.3 years ago

Hello I am working on WGBS data. I am using MethyKit package. I have created myobj as well as meth, but after that I got this error when doing clustering.

ClusterSamples(meth, dist= "correlation", method= "ward", plot=TRUE)

Error in ClusterSamples(meth, dist = "correlation", method = "ward", plot = TRUE) :

could not find function "ClusterSamples"

My question is there any other library needed to solve this? How can I solve it

Waiting for help

Regards

Shrinka Sen

Post Doctoral Fellow

BC Cancer Research

Vancouver

MethyKit package • 710 views
ADD COMMENT
0
Entering edit mode
2.3 years ago
ATpoint 64k

Lowercase c => clusterSamples

ADD COMMENT
0
Entering edit mode

I used lower case only. Thanks

ADD REPLY
0
Entering edit mode

library(methylKit) ? Or you use a methylKit version in which the function has not been included yet. Try updating. Those are the possibilities you have.

ADD REPLY
0
Entering edit mode

Thanks for your reply. How to update this?

Regards Shrinka

ADD REPLY
0
Entering edit mode
BiocManager::install("methylKit")
ADD REPLY
0
Entering edit mode

After doing this below mentioned error is showing, previously it was not there. I created myobj and with that did correlation plot. So now what to do? "NChr3-mod.txt", "KChr3-mod.txt", "CChr3-mod.txt", "EChr3-mod.txt", "PChr3-mod.txt") "SChr3-mod", "NChr3-mod", "KChr3-mod", "CChr3-mod", "EChr3-mod", "PChr3-mod"), assembly = "hg19", treatment = c(1,1,1,1,0,0,0,0), context = "CpG")

Error in data.table::fread(file = filepath, ...) :

Opened 168.6MB (176769180 bytes) file ok but could not memory map it. This is a 64bit process. There is probably not enough contiguous virtual memory available.

ADD REPLY
0
Entering edit mode

See the problem is not about the updated version of the package I think. The memory related issue what I mentioned in my last mail is not there if I use smaller file like below. Each Chr3 file was >170000 KB whereas each *Chr21 file is ~40000 KB. *Initially I used the command with Chr 3, now I ran it with Chr 21. Here I was able to create myobj "NChr21-mod.txt", "KChr21-mod.txt", "CChr21-mod.txt", "EChr21-mod.txt", "PChr21-mod.txt") "BChr21-mod", "SChr21-mod", "NChr21-mod", "KChr21-mod", "CChr21-mod", "EChr21-mod", "PChr21-mod"), assembly = "hg19", treatment = c(1,1,1,1,0,0,0,0), context = "CpG") uniting...

Error in quantile.default(sds, sd.threshold) : missing values and NaN's not allowed if 'na.rm' is FALSE

Till there is the same issue. There is no NAs in the input files, I checked. Now how can I solve this?

ADD REPLY

Login before adding your answer.

Traffic: 897 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