Preguntas etiquetadas con r

Tengo un problema al usar el 2l.normmétodo de imputación multinivel en mice. Lamentablemente, no puedo publicar un ejemplo reproducible debido al tamaño de mis datos; cuando reduzco el tamaño, el problema desaparece. Para una variable particular, miceproduce los siguientes errores y advertencias: Error in chol.default(inv.sigma2[class] * X.SS[[class]] + inv.psi) : …

11 r  missing-data  multiple-imputation  mice 

Me gustaría descomponer los siguientes datos de series temporales en componentes de temporada, tendencia y residuales. Los datos son un Perfil de Energía de Enfriamiento por hora de un edificio comercial: TotalCoolingForDecompose.ts <- ts(TotalCoolingForDecompose, start=c(2012,3,18), freq=8765.81) plot(TotalCoolingForDecompose.ts) Por lo tanto, existen obvios efectos estacionales diarios y semanales basados ​​en los …

10 r  time-series  forecasting  multiple-seasonalities  tbats 

tl; dr: comenzando con un conjunto de datos generado bajo nulo, volví a muestrear casos con reemplazo y realicé una prueba de hipótesis en cada conjunto de datos remuestreado. Estas pruebas de hipótesis rechazan el nulo más del 5% del tiempo. En la siguiente simulación, muy simple, genero conjuntos de …

10 r  bootstrap  simulation  resampling 

¿Por qué obtengo diferentes predicciones para la expansión polinómica manual y uso de la polyfunción R ? set.seed(0) x <- rnorm(10) y <- runif(10) plot(x,y,ylim=c(-0.5,1.5)) grid() # xp is a grid variable for ploting xp <- seq(-3,3,by=0.01) x_exp <- data.frame(f1=x,f2=x^2) fit <- lm(y~.-1,data=x_exp) xp_exp <- data.frame(f1=xp,f2=xp^2) yp <- predict(fit,xp_exp) lines(xp,yp) …

10 r  regression  polynomial 

Me pregunto cómo abordar adecuadamente el entrenamiento y la prueba de un modelo LASSO usando glmnet en R. Específicamente, me pregunto cómo hacerlo si la falta de un conjunto de datos de prueba externo requiere el uso de validación cruzada (u otro enfoque similar) para probar mi modelo LASSO. Déjame …

10 r  cross-validation  glmnet  lasso 

Como dice el título, estoy tratando de replicar los resultados de glmnet linear usando el optimizador LBFGS de la biblioteca lbfgs. Este optimizador nos permite agregar un término regularizador L1 sin tener que preocuparnos por la diferenciabilidad, siempre que nuestra función objetivo (sin el término regularizador L1) sea convexa. El …

10 r  regression  regularization  glmnet  elastic-net 

Entonces, antes que nada, investigué un poco en este foro, y sé que se han hecho preguntas extremadamente similares, pero por lo general no se han respondido correctamente o, a veces, las respuestas simplemente no son lo suficientemente detalladas como para que yo las entienda. Entonces esta vez mi pregunta …

10 r  regression  statistical-significance  regression-coefficients  group-differences 

Es probable que haya más de un malentendido serio en esta pregunta, pero no está destinado a hacer los cálculos correctos, sino a motivar el aprendizaje de series de tiempo con un cierto enfoque en mente. Al tratar de comprender la aplicación de series de tiempo, parece que la tendencia …

10 r  time-series  data-visualization 

Me gustaría entender cómo puedo obtener el porcentaje de varianza de un conjunto de datos, no en el espacio de coordenadas proporcionado por PCA, sino en un conjunto ligeramente diferente de vectores (rotados). set.seed(1234) xx <- rnorm(1000) yy <- xx * 0.5 + rnorm(1000, sd = 0.6) vecs <- cbind(xx, …

10 r  variance  pca  linear-algebra 

Tengo un conjunto de datos de 482 observaciones. data=Populationfull Voy a hacer un análisis de asociación de genotipo para 3 SNP. Estoy tratando de construir un modelo para mi análisis y estoy usando el aov (y ~ x, data = ...). Para un rasgo tengo varios efectos fijos y covariables …

10 r  anova 

¿Cómo llamarías a este tipo de trama, y ​​es posible crearlos en R? EDITAR: muchas gracias a todos, muy útil. Mejor título hasta ahora: ¡tramas de violín cuantizadas!

10 r  data-visualization 

Investigo sobre juegos educativos, y algunos de mis proyectos actuales implican el uso de datos de BoardGameGeek (BGG) y VideoGameGeek (VGG) para examinar las relaciones entre los elementos de diseño de los juegos (es decir, "ambientado en la Segunda Guerra Mundial", "implica lanzar dados" ) y calificaciones de jugadores de …

10 r  categorical-data  factor-analysis  binary-data 

¿Cómo puedo obtener p,d and qvalores en el ARIMA(p,d,q)modelo estimado por auto.arima(mytimeseries)? arima_model <- auto.arima (mytimeseries, ic = 'bic') Si nos fijamos en la salida de arima_model $ arma obtenemos, [1] 1 0 0 0 1 2 0 ¿Cuál es el significado de los números que aparecen en la secuencia …

10 r  arima 

Los problemas de clasificación con límites no lineales no se pueden resolver con un simple perceptrón . El siguiente código R tiene fines ilustrativos y se basa en este ejemplo en Python): nonlin <- function(x, deriv = F) { if (deriv) x*(1-x) else 1/(1+exp(-x)) } X <- matrix(c(-3,1, -2,1, -1,1, …

10 r  neural-networks  svm  kernel-trick  perceptron 

¿Existe un paquete R (o incluso una función R base) que implemente el método de Fisher o Stouffer para combinar valores p? Codificar esto debería ser casi trivial, pero prefiero usar (y citar) un paquete. Código de ejemplo en esta pregunta: Método de Fisher para peinar los valores p: ¿qué …

10 r  combining-p-values