La función R cv.glm (biblioteca: arranque) calcula el error estimado de predicción de validación cruzada K-fold para modelos lineales generalizados y devuelve delta. ¿Tiene sentido usar esta función para una regresión de lazo (biblioteca: glmnet) y, de ser así, ¿cómo se puede llevar a cabo? La...
El siguiente modelo logístico multinivel con una variable explicativa en el nivel 1 (nivel individual) y una variable explicativa en el nivel 2 (nivel grupal): π 0 j = γ 00 + γ 01 z j + u 0 j … ( 2 ) π 1 j = γ 10 + γ 11 z j + u 1 j … ( 3
¿Cómo calcular la media, la varianza, la mediana, la desviación estándar y el modo de distribución? Si genero números al azar que forman la distribución normal, he especificado la media como m=24.2desviación estándar como sd=2.2: > dist = rnorm(n=1000, m=24.2, sd=2.2) Entonces puedo hacer lo...
No estoy muy seguro de qué son exactamente VectorSource y VCorpus en el paquete 'tm'. La documentación no está clara sobre estos, ¿alguien puede hacerme entender en términos
Estoy usando lmer en R para verificar el efecto de la condición ( cond) en algún resultado. Aquí están algunos hechos de seguridad de datos, donde s es el identificador de objeto y a, by cson condiciones. library("tidyr") library("dplyr") set.seed(123) temp <- data.frame(s = paste0("S", 1:30),...
Tengo datos mensuales de 1993 a 2015 y me gustaría hacer pronósticos sobre estos datos. Utilicé el paquete tsoutliers para detectar los valores atípicos, pero no sé cómo continúo pronosticando con mi conjunto de datos. Este es mi
La arimaxfunción en el TSApaquete es, que yo sepa, el único Rpaquete que se ajustará a una función de transferencia para los modelos de intervención. Sin embargo, carece de una función de predicción que a veces es necesaria. ¿Es lo siguiente una solución alternativa para este problema,...
Problema Me gustaría hacer alguna inferencia en un sistema análogo a morir con un número desconocido de lados. El dado se tira varias veces, después de lo cual me gustaría inferir una distribución de probabilidad sobre un parámetro correspondiente al número de lados que tiene el dado,...
Estoy tratando de convertir mi columna de factor a variables ficticias: str(cards$pointsBin) # Factor w/ 5 levels ".lte100",".lte150",..: 3 2 3 1 4 4 2 2 4 4 ... labels <- model.matrix(~ pointsBin, data=cards) head(labels) # (Intercept) pointsBin.lte150 pointsBin.lte200 pointsBin.lte250...
Me han hecho una pregunta sobre un modelo lineal hecho con R lm: "¿La regresión utilizó mínimos cuadrados iterativos lineales o no lineales?" Busqué un poco y [creo que entiendo] la diferencia entre los dos, pero no pude encontrar ninguna evidencia del uso de R de mínimos cuadrados lineales...
Estoy tratando de agrupar diferentes conjuntos de datos usando algoritmos no supervisados (agrupamiento). El problema es que tengo muchas funciones (~ 500) y una pequeña cantidad de casos (200-300). Hasta ahora solía hacer solo problemas de clasificación para los cuales siempre había etiquetado...
Suponga que tiene una distribución binomial donde el previo del parámetro es uniforme. ¿Cómo puedo obtener la distribución posterior del
Estoy usando la adonis()función en el veganpaquete para determinar 1) si las especies hospedantes concurrentes varían en su comunidad microbiana en múltiples sitios, y 2) si los sitios son diferentes. He examinado todas las publicaciones sobre CV y SO, y no hay una respuesta clara sobre cómo...
Espero que alguien pueda ayudar con lo que creo que es una pregunta relativamente simple, y creo que sé la respuesta, pero sin confirmación se ha convertido en algo de lo que no puedo estar seguro. Tengo algunos datos de conteo como una variable de respuesta y quiero medir cómo cambia esa...
Si bien me gusta pensar que tengo una buena comprensión del concepto de información previa en el análisis estadístico bayesiano y la toma de decisiones, a menudo tengo problemas para comprender mi aplicación. Tengo en mente un par de situaciones que ejemplifican mis luchas, y siento que no se...
Usando R en algunos datos e intentando ver si mis datos son o no heteroscedásticos, he encontrado dos implementaciones de la prueba Breusch-Pagan, bptest (package lmtest) y ncvTest (package car). Sin embargo, estos producen resultados diferentes. ¿Cuál es la diferencia entre los dos? ¿Cuándo...
¿Cómo asigno más peso a las observaciones más recientes en R? Supongo que esto es una pregunta o deseo común, pero me resulta difícil descubrir exactamente cómo implementar esto. He intentado buscar mucho para esto, pero no puedo encontrar un buen ejemplo práctico. En mi ejemplo, tendría un gran...
He visto que hay varias preguntas sobre la interpretación y la construcción de juegos, lo que parece ilustrar la dificultad para los no estadísticos para tratar con ellos. Desafortunadamente, de ninguno de los hilos o tutoriales que leí, pude obtener una comprensión clara de cómo construir un...
Me sorprende que no pueda encontrar ningún artículo / conferencia sobre cómo uno puede incorporar Distribuciones de probabilidad de clase previa en clasificadores como Regresión logística o Bosque aleatorio. Entonces mi pregunta es: ¿Cómo se puede incorporar la distribución de probabilidad de...
Estoy usando la función randomForesten el randomForestpaquete de R para hacer una regresión. Sin embargo, cuando intento incluir un término de interacción en los siguientes códigos: library(MASS) library(randomForest) Boston_f <- within(Boston, factor(rad)) mdl <- randomForest(lstat ~ rad *...