A bug related to R factor

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…