Hi @cherman2,
I've just run the verbose flag like this:
qiime longitudinal linear-mixed-effects \
--m-metadata-file filtered-sample-metadata.tsv \
--m-metadata-file core-metrics-results/faith_pd_vector.qza \
--p-metric faith_pd \
--p-group-columns sample-type,treatment-group,administration-route,gender \
--p-state-column week \
--p-random-effects age,weight,height,bmi \
--p-individual-id-column subject-id \
--o-visualization linear-mixed-effects-faith.qzv \
--verbose
And this was the output, which to be honest I don't understand
Traceback (most recent call last):
File "/home/marina/miniconda3/envs/qiime2-2023.5/lib/python3.8/site-packages/q2cli/commands.py", line 468, in __call__
results = action(**arguments)
File "<decorator-gen-521>", line 2, in linear_mixed_effects
File "/home/marina/miniconda3/envs/qiime2-2023.5/lib/python3.8/site-packages/qiime2/sdk/action.py", line 274, in bound_callable
outputs = self._callable_executor_(
File "/home/marina/miniconda3/envs/qiime2-2023.5/lib/python3.8/site-packages/qiime2/sdk/action.py", line 558, in _callable_executor_
ret_val = self._callable(output_dir=temp_dir, **view_args)
File "/home/marina/miniconda3/envs/qiime2-2023.5/lib/python3.8/site-packages/q2_longitudinal/_longitudinal.py", line 200, in linear_mixed_effects
model_summary, model_results, model_fit, formula = _linear_effects(
File "/home/marina/miniconda3/envs/qiime2-2023.5/lib/python3.8/site-packages/q2_longitudinal/_utilities.py", line 353, in _linear_effects
mlm = mixedlm(
File "/home/marina/miniconda3/envs/qiime2-2023.5/lib/python3.8/site-packages/statsmodels/regression/mixed_linear_model.py", line 1046, in from_formula
mod = super(MixedLM, cls).from_formula(
File "/home/marina/miniconda3/envs/qiime2-2023.5/lib/python3.8/site-packages/statsmodels/base/model.py", line 229, in from_formula
mod = cls(endog, exog, *args, **kwargs)
File "/home/marina/miniconda3/envs/qiime2-2023.5/lib/python3.8/site-packages/statsmodels/regression/mixed_linear_model.py", line 741, in __init__
super(MixedLM, self).__init__(endog, exog, groups=groups,
File "/home/marina/miniconda3/envs/qiime2-2023.5/lib/python3.8/site-packages/statsmodels/base/model.py", line 270, in __init__
super().__init__(endog, exog, **kwargs)
File "/home/marina/miniconda3/envs/qiime2-2023.5/lib/python3.8/site-packages/statsmodels/base/model.py", line 95, in __init__
self.data = self._handle_data(endog, exog, missing, hasconst,
File "/home/marina/miniconda3/envs/qiime2-2023.5/lib/python3.8/site-packages/statsmodels/base/model.py", line 135, in _handle_data
data = handle_data(endog, exog, missing, hasconst, **kwargs)
File "/home/marina/miniconda3/envs/qiime2-2023.5/lib/python3.8/site-packages/statsmodels/base/data.py", line 675, in handle_data
return klass(endog, exog=exog, missing=missing, hasconst=hasconst,
File "/home/marina/miniconda3/envs/qiime2-2023.5/lib/python3.8/site-packages/statsmodels/base/data.py", line 72, in __init__
arrays, nan_idx = self.handle_missing(endog, exog, missing,
File "/home/marina/miniconda3/envs/qiime2-2023.5/lib/python3.8/site-packages/statsmodels/base/data.py", line 257, in handle_missing
raise ValueError("Shape mismatch between endog/exog "
ValueError: Shape mismatch between endog/exog and extra 2d arrays given to model.
Plugin error from longitudinal:
Shape mismatch between endog/exog and extra 2d arrays given to model.
See above for debug info.
The output for the command without the random effects parameters was the following:
Traceback (most recent call last):
File "/home/marina/miniconda3/envs/qiime2-2023.5/lib/python3.8/site-packages/q2cli/commands.py", line 468, in __call__
results = action(**arguments)
File "<decorator-gen-521>", line 2, in linear_mixed_effects
File "/home/marina/miniconda3/envs/qiime2-2023.5/lib/python3.8/site-packages/qiime2/sdk/action.py", line 274, in bound_callable
outputs = self._callable_executor_(
File "/home/marina/miniconda3/envs/qiime2-2023.5/lib/python3.8/site-packages/qiime2/sdk/action.py", line 558, in _callable_executor_
ret_val = self._callable(output_dir=temp_dir, **view_args)
File "/home/marina/miniconda3/envs/qiime2-2023.5/lib/python3.8/site-packages/q2_longitudinal/_longitudinal.py", line 200, in linear_mixed_effects
model_summary, model_results, model_fit, formula = _linear_effects(
File "/home/marina/miniconda3/envs/qiime2-2023.5/lib/python3.8/site-packages/q2_longitudinal/_utilities.py", line 353, in _linear_effects
mlm = mixedlm(
File "/home/marina/miniconda3/envs/qiime2-2023.5/lib/python3.8/site-packages/statsmodels/regression/mixed_linear_model.py", line 1046, in from_formula
mod = super(MixedLM, cls).from_formula(
File "/home/marina/miniconda3/envs/qiime2-2023.5/lib/python3.8/site-packages/statsmodels/base/model.py", line 229, in from_formula
mod = cls(endog, exog, *args, **kwargs)
File "/home/marina/miniconda3/envs/qiime2-2023.5/lib/python3.8/site-packages/statsmodels/regression/mixed_linear_model.py", line 803, in __init__
self.endog_li = self.group_list(self.endog)
File "/home/marina/miniconda3/envs/qiime2-2023.5/lib/python3.8/site-packages/statsmodels/regression/mixed_linear_model.py", line 1104, in group_list
return [np.array(array[self.row_indices[k]])
File "/home/marina/miniconda3/envs/qiime2-2023.5/lib/python3.8/site-packages/statsmodels/regression/mixed_linear_model.py", line 1104, in <listcomp>
return [np.array(array[self.row_indices[k]])
IndexError: index 71 is out of bounds for axis 0 with size 71
Plugin error from longitudinal:
index 71 is out of bounds for axis 0 with size 71
See above for debug info.