Mixed MANCOVA/ MANCOVA mit Messwiederholung

Varianzanalyse, Diskriminanzanalyse, Kontingenzanalyse, Faktorenanalyse, Clusteranalyse, MDS, ....

Moderator: EDi

Antworten
tinaR
Beiträge: 4
Registriert: Di Mai 18, 2021 7:38 pm

Mixed MANCOVA/ MANCOVA mit Messwiederholung

Beitrag von tinaR »

Hi, kann mir jemand weiterhelfen?
Ich möchte eine MANCOVA mit einem between- subjects- Faktor, zwei within- subjects- Faktoren und einer Kovariate in R berechnen. Es gibt sechs abhängige Variablen im Modell. Welche Funktion/Paket sind am Besten für mein Vorhaben geeignet? In SPSS habe ich die MANCOVA nach Anleitung gerechnet, würde jedoch lieber R als Statistikprogramm verwenden.
Ich habe es bereits mit der manova() Funktion probiert, bekomme aber falsche Werte (im Vergleich zu SPSS).

Code: Alles auswählen

manova(cbind(dv1, dv2, dv3, dv4, dv5, dv6) ~ covariate + between * within1 * within2 + Error(id/(within1*within2)), data = data)
Vielen Dank schon mal im Voraus
Benutzeravatar
EDi
Beiträge: 1599
Registriert: Sa Okt 08, 2016 3:39 pm

Re: Mixed MANCOVA/ MANCOVA mit Messwiederholung

Beitrag von EDi »

Ich habe es bereits mit der manova() Funktion probiert, bekomme aber falsche Werte (im Vergleich zu SPSS).
Ohne zu zeigen was du in SPSS gemacht hast kann niemand einschätzen was richtig oder falsch ist.

Bitte zeig uns ein reproduzierbares Beispiel, was du in SPSS gemacht hast und der erwartete output...
Bitte immer ein reproduzierbares Minimalbeispiel angeben. Meinungen gehören mir und geben nicht die meines Brötchengebers wieder.

Dieser Beitrag ist lizensiert unter einer CC BY 4.0 Lizenz
Bild.
tinaR
Beiträge: 4
Registriert: Di Mai 18, 2021 7:38 pm

Re: Mixed MANCOVA/ MANCOVA mit Messwiederholung

Beitrag von tinaR »

Vielen Dank für die schnelle Antwort!
Folgendes habe ich in SPSS versucht:
Ich habe die MANCOVA nach dieser Beschreibung gemacht http://www.statsmakemecry.com/smmcthebl ... -spss.html .
Der Between- Faktor ist zweistufig und es gibt vier Messwiederholungen, die sich auf zwei zweistufiege Within- Faktoren verteilen.
Als Output MANCOVA in SPSS erhalte ich:

Bild

Vielen Dank im Voraus!! Falls ich bei meiner Darstellung noch etwas offen gelassen habe, weist mich gerne darauf hin :)
bigben
Beiträge: 2771
Registriert: Mi Okt 12, 2016 9:09 am

Re: Mixed MANCOVA/ MANCOVA mit Messwiederholung

Beitrag von bigben »

Mit MANCoVA kenne ich mich nicht aus. Bei ANOVA-Vergleichrn zwischen R und SPSS geht es oft um sog. Typ I versus Typ III sums of squares. Dazu findet Google viel, z.B auch dir Funktion ezanova. Wie gesagt, ob sich das aufs multivariate verallgemeinern lässt weiß ich nicht, wäre aber mal eine Suche wert.
LG, Bernhard
---
Programmiere stets so, dass die Maxime Deines Programmierstils Grundlage allgemeiner Gesetzgebung sein könnte
tinaR
Beiträge: 4
Registriert: Di Mai 18, 2021 7:38 pm

Re: Mixed MANCOVA/ MANCOVA mit Messwiederholung

Beitrag von tinaR »

Danke für deinen Hinweis Bernhard, ich lese mich gleich morgen mal ein.
Für ein besseres verständnis noch mal der R Output: (von

Code: Alles auswählen

summary(mancovaModell, test="Pillai")
)

Code: Alles auswählen

