En modelos lineales y particularmente en ANOVA, un contraste es una combinación lineal de parámetros con coeficientes que suman cero. Se utiliza para probar la hipótesis nula correspondiente. Los contrastes se utilizan especialmente con predictores categóricos (factores) para hacer comparaciones entre los grupos (categorías). [Consulte también la etiqueta "codificación-categórica"]
Proporcione un código R que le permita a uno realizar un ANOVA entre sujetos con -3, -1, 1, 3 contrastes. Entiendo que hay un debate con respecto al tipo apropiado de Suma de cuadrados (SS) para dicho análisis. Sin embargo, como el tipo predeterminado de SS utilizado en SAS y SPSS (Tipo III) se...
He realizado un ANOVA de medidas repetidas en R, de la siguiente manera: aov_velocity = aov(Velocity ~ Material + Error(Subject/(Material)), data=scrd) summary(aov_velocity) ¿Qué sintaxis en R puede usarse para realizar una prueba post hoc después de un ANOVA con medidas repetidas? ¿Sería...
No puedo entender el uso de contrastes polinómicos en el ajuste de regresión. En particular, me refiero a una codificación utilizada por Rpara expresar una variable de intervalo (variable ordinal con niveles igualmente espaciados), descrita en esta página . En el ejemplo de esa página , si entendí...
Digamos que tengo un simple experimento factorial 2x2 en el que quiero hacer ANOVA. Así, por ejemplo: d <- data.frame(a=factor(sample(c('a1','a2'), 100, rep=T)), b=factor(sample(c('b1','b2'), 100, rep=T))); d$y <- as.numeric(d$a)*rnorm(100, mean=.75, sd=1) + as.numeric(d$b)*rnorm(100,...
Si se ajusta a un modelo lineal o un modelo mixto, hay diferentes tipos de codificaciones disponibles para transformar una variable categórica o nominal en una serie de variables para las que se estiman los parámetros, como la codificación ficticia (el valor predeterminado de R) y la codificación...
Estoy tratando de ejecutar medidas repetidas Anova en R seguido de algunos contrastes específicos en ese conjunto de datos. Creo que el enfoque correcto sería usar Anova()el paquete del automóvil. Vamos a ilustrar mi pregunta con el ejemplo tomado del ?Anovauso de los OBrienKaiserdatos (Nota:...
Publiqué esto en mathoverflow y nadie respondió: El método de Scheffé para identificar contrastes estadísticamente significativos es ampliamente conocido. Un contraste entre las medias , de poblaciones es una combinación lineal en la que , y un múltiplo escalar de un contraste es esencialmente el...
Usuario principiante de R aquí luchando con un ANOVA de medidas repetidas. Tengo un conjunto de datos que consiste en uno entre el factor de sujetos con 4 niveles (codificado en una sola variable llamada 'grupos') y uno dentro del factor de sujetos con 4 niveles (codificado en cuatro variables...
Dame una idea de cómo recodificar eficientemente una variable categórica (factor) en el conjunto de variables de contraste polinomiales ortogonales. Para muchos tipos de variables de contraste (por ejemplo, desviación, simple, Helmert, etc.) el pase es: Componga la matriz de coeficientes de...
Estoy haciendo un ANOVA unidireccional (por especie) con contrastes personalizados. [,1] [,2] [,3] [,4] 0.5 -1 0 0 0 5 1 -1 0 0 12.5 0 1 -1 0 25 0 0 1 -1 50 0 0 0 1 donde comparo la intensidad 0.5 contra 5, 5 contra 12.5 y así sucesivamente. Estos son los datos en los que estoy...
Usando R, he ajustado un modelo lineal para una variable de respuesta única a partir de una combinación de predictores continuos y discretos. Esto es súper básico, pero tengo problemas para comprender cómo funciona un coeficiente para un factor discreto. Concepto: Obviamente, el coeficiente de la...
Tengo un modelo de regresión que se ve así:Y=β0+β1X1+β2X2+β3X3+β12X1X2+β13X1X3+β123X1X2X3Y=β0+β1X1+β2X2+β3X3+β12X1X2+β13X1X3+β123X1X2X3Y = \beta_0+\beta_1X_1 + \beta_2X_2 + \beta_3X_3 +\beta_{12}X_1X_2+\beta_{13}X_1X_3+\beta_{123}X_1X_2X_3 ... o en notación R: y ~ x1 + x2 + x3 + x1:x2 + x1:x3 +...
Estoy usando lmer en R para verificar el efecto de la condición ( cond) en algún resultado. Aquí están algunos hechos de seguridad de datos, donde s es el identificador de objeto y a, by cson condiciones. library("tidyr") library("dplyr") set.seed(123) temp <- data.frame(s = paste0("S", 1:30),...
Tengo 2 factores Ay B(5 × 3) y una covariable Xen un diseño dentro de la asignatura. Así es como especifico mi modelo general: lme.out = lme(y~ A*B*X, random=~1|Subject, data=mydata) Mi interpretación es que estoy mirando un gráfico y~x, donde la pendiente cambia debido a la covariable, y la...
La subitización es la enumeración rápida y precisa de las pantallas de baja numerosidad, que se distingue del conteo por una no linealidad aguda en la gráfica de tiempos de respuesta. A continuación se muestra una trama representativa, de Watson, DG, Maylor, EA y Bruce, LAM (2007). Observe que los...
Supongamos que tengo un resultado continuo yy dos predictores factoriales, cada uno con dos niveles. Uno de mis predictores categóricos drug, puede tener dos niveles ("A" o "B"), el otro es smokeYes. Cuando ejecuto un modelo de regresión, puedo elegir la línea base o el nivel de referencia drugpara...