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