...
Error: Within
                                             Df   Pillai approx F num Df den Df    Pr(>F)    
data$covariate                                1 0.055056   15.673      6   1614 < 2.2e-16 ***
data$between                           	1 0.085742   25.228      6   1614 < 2.2e-16 ***
data$within1                                  1 0.161122   51.666      6   1614 < 2.2e-16 ***
data$within2                                  1 0.029710    8.237      6   1614 8.411e-09 ***
data$between:data$within1                     1 0.009592    2.605      6   1614   0.01622 *  
data$between:data$within2                     1 0.004653    1.258      6   1614   0.27402    
data$within1:data$within2                 1 0.010176    2.765      6   1614   0.01118 *  
data$between:data$within1:data$within2        1 0.002321    0.626      6   1614   0.70977    
Residuals                                  1619                                              
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Ich habe auch schon verschiedene Funktionen aus afex und jmv probiert, aber nie die selben Ergebnisse für eine MANCOVA erhalten. Mit jmv kann ich eine ANCOVA rechnen die ich auch so in SPSS reproduzieren kann. Stimmt den mein Syntax für die manova() Funktion?
Danke im Voraus!!
Liebe Grüße Tina
Benutzeravatar
EDi
Beiträge: 1599
Registriert: Sa Okt 08, 2016 3:39 pm

Re: Mixed MANCOVA/ MANCOVA mit Messwiederholung

Beitrag von EDi »

Ohne mich mit SPSS auszukennen, sieht man dass du dort zwei Dreifachinteraktionen hast, im Model in R aber nur eine! Das sind also zwei verschiedene Modelle.

Wie sieht denn das Model aus, dass du fitten möchtest? Kannst du es uns in statistischer Model Notation hier aufschreiben?
Bitte immer ein reproduzierbares Minimalbeispiel angeben. Meinungen gehören mir und geben nicht die meines Brötchengebers wieder.

Dieser Beitrag ist lizensiert unter einer CC BY 4.0 Lizenz
Bild.
tinaR
Beiträge: 4
Registriert: Di Mai 18, 2021 7:38 pm

Re: Mixed MANCOVA/ MANCOVA mit Messwiederholung

Beitrag von tinaR »

Hey EDi,

danke für den Hinweis!

Das Model in R hab ich angepasst, sodass ich die selben Interaktionsterme habe, bekomme aber immer noch unterschiedliche Ergebnisse. :(

Code: Alles auswählen

mancovaModell <- manova(avs ~ covariate *  within1 * within2 + between * within1 * within2 + Error(data$ID/(within1*within2)), data = data)
summary(mancovaModell, test="Pillai")

...
Error: Within
                                             Df   Pillai approx F num Df den Df    Pr(>F)    
covariate                                    1 0.055627   15.816      6   1611 < 2.2e-16 ***
within1                                      1 0.161601   51.753      6   1611 < 2.2e-16 ***
within2                                      1 0.029786    8.243      6   1611 8.271e-09 ***
between                                      1 0.085998   25.263      6   1611 < 2.2e-16 ***
covariate:within1                            1 0.020194    5.534      6   1611 1.093e-05 ***
covariate:within2                            1 0.001368    0.368      6   1611   0.89956    
within1:within2                              1 0.010186    2.763      6   1611   0.01124 *  
within1:between                              1 0.007560    2.045      6   1611   0.05682 .  
within2:between                              1 0.004513    1.217      6   1611   0.29437    
covariate:within1:within2                    1 0.001195    0.321      6   1611   0.92612    
within1:within2:between                      1 0.002267    0.610      6   1611   0.72252    
Residuals                                    1616                                              
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Ich bin leider keine Expertin auf diesem Gebiet, weshalb ich kurz nachfrage, was genau du damit meinst bei deiner Frage welches Model ich denn fitten möchte? Gerne würde ich die statistische Model Notation schreiben, wenn ich genau wüsste wonach gefragt wird :)

Ich habe auch noch diese herangehensweise über die lm Funktion gefunden https://stackoverflow.com/questions/630 ... ova-with-r allerdings ist hier die Frage, wie man die Kovariate integrieren würde.

Danke im Voraus!
LG Tina
Antworten