Displaying 1 result from an estimated 1 matches for "clmatch1".
2006 Jul 29
0
Help with clogit in survival - conditional logistic regression
...,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
ii <- order(clmatch1$age,clmatch1$case)
clmatch2 <- clmatch1[ii,]
clmatch2
clogit(case~scc+agescc+strata(age),data=clmatch2)
-----------------______...