Re: [R] prediction based on conditional logistic regression, clogit

2014-06-18 Thread peter dalgaard
On 17 Jun 2014, at 16:38 , Therneau, Terry M., Ph.D. wrote: > As Peter D said, the clogit function simply sets up a special data set and > then calls coxph, and is based on an identity that the likelihood for the > conditional logistic is identical to the likelihood of a Cox model for a > spe

Re: [R] prediction based on conditional logistic regression, clogit

2014-06-17 Thread Therneau, Terry M., Ph.D.
I ususally scan the digest for "surv", so missed your question on the first round. You caught predict with a case that I never thought of; I'll look into making it smarter. As Peter D said, the clogit function simply sets up a special data set and then calls coxph, and is based on an identity

Re: [R] prediction based on conditional logistic regression clogit

2014-06-16 Thread Charles Berry
peter dalgaard gmail.com> writes: > > > On 16 Jun 2014, at 05:22 , array chip yahoo.com> wrote: > > > Hi, I am using clogit() from survival package to do conditional > > logistic regression. I also need to make prediction on an > > independent dataset to calculate predicted probability. Here

Re: [R] prediction based on conditional logistic regression clogit

2014-06-16 Thread array chip
Thank you Peter. Any other suggestions are absolutely welcome!! John From: peter dalgaard Cc: "r-help@r-project.org" Sent: Monday, June 16, 2014 2:22 AM Subject: Re: [R] prediction based on conditional logistic regression clogit > H

Re: [R] prediction based on conditional logistic regression clogit

2014-06-16 Thread peter dalgaard
On 16 Jun 2014, at 05:22 , array chip wrote: > Hi, I am using clogit() from survival package to do conditional logistic > regression. I also need to make prediction on an independent dataset to > calculate predicted probability. Here is an example: > > >> dat <- data.frame(set=rep(1:50,each=

[R] prediction based on conditional logistic regression clogit

2014-06-15 Thread array chip
Hi, I am using clogit() from survival package to do conditional logistic regression. I also need to make prediction on an independent dataset to calculate predicted probability. Here is an example: > dat <- data.frame(set=rep(1:50,each=3), status=rep(c(1,0,0),50), > x1=rnorm(150,5,1), x2=rnorm