Preguntas etiquetadas con r

Dado que el error estándar de una regresión lineal generalmente se da para la variable de respuesta, me pregunto cómo obtener intervalos de confianza en la otra dirección, por ejemplo, para una intersección x. Puedo visualizar lo que podría ser, pero estoy seguro de que debe haber una forma directa …

9 r  regression  confidence-interval  bootstrap 

¿Alguien podría decirme por qué obtengo resultados diferentes de Rmínimos cuadrados ponderados y solución manual por operación de matriz ? Específicamente, estoy tratando de resolver manualmente , donde es la matriz diagonal en los pesos, es la matriz de datos, es la respuesta vector. W A bW A x = …

9 r  regression  least-squares  weighted-regression  weighted-data 

En resumen, mi pregunta: ¿existen métodos para mejorar el tiempo de ejecución de R MICE (imputación de datos)? Estoy trabajando con un conjunto de datos (30 variables, 1.3 millones de filas) que contiene (bastante al azar) datos faltantes. Alrededor del 8% de las observaciones en aproximadamente 15 de las 30 …

9 r  multiple-imputation  mice 

¿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 conjunto …

9 r  time-series  regression  random-forest 

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 debería elegir usar uno u otro? …

9 r  regression  heteroscedasticity  assumptions  breusch-pagan 

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 …

9 r  binomial  proportion  glmm  lme4-nlme 

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 los datos como conjuntos de entrenamiento. Allí …

9 r  clustering  feature-selection  unsupervised-learning 

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 pointsBin.lte300 # 741 1 0 0 0 0 # 407 …

9 r  categorical-data  categorical-encoding 

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, θ. Intuición Si …

9 r  probability  bayesian  jags 

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), a = rnorm(30, -2, …

9 r  lme4-nlme  contrasts 

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 código: product.outlier<-tso(product,types=c("AO","LS","TC")) plot(product.outlier) Esta es mi salida del paquete tsoutliers ARIMA(0,1,0)(0,0,1)[12] Coefficients: …

9 r  time-series  forecasting  arima  outliers 

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 simples?

9 r  text-mining 

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 biblioteca glmnet utiliza una …

9 r  regression  cross-validation  lasso  glmnet 

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 …

9 r  logistic  generalized-linear-model  simulation  multilevel-analysis 

¿Cómo debo entender el anovaresultado al comparar dos modelos? Ejemplo: Res.Df RSS Df Sum of Sq F Pr(>F) 1 9 54.032 2 7 4.632 2 49.4 37.329 0.0001844 *** La página de manual dice: "Calcular las tablas de análisis de varianza (o desviación) para uno o más objetos de modelo …

9 r  regression  anova