Note a bug in my code today. Sometimes you need to put a certain level (e.g. healthy control) in the first position for your covariance.┬áHere is my old code:dds[[variable]]=factor(dds[[variable]])levels(dds[[variable]])= union(variable_REF, levels…

Continue reading: A bug related to R …