similar to: MLG con Binomial Negativa

Estimados, Mi problema consiste en probar efectos de dos sistemas de siembra en la abundancia de insectos (variable de conteo). Tengo un diseño en bloques con el sistema de siembra, con muestreo semanales de insectos (10 semanas) en cada unidad experimental. Quiero aplicar un modelo lineal generalizado con distribución de Poisson, con medidas repetidas en el tiempo. Agradecería ayuda con la

Estimados, Necesito estudiar los efectos de dos factores sobre una variable de conteo con un rango estrecho de variación (es decir, probablemente sin distribución Normal). El experimento tiene más de 10 (fechas de muestreo) observaciones repetidas en el tiempo. En primer lugar, considero que debería aplicar un modelo lineal generalizado (¿Mixto?) usando Poisson como distribución de referencia.

Hola!! Soy novato en esto del R. Consulto por que ya estoy agotando los recursos para resolver un problema con tutoriales, consultas que se encuentran en google, etc. Si bien el problema parece ser simple no le encuentro la vuelta. Tengo un data.frame que armé importando un csv. Algunas de las columnas son mayoritariamente compuestas por valores números pero class character ya que al contener

Hola a todos, Necesito estimar dos parametros utilizando la función nlm; fit<-nlm(hood2par,c(x01[i],x02[j]),iterlim=300, catch=x[,c(3,4,5)],sp=.5) donde hood2par es una logística modificada. Pero en mi caso, la convergencia de nlm depende de los valores iniciales de dichos parámetros. Para buscar dichos valores iniciales de manera automática, genero dos vectores de valores iniciales

Hola, yo quiero construir un data.table donde una columna (Parametros) son caracteres y otra el resultado de la función information.gain, que devuelve un data.frame. El código que he usado es este, pero me da error PesosParam <- data.table(,.(Parametros, Peso:= information.gain(In.hospital_death~., ParamCol))) Es posible hacer lo que digo? o debo hacer una transformación del data.frame a

Just curious as to which font the logo on the site is using. Looks a lot like Kaapeli Bold, but when I zoom in closely there appear to be subtle differences (perfectly round hole in the a, slightly more slanted l). Thanks, MLG

_____ Von: Pascal Speck [mailto:p.speck@ewersbach.net] Gesendet: Montag, 5. September 2005 10:37 An: 'asterisk-users@lists.digium.org' Betreff: Timeout when Dialing - HELP When i try to do a call I get this message after a few seconds: I IND :TIMEOUT pid:1 mode:NT addr:51400102 port:2 --> l3id:10040 cause:16 dad:800759 oad:20 channel:1 port:2 --> lib: prim 34582

Hola, ¿qué tal? Voy a plantear una pregunta que no creo que sea enteramente relacionada con R. Incluso pudiera ser una pregunta que tuviera que ver más con Windows y su gestión de memoria. Tengo una función foo que toma un argumento numérico y devuelve otro número. Durante el proceso crea un objeto muy grande. La llamo muchas veces haciendo algo así como: resultados <- sapply( entradas, foo

Buenos días a todos! Quisiera saber si me pueden ayudar con un Script. Necesito hacer un ANOVA para medidas repetidas. Sé que el modelo es el *lme,* pero no se indicarle las medidas repetidas. El conjunto de datos tiene como variable respuesta volumen ingresado y como factores tiempo y tratamientos. El tiempo se midió a las 2,4,6 y 8 min para los distintos tratamientos, son medidas repetidas

Dear listmembers, I've adjusted a mixed model with glmmPQL: nulo<-glmmPQL(POS~1,random=~1|GRUPO, family=binomial, data=new) and I've reched the following results: Linear mixed-effects model fit by maximum likelihood Data: new AIC BIC logLik 53238.5 53260.74 -26616.25 Random effects: Formula: ~1 | GRUPO (Intercept) Residual StdDev: 0.3402137 0.9952645

Una distribución binomial negativa podría ser adecuada. De hecho, una extensión de la binomial negativa, como la distribución de Waring generalizada, termina convergiendo a una distribución binomial negativa. Puedes ver que la mejora en la bondad del ajuste de la binomial negativa con respecto a la Poisson es sustancial (te muestro también las salidas.). datos<-c(1280, 1262, 1290, 1321,

Gracias a todos por las recomendaciones. Ya me he puesto en contacto con el autor y le explicado donde esta el error. Es tan fácil como que han cambiado la estructura de los objetos tipos lmer desde la versión 2.15.3 con lo que su paquete dejo de funcionar. Por suerte la nueva estructura conserva por lo menos la parte necesaria para hacer sus funciones. Se localizar el fallo y como arreglarlo, lo

Hace unos meses os escribir para comunicaros que había un fallo en esta función. Como os prometí os comento la respuesta por si alguno está interesado en utilizar el paquete LMERconvenientsfucntions Dear Javier, The package has been updated and should work for you fine now. Note that function mcp.fnc does not return the fourth plot (dffits) anymore. We still have to figure out a way to compute

Estimados quisiera poder salvar a csv y a excell la base de datos titanic data("titanic", package = "prLogistic") saludos José

Hi. I need to draw a circle whit center (a,b) and radio r. So I use the R code below a<-1.975 # valore x del centro b<-1.215 # valores y del centro r<-1.46 # radio x1<-seq(a-r,a+r,by=0.01); #los valores de x yp<-sqrt(r^2-(x1-a)^2)+b; # los valores y a partir de la ra??z positiva yn<-(-1)*sqrt(r^2-(x1-a)^2)+b; # los valores y a partir de la ra??z negativa x<-c(x1,x1);

Gracias Miguel Angel, me queda todo mucho mas claro, es que con el tema legal no tengo mucho conocimiento y me invento (o me hacen inventar) algunas cosas, con tu explicación y con la de los anteriores me queda claro, es que alguna cosa que he escuchado por allí y que me daba un poco de miedo ya que alguno ha argumentado en alguna reunion para criticar el software libre y es el tema de la licencia

Estimado Luisfo Chiroque Hay aristas múltiples, eso es por la preparación de los datos, es que separe por componentes, no me refiero al término componente estadístico, sino a descomposición de algo en los componentes que lo construyen, algunos son comunes y otros no, cuándo son comunes hay una relación en la red. Voy a probar con simplify, cuándo lo leí yo lo descarté, lo releeré, ¿Qué opina de

Estimados miembros de la lista. Estoy usando 2 funciones del paquete vcd y obtengo resultados que parecen discrepar entre ambas funciones. > dput(peces) c(26L, 63L, 1L, 21L, 20L, 50L, 0L, 104L, 19L, 0L, 8L, 32L, 0L, 24L, 5L, 2L, 0L, 36L, 0L, 5L, 27L, 1L, 13L, 24L, 15L, 25L, 23L, 5L, 37L, 26L, 29L, 14L, 15L) La función Ord_plot me dice que la mejor distribución para los datos es binomial

building on ubuntu with wx2.6 gives only one error: an argument type has changed from long to wxBitmapType in a constructor of the wxIcon class. These changes were needed to fix the build. This brings up the question what is the status of the swig/classes/include files which were generated using the xml extract script which is no longer available in CVS. Should that script be rerun using wx2.6

Hello, The following piece of code works fine in R.2.10 (ubuntu): switch(distr,         normal    = {if (is.infinite(param["desv"]))                      n <- c(n,"La desv. estándar no puede ser Inf.")                    if (param["desv"]<0)                      n <- c(n,"La desv. estándar no puede ser <0.")                    },