Multivariate Statistik - Aufgabe 3 - Multivariate Anova

From StatWiki

Jump to: navigation, search

Contents


Beispiel 1

Man konstruiere einen Datensatz mit 2 stetigen Variablen und einer Gruppenvariable mit 3 Ausprägungen und 10 Beobachtungen pro Gruppe (also eine 15x3-Tabelle), bei dem beiden univariaten Varianzanalysen nicht signifikant sind, die multivariate jedoch schon (jeweils %5=α), und zwar so, dass die stetigen Variablen eine Korrelation von

  1. rXY=0,9 bzw.
  2. rXY=0,3
  3. rXY=0

aufweisen.

<Rworkspace> set.seed(123) groupData = function(mu,sd=1) { as.vector(sapply(mu, function(mu) { rnorm(10, mu, sd) })) } g=c(rep(1,10),rep(2,10),rep(3,10)) x=groupData(c(1,1.2,1.4))

library(foreign) sids = read.spss(paste(rfiles,"MANOVA_Diskriminanzanalyse_SIDS.raw",sep="/"),to.data.frame=TRUE) iris = read.spss(paste(rfiles,"iris.sav",sep="/"),to.data.frame=TRUE) </Rworkspace>

Beispiel 1.1

\Sigma = \begin{pmatrix} \sigma_1^2 & \sigma_1 \cdot \sigma_2 \cdot 0.9 \\ \sigma_1 \cdot \sigma_2 \cdot 0.9 & \sigma_2^2 \end{pmatrix}

\begin{pmatrix} \sigma_1 & 0 \\ 0 & \sigma_2 \end{pmatrix} \cdot R \cdot \begin{pmatrix} \sigma_1 & 0 \\ 0 & \sigma_2 \end{pmatrix}

\rho = \frac{\mathrm{cov}(x,y)}{\sqrt{\mathrm{var}(x) \mathrm{var}(y)}}


1. Generiere Normalverteilte Beobachtungen ~N(0,1)

G1 G2 G3
X1 X2 X3
Y1 Y2 Y3

2. Generiere Bivariate Normalverteilungen mit Korrelationsmatrix

R = \begin{pmatrix} 1 & 0.9 \\ 0.9 & 1 \end{pmatrix}

R^{\frac{1}{2}}(\vec x,\vec y)


3. Parameter \vec \mu_1, \vec \mu_2, \vec \mu_3

G_1: \,

\vec \mu_1 + R^{\frac{1}{2}} \cdot (\vec x_1, \vec y_1)

\vec \mu_2 + R^{\frac{1}{2}} \cdot (\vec x_2, \vec y_2)

\vec \mu_3 + R^{\frac{1}{2}} \cdot (\vec x_3, \vec y_3)

Wähle μ1,μ2,μ3 derart, dass Unterschiede in 1-Achse

4. Varianzanalyse signifikant, falls \frac{\mathrm{MSB}}{\mathrm{MSR}}

MSR groß, falls σ groß.

\cos(\alpha) = \rho \cdot \frac{\sigma y}{\sigma x}

Wähle \sigma x \, und \sigma y \, derarts, dass Unterschiede in 2. Achse signifikant

c) Geht nicht, kein Unterschied! 


Error in loadNamespace(name) : there is no package called 'distrEx'
Calls: sys.load.image ... tryCatch -> tryCatchList -> tryCatchOne -> 
Execution halted

> sys.load.image("/var/www/localhost/htdocs/StatWiki/Rfiles/true", TRUE)
Error in loadNamespace(name) : there is no package called 'distrEx'
Calls: sys.load.image ... tryCatch -> tryCatchList -> tryCatchOne -> 
Execution halted

> sys.load.image("/var/www/localhost/htdocs/StatWiki/Rfiles/true", TRUE)
Error in loadNamespace(name) : there is no package called 'distrEx'
Calls: sys.load.image ... tryCatch -> tryCatchList -> tryCatchOne -> 
Execution halted

Beispiel 1.2 


Error in loadNamespace(name) : there is no package called 'distrEx'
Calls: sys.load.image ... tryCatch -> tryCatchList -> tryCatchOne -> 
Execution halted

