Post by amy on Jun 13, 2018 19:16:16 GMT
I'm having trouble getting gcta to match the list of variants from my conditional list to the genotype and summary file, although it looks like the SNPs are written the same way in each of the files.
As an example, here is one variant "1:25768937" which is listed in the genotype, the summary, and the conditional SNP list:
$ grep 1:25768937 AF5.chr1.bim
1 1:25768937 0 25768937 A G
$ grep 1:25768937 AA.HDL.CurSmk.Int.1.txt
1:25768937 A G 0.2875 -0.0036 0.0073 0.6239 23748
$ grep 1:25768937 Known.Conditional.txt
1:25768937
The output that I receive is pasted below:
Options:
--bfile AF5.chr1
--extract AA.HDL.CS.1.AAref.candidate.txt
--cojo-file AA.HDL.CurSmk.Int.1.txt
--cojo-cond Known.Conditional.txt
--out OUTPUT/AA.HDL.CS.Int.Conditioned
Reading PLINK FAM file from [AF5.chr1.fam].
8425 individuals to be included from [AF5.chr1.fam].
Reading PLINK BIM file from [AF5.chr1.bim].
1573279 SNPs to be included from [AF5.chr1.bim].
Reading a list of SNPs from [AA.HDL.CS.1.AAref.candidate.txt].
1 SNPs are extracted from [AA.HDL.CS.1.AAref.candidate.txt].
Reading PLINK BED file from [AF5.chr1.bed] in SNP-major format ...
Genotype data for 8425 individuals and 1 SNPs to be included from [AF5.chr1.bed].
Reading GWAS summary-level statistics from [AA.HDL.CurSmk.Int.1.txt] ...
GWAS summary statistics of 14607517 SNPs read from [AA.HDL.CurSmk.Int.1.txt].
Phenotypic variance estimated from summary statistics of all 14607517 SNPs: 0.567725 (variance of logit for case-control studies).
Matching the GWAS meta-analysis results to the genotype data ...
1 SNPs are matched to the genotype data.
Calculating allele frequencies ...
Calculating the variance of SNP genotypes ...
Reading a list of given SNPs from [Known.Conditional.txt].
Error: none of the given SNPs can be matched to the genotype and summary data.
I can't think of what else to check to see why these aren't matching. Any ideas?
As an example, here is one variant "1:25768937" which is listed in the genotype, the summary, and the conditional SNP list:
$ grep 1:25768937 AF5.chr1.bim
1 1:25768937 0 25768937 A G
$ grep 1:25768937 AA.HDL.CurSmk.Int.1.txt
1:25768937 A G 0.2875 -0.0036 0.0073 0.6239 23748
$ grep 1:25768937 Known.Conditional.txt
1:25768937
The output that I receive is pasted below:
Options:
--bfile AF5.chr1
--extract AA.HDL.CS.1.AAref.candidate.txt
--cojo-file AA.HDL.CurSmk.Int.1.txt
--cojo-cond Known.Conditional.txt
--out OUTPUT/AA.HDL.CS.Int.Conditioned
Reading PLINK FAM file from [AF5.chr1.fam].
8425 individuals to be included from [AF5.chr1.fam].
Reading PLINK BIM file from [AF5.chr1.bim].
1573279 SNPs to be included from [AF5.chr1.bim].
Reading a list of SNPs from [AA.HDL.CS.1.AAref.candidate.txt].
1 SNPs are extracted from [AA.HDL.CS.1.AAref.candidate.txt].
Reading PLINK BED file from [AF5.chr1.bed] in SNP-major format ...
Genotype data for 8425 individuals and 1 SNPs to be included from [AF5.chr1.bed].
Reading GWAS summary-level statistics from [AA.HDL.CurSmk.Int.1.txt] ...
GWAS summary statistics of 14607517 SNPs read from [AA.HDL.CurSmk.Int.1.txt].
Phenotypic variance estimated from summary statistics of all 14607517 SNPs: 0.567725 (variance of logit for case-control studies).
Matching the GWAS meta-analysis results to the genotype data ...
1 SNPs are matched to the genotype data.
Calculating allele frequencies ...
Calculating the variance of SNP genotypes ...
Reading a list of given SNPs from [Known.Conditional.txt].
Error: none of the given SNPs can be matched to the genotype and summary data.
I can't think of what else to check to see why these aren't matching. Any ideas?