ggloor / ALDEx2_dev

ALDEx tool to examine compositional high-throughput sequence data with Welch's t-test
GNU Affero General Public License v3.0
12 stars 6 forks source link

Does Aldex2 supports confounding covariate adjustment? #22

Open Ivy-ops opened 3 years ago

Ivy-ops commented 3 years ago

Hi, May I know does Aldex2 supports covariates adjustment? I see there is a function aldex.glm, is it the correct function to achieve this goal? If yes, which column is the correct pvalue to be extracted? Thanks!

(Intercept) Estimate (Intercept) Std. Error (Intercept) t value (Intercept) Pr(>|t|) model.A Estimate S:D:A:D 2.32201503 1.6048552 1.44686884 0.17582287 -0.8645702 S:D:A:E 2.34043611 0.9274932 2.52339970 0.02830383 -0.6441754 S:E:A:D 1.03699020 0.9955797 1.04159433 0.31995195 1.5643757 S:E:A:E 2.67510968 1.0422060 2.56677633 0.02620043 -1.1223356 S:D:C:D -0.01981533 0.7282679 -0.02720884 0.97878052 1.0472986 S:D:C:E 1.02661028 0.7654643 1.34116026 0.20690355 -0.3351250 model.A Std. Error model.A t value model.A Pr(>|t|) model.B Estimate model.B Std. Error model.B t value S:D:A:D 1.8223933 -0.4744147 0.6444825 0.32802184 1.8037016 0.1818604 S:D:A:E 1.0532149 -0.6116277 0.5532148 2.91606412 1.0424124 2.7974189 S:E:A:D 1.1305305 1.3837537 0.1938683 3.78552532 1.1189350 3.3831503 S:E:A:E 1.1834770 -0.9483375 0.3633269 5.48445133 1.1713384 4.6822090 S:D:C:D 0.8269846 1.2664064 0.2315279 -0.09481651 0.8185025 -0.1158414 S:D:C:E 0.8692229 -0.3855455 0.7071849 0.06993148 0.8603076 0.0812866 model.B Pr(>|t|) (Intercept) Pr(>|t|).BH model.A Pr(>|t|).BH model.B Pr(>|t|).BH S:D:A:D 0.8590002535 1 1 1.0000000 S:D:A:E 0.0173547079 1 1 1.0000000 S:E:A:D 0.0061084312 1 1 1.0000000 S:E:A:E 0.0006691559 1 1 0.2542792 S:D:C:D 0.9098659327 1 1 1.0000000 S:D:C:E 0.9366741895 1 1 1.0000000