


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]]

