Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Using --mmscore and --interaction will not output a column with cov_SNP_int_SNP_covariate #2

Open
lckarssen opened this issue Nov 4, 2015 · 0 comments
Milestone

Comments

@lckarssen
Copy link
Member

This is a 'liftover' of bug #6045 from the R-forge bug tracker.

This was reported several times in the forum:

http://forum.genabel.org/viewtopic.php?f=10&t=830&p=1662
http://forum.genabel.org/viewtopic.php?f=10&t=896
http://forum.genabel.org/viewtopic.php?f=10&t=885&p=1881

Some test output:

../src/palinear -p testphenommscore.phe -d inputfiles/mmscore_gen.dose.fvi -i inputfiles/mmscore_gen.mlinfo -o lck2 --mmscore=inputfiles/mmscore_InvSigma_aj.sex.age.dat --interaction=1
probabel v. 0.4.4
(C) Yurii Aulchenko, Lennart C. Karssen, Maarten Kooyman, Maksim Struchalin, The GenABEL team, EMC Rotterdam

Using EIGEN version 3.2.0 for matrix operations
Options in effect:
--pheno = testphenommscore.phe
--info = inputfiles/mmscore_gen.mlinfo
--dose = inputfiles/mmscore_gen.dose.fvi
(using FVF data)
--ntraits = 1
--ngpreds = 1
--interaction = 1
--interaction_only = 0
--mmscore = inputfiles/mmscore_InvSigma_aj.sex.age.dat
--map = not in output
--nids = estimated from data
--chrom = not in output
--out = lck2
--skipd = 2
--separat = " "
--score = OFF
--nohead = OFF
--allcov = OFF
--robust = OFF

Reading info data...
Number of SNPs = 5
Reading phenotype data...
Actual number of people in phenofile = 500; using all of these
Linear model: ( height_residuals ) ~ mu + cov1 + cov2 + SNP_A1 + cov1*SNP_A1
You are running mmscore...
loaded InvSigma...
Reading genotype data... done.
loaded null data... estimated null model... formed regression object...
Analysis: 100.00%... done

Gives:

name A1 A2 Freq1 MAF Quality Rsq n Mean_predictor_allele beta_SNP_addA1 sebeta_SNP_addA1 beta_SNP_cov1 sebeta_SNP_cov1 chi2_SNP_add
rs735579 AC B 0.5847 0.5847 0.5847 0.5847 491 0.534623 1.73386 0.909754 0.296646 1.59171 3.63231
rs9088604 A CTG 0.5847 0.5847 0.5847 0.5847 497 0.435614 -2.07862 0.961285 1.87346 1.6741 4.67569
rs1413801 A B 0.5847 0.5847 0.5847 0.5847 492 0.311992 -2.18676 1.12553 1.20374 1.87225 3.77478
rs4911638 A B 0.5847 0.5847 0.5847 0.5847 496 0.150202 0.675667 1.30507 0.75037 2.25757 0.268039
rs70099 A B 0.5847 0.5847 0.5847 1.2e-17 493 0.0922921 NaN NaN NaN NaN NaN
@lckarssen lckarssen added the bug label Nov 4, 2015
@lckarssen lckarssen added this to the v0.5.0 milestone Nov 4, 2015
@lckarssen lckarssen changed the title Using --mmscore and --interaction will not output a column with cov_SNP_int_SNP_covariate Using --mmscore and --interaction will not output a column with cov_SNP_int_SNP_covariate Nov 4, 2015
@lckarssen lckarssen modified the milestones: v1.0.0, v0.5.0 Feb 8, 2016
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Projects
None yet
Development

No branches or pull requests

1 participant