Displaying 1 result from an estimated 1 matches for "agescc".
Did you mean:
ageocc
2006 Jul 29
0
Help with clogit in survival - conditional logistic regression
...am not able to run clogit
on my dataset. I appreciate any input.
Regards,
Ashraf
age <- c(rep(0,8),rep(1,8))
scc <- c(1,1,1,0,0,1,0,0,1,1,1,0,0,1,0,0)
case <- rep(c(1,0),8)
count <- c(5,5,216,216,110,110,40,40,5,5,308,308,212,212,21,21)
#count <- c(2,5,3,5,4,6,2,4,2,5,1,3,1,4,2,5)
agescc <- age*scc
# Grouped data
clmatch <- data.frame(age,case,scc,agescc,count)
# Ungrouping the data as clogit as I can not use weights=count
clmatch1 <- data.frame(lapply(clmatch, function(x)
rep(x,clmatch$count)))[,1:4]
# Sorting by age and case for clogit - I think it is not necessary
i...