Estimado Matias Ledesma
Posiblemente el modelo esta mal escrito, no porque el razonamiento fuese
incorrecto desde el punto de vista científico (desconozco las
particularidades de su caso), pero me paso al realizar un cbind, y luego
ingresar este al modelo, no todos los paquetes o librerías trabajan de la
misma forma, y al utilizar esa alternativa hay inconvenientes.
Lógicamente depende de como se pueden preparar los datos, yo lo solucione de
la siguiente forma.
Supongamos que el modelo es cabeza + pie = altura, peso ....
Lo que usted realizó es cdind (cabeza, pie)
Lo que habría que buscar es parte = código, algura, peso
donde código es 1 para cabeza y dos para pie
Pero, claro, es medio complicado, hay que pensar mucho y no en todos los
casos es posible, desconozco si para su tipo de análisis esto es adecuado o
técnicamente posible, porque hay que preparar los datos, recodificarlos,
aunque, posiblemente hay ejemplos adecuados para su trabajo, alguna
documentación que facilite su labor.
Javier Marcuzzi
-----Mensaje original-----
From: Matias Ledesma
Sent: Friday, April 26, 2013 8:12 AM
To: R
Subject: [R-es] duda sobre error en modelo de hurdle
Hola a todos,
Tengo una
pregunta sobre el modelo de Hurdle, estoy tratando de analizar una
variable la cual presenta las siguientes caracteristicas: gran variacion,
muchos ceros , con asimetria
positive y es una fraction (embryones malformados/embryones
sanos + embryones malformados por individuo).
La data
proviene del monitoreo anual que se
realiza en el Mar Baltico, donde se utiliza un amphipodo (monoporeia
affinis)
como bioindicador de contaminacion entre otras cosas. Cada hembra capturada
tiene un marsupio con embryones (entre unos pocos a 30 o 40 embryones).
El muestreo
se realiza en estaciones de referencia por lo cual los embryones malformados
son un porcentaje muy bajo del total , ademas se cree que ahi un factor
genetico lo cual hace que los
embryones no se puedan tomar como independientes.
Eh
intentado distintos caminos y actualmente estoy probando el modelo de
hurdle,
probando como varia las malformaciones con las estaciones y años.
dat<-read.delim("HurdleR1213.txt")
malform<-cbind(dat$Tot.dameggs,dat$FECUND-dat$Tot.dameggs)
library(pscl)
hurdle(malform~dat$STATN*as.factor(dat$YEAR),dist=c("negbin"),zero.dist=c("binomial"))
El problema
que tengo es que me sale el siguiente error y no se como solucionarlo,
alguien
tiene alguna idea?
Error in
glm.fit(X, Y, family = poisson(), weights = weights, offset = offsetx) :
(subscript) logical subscript too long
Saludos,
/Matias
[[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