[R-es] Una pregunta de estadística (marginalmente relacionada con R)

[Carlos J. Gil Bellosta]

Muchas gracias por la contestación (y bienvenido a la lista).

Pensé en utilizar AIC pero me da algo de miedo cuando los modelos no
están anidados (la exponencial lo está, en cierto modo, en la
weibull), pero, en general...

¿Has probado a experimentar con tu código y a probar distintos valores
del número de datos y de distribuciones con los que los generas?

Un saludo,

Carlos J. Gil Bellosta
http://www.datanalytics.com

El día 30 de abril de 2009 16:23, Pablo Emilio Verde
<PabloEmilio.Verde en uni-duesseldorf.de> escribió:
> Este ejemplo te puede servir. Utilizo las distribuciones que estan la
> funcion survreg() del paquete survival y extraigo en AIC con la
> funcion extractAIC() del paquete MASS.
>
> #################################################
> set.seed(123)
> x <- rweibull(100, shape = 2, scale = 10)
>
> library(MASS) # para aplicar extractAIC
> library(survival) # para survreg
>
> distrib <-c("weibull", "exponential", "gaussian", "logistic",
>  "lognormal", "loglogistic")
>
> for( Dis in distrib){ tmpf <-extractAIC(survreg(Surv(x)~1, dist=Dis))
>                             cat(Dis, ", AIC = ",tmpf[2], "\n")
>               }
> #################################################
>
> El menor valor de AIC indica la distribucion de probabiliad que
> potencialmente genero los datos.
>
> Pablo