This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Error While Using BSmooth in bsseq

I am using the BSmooth.tstat() function in the bsseq R package but I get the following error. I can use the test data without error but not my own data.

  • meth_fit = BSmooth(BSseq = meth, verbose = TRUE)
  • meth_stat = BSmooth.tstat(meth_fit, group1 = c(1, 2, 3), group2 = c(4, 5, 6))

[BSmooth.tstat] preprocessing ... done in 5.7 sec [BSmooth.tstat] computing stats within groups ... done in 3.0 sec [BSmooth.tstat] computing stats across groups ... Error in h(simpleError(msg, call)) : error in evaluating the argument 'args' in selecting a method for function 'do.call': need at least two non-NA values to interpolate

enter image description here

error bsmooth wgbs bsseq methylation
• 1,076 views • link
written by Burak Özyürek • Burak Özyürek

1 answer

Solved that very day by filtering pseudo-chromosomes with chrSelectBSseq() function :)

• 0 views • link
written by Burak Özyürek • Burak Özyürek

Log in to answer this question.

More posts like this
  • Something wrong with my code concerning bsseq visualization
    written by diqixiaoyaoer • diqixiaoyaoer
    Here is my code: #### library(bsseq) library(bsseqData) ## ----showData----------------------------------------------------------------- data(BS.cancer.ex) BS.cancer.ex <- updateObject(BS.cancer.ex) BS.cancer.ex pData(BS.cancer.ex) ## ----smooth,eval=FALSE-------------------------------------------------------- # BS.cancer.ex.fit <- BSmooth( # BSseq = BS.cancer.ex, …
  • BSmooth.tstat error due to NA
    written by ttsutsui1028 • ttsutsui1028
    Dear All, I use bsseq package for WGBS data analysis. After smoothing, I had tried to BSmooth.tstat for computing t-statistics. But BSmooth.tstat gave me following …