Entering edit mode
Hello. I ran into a discrepancy between the manual for EBseq and the output genarated. When I run the tutorial code, instead of a list of the parameters used for alpha, beta and P I get:
Alpha: A single value Beta: a large list of values that does not matches the number of iterations P: two values this is the code, straight from the manual:
data(GeneMat)
str(GeneMat)
Sizes=MedianNorm(GeneMat)
Conditions=as.factor(rep(c("C1","C2"),each=5))
EBOut=EBTest(Data=GeneMat, Conditions=as.factor(rep(c("C1","C2"),each=5)),sizeFactors=Sizes, maxround=5)
EBOut$Alpha
EBOut$Beta
EBOut$P
Otherwisze the script runs as expected. Is it something wrong or am I missing something? Thanks!