similar to: modelos mixtos con familia quasibinomial

Hola a todos, Me gustaria saber si en un modelo mixto se puede usar el tamaño de muestra como peso o es incorrecto hacerlo. Por ejemplo en el comando 'lmer' del paquete lme4 hay una opcion 'weights' (igual que en 'lm' de stats). Si tengo datos de una medida (por ejemplo el peso) de 10 especies de aves (rango, de 50 a 3000 datos segun la especie) ¿puedo usar esta n en

Hola a todos. Estoy utilizando la función lmer del paquete lme4 para ajustar un modelo mixto. Tengo varias variables en mi data.frame, unas son a nivel individual y otras a nivel de comarcas. Listo algunas. ingre_6 : Ingresos (nivel individual) iscs_a : un indicador sintético resumen de otras variables, calculado mediante componentes principales. sau_com : superficie agraria útil de cada una

Hola a todos, Estoy trabajando con modelos lineales generalizados (GLM), en particular con las funciones glm y anova.glm de R. Tengo una pregunta que es más bien técnica en el sentido estadístico (y un poco Off Topic y probablemente un tanto naive). No tengo claro si es correcto decir análisis de varianza o de devianza al utilizar GLMs, veo que en el help de anova.glm dice análisis de

Is this a bug? seasonal.dummies <- function(x, contr=NULL) { # takes a time series and returns a matrix of seasonal dummies for # x. This is almost cycle(x), we only have to make it into a factor # and add suitable level names. # return a matrix which includes a constant! # level names here assumes frequency is 12! cyc <- factor( cycle(x), labels=c("ene","feb",

Dispongo de un experimento en el que cinco tratamientos ha sido aplicados a cinco grupos de voluntarios. En cada grupo había tres personas y a cada persona se le tomaron 3 medidas, En total dispongo de 45 medidas, pero evidentemente no son independientes entre sí. Si no tomo en cuenta que las medidas de la misma persona son más parecidas entre sí (bloques anidados) estaría incurriendo en

Buenas a todos. Desde hace un tiempo estoy tratando de aplicar las funciones de la familia *pply en todo lo que puedo, pero todavía no es algo que me surja tan rápidamente o naturalmente al momento de los loops como usar for(). Conozco las ventajas de usar estas funciones y por eso mi intento de hacerme de ellas. Por ejemplo en este problema: data=matrix(rnorm(100*20),20,100)

Hola Jorge, muchas gracias por tu pronta respuesta, no me di cuenta que el formateo podría causar problemas, envío de nuevo el código sin formatos. La idea básica es para un set de números de columnas (desordenados) y un set de numeros de fila el loop lo que hace es ir a la fila y columna correspondiente de data, tomar el valor y luego hacer la media sobre esos.

Hi, How do I build a table from a regression model adjusted using xtable? Commands are: modelo1 = lm(Y~X1 + X2) influencia = influence.measures(modelo1) require(xtable) xtable(influencia) but it isn't work. Thanks, -------------------------------------- Silvano Cesar da Costa Departamento de Estat?stica Universidade Estadual de Londrina Fone: 3371-4346

Javier, es necesario para así saber el valor de la masa seca, la planta no se va a ver muy afectada en terminos de su fisiología, solo se toma una hoja por individuo, el problema aquí es calcular el área defoliada, si lo hago con una foto va a ser un poco subjetivo, porque hay hojas que están plegadas y eso puede afectar la estimación, la extrapolación a otras hojas perdería precisión. Este es

Muchas gracias a ambos Carlos y Jorge por las respuestas. Pido disculpas en la demora de respuesta, pero estuvo complicada la semana. La pregunta era un ejercicio de ejemplo para poder entender mejor los usos, creo que me armaré una guía en markdown con ejemplos varios para ir consultando cuando me salgan dudas de como usarlos. En realidad no importaba tanto si mejorara demasiado los tiempos

Hola a todos. ¿Cómo puedo ajustar un modelo de regresión logística en el que la varianza de los errores no son iguales, sino que esta puede ser modelada por otras variables ? ¿Utilizando modelos mixtos o voy desencaminado? Gracias..

An embedded and charset-unspecified text was scrubbed... Name: not available URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20091125/3f0e364b/attachment-0001.pl>

Muchas gracias por sus respuestas Manuel y Javier. El problema al que me enfrento es que a veces estas afectaciones por herbivoría dejan cicatrices que pueden deformar la hoja, a veces se pliegan sobre sí mismas y otras veces se pierde bastante la configuración de la hoja. Por el momento solo me interesa tratar de reconstruir la hoja para poder obtener el área que tendría esa hoja si no se hubiese

Estimado Quizás, por lo que dice Carlos Ortega, antes de pasarlo a escala de grises, podría pasarlo a escala RGB, jugar un poco, el entrenamiento no distinguirá si es escala de gris, como mezcla de negro y blanco, puede intentar rojo y blanco, rojo y verde, con la idea de detectar cambios en el color relacionado a posibles cambios en la humedad o ?madurez? de las hojas, por ahí, con muy poco

Muchas gracias por todos sus comentarios. Muy valiosa la información, hay algunas cosas que no entiendo muy bien. Si tuvieran algunos documentos que pudiera leer les agradecería. Yo he trabajado con las relaciones entre largo y ancho de la hoja, esas relaciones son muy precisas, sin embargo, el problema está cuando las hojas son afectadas, se pierde bastante la estructura. Al yo obtener el área

Estimado Jimmy Nooo ¿como la cortas? Entiendo que no sirve más. Creo que el fondo blanco te perjudica. Mi experiencia, realicé una reconstrucción 3d, foto, foto con fotógrafo que eliminó el resto, filmación pasada a fotos. Lo peor, la foto con fotógrafo, la de mejor calidad. Quizás deberías estudiar todo nuevamente, tomar los errores, buscar soluciónes y rehacer el trabajo. Javier Marcuzzi

Hola, ¿qué tal? Voy a ejercer de portavoz oficioso del "petit comité" que inició esta lista de correo y pedir un favor a los actuales miembros de la lista. Como más de uno ya sabe, crearla, aparte del interés "per se", tenía un objetivo adicional y de mayor alcance: reunir a la comunidad hispanohablante de usuarios de R y, en particular, a la española para intercambiar

Hi, I have tried to run the following (I know it's a huge data set but I tried to perform it with a 1 GB RAM computer): library(foreign) library(MASS) library(nlme) datos<-read.spss(file="c:\\Documents and Settings\\Administrador\\Escritorio\\datosfin.sav",to.data.frame=TRUE) str(datos) `data.frame': 1414 obs. of 5 variables: $ POB : Factor w/ 6 levels

Hola, Si, lo eh pensado pero no eh tenido tiempo de ondar y tengo que profundisar en la biologia y en el procediminto del monitoreo. Dado que todos los años se muestera la misma estacion talvez habria que utilizar un glmm en vez de glm... En realidad lo que estoy aplicando es el modelo de Hurdle, pero haciendolo en dos partes separadas ya que primero analizo la possiblidad que hay de encontrar

Gracias! El 15 de junio de 2015, 16:54, Freddy Omar López Quintero < freddy.vate01 en gmail.com> escribió: > ?Holap.? > > ran out of iterations and failed to converge > > > ?Prueba aumentando el número de iteraciones, con el argumento maxit: > > ?GLM <- bigglm(In.hospital_death ~ GCS + BUN, data = DatosGLM, family = >> binomial(logit), maxit=1000)? >