Multivariate Statistik - Abschluss Aufgabe

From StatWiki
Jump to: navigation, search

Anova für Faktoren sex und Unique national programme code, abhängige Variablen IC5a - IC5l

Diskriminanzanalyse für die erzielten Werte in den beiden Tests, wie gut kann zwischen den einzelnen Schultypen diskriminiert werden.

MAnova Unterschiede zwischen den Faktoren SEX und Schultyp in Bezug auf die Testwerte

Anova Unterschiede im Faktor Geschlecht in Bezug auf die Testwerte über die unterschiedlichen Schultypen


Mathematik-Leistung vs. Computernutzung

Präsentation ist am 27.7.2006 (max. 20min pro Gruppe)

Abgabetermin ist am 12.7.2006


Engine.php: found forbidden command "read.csv("
in
  1. Das hier ist R-code

library(xtable)

  1. load data

pisa2 = read.csv(paste(rfiles,"8_PISA2.csv",sep="/"))

  1. fix factor
  2. rawLabels = attr(raw,"label.table")$PROGN

raw = read.csv(paste(rfiles,"8_PISA2.raw.csv",sep="/")) rawLabels = read.csv(paste(rfiles,"8_PISA2.raw.label.csv",sep="/"))

  1. colnames(pisa2)[1] = "sex"
  2. colnames(pisa2)[2] = "school"
  1. pisa2$sex = factor(as.integer(pisa2$sex),labels=c("Female","Male"))
  1. get proper labels and cut off AUT: and the english translation
  2. pisa2$school = as.factor(sub("\\s*\\(.*\\)","",sub("AUT: ","",names(rawLabels)[sapply(raw$PROGN, grep, rawLabels, grep)])))

Aufgabe 2

Anova

Engine.php: unknown attribute(s) "direct"
in

dp=data.frame(t(sapply(colnames(pisa2)[15:22], function(i) { f=formula(paste(i,"~sex*school"))

# test auf homogenität (bartlett.test) # test auf normalvt über alle elemente # test auf normalvt für jede gruppe c(

       bartlett.test(pisa2[,i],pisa2$sex)$p.value,
       bartlett.test(pisa2[,i],pisa2$school)$p.value,

shapiro.test(pisa2[,i])$p.value, by(pisa2[,i], pisa2$sex, function(x) { shapiro.test(x)$p.value }), by(pisa2[,i], pisa2$school, function(x) { shapiro.test(x)$p.value }), summary(aov(f,data=pisa2))[[1]][["Pr(>F)"]][-4]

       )

})))

colnames(dp)[1:3] = c("Bartlett.sex","Bartlett.school","Shapiro") colnames(dp)[12:14] = c("F.sex","F.school","F.sex*school")

print(xtable(dp,digits=c(10,2,4,2,rep(3,8),3,-4,2)),type="html")

Die Spalten Female bist Polytechnische.Schule sind p-Werte eines Shapiro Wilk Testest auf Normalverteilung. Einzig PV1READ macht Probleme.

Box Plot

REngine.php: > rpdf<-'/var/www/localhost/htdocs/StatWiki/Rfiles/R/ef386787a10017c719520b39f307c02f6ca9ae87_%i.pdf'
> rpdfno<-0
> rhtml<-''
> rfiles<-'/var/www/localhost/htdocs/StatWiki/Rfiles/R/'
> source('/var/www/localhost/htdocs/StatWiki/Rfiles/R/@.R')
> rout<-'display'
> cat('<!--- Start of program --->\n')
<!--- Start of program --->
> pdf(rpdf)
>  par(mar=c(3,5,4,2)+0.1)
>  boxplot(PV1MATH~sex,data=pisa2,horizontal=T,las=1,main="PV1MATH")
Error in eval(expr, envir, enclos) : object 'pisa2' not found
Calls: boxplot -> boxplot.formula -> eval -> eval
Execution halted
in
 
pdf(rpdf)
par(mar=c(3,5,4,2)+0.1)
boxplot(PV1MATH~sex,data=pisa2,horizontal=T,las=1,main="PV1MATH")
REngine.php: > rpdf<-'/var/www/localhost/htdocs/StatWiki/Rfiles/R/a3b50897a6a2236a62a1019c3aaabda680d684a3_%i.pdf'
> rpdfno<-0
> rhtml<-''
> rfiles<-'/var/www/localhost/htdocs/StatWiki/Rfiles/R/'
> source('/var/www/localhost/htdocs/StatWiki/Rfiles/R/@.R')
> rout<-'display'
> cat('<!--- Start of program --->\n')
<!--- Start of program --->
> pdf(rpdf)
>  par(mar=c(3,12,4,2)+0.1)
>  boxplot(PV1MATH2~school,data=pisa2,horizontal=T,las=1,main="PV1MATH2")
Error in eval(expr, envir, enclos) : object 'pisa2' not found
Calls: boxplot -> boxplot.formula -> eval -> eval
Execution halted
in
 
pdf(rpdf) 
par(mar=c(3,12,4,2)+0.1)
boxplot(PV1MATH2~school,data=pisa2,horizontal=T,las=1,main="PV1MATH2")
REngine.php: > rpdf<-'/var/www/localhost/htdocs/StatWiki/Rfiles/R/ea6531308b5fccd71ef0a975120116c231b50a05_%i.pdf'
> rpdfno<-0
> rhtml<-''
> rfiles<-'/var/www/localhost/htdocs/StatWiki/Rfiles/R/'
> source('/var/www/localhost/htdocs/StatWiki/Rfiles/R/@.R')
> rout<-'display'
> cat('<!--- Start of program --->\n')
<!--- Start of program --->
> pdf(rpdf)
>  par(mar=c(3,14,4,2)+0.1)
>  boxplot(PV1MATH2~sex+school,data=pisa2,horizontal=T,las=1,main="PV1MATH2")
Error in eval(expr, envir, enclos) : object 'pisa2' not found
Calls: boxplot -> boxplot.formula -> eval -> eval
Execution halted
in
 
pdf(rpdf) 
par(mar=c(3,14,4,2)+0.1)
boxplot(PV1MATH2~sex+school,data=pisa2,horizontal=T,las=1,main="PV1MATH2") 

Pairwise T-Test

REngine.php: <!--- Start of program --->
Error in factor(g) : object 'pisa2' not found
Calls: pairwise.t.test -> factor
Execution halted
in
  1. das alternative geht nicht so recht

pairwise.t.test(pisa2$PV1MATH2,pisa2$school,alternative="less")