PLINK multivariable GLM output: p-value without BETA and SE
0
1
Entering edit mode
22 days ago

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 

pLINK output

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:

enter image description here

Is anyone familiar with this output and can confirm what it means?

Thanks much in advance!

linear-model glm plink • 471 views
ADD COMMENT
0
Entering edit mode
  1. Yes, USER_10DF refers to the joint test you requested with --tests. I will add a note about this to the online documentation.
  2. https://www.cog-genomics.org/plink/2.0/formats#adjusted describes the .adjusted file format.
ADD REPLY

Login before adding your answer.

Traffic: 1690 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6