He utilizado rpart.controlpara minsplit=2, y obtuvo los siguientes resultados de rpart()la función. Para evitar sobreajustar los datos, ¿necesito usar splits 3 o splits 7? ¿No debería usar splits 7? Por favor hagamelo saber. Variables realmente utilizadas en la construcción de árboles: [1] ct_a...
He estado usando R durante varios años, pero siempre he usado una GUI como Tinn-R, JGR y, más recientemente, R-Studio. Aunque me gusta la interfaz de R-Studio, creo que para programas más largos con pocos / sin gráficos, podría codificar de manera más eficiente aprovechando los accesos directos...
Necesito dibujar 20 distribuciones en un solo gráfico en R, y no se ve bien (desordenado) con diagrama de caja regular (20 cuadros) incluso con boxwex = 0.3. ¿Podría sugerirme cómo puedo trazar un tipo de diagrama de caja en R para las 20 distribuciones, con puntos para la mediana y solo una línea...
Tengo la siguiente serie de tiempo obtenido utilizando los datos publicados a continuación. Para un tamaño de ventana deslizante de 10, estoy tratando de calcular la divergencia KL entre el PMF de los valores dentro de la ventana deslizante actual y el PMF del historial con el objetivo final de...
Necesito hacer inferencia sobre un parámetro positivo . Para acomodar la positividad reparametrice . Usando la rutina MLE, calculé la estimación puntual y se para q . La propiedad de invariancia del MLE me da directamente una estimación puntual para p , pero no estoy seguro de cómo calcular se para...
Ejecuté un ANOVA encontrando, por ejemplo, una interacción entre género y grado de la que quiero saber en qué grados difieren los niños y las niñas, pero en muchos casos encuentro valores p (ajustados) de 0 y 1. ¿Cómo / por qué es esto posible? No parece correcto ... as.factor(gender) 1 16 16.2...
Bloqueado . Esta pregunta y sus respuestas están bloqueadas porque la pregunta está fuera de tema pero tiene un significado histórico. Actualmente no acepta nuevas respuestas o interacciones. Las funciones con las que estoy familiarizado incluyen escalar desde la base...
Entiendo que no deberíamos permitir que el mismo conjunto de datos que estamos analizando conduzca / defina cómo se verían las distribuciones anteriores en un análisis bayesiano. Específicamente, no es apropiado definir distribuciones previas para un análisis bayesiano basado en estadísticas...
Mientras comprobaba si dos variables estaban correlacionadas, observé que la aplicación de la correlación de Pearson arrojó números tan bajos como 0.1, lo que indica que no hay correlación. ¿Hay algo que pueda hacer para fortalecer esta afirmación? El conjunto de datos (subconjunto debido a los...
Estoy revisando el libro 'Introductory Time Series with R' de Cowpertwait y Metcalfe. En la página 36, dice que las líneas están en: . He leído aquí el foro R que las líneas están en . - 1 / n ± 2 / n--√-1/ /norte±2/ /norte-1/n \pm 2/\sqrt{n}± 1.96 / n--√±1,96/ /norte\pm 1.96/\sqrt{n} Ejecuté el...
Wikipedia sugiere que una forma de ver la confiabilidad entre evaluadores es usar un modelo de efectos aleatorios para calcular la correlación intraclase . El ejemplo de correlación intraclase habla de mirar σ2ασ2α+ σ2ϵσα2σα2+σϵ2\frac{\sigma_\alpha^2}{\sigma_\alpha^2+\sigma_\epsilon^2} de un...
Tengo el vector x <- c(1,2,3,4,5,5,5,6,6,6,6, 7,7,7,7,7,7,7,7,7,7,7,7,7,7,7,7,7,7,7, 7,7,7,7,7,7,7,7,8,8,8,8,9,9,9,10) (mi vector real tiene una longitud de> 10,000) y me gustaría encontrar los intervalos donde se encuentra el 90% de la densidad. ¿Es quantile(x, probs=c(0.05,0.95),...
Tengo el siguiente tipo de datos. He evaluado a 10 personas cada una repetidas 10 veces. Tengo una matriz de relación de 10x10 (relación entre todas las combinaciones de los individuos). set.seed(1234) mydata <- data.frame (gen = factor(rep(1:10, each = 10)), repl = factor(rep(1:10, 10)), yld...
¿Hay alguna diferencia entre los siguientes términos o son los mismos? Parcialidad Sesgo sistemático Errores sistemáticos Si existen algunas diferencias, explíquelas. ¿Se pueden reducir estos errores cuando se aumenta el tamaño de la muestra? ACTUALIZACIÓN: Mi campo de interés es la...
En R, ¿cómo especifico el modelo lmer sin efecto fijo global? Por ejemplo, si digo algo como lmer(y ~ (1 | group) + (0 + x | group), data = my_df) el modelo ajustado será yij=a+αi+βixijyij=a+αi+βixijy_{ij} = a + \alpha_{i} + \beta_{i} x_{ij} ¿Cómo encajo modelo? yij=αi+βixijyij=αi+βixijy_{ij}...
Estoy tratando de visualizar algunos datos del consumidor, que tiene 4 categorías. Los usuarios son libres de cambiar entre diferentes categorías. Me gustaría visualizar los últimos tres o cuatro interruptores para cada individuo. Entonces, comenzaríamos con un gráfico con una columna con 4...
Tengo parcelas de 13 grupos que muestro en una parcela. Los grupos tienen poblaciones desequilibradas y normalmente no están distribuidos. Quiero mostrar qué pares son estadísticamente similares (es decir, tienen kruskal.test p-value <0.05) colocando a, b, c, etc. en la parte superior de los...
Si tengo una variable con 4 niveles, en teoría necesito usar 3 variables ficticias. En la práctica, ¿cómo se lleva a cabo esto realmente? ¿Uso 0-3, uso 1-3 y dejo en blanco los 4? ¿Alguna sugerencia? NOTA: voy a estar trabajando en R. ACTUALIZACIÓN: ¿Qué sucedería si solo uso una columna que usa...
Digamos que genero la probabilidad de un resultado basado en cierto factor y trazo la curva de ese resultado. ¿Hay alguna forma de extraer la ecuación para esa curva de R? > mod = glm(winner~our_bid, data=mydat, family=binomial(link="logit")) > summary(mod) Call: glm(formula = winner ~...
Tengo algunos datos de series temporales en los que la variable medida es enteros positivos discretos (recuentos). Quiero probar si hay una tendencia al alza con el tiempo (o no). La variable independiente (x) está en el rango 0-500 y la variable dependiente (y) está en el rango 0-8. Pensé que...