[R-es] duda sobre error en modelo de hurdle

[Marcuzzi, Javier Rubén]

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