Hi all,
I have problem with quantitative trait association (plink).I used the following command to run this test:
./plink --file file_1 --assoc
I am getting sample.assoc file which is not the output. output file looks something like this:
CHR SNP BP A1 F_A F_U A2 CHISQ P OR
1 rs10458597 564621 T 0.012 NA C NA NA NA
1 rs12565286 721290 C 0.02846 NA G NA NA NA
1 rs12082473 740857 A 0.01181 NA G NA NA NA
in place of NA there should be numerical values which is not in this case.
Thanks in advance
1 answer
Bases on your output, your command is doing a basic case-control test, where the case/control status is defined in the 6th column of your ped file (taking values 0,1 or 2). In your file, the 6th column probably only contain the value 2, explaining why only F_A contains a non missing value. See this link for case/control tests:
http://pngu.mgh.harvard.edu/~purcell/plink/anal.shtml#cc
You want to do a quantitative trait, make sure the 6th column contains values other than 0,1,2, or else specify a .phen file. See this link:
Log in to answer this question.
-
The result of plink --freq is filled with NA
written by yoshifumimiya •
-
PLINK : Problem with .assoc file
written by cdg •
-
Odd Ratio calculation in plink
written by dandanli0365 •
-
Association test in Plink
written by biogirl
-
CrossMap liftover: REF alleles in VCF files overwritten by reference genome file.
written by freezingsleet97 •
-
Obtain A1 and A2 alleles from plink output
written by cdeboever3 •
-
Plot with locuszoom and chr:pos SNPs
written by alesssia
-
Permutation using Plink
written by Molly_K
-
Problems with sex determination with genotyping data Illumina
written by vicbp1 •
-
PLINK returns NA values for all minor allele frequencies in my data
written by yorgos.athanasiadis •
Could you please paste few lines from your input here?
input file test.map looks like this:
I need few lines from test.map (already provided by you, above) and test.ped (for the snps in example map) so that I can reproduce the error. There are multiple columns with NA. Which column are you expecting to have values?
Thanks. I guess your question is addressed above.