Realicé un experimento donde crié diferentes familias provenientes de dos poblaciones de origen diferentes. A cada familia se le asignó uno de dos tratamientos. Después del experimento medí varios rasgos en cada individuo. Para probar un efecto del tratamiento o la fuente, así como su interacción,...
Tengo un conjunto de datos de casos incidentes por temporada de una enfermedad rara. Por ejemplo, digamos que hubo 180 casos en la primavera, 90 en el verano, 45 en el otoño y 210 en el invierno. Me cuesta saber si es apropiado adjuntar errores estándar a estos números. Los objetivos de la...
Tengo algunos datos dicotómicos, solo variables binarias, y mi jefe me pidió que realizara un análisis factorial utilizando la matriz de correlaciones tetracóricas. Anteriormente, he podido enseñarme a mí mismo cómo ejecutar diferentes análisis basados en los ejemplos aquí y en el sitio de...
Considere un problema de análisis de datos clásico en el que tiene un resultado y cómo se relaciona con varios predictores . El tipo básico de aplicación en mente aquí es que X i 1 , . . . , X i pYiYiY_{i}Xi1,...,XipXi1,...,XipX_{i1}, ..., X_{ip} YiYiY_{i} es un resultado a nivel de grupo, como...
Tengo los siguientes datos y me gustaría ajustarle un modelo de crecimiento exponencial negativo: Days <- c( 1,5,12,16,22,27,36,43) Emissions <- c( 936.76, 1458.68, 1787.23, 1840.04, 1928.97, 1963.63, 1965.37, 1985.71) plot(Days, Emissions) fit <- nls(Emissions ~ a* (1-exp(-b*Days)),...
Estoy usando el paquete nnet en R para intentar construir un ANN para predecir los precios inmobiliarios de condominios (proyecto personal). Soy nuevo en esto y no tengo experiencia en matemáticas, así que por favor, descuida conmigo. Tengo variables de entrada que son binarias y continuas. Por...
Estaba leyendo sobre el cálculo de la estimación imparcial de la desviación estándar y la fuente que leí decía (...) excepto en algunas situaciones importantes, la tarea tiene poca relevancia para las aplicaciones de estadísticas, ya que su necesidad se evita mediante procedimientos estándar,...
Me preguntaba qué paquete estadístico de software recomiendan ustedes para realizar la Inferencia Bayesiana. Por ejemplo, sé que puede ejecutar openBUGS o winBUGS como independientes o también puede llamarlos desde R. Pero R también tiene varios de sus propios paquetes (MCMCPack, BACCO) que pueden...
Lectura mgcv::gamde la página de ayuda: los intervalos de confianza / credibilidad están disponibles para cualquier cantidad predicha usando un modelo ajustado Sin embargo, no puedo encontrar una manera de conseguir uno. Pensé predict.gamque tendría un type=confidencey un levelparámetro, pero...
Necesito estimar la función de densidad basada en un conjunto de observaciones usando el estimador de densidad del núcleo. Basado en el mismo conjunto de observaciones, también necesito estimar la primera y segunda derivadas de la densidad usando las derivadas del estimador de densidad del núcleo....
Estoy tratando de trazar un diagrama QQ con dos conjuntos de datos de aproximadamente 1,2 millones de puntos, en R (usando qqplot y alimentando los datos en ggplot2). El cálculo es bastante fácil, pero el gráfico resultante es muy lento de cargar, porque hay muchos puntos. He intentado una...
Tengo un montón de datos de dos muestras (control y tratamiento), cada una con varios miles de valores que se someterán a pruebas de significación en R. Teóricamente, los valores deberían ser continuos, pero debido al redondeo realizado por el software de medición no son t y tienen lazos. Las...
Recientemente implementé un filtro de Kalman en el ejemplo simple de medir la posición de una partícula con una velocidad y aceleración aleatorias. Descubrí que el filtro Kalman funcionaba bien, pero luego me pregunté cuál era la diferencia entre esto y simplemente hacer un promedio móvil. Descubrí...
Acabo de empezar a aprender a usar Stan y rstan. A menos que siempre haya estado confundido acerca de cómo funcionaba JAGS / BUGS, pensé que siempre tenía que definir una distribución previa de algún tipo para cada parámetro del modelo que se extraería. Sin embargo, parece que no tiene que hacer...
Tengo un vector de números que he subido aquí (... / code / MyData.Rdata) usando dput. Me gustaría obtener el bca ci, así que escribí este código: my.mean <- function(dat, idx){ return (mean(dat[idx], na.rm = TRUE)) } boot.out<-boot(data=my.data, statistic = my.mean, R=1000) Pero...
Utilicé randomForest para clasificar 6 comportamientos de animales (p. Ej., Estar de pie, caminar, nadar, etc.) en base a 8 variables (diferentes posturas corporales y movimiento). El MDSplot en el paquete randomForest me da este resultado y tengo problemas para interpretar el resultado. Hice un...
Estoy haciendo una validación cruzada usando el método de dejar uno afuera. Tengo una respuesta binaria y estoy usando el paquete de arranque para R y la función cv.glm . Mi problema es que no entiendo completamente la parte del "costo" en esta función. Por lo que puedo entender, esta es la función...
Cerrado. Esta pregunta está fuera de tema . Actualmente no está aceptando respuestas. ¿Quieres mejorar esta pregunta? Actualice la pregunta para que esté en el tema de Cross Validated. Cerrado hace 2 años . Tengo un conjunto de datos bastante simple que...
Actualmente estoy trabajando en un modelo de regresión donde solo tengo variables categóricas / factoriales como variables independientes. Mi variable dependiente es una relación transformada logit. Es bastante fácil ejecutar una regresión normal en R, ya que R sabe automáticamente cómo codificar...
Me han dado datos x = c(21,34,6,47,10,49,23,32,12,16,29,49,28,8,57,9,31,10,21,26,31,52,21,8,18,5,18,26,27,26,32,2,59,58,19,14,16,9,23,28,34,70,69,54,39,9,21,54,26) y = c(47,76,33,78,62,78,33,64,83,67,61,85,46,53,55,71,59,41,82,56,39,89,31,43,29,55,