I am running fitZig for differential abundance testing. When I place:
mod = model.matrix (Treatment + Before + normFactor), the answer is:
Error: $ operator is invalid for atomic vectors
Based on internet there is an error when applying
$ to a vector. But how can I fix the problem in the script?