Linear mixed effects with individual features

Hi,

I am encountering a strange error when running qiime longitudinal linear-mixed-effects with individual feature IDs. For all feature IDs beginning with a letter the command runs fine, however if the feature ID begins with a number I get the following error:

Plugin error from longitudinal:

invalid syntax (, line 1)
Debug info has been saved to /tmp/qiime2-q2cli-err-lrhc59gy.log (or other file)

I have attached some of the debug files incase they help.
qiime2-q2cli-err-lrhc59gy.log.txt (3.0 KB)
qiime2-q2cli-err-32rrmec8.log.txt (3.0 KB)
qiime2-q2cli-err-lrhc59gy.log.txt (3.0 KB)

I have also pasted the loop I am running at the bottom of the post, however as I mentioned before all features beginning with a letter worked fine and produced a correct output qzv, whilst all beginning with a number failed.

Thank you,
Jen

for i in {'3de5f525094ec393ec1767ae0d16de8d','6e5ebd83a964b0bee237cd9b221aedce','f22088afb95783a66916935ad3be16ee','18d1d3003cc423b568ba5b5ccb433023','f257f8bb0d9b4aabfa6ccd6ec84b256d','7b58a25f42f2aa49626ddffe0e69ae7b','a2ed5d5b89e82d3e6289e05cbbdf31a8','f3028e6e91b09bef83702f3ed00265b1','6089e3455fdfb31efad4063931738ce7','3bb9edc70f70cd01c186d491c64be80f','d7f3bdeea12378c23ab4dfa79d47aca5','2d2ba432de65fdeac8600f615ad13175','1ad106850823904ec6b0677ec26a6bf4','e8c6023a59abc1f9a47712651fecf706','2e749ec8309ac6c24c13f39cdbc9e3ad','2bece401405f1792df34008de663fcfc','fa8ad273584a7d5f644b2aaa70acf1fa','64c48d34bfc2889dfe1788d653ecda6d','a81eeef79efbeed60fe844aa4e71619b','dae0dac699437281dbcdf38b98735d0f','91dd9ad005a3e7c18545584be1a3ba33','e85b9d0833bdc9dbc76e217b4e40fca'}

do
qiime longitudinal linear-mixed-effects
--m-metadata-file Metadata/VMR_sample_metadata.txt
--i-table QIIME2_analysis/dada2_output/table_filtered_MC_relative_frequency.qza
--p-metric "$i"
--p-state-column Week_number
--p-individual-id-column Calf_ID
--p-group-columns Milk_replacer_ID
--o-visualization QIIME2_analysis/GNEISS_analysis/sig_balance_taxonomies/lme_y51_features/"$i"_week_number_milk_replacer.qzv;
done

2 Likes

Nice bug catch, @JenKelly! :bug:

And thank you for sharing your code and error logs. Nice detective work on figuring out that it is just numeric feature names! Your help should make this quick and easy to fix.:mag:

I have raised an issue here to track this. Stay tuned for more details! :radio:

Thanks

1 Like

This topic was automatically closed 31 days after the last reply. New replies are no longer allowed.