[R] Model Comparision for case control studies in R

anteneh asmare h@n@tezer@ @end|ng |rom gm@||@com
Wed Jun 15 19:09:53 CEST 2022

Dear Tim, Thanks. the first vector
y<-c(0,1,1,0,0,1,0,0,1,1,1,0,1,1,1,0,0,0,0,1) is the disease status y=
(1=Case,0=Control). The covariate age, smoking status and hypertension
are independent(uncorrelated). The logistic regression (unconditional)
will used. But I need to compare other models with logistic regression
instead of fitting it directly to logistic regression.
There is no matching on the data to use conditional logistics regression.
On 6/15/22, Ebert,Timothy Aaron <tebert using ufl.edu> wrote:
> Disease status is missing from the sample data.
> Are age, disease, smoking, and/or hypertension correlated in any way or are
> they independent (correlation=0)?
> Are the correlations large enough to adversely influence your model?
> Tim
> -----Original Message-----
> From: R-help <r-help-bounces using r-project.org> On Behalf Of anteneh asmare
> Sent: Wednesday, June 15, 2022 7:29 AM
> To: r-help using r-project.org
> Subject: [R] Model Comparision for case control studies in R
> [External Email]
> y<-c(0,1,1,0,0,1,0,0,1,1,1,0,1,1,1,0,0,0,0,1)
> age<-c(45,23,56,67,23,23,28,56,45,47,36,37,33,35,38,39,43,28,39,41)
> smoking<-c(0,1,1,1,0,0,0,0,0,1,1,0,0,1,0,1,1,1,0,1)
> hypertension<-c(1,1,0,1,0,1,0,1,1,0,1,1,1,1,1,1,0,0,1,0)
> data<-data.frame(y,age,smoking,hypertension)
> data
> model<-glm(y~age+factor(smoking)+factor(hypertension), data, family =
> binomial(link = "logit"),na.action = na.omit)
> summary(model)
> from above sample data I want to study a case-control study on male
> individuals with my response variable y, disease status (1=Case,
> 0=Control) with covariates age, smoking status(1=Yes, 0=No)  and
> hypertension, hypertensive (1=Yes, 0=No). I want to fit the model to predict
> the disease status using at least two different methods. And to make model
> comparisons. I think logistic regression will be the best fit for this case
> control study. Do we have other options in addition to logistic regression?
> My objective is to fit the model to predict the disease status using at
> least two different methods.
> Kind regards,
> Hana
> ______________________________________________
> R-help using r-project.org mailing list -- To UNSUBSCRIBE and more, see
> https://urldefense.proofpoint.com/v2/url?u=https-3A__stat.ethz.ch_mailman_listinfo_r-2Dhelp&d=DwICAg&c=sJ6xIWYx-zLMB3EPkvcnVg&r=9PEhQh2kVeAsRzsn7AkP-g&m=l7afPQ_gGAoV2EsNoYSYul0qAISEiXLmTmu0IQ03nZO4rcAi9xHZGsWwwig4oYOB&s=ztyDthknydhlcM49F33Gz6xRl6G7U9s8aIhB1VN-EKY&e=
> PLEASE do read the posting guide
> https://urldefense.proofpoint.com/v2/url?u=http-3A__www.R-2Dproject.org_posting-2Dguide.html&d=DwICAg&c=sJ6xIWYx-zLMB3EPkvcnVg&r=9PEhQh2kVeAsRzsn7AkP-g&m=l7afPQ_gGAoV2EsNoYSYul0qAISEiXLmTmu0IQ03nZO4rcAi9xHZGsWwwig4oYOB&s=tcsGkhvtVvoVvb1Ehah-vLRC6an40rJXQXqqfX2f0gI&e=
> and provide commented, minimal, self-contained, reproducible code.

More information about the R-help mailing list