I processed my data along the new tutorial "Performing longitudinal and paired sample comparisons with q2-longitudinal". While pairwise difference comparisons and pairwise distance comparisons worked out well on my data I encountered the following errors with longitudinal linear-mixed-effects:
either
Plugin error from longitudinal:
Singular matrix
Debug info has been saved to /tmp/qiime2-q2cli-err-2t5n23zo.log
or
Plugin error from longitudinal:
'AxesSubplot' object does not support indexing
Debug info has been saved to /tmp/qiime2-q2cli-err-83vjyolh.log
depending on my input for --p-group-categories --p-state-column and --p-individual-id-column.
I would prefer to send you the files directly by email. However I could not find your email contact at Northern Arizona University. Where can I find your contact?
@Mechah thanks for raising these two issues and sending your data to help me debug.
The first issue appears to be related to the independent variables that you are providing to LME. They are not actually independent, but appear to be covariates (e.g., most Descriptions only use one type of ProductTreatment), resulting in this error. If instead you provide --p-group-categories CopiesMW,ProductTreatment, LME will run fine (though I recognize that this is probably biologically meaningless, it's useful for diagnosing which inputs cause the error). Thank you for finding this issue: q2-longitudinal should provide more information on why this error is occurring, and I have opened this ticket to improve in the next release. For now, you should try to identify other, non-covariate factors to use as group-categories.
This issue is a minor bug in the way plots are created/indexed, so thank you for detecting this. The code was written assuming that 2+ group-categories would be used; a ticket for this issue has been raised and will be solved in the next release (next month). For now, LME will only work when 2 or more group-categories are provided.
Thank you for your help in debugging linear-mixed-effects!
Hi @Nicholas_Bokulich,
wow thanks for your fast support!
I played around with other variables and it worked. Nevertheless, Iâm looking forward to the revised version of LME in an upcoming release.
Cheers!
I had the same exact problem! I was receiving the error ââAxesSubplotâ object does not support indexingâ and could not figure out what I was doing wrong. It is working now if I supply two group categories.
I am wondering if is possible to use linear-mixed-effects with a unifrac distance matrix. I saw in the documentation that it is possible to determine the relationship between a specific taxa in a feature table and several variable. If I supplied a unifrac distance matrix as --i-table what would I use for 'âp-metric` and is this allowable?
Hi @Maggi_Mars, linear-mixed-effects cannot accept a DistanceMatrix as input â the table artifact only accepts FeatureTable[RelativeFrequency] artifacts (an OTU table).
linear-mixed-effects does accept pcoa results as input, so it would be possible to test, for example, how PC1 changes over time in relation to your selected variables.
If you are specifically interested in distance between samples, rather than principal coordinates, you could use the pairwise-distances visualization â though that action only examines distance between paired samples at two time points in relation to one variable.
I hope that helps. If not, let us know â q2-diversity or another plugin may be able to accomplish what you want to do.
Both of the q2-longitudinal issues described in this topic (support for single variables, and better error message with invalid input) have been fixed in the 2017.9 release!