From Genome Analysis Wiki
Jump to navigationJump to search
48 bytes added
, 16:05, 5 July 2012
Line 38: |
Line 38: |
| % emmax -v -d 10 -t [tped_prefix] -p [pheno_file] -k [kin_file] -o [out_prefix] | | % emmax -v -d 10 -t [tped_prefix] -p [pheno_file] -k [kin_file] -o [out_prefix] |
| | | |
− | This will generate the following files:
* [out_prefix].reml : REML output with 6 lines, where each line represents (1) Log-likelihood with variance component (2) Log-likelihood without variance component, (3) \delta = \sigma_e^2 / \sigma_g^2 (Ratio between variance parameters) (4) \sigma_g^2 (genetic variance parameter), and (5) \sigma_e^2 (residual variance parameter), and (6) The pseudo-heritability estimates
* [out_prefix].ps : Each line consist of [SNP ID], [beta], [SE-beta], [p-value]. | + | This will generate the following files: |
| + | * [out_prefix].reml : REML output with 6 lines, where each line represents |
| + | ##Log-likelihood with variance component |
| + | ##Log-likelihood without variance component |
| + | ## \delta = \sigma_e^2 / \sigma_g^2 (Ratio between variance parameters) |
| + | ## \sigma_g^2 (genetic variance parameter) |
| + | ## sigma_e^2 (residual variance parameter) |
| + | ## The pseudo-heritability estimates
. (Explained variance by the kinship matrix) |
| + | |
| + | * [out_prefix].ps : Each line consist of |
| + | ## SNP ID |
| + | ## Beta (1 is effect allele) |
| + | ## SE(beta) |
| + | ## p-value. |
| | | |
| === Incorporating Covariates === | | === Incorporating Covariates === |