Post by luba on May 6, 2019 12:53:58 GMT
Dear all,
I am running a conditional GWAS using actual genotypes.
The program seems to work but it does not give me any output.
I am uSing the latest version of gcta and commands"
--bfile /home/xx/xx/Plinkbedfiles ###
--maf 0.03
--cojo-file METAANALYIS.txt
--cojo-slct
--cojo-actual-geno
--cojo-p 0.0005
--out test
--thread-num 6
--maf 0.03
--cojo-file METAANALYIS.txt
--cojo-slct
--cojo-actual-geno
--cojo-p 0.0005
--out test
--thread-num 6
I get some log but not output. Am I doing something wrong?
REGARDS,
lUBA
HERE THE LOG i GET
Reading GWAS summary-level statistics from [METAANALYIS.txt] ...
GWAS summary statistics of 8024069 SNPs read from [ METAANALYIS.txt].
Phenotypic variance estimated from summary statistics of all 8024069 SNPs: 256.766 (variance of logit for case-control studies).
Matching the GWAS meta-analysis results to the genotype data ...
Warning: can't match the reference alleles or of 4119 SNPs to those in the genotype data. These SNPs have been saved in [test.badsnps].
5037535 SNP(s) have large difference of allele frequency among the GWAS summary data and the reference sample. These SNPs have been saved in [test.freq.badsnps].
861871 SNPs are matched to the genotype data.
Calculating the variance of SNP genotypes ...
Performing stepwise model selection on 861871 SNPs to select association signals ... (p cutoff = 0.0005; collinearity cutoff = 0.9)
5 associated SNPs have been selected.
10 associated SNPs have been selected.
15 associated SNPs have been selected.
20 associated SNPs have been selected.
25 associated SNPs have been selected.
30 associated SNPs have been selected.
35 associated SNPs have been selected.
40 associated SNPs have been selected.
45 associated SNPs have been selected.
50 associated SNPs have been selected.
GWAS summary statistics of 8024069 SNPs read from [ METAANALYIS.txt].
Phenotypic variance estimated from summary statistics of all 8024069 SNPs: 256.766 (variance of logit for case-control studies).
Matching the GWAS meta-analysis results to the genotype data ...
Warning: can't match the reference alleles or of 4119 SNPs to those in the genotype data. These SNPs have been saved in [test.badsnps].
5037535 SNP(s) have large difference of allele frequency among the GWAS summary data and the reference sample. These SNPs have been saved in [test.freq.badsnps].
861871 SNPs are matched to the genotype data.
Calculating the variance of SNP genotypes ...
Performing stepwise model selection on 861871 SNPs to select association signals ... (p cutoff = 0.0005; collinearity cutoff = 0.9)
5 associated SNPs have been selected.
10 associated SNPs have been selected.
15 associated SNPs have been selected.
20 associated SNPs have been selected.
25 associated SNPs have been selected.
30 associated SNPs have been selected.
35 associated SNPs have been selected.
40 associated SNPs have been selected.
45 associated SNPs have been selected.
50 associated SNPs have been selected.