Displaying 2 results from an estimated 2 matches for "m1_null".
Did you mean:
dt_null
2010 Dec 10
1
survreg vs. aftreg (eha) - the relationship between fitted coefficients?
...ata for simplicity (the principle should be the same with left truncation I hope)
library(survival)
library(eha)
# COMPARE coefs between survreg ('survival' pkg) and aftreg ('eha' pkg)
#Fitting NULL models (no covariates) results in (approximately) the same coefs (which is good!)
m1_NULL=survreg(Surv(futime/365, status==1) ~ 1, data=pbcseq)
m2_NULL=aftreg(Surv(futime/365, status==1) ~ 1, data=pbcseq)
c(m1_NULL$coef, 1/m1_NULL$scale) #--> intercept= 3.878656 , shape = 1.478177
c(m2_NULL$coef[1], exp(m2_NULL$coef[2])) #--> intercept= 3.878859 , shape=1.478150
# NOW I adju...
2010 Dec 11
0
is there a packge or code to generate markov chains in R
...ata for simplicity (the principle should be the same with left truncation I hope)
library(survival)
library(eha)
# COMPARE coefs between survreg ('survival' pkg) and aftreg ('eha' pkg)
#Fitting NULL models (no covariates) results in (approximately) the same coefs (which is good!)
m1_NULL=survreg(Surv(futime/365, status==1) ~ 1, data=pbcseq)
m2_NULL=aftreg(Surv(futime/365, status==1) ~ 1, data=pbcseq)
c(m1_NULL$coef, 1/m1_NULL$scale) #--> intercept= 3.878656 , shape = 1.478177
c(m2_NULL$coef[1], exp(m2_NULL$coef[2])) #--> intercept= 3.878859 , shape=1.478150
# NOW I adju...