Hola Lorena,
Creo que el error viene de que no has escrito bien la función
"LogLike()".
Es su lugar has escrito "logLik()" que no existe en el paquete que has
utilizado "spatcounts".
Si estás utilizando otro paquete para este test, ayuda mucho el que lo
indiques en la pregunta.
El que hayas incluido un ejemplo breve está muy bien, pero si incluso
añades un conjundo de datos con el que se pueda reproducir el error, pues
mucho mejor.
Saludos,
Carlos Ortega
www.qualityexcellence.es
El 27 de octubre de 2014, 18:58, Lorena Tudela Marco <
lorenatudelamarco en gmail.com> escribió:
> Hola a todxs,
>
> Os escribo para compartir una duda sobre Test de Clarke:
>
> *Clarketest(LogLike1, LogLike2, alpha = 0.05, p = NULL, q = NULL, *
> *correction = TRUE)*
>
> Estoy intentando comparar dos modelos (non nested models comparison) : Zero
> inflated Model vs Negative Bonimial Model. He aplicado el Test de Vuong
> (que se decanta por el NBM) y los criterios BIC y AIC (que apuestan por el
> ZIM). Por tener un estadístico mas, quería aplicar el test de Clarke:
>
> 1º Estimo los "loglik":
>
> NB<-logLik(NegativeB, df=TRUE, nobs=TRUE)
> ZI<-logLik(ZIM, df=TRUE, nobs=TRUE)
>
> 2º Aplico el test:
>
> Clarketest(NB, ZI, alpha = 0.05, p = NULL, q = NULL, correction = TRUE)
>
> Me da el siguiente error en la consola:
> Error in LogLike1$Coef : $ operator is invalid for atomic vectors
>
> ¿Sabéis porque me da ese error? ¿Estoy introduciendo mal el Loglike?
> Creo...que me falta algo por completar...
>
> Sin duda, cualquier aportación o sugerencia es bienvenida! muchas gracias
>
> y buen dia,
>
> Lorena
>
> [[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
>
--
Saludos,
Carlos Ortega
www.qualityexcellence.es
[[alternative HTML version deleted]]