similar to: Una pregunta de estadística (marginalmente relacionada con R)

Se ha borrado un adjunto en formato HTML... URL: <https://stat.ethz.ch/pipermail/r-help-es/attachments/20150128/9e23caec/attachment.html>

Hola, LLevo buscando desde hace tiempo como hacer el Goodness of fit test en R. Es decir, me explico, intento hacer una cosa parecida que se hace en Minitab, por ejemplo, yo tengo un conjunto de datos, y lo que quiero es sabes que tipo de distibución es, en minitab se hace un histograma para ver si se ajusta bien o no a la campana de Gauss, luego vemos si aproximar la distribución de la muestra

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

Estimados Escribo para consultar sobre el uso de modelos mixtos anidados. Los datos que estoy analizando provienen de censos de malezas en cuatro tipos de paisajes de la región pampeana, en los que seleccioné al azar igual número de lotes agrícolas cultivados con tres cultivos (maíz, soja y trigo-soja). En cada lote censé el número de especies de malezas en tres posiciones: el alambrado, el borde

Se ha borrado un adjunto en formato HTML... URL: <https://stat.ethz.ch/pipermail/r-help-es/attachments/20150127/8c94355b/attachment.html>

Buenas tardes, Estoy interesado en generar observaciones de una distribucion binomial bivariada en la que hay _cierto_ grado de correlacion (denotemoslo rho). Podria por favor alguien indicarme como hacerlo en R? Este es el contexto. Supongamos que se tienen dos experimentos en los que la variable respuesta _sigue_ una distribucion binomial, i.e., X_i ~Binomial(n_i, p_i), i=1,2 y que, por ahora,

Hola a tod en s, tengo una duda que la comunidad R me puede ayudar. Estoy trabajando con gamlss, porque tengo una variable respuesta con valores entre 0 y 1 e incluidos estos. La distribución que utilizo com gamlss para este caso es "beta inflated" (Stasinopulos and Rigby 2007. Journal of Statistical Software 23(7)). El modelo que intento correr es: m1<-gamlss(Teleosteos ~

Cordial saludo. Estoy trabajando con una densidad conjunta algo elaborada (suma y cociente de dos variables aleatorias positivas usando la cópula Frank) y requiero encontrar la distribución marginal de cada una de las variables. A manera de ejemplo (no es el caso específico, que tiene mucho más código), si partimos de la distribución exponencial bivariada #----- theta<-0.5

Hola, ¿Alguien sabe cómo lidia la función aov con diseño no balanceado? ¿Utiliza el las funciones de suma de cuadrados type II ? Gracias!! [[alternative HTML version deleted]]

Alguien tiene experiencia en escribir una pdMat. Para aquellos que no lo recuerden son las matrices de covarianzas de los efectos aleatorios que ajusta la función lme de la librería nlme Estas matrices tiene especial importancia en aplicaciones de genética de poblaciones y en particular en mapeo de asociación. Pinheiro y Bates dicen que el usuario puede crear sus propias pdMat y sugiere como

Lo he buscado y no encuentro la forma de aplicar una restricción lineal a los parámetros de una regresión. Podrían ayudarme?? Quiero una regresión en la que las suma de los coeficientes de igual a 1. Saludos -- Deybi Morales León

Se ha borrado un adjunto en formato HTML... URL: <https://stat.ethz.ch/pipermail/r-help-es/attachments/20120425/dafc9a68/attachment.html>

Buenos días para cada una/o. Acabo de suscribirme. Tengo una inquietud que vi tratada por partes (y usando distintos paquetes) en las consultas anteriores. ¿Cuál es la mejor forma de ajustar distribuciones datos a distribuciones de probabilidad teóricas? ¿Cuáles son los paquetes para hacerlo en forma ágil? Gracias por su atención. Saludos. César Escalante C.

Colegas: Necesito hacer una serie de ajustes de un modelo de decaimiento exponencial a unos datos de concentración de compuestos fluorescentes contra el tiempo. Para la mayoría de los experimentos tengo tres réplicas por cada tiempo, y haciendo los gráficos correspondientes parece haber diferencias muy grandes en la dispersión de los datos, siendo generalmente mas grandes al principio del

Hola compañeros de la lista. Tengo el siguiente set de datos: Repeticiones <- c(rep("RI", 14), rep("RII", 14), rep("RIII", 14)) Tiempo <- rep(c(0, 2, 4, 6, 8, 10, 12, 24, 36, 48, 60, 72, 96, 120), 3) Concentracion_celular <- c(0.4862, 0.5375, 0.4309, 0.4390, 0.4603, 0.4733, 0.3936, 0.9085, 0.5838, 0.5477, 0.6331, 0.8693, 1.0092, 0.6341, 0.5350,

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

Hola. Como hago para calcular la "Maxima Verosimilitud" de unos datos con distribucion uniforme? saludos. -- Leo [[alternative HTML version deleted]]

Buenas tardes, ¿Cómo interpretarías el intercepto que da R en un modelo de ceros inflados? Por un lado en la parte de conteo tengo un intercepto de -4.2 y en la parte de ceros de 102, ambos salen significativos (***). ¿Qué me dirían? Gracias

Hola a tod@s, Hace poco en la lista se discutieron algunas aproximaciones para determinar la distribución de probabilidad que potencialmente podrían haber generado los datos (ver [1]) y una de ellas fué el AIC. Haciendo uso del programa enviado por Pablo Verde (ver [1]), estos son los resultados para mis datos: weibull , AIC = 69839.44 exponential , AIC = 79488.77 gaussian , AIC = 69413.03

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