-
Notifications
You must be signed in to change notification settings - Fork 12
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
gwasurvivr and R "survival" package give different results #7
Comments
Please send a reproducible example and I can look at it with you. In theory they shouldn't be off by an order of magnitude unless the genotyping (or imputation) quality is poor. We are using the same |
Thanks for your reply. I am sorry the data is in the server I cannot copy it. Below are codes of two methods
bed.file <- "data_clean/Two/cov_id.bed" plinkCoxSurv(bed.file=bed.file, I transferred bed file to 0/1/2 and it is contained in cov_match.txt used by cox.ph() Also, I checked the source code of gwasurvivr and I did find that it uses coxph.fit() rather than coxph. Is there any difference? I find few manual about coxph.fit() At last, you mentioned "The only difference is that in gwasurvivr we put in initial estimates that we obtain by fitting the model just the covariates first." When I use coxph() in R survial, is there anything I need to modify? Thanks for your help! |
Hi and thanks for your time! Today I checked the sample data within the gwasurvivr package. The results of At last, returning to my original problem, my data is case-only, that is, all subjects experienced the event. Is this a possible reason that two methods gave me different coefs? Thank you! |
Hi,
I find that given the same genotype data and phenotype data. gwasurvivr and R survival package cox.ph() give different p-value (both significant, but p-values are quite different, for example, one might be 1e-6 and the other 1e-7)and coefs.
Could you please give me any help? Is there any possible reason?
Thanks!
The text was updated successfully, but these errors were encountered: