cozygene / FEAST

Fast expectation maximization for microbial source tracking
Other
118 stars 60 forks source link

Error in if (sum(x[i, ]) > maxdepth) { : missing value where TRUE/FALSE needed #51

Closed rubmc97 closed 1 year ago

rubmc97 commented 1 year ago

Hello,

I'm trying to run FEAST with a 16S amplicon sequencing dataset with 2 sinks and 2 sources. Looks like I messed something up with the "C" and "metadata" as I keep getting the error message "Error in if (sum(x[i, ]) > maxdepth) { : missing value where TRUE/FALSE needed" when running "feast_out = FEAST::FEAST(C=otus, metadata = meta_feast, outfile = "feast_out", different_sources_flag = 1)"

Any tips or help would be very much appreciated,

Best, Rubén

liashenhav commented 1 year ago

Hi Rubén,

Can you share your metadata and otu matrix? That would help me debug it.

Thanks, Liat

rubmc97 commented 1 year ago

Hi Liat,

Thanks a lot for the quick response. Looks like it was my fault for not understanding how to correctly label the samples with the "id" number. Did not realise that several sinks cannot have the same number. Your response to another user was quite helpful. I'll close the issue.

Best, Rubén