Preguntas etiquetadas con r

Quiero comparar 2 medias de muestra para devoluciones de existencias de 1 minuto. Supongo que están distribuidos en Laplace (ya marcados) y dividí los retornos en 2 grupos. ¿Cómo puedo verificar si son significativamente diferentes? Creo que no puedo tratarlos como una distribución Normal, porque a pesar de que son …

10 r  confidence-interval  laplace-distribution 

En la lmerfunción dentro lme4de Rhay una llamada para la construcción de una matriz de modelo de efectos aleatorios, , como se explica aquí , páginas 7 - 9.ZZZ El cálculo de implica productos de KhatriRao y / o Kronecker de dos matrices, y . ZZZJyoJyoJ_iXyoXyoX_i La matriz es un …

10 r  mixed-model  lme4-nlme  matrix 

Tengo los siguientes vectores X e Y simples: > X [1] 1.000 0.063 0.031 0.012 0.005 0.000 > Y [1] 1.000 1.000 1.000 0.961 0.884 0.000 > > plot(X,Y) Quiero hacer una regresión usando el registro de X. Para evitar obtener el registro (0), trato de poner +1 o +0.1 …

10 r  regression  lognormal 

Imagine una situación: tenemos registros históricos (20 años) de tres minas. ¿La presencia de plata aumenta la probabilidad de encontrar oro en el próximo año? ¿Cómo probar tal pregunta? Aquí hay datos de ejemplo: mine_A <- c("silver","rock","gold","gold","gold","gold","gold", "rock","rock","rock","rock","silver","rock","rock", "rock","rock","rock","silver","rock","rock") mine_B <- c("rock","rock","rock","rock","silver","rock","rock", "silver","gold","gold","gold","gold","gold","rock", "silver","rock","rock","rock","rock","rock") mine_C <- c("rock","rock","silver","rock","rock","rock","rock", "rock","silver","rock","rock","rock","rock","silver", "gold","gold","gold","gold","gold","gold") time …

10 r  time-series  hypothesis-testing  stochastic-processes 

Dadas tres variables, yy x, que son positivas continuas, y z, que es categórica, tengo dos modelos candidatos dados por: fit.me <- lmer( y ~ 1 + x + ( 1 + x | factor(z) ) ) y fit.fe <- lm( y ~ 1 + x ) Espero comparar estos …

10 r  hypothesis-testing  mixed-model  lme4-nlme 

Estoy buscando agrupar / fusionar nodos en un gráfico usando la agrupación de gráficos en 'r'. Aquí hay una variación asombrosamente juguetona de mi problema. Hay dos "grupos" Hay un "puente" que conecta los grupos Aquí hay una red de candidatos: Cuando miro la distancia de conexión, el "conteo", si …

10 r  clustering  data-visualization  numerics 

Tengo algunas variables y estoy interesado en encontrar relaciones no lineales entre ellas. Así que decidí ajustar algunas spline o loess, e imprimir bonitos gráficos (ver el código a continuación). Pero también quiero tener algunas estadísticas que me den una idea de la probabilidad de que la relación sea una …

10 r  regression  splines  loess 

Hice un arranque con un modelo mixto (varias variables con interacción y una variable aleatoria). Obtuve este resultado (solo parcial): > boot_out ORDINARY NONPARAMETRIC BOOTSTRAP Call: boot(data = a001a1, statistic = bootReg, R = 1000) Bootstrap Statistics : original bias std. error t1* 4.887383e+01 -1.677061e+00 4.362948e-01 t2* 3.066825e+01 1.264024e+00 5.328387e-01 …

10 r  confidence-interval  bootstrap  bias 

Estoy ejecutando una regresión OLS agrupada usando el paquete plm en R. Sin embargo, mi pregunta es más sobre estadísticas básicas, por lo que intento publicarlo aquí primero;) Dado que mis resultados de regresión producen residuos heteroscedasticos, me gustaría intentar usar errores estándar robustos de heteroscedasticidad. Como resultado coeftest(mod, vcov.=vcovHC(mod, …

10 r  anova  multiple-regression  heteroscedasticity  robust-standard-error 

Estoy ejecutando glms en R (modelos lineales generalizados). Pensé que conocía los valores, hasta que vi que llamar a un resumen para un glm no le da un valor de referencia superior del modelo en su conjunto, al menos no en el lugar donde lo hacen los modelos lineales. Me …

10 r  statistical-significance  generalized-linear-model  p-value  descriptive-statistics 

Estoy jugando con el conjunto de datos de cáncer de mama y creé un diagrama de dispersión de todos los atributos para tener una idea de cuáles tienen el mayor efecto en la predicción de la clase malignant(azul) de benign(rojo). Entiendo que la fila representa el eje xy la columna …

10 r  data-visualization  interpretation  scatterplot 

Toda la literatura de Modelado de Distribución de Especies sugiere que cuando se predice la presencia / ausencia de una especie utilizando un modelo que genera probabilidades (por ejemplo, RandomForests), la elección del umbral de probabilidad para clasificar una especie como presencia o ausencia es importante y uno debe No …

10 r  classification  random-forest  threshold 

Tengo problemas para entender cómo varImpfunciona la función para un modelo randomForest con el caretpaquete. En el ejemplo a continuación, la característica var3 tiene cero importancia usando la varImpfunción de caret , pero el modelo final randomForest subyacente tiene una importancia distinta de cero para la característica var3. ¿Por qué …

10 r  caret  random-forest 

Estoy trabajando con algunos datos del mundo real y los modelos de regresión están dando algunos resultados contraintuitivos. Normalmente confío en las estadísticas, pero en realidad algunas de estas cosas no pueden ser ciertas. El principal problema que estoy viendo es que un aumento en una variable está causando un …

10 r  regression  regression-coefficients  computational-statistics 

Soy bastante nuevo en estadísticas y R. Me gustaría conocer el proceso para determinar los parámetros ARIMA para mi conjunto de datos. ¿Puedes ayudarme a descubrir lo mismo usando R y teóricamente (si es posible)? Los datos van del 12 de enero al 14 de marzo y muestran las ventas …

10 r  arima  box-jenkins