hotfix ANOVA interactions

ANOVA interaction term now included in outputs

Author: grlloyd

DESCRIPTION

@@ -1,7 +1,7 @@
 Package: structToolbox
 Type: Package
 Title: Data processing & analysis tools for Metabolomics and other omics
-Version: 1.19.0
+Version: 1.19.1
 Authors@R: c(
     person(
         c("Gavin","Rhys"),

R/anova_class.R

@@ -144,20 +144,23 @@ setMethod(f="model_apply",
             return(A)
         })
 
+        out_cols=terms(M$formula)
+        out_cols=attributes(out_cols)$term.labels
+        
         f_statistic=sapply(output,function(x){
             x$`F value`
         })
         f_statistic=as.data.frame(t(f_statistic))
         colnames(f_statistic)=rownames(output[[1]])
-        f_statistic=f_statistic[,colnames(y),drop=FALSE]
+        f_statistic=f_statistic[,out_cols,drop=FALSE]
 
 
         p_value=sapply(output,function(x){
             x$`Pr(>F)`
         })
         p_value=as.data.frame(t(p_value))
         colnames(p_value)=rownames(output[[1]])
-        p_value=p_value[,colnames(y),drop=FALSE]
+        p_value=p_value[,out_cols,drop=FALSE]
 
         # fdr correct the p.values
         for (k in 1:ncol(p_value)) {