Preguntas etiquetadas con r

En una pregunta en otra parte de este sitio, varias respuestas mencionaron que el AIC es equivalente a la validación cruzada de dejar uno fuera (LOO) y que el BIC es equivalente a la validación cruzada K-fold. ¿Hay alguna manera de demostrar esto empíricamente en R de modo que las …

26 r  aic  cross-validation  bic 

Proporcione un código R que le permita a uno realizar un ANOVA entre sujetos con -3, -1, 1, 3 contrastes. Entiendo que hay un debate con respecto al tipo apropiado de Suma de cuadrados (SS) para dicho análisis. Sin embargo, como el tipo predeterminado de SS utilizado en SAS y …

26 r  anova  contrasts  sums-of-squares 

Soy nuevo en R, ordené la regresión logística y polr. La sección "Ejemplos" en la parte inferior de la página de ayuda para polr (que ajusta un modelo de regresión logística o probit a una respuesta de factor ordenado) muestra options(contrasts = c("contr.treatment", "contr.poly")) house.plr <- polr(Sat ~ Infl + …

26 r  logistic 

Estoy ejecutando modelos de regresión LOESS en R, y quiero comparar los resultados de 12 modelos diferentes con diferentes tamaños de muestra. Puedo describir los modelos reales con más detalles si me ayuda a responder la pregunta. Aquí están los tamaños de muestra: Fastballs vs RHH 2008-09: 2002 Fastballs vs …

26 r  regression  loess 

Tengo una matriz de correlación de retornos de seguridad cuyo determinante es cero. (Esto es un poco sorprendente ya que la matriz de correlación de la muestra y la matriz de covarianza correspondiente deberían ser teóricamente definitivas positivas). Mi hipótesis es que al menos un valor depende linealmente de otros …

26 r  correlation  pca  linear-model  svd 

¿Es posible encontrar el valor p en la correlación de pearson en R? Para encontrar la correlación de Pearson, generalmente hago esto col1 = c(1,2,3,4) col2 = c(1,4,3,5) cor(col1,col2) # [1] 0.8315218 Pero, ¿cómo puedo encontrar el valor p de esto?

26 r  correlation  p-value  pearson-r 

Actualmente estoy ejecutando algunos modelos lineales de efectos mixtos. Estoy usando el paquete "lme4" en R. Mis modelos toman la forma: model <- lmer(response ~ predictor1 + predictor2 + (1 | random effect)) Antes de ejecutar mis modelos, verifiqué la posible multicolinealidad entre los predictores. Hice esto por: Hacer un …

26 r  correlation  mixed-model  lme4-nlme  multicollinearity 

He explorado muchos sitios de ayuda y todavía estoy confundido acerca de cómo especificar términos anidados más complicados en un modelo mixto también. También estoy confundido en cuanto al uso de :y /y |en la especificación de interacciones y anidando con factores aleatorios usando lmer()en el lme4paquete de R. A …

26 r  mixed-model  lme4-nlme 

Esto es similar a los métodos de muestreo de Caret , aunque eso realmente nunca respondió a esta parte de la pregunta de una manera acordada. la función de tren de caret ofrece cvy repeatedcv. ¿Cuál es la diferencia en decir haciendo: MyTrainControl=trainControl( method = "cv", number=5, repeats=5 ) vs …

26 r  machine-learning  caret 

Inspirado por la charla de Peter Donnelly en TED , en la que analiza cuánto tiempo le tomaría a un cierto patrón aparecer en una serie de lanzamientos de monedas, creé el siguiente guión en R. Dados dos patrones 'hth' y 'htt', calcula cuánto tiempo tarda (es decir, cuántos lanzamientos …

26 r  probability  stochastic-processes 

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 6 años . Me gustaría realizar una normalización en columna de una matriz en R. Dada una matriz m, …

26 r  data-transformation  normalization  matrix 

Anteriormente utilicé Forecast Pro para pronosticar series de tiempo univariantes, pero estoy cambiando mi flujo de trabajo a R. El paquete de pronósticos para R contiene muchas funciones útiles, pero una cosa que no hace es ningún tipo de transformación de datos antes de ejecutar Auto .arima (). En algunos …

26 r  time-series  data-transformation  forecasting  arima 

Considere un modelo de obstáculo que predice datos de conteo yde un predictor normal x: set.seed(1839) # simulate poisson with many zeros x <- rnorm(100) e <- rnorm(100) y <- rpois(100, exp(-1.5 + x + e)) # how many zeroes? table(y == 0) FALSE TRUE 31 69 En este caso, …

25 r  count-data  zero-inflation 

La siguiente cita de bloque, de los líderes en el campo del modelado de efectos mixtos, afirma que los cambios de coordenadas en modelos con correlación cero entre efectos aleatorios (modelos 'ZCP') cambian las predicciones del modelo. Pero, ¿alguien puede elaborar o justificar aún más sus afirmaciones? Las declaraciones en …

25 r  mixed-model  lme4-nlme 

Me gustaría mostrar cómo los valores de ciertas variables (~ 15) cambian con el tiempo, pero también me gustaría mostrar cómo las variables difieren entre sí en cada año. Entonces creé esta trama: Pero incluso al cambiar el esquema de color o agregar diferentes tipos de línea / forma, esto …

25 r  data-visualization