Error in 'differential_abundance'
0
0
Entering edit mode
strep • 0
@688ee615
Last seen 2.1 years ago
United Kingdom

Hi there - asking for some help with an error that has come up using the 'differential_abundance' command. (nb: I am a beginner)

I have my code

differential_abundance(physeqDE, grouping_column = "Compartment", pvalue.threshold = 0.05, lfc.threshold = 0, filename = F, output_norm = "log-relative")

And I get the output error

converting counts to integer mode
estimating size factors
estimating dispersions
gene-wise dispersion estimates
mean-dispersion relationship
final dispersion estimates
fitting model and testing
-- replacing outliers and refitting for 4 genes
-- DESeq argument 'minReplicatesForReplace' = 7 
-- original counts are preserved in counts(dds)
estimating dispersions
fitting model and testing
Error in y - ymean : non-numeric argument to binary operator
In addition: Warning messages:
1: In DESeqDataSet(se, design = design, ignoreRank) :
  some variables in design formula are characters, converting to factors
2: In randomForest.default(m, y, ...) :
  The response has five or fewer unique values.  Are you sure you want to do regression?
3: In mean.default(y) : argument is not numeric or logical: returning NA

I have a colleague who has used this command in the past and it has worked for them but they get the same error as well when they run their code, whilst still on their original environment too.

Any help would be appreciated :)

phyloseq devtools • 677 views
ADD COMMENT
0
Entering edit mode

y - ymean is not coming from anywhere in DESeq2 codebase, so I can't say the origin of the error.

ADD REPLY

Login before adding your answer.

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