Dans les modèles linéaires et en particulier dans l'ANOVA, un contraste est une combinaison linéaire de paramètres avec des coefficients totalisant zéro. Il est utilisé pour tester l'hypothèse nulle correspondante. Les contrastes sont particulièrement souvent utilisés avec des prédicteurs catégoriels (facteurs) pour faire des comparaisons entre les groupes (catégories). [Voir aussi la balise 'categorical-encoding']
Veuillez fournir le code R qui permet d'effectuer une ANOVA entre sujets avec -3, -1, 1, 3 contrastes. Je comprends qu'il y a un débat concernant le type de somme de carrés (SS) approprié pour une telle analyse. Cependant, comme le type par défaut de SS utilisé dans SAS et SPSS (Type III) est...
J'ai effectué une ANOVA de mesures répétées en R, comme suit: aov_velocity = aov(Velocity ~ Material + Error(Subject/(Material)), data=scrd) summary(aov_velocity) Quelle syntaxe en R peut être utilisée pour effectuer un test post hoc après une ANOVA avec des mesures répétées? Le test de Tukey avec...
Je ne peux pas comprendre l'utilisation des contrastes polynomiaux dans l'ajustement de régression. En particulier, je me réfère à un codage utilisé par Rafin d'exprimer une variable d'intervalle (variable ordinale avec des niveaux également espacés), décrit sur cette page . Dans l'exemple de cette...
Disons que j'ai une simple expérience factorielle 2x2 sur laquelle je veux faire de l'ANOVA. Comme ça, par exemple: 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) +
Si vous ajustez un modèle linéaire ou un modèle mixte, il existe différents types de codages disponibles pour transformer une varibale catégorielle ou nominale en un certain nombre de variables pour lesquelles les paramètres sont estimés, tels que le conditionnement fictif (la valeur par défaut R)...
J'essaie d'exécuter une mesure répétée Anova dans R suivie de quelques contrastes spécifiques sur cet ensemble de données. Je pense que la bonne approche serait d'utiliser à Anova()partir du package de voiture. Permet d'illustrer ma question avec l'exemple tiré de l' ?Anovautilisation des...
J'ai posté ceci sur mathoverflow et personne ne répond: La méthode de Scheffé pour identifier les contrastes statistiquement significatifs est largement connue. Un contraste entre les moyennes , de populations est une combinaison linéaire dans laquelle , et un multiple scalaire d'un contraste est...
Utilisateur débutant de R aux prises avec une ANOVA de mesures répétées. J'ai un ensemble de données qui comprend un facteur entre les sujets avec 4 niveaux (codé dans une seule variable appelée «groupes») et un facteur entre les sujets avec 4 niveaux (codé en quatre variables distinctes «DV1»,...
Veuillez me donner une idée de la façon de recoder efficacement une variable catégorielle (facteur) dans l'ensemble de variables de contraste polynomiales orthogonales. Pour de nombreux types de variables de contraste (par exemple, déviation, simple, Helmert, etc.), la réussite est la suivante:...
Je fais une ANOVA à sens unique (par espèce) avec des contrastes personnalisés. [,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 où je compare l'intensité 0,5 contre 5, 5 contre 12,5 et ainsi de suite. Ce sont les données sur lesquelles je travaille avec les...
En utilisant R, j'ai ajusté un modèle linéaire pour une variable de réponse unique à partir d'un mélange de prédicteurs continus et discrets. C'est uber-basique, mais j'ai du mal à comprendre comment fonctionne un coefficient pour un facteur discret. Concept: De toute évidence, le coefficient de la...
J'ai un modèle de régression qui ressemble à ceci:Oui= β0+ β1X1+ β2X2+ β3X3+ β12X1X2+ β13X1X3+ β123X1X2X3Oui=β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 en notation R: y ~ x1 + x2 + x3 +...
J'utilise lmer dans R pour vérifier l'effet de condition ( cond) sur un résultat. Voici quelques données composées, où s est l'identifiant du sujet et a, bet csont des conditions. library("tidyr") library("dplyr") set.seed(123) temp <- data.frame(s = paste0("S", 1:30), a = rnorm(30, -2, 1), b =...
J'ai 2 facteurs Aet B(5 × 3) et une covariable Xdans une conception intra-sujet. Voici comment je spécifie mon modèle global: lme.out = lme(y~ A*B*X, random=~1|Subject, data=mydata) Mon interprétation est que je regarde un graphique y~x, où la pente change en raison de la covariable, et la ligne se...
La subitisation est l'énumération rapide et précise des affichages à faible numérosité, qui se distingue du comptage par une forte non-linéarité dans le tracé des temps de réponse. Vous trouverez ci-dessous un graphique représentatif de Watson, DG, Maylor, EA et Bruce, LAM (2007). Notez que les...
Supposons que j'ai un résultat continu yet deux prédicteurs factoriels, chacun avec deux niveaux. L'un de mes prédicteurs catégoriques drug, peut avoir deux niveaux ("A" ou "B"), l'autre l'est smokeYes. Lorsque j'exécute un modèle de régression, je peux choisir la ligne de base ou le niveau de...