Preguntas etiquetadas con r

Tengo una pregunta sobre la interpretación de la llamada tsboot en R. Revisé la documentación de Kendall y el paquete de arranque, pero no soy más inteligente que antes. Cuando ejecuto un bootstrap usando, por ejemplo, el ejemplo en el paquete Kendall, donde la estadística de prueba es la tau …

11 r  bootstrap 

Usé el siguiente código R para ajustar un modelo probit: p1 <- glm(natijeh ~ ., family=binomial(probit), data=data1) stepwise(p1, direction='backward/forward', criterion='BIC') ¿Quiero saber qué hace stepwisey qué hace backward/forwardexactamente y cómo seleccionar las variables?

11 r  probit  stepwise-regression 

Usando Amelia en R, obtuve múltiples conjuntos de datos imputados. Después de eso, realicé una prueba de medidas repetidas en SPSS. Ahora, quiero agrupar los resultados de las pruebas. Sé que puedo usar las reglas de Rubin (implementadas a través de cualquier paquete de imputación múltiple en R) para agrupar …

11 r  spss  p-value  multiple-imputation  pooling 

¿Existe una estrategia para elegir la cantidad de árboles en un GBM? Específicamente, el ntreesargumento en Rla gbmfunción de. No veo por qué no debe establecer ntreesel valor razonable más alto. Me di cuenta de que una mayor cantidad de árboles reduce claramente la variabilidad de los resultados de múltiples …

11 r  machine-learning  gbm 

Estoy tratando de usar el neuralnetpaquete de R (documentación aquí ) para la predicción. Aquí lo que estoy tratando de hacer: library(neuralnet) x <- cbind(runif(50, min=1, max=500), runif(50, min=1, max=500)) y <- x[, 1] * x[, 2] train <- data.frame(x, y) n <- names(train) f <- as.formula(paste('y ~', paste(n[!n %in% …

11 r  neural-networks  prediction  error 

Tengo un conjunto de datos de 20 años de un recuento anual de abundancia de especies para un conjunto de polígonos (~ 200 polígonos continuos de forma irregular). He estado utilizando análisis de regresión para inferir tendencias (cambio en el conteo por año) para cada polígono, así como agregaciones de …

11 r  time-series  correlation  spatial  spatio-temporal 

Considere el siguiente código y salida: par(mfrow=c(3,2)) # generate random data from weibull distribution x = rweibull(20, 8, 2) # Quantile-Quantile Plot for different distributions qqPlot(x, "log-normal") qqPlot(x, "normal") qqPlot(x, "exponential", DB = TRUE) qqPlot(x, "cauchy") qqPlot(x, "weibull") qqPlot(x, "logistic") Parece que el gráfico QQ para log-normal es casi el …

11 r  data-visualization  interpretation  qq-plot 

Para los datos de recuento que he recopilado, uso la regresión de Poisson para construir modelos. Hago esto usando la glmfunción en R, donde uso family = "poisson". Para evaluar posibles modelos (tengo varios predictores) utilizo el AIC. Hasta aquí todo bien. Ahora quiero realizar una validación cruzada. Ya logré …

11 r  generalized-linear-model  poisson-distribution 

Tengo un conjunto de datos que contiene la cantidad de acciones realizadas por individuos en el transcurso de 7 días. La acción específica no debería ser relevante para esta pregunta. Aquí hay algunas estadísticas descriptivas para el conjunto de datos: RangoMediaDiferenciaNumero de observaciones0 - 77218,22791696Range0−772Mean18.2Variance2791Number of observations696 \begin{array}{|c|c|} \hline \text{Range} …

11 r  distributions  poisson-distribution  mean  sample 

La implementación randomForest no permite el muestreo más allá del número de observaciones, incluso cuando se realiza un muestreo con reemplazo. ¿Por qué es esto? Funciona bien: rf <- randomForest(Species ~ ., iris, sampsize=c(1, 1, 1), replace=TRUE) rf <- randomForest(Species ~ ., iris, sampsize=3, replace=TRUE) Lo que quiero hacer: rf …

11 r  sampling  random-forest  stratification  oversampling 

Para analizar los recuentos de aves inflados a cero, me gustaría aplicar modelos de recuento inflados a cero usando el paquete R pscl . Sin embargo, al observar el ejemplo proporcionado en la documentación para una de las funciones principales ( ? Zeroinfl ), empiezo a dudar de cuál es …

11 r  poisson-distribution  zero-inflation 

Me gustaría usar lme4para ajustar una regresión de efectos mixtos y multcompcalcular las comparaciones por pares. Tengo un conjunto de datos complejo con múltiples predictores continuos y categóricos, pero mi pregunta se puede demostrar utilizando el ChickWeightconjunto de datos incorporado como ejemplo: m <- lmer(weight ~ Time * Diet + …

11 r  mixed-model  multiple-comparisons 

Me gustaría probar en qué regresión se ajusta mejor a mis datos. Mi variable dependiente es un recuento y tiene muchos ceros. Y necesitaría ayuda para determinar qué modelo y familia usar (poisson o cuasipoisson, o regresión de poisson inflada a cero) y cómo probar los supuestos. Regresión de Poisson: …

11 r  count-data  zero-inflation  poisson-regression 

Quiero visualizar los resultados de una agrupación (producida con protoclust{protoclust}) creando gráficos de dispersión para cada par de variables utilizadas para clasificar mis datos, colorear por clases y superponer las elipses para el intervalo de confianza del 95% para cada una de las clases (para inspeccionar cuáles las clases de …

11 r  confidence-interval  ggplot2  scatterplot 

Hola, estoy tomando un curso de postgrado en Estadística y hemos estado cubriendo estadísticas de exámenes y otros conceptos. Sin embargo, a menudo puedo aplicar las fórmulas y desarrollar una especie de intuición sobre cómo funcionan las cosas, pero a menudo me queda la sensación de que tal vez si …

11 r  hypothesis-testing  sas  simulation  computational-statistics