Aldex2 error following tutorial - complex study design
0
2
Entering edit mode
@mstagliamonte-18924
Last seen 4.7 years ago

Hi,

I am using Aldex2 Version: 1.17.0 (Date: 2019-2-28) and I am having some issues following the tutorial I found here:

https://bioconductor.org/packages/devel/bioc/vignettes/ALDEx2/inst/doc/ALDEx2_vignette.pdf (check page 7 for complex study design and glm)

This is my command line:

covariates <- data.frame("A" = sample(0:1, 14, replace = TRUE),
                         "B" = c(rep(0, 7), rep(1, 7)))
mm <- model.matrix(~ A + B, covariates)
selex.sub<- selex[1:400,]

Xtestclr <- aldex.clr(selex.sub, mm, mc.samples=8, denom="all")
glmTest <- aldex.glm(Xtestclr, mm)
#
#Features that are differential under any of the variables
#should be validated using aldex.effect:
XtestEffect<- aldex.effect(clr = Xtestclr)

I am receiving the following error:

#Error in aldex.effect(clr = Xtestclr) : 
 # mismatch btw 'length(conditions)' and 'ncol(reads)'

#

I guess the problem is the clr output, which I should probably re-create it for individual variables ? I am not sure about it though., I would appreciate some advice on how to proceed.

Thank you for your kind attention, Max

software error aldex2 • 1.0k views
ADD COMMENT

Login before adding your answer.

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