[En preguntas recientes, estaba buscando generar vectores aleatorios en R , y quería compartir esa "investigación" como un Q&A independiente sobre un punto específico.] La generación de datos aleatorios con correlación se puede hacer usando la descomposición de Cholesky de la matriz de...
Tengo un conjunto de datos de 17 personas, clasificando 77 declaraciones. Quiero extraer componentes principales en una matriz de correlación transpuesta de correlaciones entre personas (como variables) a través de declaraciones (como casos). Yo sé, es raro, se llama Metodología Q . Quiero...
Tengo dos muestras bastante pequeñas ( para cada una) y obtengo una potencia de prueba t de dos muestras para ellas 0.49, que es relativamente pequeña (calculada por R ). Sin embargo, la prueba t de Welch Two Sample da un valor p de aquí, por lo que la hipótesis nula puede ser rechazada.n =...
Tengo la siguiente tabla de valores: 25 75 38 162 La razón de posibilidades es 0.7037 y el registro (OR) es -0.3514. Para una tabla de contingencia con valores a, b, c y d, la varianza de log (OR) viene dada por (1/a + 1/b + 1/c + 1/d) ¿Cómo puedo calcular el p.value de log (OR) a partir de...
La asimetría y la curtosis se definen como: ζ3=mi[ ( X- μ)3]mi[ ( X- μ)2]3 / 2=μ3σ3ζ3=mi[(X-μ)3]mi[(X-μ)2]3/ /2=μ3σ3\zeta_3 = \frac{E[(X-\mu)^3]}{E[(X-\mu)^2]^{3/2}} = \frac{\mu_3}{\sigma^3} ζ4 4=mi[ ( X- μ)4 4]mi[ ( X- μ)2]2=μ4 4σ4 4ζ4 4=mi[(X-μ)4 4]mi[(X-μ)2]2=μ4 4σ4 4\zeta_4 =...
Estoy realizando comparaciones post hoc después de una prueba de Kruskal-Wallis. Estoy usando el paquete PMCMR . > posthoc.kruskal.nemenyi.test( preference ~ instrument) Pairwise comparisons using Tukey and Kramer (Nemenyi) test with Tukey-Dist approximation for independent samples data:...
Tengo dos series de tiempo (parámetros de un modelo para hombres y mujeres) y mi objetivo es identificar un modelo ARIMA apropiado para hacer pronósticos. Mi serie de tiempo se ve así: El gráfico y el ACF muestran no estacionario (los picos del ACF se cortan muy lentamente). Por lo tanto, uso la...
Estoy tratando de determinar qué alfa usar en mi glmnetfunción, pero el archivo de ayuda me dice: Tenga en cuenta que cv.glmnet NO busca valores para alfa. Se debe proporcionar un valor específico; de lo contrario, se supone alfa = 1 de forma predeterminada. Si los usuarios desean validar...
Supongamos que tengo un conjunto de datos en el que la incertidumbre en las mediciones (que provienen de la propagación de errores sistemáticos del aparato de medición) es diferente para cada punto. Si hago una regresión lineal en el conjunto de datos, ¿cómo calculo la incertidumbre en la...
Mi pregunta trata sobre los árboles de clasificación . Considere el siguiente ejemplo del conjunto de datos de Iris: Quiero seleccionar manualmente el mejor predictor para la primera división. Según el algoritmo CART, la mejor característica para hacer una división es la que maximiza la...
Estoy tratando de entender el Alfa de Cronbach intuitivamente. ¿Cuál es la idea general detrás de esta construcción? ¿Qué propiedades estaban tratando de asegurar que
Estoy ejecutando algunas optimizaciones con la implementación de BFGS de optim. La función objetivo es en realidad un algoritmo computacional, no solo matemático. Descubrí que cuando agrego una penalización L1, las cosas se ralentizan bastante. ¿Por qué podría ser esto? ¿Hay algo en L1 que...
Tengo un conjunto de datos x, y que estoy usando para construir un bosque aleatorio. Los datos x son un vector de valores que incluye algunos NA. Así que solía rfImputemanejar los datos faltantes y crear un bosque aleatorio. Ahora tengo una nueva observación invisible x (con un NA) y quiero...
Estoy usando R para realizar una regresión lineal. He visto formas de calcular los intervalos de predicción, pero estos dependen de datos homoscedastic. ¿Hay alguna forma de calcular los intervalos de predicción con datos
Estoy trabajando con el conjunto de datos "géiser" del paquete MASS y comparando las estimaciones de densidad del núcleo del paquete np. Mi problema es comprender la estimación de densidad utilizando la validación cruzada de mínimos cuadrados y el núcleo
Comunidad SE, espero obtener algunas ideas sobre el siguiente problema. Dado un modelo de regresión lineal simple Bajo una función de probabilidad gaussiana con términos de error homoscedastic, la distribución condicional de la variable dependiente toma la forma Asigno un conjugado condicional (no...
Cerrado. Esta pregunta está fuera de tema . Actualmente no está aceptando respuestas. ¿Quieres mejorar esta pregunta? Actualice la pregunta para que esté en el tema de Cross Validated. Cerrado hace 4 años . Estoy trabajando en algunos puntos de...
Tengo una regresión logística binaria con solo un predictor binario de factor fijo. La razón por la que no lo hago como un Chi cuadrado o la prueba exacta de Fisher es que también tengo una serie de factores aleatorios (hay múltiples puntos de datos por individuo y los individuos están en grupos,...
Estoy usando la función Deming proporcionada por Terry T. en este hilo archivado de r-help . Estoy comparando dos métodos, por lo que tengo datos que se ven así: y x stdy stdx 1 1.2 0.23 0.67 2 1.8 0.05 0.89 4 7.5 1.13 0.44 ... ... ... ... He hecho mi regresión de Deming (también llamada...
¿Cuál es la diferencia entre (1|DNA.concentration/mouse.id)y (DNA.concentration|mouse.id)? ¿Qué significan los símbolos |y /significan dentro de la sintaxis para el efecto aleatorio?