Preguntas etiquetadas con r

¿Por qué el promedio del valor más alto de 100 sorteos de una distribución normal es diferente del percentil 98% de la distribución normal? Parece que, por definición, deberían ser lo mismo. Pero... Código en R: NSIM <- 10000 x <- rep(NA,NSIM) for (i in 1:NSIM) { x[i] <- max(rnorm(100)) …

8 r  distributions  maximum 

Ayudo a un colega a poner en marcha un modelo de meta-análisis de efectos mixtos utilizando el marco del paquete metafor R escrito por @Wolfgang. Interesante y preocupante, para uno de los coeficientes del modelo obtengo una distribución bimodal al arrancar (ver el panel inferior derecho de la figura a …

8 r  mixed-model  bootstrap  meta-analysis  bimodal 

La siguiente pregunta se basa en la discusión que se encuentra en esta página . Dada una variable de respuesta y, una variable explicativa continua xy un factor fac, es posible definir un Modelo Aditivo General (GAM) con una interacción entre xy facusando el argumento by=. De acuerdo con el …

8 r  interaction  gam  mgcv 

Me está costando llegar a la misma página que mi supervisor cuando se trata de validar mi modelo. He analizado los residuos (observados contra los valores ajustados) y utilicé esto como argumento para discutir los resultados obtenidos por mi modelo, sin embargo, mi supervisor insiste en que la única forma …

8 r  references  glmm  validation  glmmtmb 

He estimado algunas medidas repetidas de modelos de efectos fijos, con un componente de error anidado, basado en variables de agrupación, es decir, modelos no anidados, utilizando plm. Ahora estoy interesado en pruebe si los modelos completos son significativamente diferentes, es decir, donde es el modelo completo para y es …

8 r  hypothesis-testing  repeated-measures  nested-data  plm 

Tengo un conjunto de datos con 5.818.446 líneas y 51 columnas, de las cuales 50 son predictores. Mi respuesta es cuantitativa, así que estoy interesado en un modelo de regresión. Estoy tratando de ajustar un bosque aleatorio a mis datos usando el paquete caret. Sin embargo, no tengo suficiente RAM …

8 r  random-forest  large-data 

Estoy viendo cómo la distancia euclidiana mínima esperada entre puntos aleatorios uniformes y el origen cambia a medida que aumentamos la densidad de puntos aleatorios ( puntos por unidad de cuadrado ) alrededor del origen. He logrado llegar a una relación entre los dos descritos como tales: Expected Min Distance=12Density−−−−−−√Expected …

8 r  expected-value  monte-carlo  uniform  minimum 

Tengo un conjunto de datos con 338 predictores y 570 instancias (desafortunadamente no se puede cargar) en el que estoy usando el Lazo para realizar la selección de funciones. En particular, estoy usando la cv.glmnetfunción de la glmnetsiguiente manera, donde mydata_matrixhay una matriz binaria de 570 x 339 y la …

8 r  classification  lasso  glmnet 

Podemos simular una regresión lineal sin aleatoriedad, lo que significa que hacemos lugar de . Entonces, si ajustamos un modelo lineal, los coeficientes serán idénticos a la "verdad fundamental". Aquí hay un ejemplo.y= Xβy=Xβy=X\betay= Xβ+ ϵy=Xβ+ϵy=X\beta+\epsilon set.seed(0) n <- 1e5 p <- 3 X <- matrix(rnorm(n*p), ncol=p) beta <- runif(p) …

8 r  logistic  regularization  linear  separation 

Considere medidas repetidas ANOVA (RM-ANOVA) con un factor dentro del sujeto Ay varias mediciones por sujeto para cada nivel de A. Está estrechamente relacionado con ANOVA de dos vías con dos factores: Ay subject. Utilizan la descomposición idéntica de la suma de los cuadrados en cuatro partes: A, subject, A⋅subject, …

8 r  anova  repeated-measures 

Tengo un conjunto de datos de 1 dimensión y uso la boxplotfunción para hacer un diagrama de caja. Entonces puedo ver que tengo algunos valores atípicos. ¿Cuentan los valores atípicos cuando se determinan los cuantiles? ¿Hay una forma correcta / incorrecta o ambas son correctas siempre que tengamos claro qué …

8 r  data-visualization  boxplot 

Quería estimar el intervalo de confianza para la desviación estándar de algunos datos. El código R se ve así: library(boot) sd_boot <- function (x, ind) { res <- sd(x$ReadyChange[ind], na.rm = TRUE) return(res) } data_boot <- boot::boot(data, statistic = sd_boot, R = 10000) plot(data_boot) Y tengo la siguiente trama: Estoy …

8 r  confidence-interval  standard-deviation  bootstrap 

Realicé el lazo y luego dejé una validación cruzada cv<-cv.glmnet(df, df$Price, nfolds = 1500) Cuando trazo cv me sale lo siguiente: También noté que obtengo 2 lambdas diferentes: lambda.minylambda.1se ¿Cuál es la diferencia entre estas lambdas? ¿Qué puedo entender de la gráfica anterior en general (de qué se tratan estos …

8 r  cross-validation  interpretation  lasso 

Tengo datos de recolección a largo plazo y me gustaría probar si el número de animales recolectados está influenciado por los efectos del clima. Mi modelo se ve a continuación: glmer(SumOfCatch ~ I(pc.act.1^2) +I(pc.act.2^2) + I(pc.may.1^2) + I(pc.may.2^2) + SampSize + as.factor(samp.prog) + (1|year/month), control=glmerControl(optimizer="bobyqa", optCtrl=list(maxfun=1e9,npt=5)), family="poisson", data=a2) Explicación de …

8 r  poisson-distribution  heteroscedasticity  glmm 

Estoy tratando de determinar si las orugas que comen una dieta natural (monkeyflower) son más resistentes a los depredadores (hormigas) que las orugas que comen una dieta artificial (una mezcla de germen de trigo y vitaminas). Hice un estudio de prueba con un tamaño de muestra pequeño (20 orugas; 10 …

8 r  anova