Hi, I am using PLINK 2 GLM to associate the genetic variants that I have found with a phenotype in the tested population. This is the PLINK command I use:
./plink2 --pfile PGEN_FILE --glm --covar COVARIATES.TAB --tests all --adjust --out OUTPUT_FILE
This is the output in the ".linear" file, listing every covariate in a different line.
I understand most of the output but I'm unsure about the last line, containing "USER_10DF". There is no BETA or SE in that line but a p-value. I assume it is some kind of overall model-based test and that the 10DF stands for 10 degrees of freedom but I can't find clear documentation in the pLINK manual or anywhere else.
The p-value given in the "USER_10DF" line is also the one that is saved in the "linear.adjusted" file under the same variant name (i.e., ID column). Here's the corresponding output from that file:
Is anyone familiar with this output and can confirm what it means?
Thanks much in advance!
- Yes, USER_10DF refers to the joint test you requested with --tests. I will add a note about this to the online documentation.
- https://www.cog-genomics.org/plink/2.0/formats#adjusted describes the .adjusted file format.