Hi all,
I'm just starting to work with baySeq (switching from DESeq2 to see what's different), and I've run into a problem pretty early on with getPriors
> head(all)
S1IP_1 S1IP_2 S1IP_3 S2IP_1 S2IP_2 S2IP_3
[1,] 93 35 92 271 231 350
[2,] 2 3 7 9 7 10
[3,] 0 5 1 0 0 1
[4,] 1 14 9 3 5 0
[5,] 0 0 0 0 1 0
[6,] 5 98 15 17 19 48
> seglens <- input [,5]
> replicates <-c("S1","S1","S1","S2","S2","S2")
> groups <- list(NDE = c(1,1,1,1,1,1), DE = c(1,1,1,2,2,2))
> CD <- new("countData", data = all, seglens = seglens, replicates = replicates, groups = groups)
> CD@annotation <- as.data.frame(cname)
> libsizes(CD) <- getLibsizes(CD)
>CDpriors <- getPriors.NB(CD, samplesize = 10^5, estimation = "QL", cl=cl)
Finding priors...Error in checkForRemoteErrors(val) :
one node produced an error: 'xmin' not less than 'xmax'
Any ideas about what's going on and how I can fix it?
I'm just starting to work with baySeq (switching from DESeq2 to see what's different), and I've run into a problem pretty early on with getPriors
> head(all)
S1IP_1 S1IP_2 S1IP_3 S2IP_1 S2IP_2 S2IP_3
[1,] 93 35 92 271 231 350
[2,] 2 3 7 9 7 10
[3,] 0 5 1 0 0 1
[4,] 1 14 9 3 5 0
[5,] 0 0 0 0 1 0
[6,] 5 98 15 17 19 48
> seglens <- input [,5]
> replicates <-c("S1","S1","S1","S2","S2","S2")
> groups <- list(NDE = c(1,1,1,1,1,1), DE = c(1,1,1,2,2,2))
> CD <- new("countData", data = all, seglens = seglens, replicates = replicates, groups = groups)
> CD@annotation <- as.data.frame(cname)
> libsizes(CD) <- getLibsizes(CD)
>CDpriors <- getPriors.NB(CD, samplesize = 10^5, estimation = "QL", cl=cl)
Finding priors...Error in checkForRemoteErrors(val) :
one node produced an error: 'xmin' not less than 'xmax'
Any ideas about what's going on and how I can fix it?