
Hi there,

I perform a stratified cox and then I need the strata as a  numeric array "straft.ln"
ft.ln <- coxph(Surv(times,deaths)~ages+chemos+chemos:f1+chemos:f2+horms+horms:f1+horms:f2+grades+grades:f1+grades:f2+positives+positives:f1+positives:f2+sizes+sizes:f1+sizes:f2+strata(stra),data=ddd)

basehazzft.ln=basehaz(ft.ln,centered=FALSE)
H0ft.ln=c(basehazzft.ln[,1])
timeft.ln=c(basehazzft.ln[,2])
straft.ln=c(basehazzft.ln[,3])

I notice that basehazzft.ln[,3] is not the same with straft.ln with respect to the values i.e.:

(basehazzft.ln$stra[1])
[1] stra=1
134 Levels: stra=1 stra=10 stra=100 stra=101 stra=102 ... stra=99
> c(basehazzft.ln$stra[1])
[1] 1


so far so good but:

(basehazzft.ln$stra[285])
[1] stra=2
134 Levels: stra=1 stra=10 stra=100 stra=101 stra=102 ... stra=99
> c(basehazzft.ln$stra[285])
[1] 47

while the desired value is 2, I get a 47. What am I doing wrong? I tried the as.numeric function but I have the same problem..

I just want instead of the [ "stra=1" "stra=1" "stra=1" etc]  the [1,1,1 etc]

Thanx in advance for any answers!




      
___________________________________________________________ 
Χρησιμοποιείτε Yahoo!; 
Βαρεθήκατε τα ενοχλητικά μηνύματα (spam); Το Yahoo! Mail 
διαθέτει την καλύτερη δυνατή προστασία κατά των ενοχλητικών 
μηνυμάτων http://login.yahoo.com/config/mail?.intl=gr

	[[alternative HTML version deleted]]