> sys.load.image("/var/www/localhost/htdocs/StatWiki/Rfiles/true", TRUE)
Error in loadNamespace(name) : there is no package called 'distrEx'
Calls: sys.load.image ... tryCatch -> tryCatchList -> tryCatchOne -> 
Execution halted

> sys.load.image("/var/www/localhost/htdocs/StatWiki/Rfiles/true", TRUE)
Error in loadNamespace(name) : there is no package called 'distrEx'
Calls: sys.load.image ... tryCatch -> tryCatchList -> tryCatchOne -> 
Execution halted

Beispiel 2

File: MANOVA_Diskriminanzanalyse_SIDS.raw. Man prüfe sowohl mit univariater Varianzanalyse (auch graphisch), als auch mit multivariater Varianzanalyse, ob Puls, Geburtgewicht, Herz-und Lungenfunktion, bzw. Gestationsalter bei Säuglingen, die dem plötzlichen Säuglingstod (SIDS) erliegen, signifikant anders sind (%5=α) als bei gesunden. Man prüfe auch die Modellvoraussetzungen für sämtliche verwendete Verfahren (Varianzhomogenität, Normalverteilung).

Die Voraussetzungen sollten eigentlich auf die einzelnen Gruppen angewendet werden. 


Error in loadNamespace(name) : there is no package called 'distrEx'
Calls: sys.load.image ... tryCatch -> tryCatchList -> tryCatchOne -> 
Execution halted

> sys.load.image("/var/www/localhost/htdocs/StatWiki/Rfiles/true", TRUE)
Error in loadNamespace(name) : there is no package called 'distrEx'
Calls: sys.load.image ... tryCatch -> tryCatchList -> tryCatchOne -> 
Execution halted

> sys.load.image("/var/www/localhost/htdocs/StatWiki/Rfiles/true", TRUE)
Error in loadNamespace(name) : there is no package called 'distrEx'
Calls: sys.load.image ... tryCatch -> tryCatchList -> tryCatchOne -> 
Execution halted

> sys.load.image("/var/www/localhost/htdocs/StatWiki/Rfiles/true", TRUE)
Error in loadNamespace(name) : there is no package called 'distrEx'
Calls: sys.load.image ... tryCatch -> tryCatchList -> tryCatchOne -> 
Execution halted

> sys.load.image("/var/www/localhost/htdocs/StatWiki/Rfiles/true", TRUE)
Error in loadNamespace(name) : there is no package called 'distrEx'
Calls: sys.load.image ... tryCatch -> tryCatchList -> tryCatchOne -> 
Execution halted


Error in loadNamespace(name) : there is no package called 'distrEx'
Calls: sys.load.image ... tryCatch -> tryCatchList -> tryCatchOne -> 
Execution halted

Beispiel 3

File: iris.sav. Man prüfe mittels multivariater Varianzanalyse, ob es einen Unterschied in Kelchblatt (sepal) und Blütenblatt (petal) bei den drei verschiedenen Arten der Iris gibt (%5=α). Falls ja, prüfe man mit univariaten Varianzanalysen, wo es Unterschiede gibt (in welchen Variablen bzw. zwischen welchen Gruppen) und stelle Unterschiede auch grafisch dar. Die Voraussetzungen der verwendeten Verfahren sind ebenfalls zu prüfen.

Vorraussetzungen nicht gegeben.

Daher Nicht-Parametrische Statistiken verwenden (z.b. Kurskal-Wallis). Das Problem bei Anova: Die Zerlegung stimmt trotzdem, aber die Annahmen für die F-Verteilungen stimmen nicht mehr und daher eine falsche Aussage liefern. 


Error in loadNamespace(name) : there is no package called 'distrEx'
Calls: sys.load.image ... tryCatch -> tryCatchList -> tryCatchOne -> 
Execution halted
Retrieved from "http://www.eisber.net/StatWiki/index.php/Multivariate_Statistik_-_Aufgabe_3_-_Multivariate_Anova"
Personal tools