Em modelos lineares e particularmente em ANOVA, um contraste é uma combinação linear de parâmetros com coeficientes que somam zero. É usado para testar a hipótese nula correspondente. Os contrastes são especialmente usados com preditores categóricos (fatores) para fazer comparações entre os grupos (categorias). [Consulte também a tag 'codificação categórica']
Forneça um código R que permita realizar uma ANOVA entre sujeitos com -3, -1, 1, 3 contrastes. Entendo que há um debate sobre o tipo apropriado de soma de quadrados (SS) para essa análise. No entanto, como o tipo padrão de SS usado no SAS e SPSS (Tipo III) é considerado o padrão na minha área....
Realizei ANOVA de medidas repetidas em R, da seguinte maneira: aov_velocity = aov(Velocity ~ Material + Error(Subject/(Material)), data=scrd) summary(aov_velocity) Qual sintaxe em R pode ser usada para realizar um teste post hoc após uma ANOVA com medidas repetidas? O teste de Tukey com...
Não consigo entender o uso de contrastes polinomiais no ajuste de regressão. Refiro-me, em particular, a uma codificação usada Rpara expressar uma variável de intervalo (variável ordinal com níveis igualmente espaçados), descrita nesta página . No exemplo dessa página , se eu entendi corretamente,...
Digamos que eu tenha um experimento fatorial 2x2 simples no qual eu quero fazer a ANOVA. Assim, por exemplo: 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) +
Se você ajustar um modelo linear ou um modelo misto, existem diferentes tipos de codificação disponíveis para transformar uma variedade categórica ou nominal em um número de variáveis para as quais os parâmetros são estimados, como a dummy conding (o padrão R) e a codificação de efeitos. Ouvi...
Estou tentando executar uma Anova de medidas repetidas em R, seguida de alguns contrastes específicos nesse conjunto de dados. Eu acho que a abordagem correta seria usar Anova()o pacote do carro. Vamos ilustrar minha pergunta com o exemplo retirado do ?Anovauso dos OBrienKaiserdados (Nota: omiti...
Eu postei isso no mathoverflow e ninguém está respondendo: O método de Scheffé para identificar contrastes estatisticamente significativos é amplamente conhecido. Um contraste entre as médias , i = 1 , … , r de r populações é uma combinação linear ∑ r i = 1 c i μ i na qual ∑ r i = 1 c i =...
Usuário iniciante de R aqui lutando com uma ANOVA de medidas repetidas. Eu tenho um conjunto de dados que consiste em um fator entre sujeitos com 4 níveis (codificado em uma única variável chamada 'grupos') e um fator dentro de sujeitos com 4 níveis (codificado em quatro variáveis separadas...
Por favor, me dê uma idéia de como recodificar com eficiência uma variável categórica (fator) no conjunto de variáveis de contraste polinomial ortogonais. Para muitos tipos de variáveis de contraste (por exemplo, desvio, simples, Helmert etc.), o passe é: Componha a matriz dos coeficientes...
Estou fazendo uma ANOVA de sentido único (por espécie) com 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 onde eu comparo a intensidade 0,5 contra 5, 5 contra 12,5 e assim por diante. Estes são os dados em que estou trabalhando com...
Usando R, ajustei um modelo linear para uma variável de resposta única a partir de uma mistura de preditores contínuos e discretos. Isso é super básico, mas estou tendo problemas para entender como funciona um coeficiente para um fator discreto. Conceito: Obviamente, o coeficiente da variável...
Eu tenho um modelo de regressão parecido com este: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 ... ou na notação R: y ~ x1 + x2 + x3 + x1:x2 +...
Estou usando o lmer no R para verificar o efeito da condição ( cond) em algum resultado. Aqui estão alguns dados compostos, onde s é o identificador de sujeito e a, be csão condições. library("tidyr") library("dplyr") set.seed(123) temp <- data.frame(s = paste0("S", 1:30), a = rnorm(30, -2,...
Eu tenho 2 fatores Ae B(5 × 3) e uma covariável Xem um design de assunto interno. Veja como eu especifico meu modelo geral: lme.out = lme(y~ A*B*X, random=~1|Subject, data=mydata) Minha interpretação é que estou olhando para um gráfico y~x, onde a inclinação muda devido à covariável e a linha...
Subitização é a enumeração rápida e precisa de telas de baixa numerosidade, diferenciada da contagem por uma não linearidade acentuada no gráfico dos tempos de resposta. Abaixo está um gráfico representativo, de Watson, DG, Maylor, EA e Bruce, LAM (2007). Observe que o tempo médio de enumeração...
Suponha que eu tenha um resultado contínuo ye dois fatores preditores fatoriais, cada um com dois níveis. Um dos meus preditores categóricos,, drugpode ter dois níveis ("A" ou "B"), o outro é smokeYes. Quando executo um modelo de regressão, posso escolher a linha de base ou o nível de referência...