Detailed description |
See also the forum post linked above.
From the forum post:
When our phenotype file includes four columns (ID, Age, Event, cov1) and the interaction parameter is specified as interaction==1, the model given is incorrect (Coxph model: ( Age , Event ) ~ mu + cov1 + SNP_A1 + mu*SNP_A1).
When the interaction parameter is specified as interaction==2, the model shown is now correct, Coxph model: ( Age , Event ) ~ mu + cov1 + SNP_A1 + cov1*SNP_A1, but the analysis doesn't run, giving the following error
error: Interaction parameter is out of range (interaction=2)
When we include a phenotype file with an additional covariate column (ID, Age, Event, cov1, cov2) and specify interaction==2, the model looks ok (with the exception of a covariate we didn't want to include), Coxph model: ( Age , Event ) ~ mu + cov1 + cov2 + SNP_A1 + cov1*SNP_A1, and the analysis seems to run, but it's not clear what analysis is being run to give the results; we can't replicate any of the results shown in R.
|
|