[R-es] Sobre survival analysis

[Francisco Viciana]

El scrip:

modelo<-coxp(surv(deaths,day)~size+strata(block))

resuelve un modelo de Cox, distinto para cada bloque, pero no suministra 
un efecto size genérico. Si introduces el" block" como efecto fijo, si 
tendrás un efecto "size" generico, pero con un  estimación de la 
varianza subestimada.

Si quieres introducir los bloques como efectos aleatorios en el modelo 
de cox, en R creo que tienes dos alternativas  bien con la  funcion 
"frailty" de "sulvival", con los cual el scrip quedaria algo asi

modelo<-coxp(surv(deaths,day)~size+frailty(block))

O bien, lo que parece mas recomendable usar el paquete 'coxme'  que 
permite plantear modelos multinivel con regresión de cox, que tu scrip 
quedaria algo asi:

require(coxme)

modelo<-coxme(surv(deaths,day)~size+(1|block))



El 28/06/2012 15:35, Flores-de-Gracia, Eric escribió:
> Hola
>
>
>
> Estoy tratando de correr un survival analysis usando un Cox regression model.
>
> Tengo una duda respecto a la organizacion del script. Tengo una variable que es -tamano del individuo- y quiero ver si hay diferencia en sobrevivencia respecto a tamano. Como diseno de campo los tamanos fueron ubicados de forma aleatoria en bloques al azar.
>
> Cuado planteo el script tengo algo como:
>
>
>
> modelo<-coxp(surv(deaths,day)~size+strata(block))
>
>
>
> El factor bloque es random pero no se si esta es la nomenclatura correcta para testar mi hipotesis, o simplemente debo omitir el termino -strata-
>
>
>
> Gracias,saludos
>
>
>
>
>
> Eric Flores De Gracia
> School of Biosciences
> University of Exeter, Cornwall Campus
> Penryn, Cornwall
> TR10 9EZ
> United Kingdom
> Mobilephone: +044 07578724705
> http://biosciences.exeter.ac.uk/cec/staff/postgradresearch/index.php?web_id=Eric_Flores_de_Gracia&tab=profile
>
> 	[[alternative HTML version deleted]]
>
> _______________________________________________
> R-help-es mailing list
> R-help-es en r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>


-- 
+--------------------------------------------------------------
| Francisco J. Viciana Fernández
| Coordinador del Registro de Población
| Servicio de Estadísticas Demográficas y Sociales
| Instituto de Estadística y Cartografía de Andalucía
| Leonardo Da Vinci, nº 21. Isla de La Cartuja.
| 41071 SEVILLA.
| franciscoj.viciana en juntadeandalucia.es
| Tlf.: +(34) 95 503 38 21
+--------------------------------------------------------------