par5 <- '1'
par4 <- '2'
par3 <- '4'
par2 <- '3'
par1 <- '5'
cat1 <- as.numeric(par1) #
cat2<- as.numeric(par2) #
cat3 <- as.numeric(par3)
cat4 <-as.numeric(par4)
cat5 <-as.numeric(par5)
x <- t(x)
x1<-as.numeric(x[,cat1])
wf1<-as.character(x[,cat2])
wf2 <- as.character(x[,cat3])
bf1 <- as.character(x[,cat4])
sid<- as.character(x[,cat5]) # author of ez changed within subjects variable name from sid to wid
xdf<-data.frame(x1,wf1, wf2, bf1, sid)
(V1<-dimnames(y)[[1]][cat1])
(V2<-dimnames(y)[[1]][cat2])
(V3 <-dimnames(y)[[1]][cat3])
(V4 <-dimnames(y)[[1]][cat4])
(V5 <-dimnames(y)[[1]][cat5])
names(xdf)<-c(V1, V2, V3, V4, V5)
library(ez)
library(Cairo)
(ezout <- ezANOVA(data=xdf, dv=.(mean_rt), wid=.(sid), within=.(cue, flanker), between=.(group) ) )
load(file='createtable')
a<-table.start()
nr <- nrow(ezout$ANOVA)
nc <- ncol(ezout$ANOVA)
a<-table.row.start(a)
a<-table.element(a,'Repeated Measures ANOVA', nc+1,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'effect', 1,TRUE)
a<-table.element(a,'Dfn',1,TRUE)
a<-table.element(a,'DFd', 1,TRUE)
a<-table.element(a, 'F', 1,TRUE)
a<-table.element(a,'p', 1,TRUE)
a<-table.element(a,'p<0.05', 1,TRUE)
a<-table.element(a, 'ges', 1,TRUE) # generalized eta-sq - was partial eta-sq in earlier version
a<-table.row.end(a)
for ( i in 1:nr){
a<-table.row.start(a)
a<-table.element(a,ezout$ANOVA$Effect[i], 1, TRUE)
for(j in 2:nc){
if ( j != 6) # author of ez reduced number of columns in output from 8
a<-table.element(a,round(ezout$ANOVA[[j]][i], digits=3), 1, FALSE)
else a<-table.element(a, ezout$ANOVA[[j]][i], 1, FALSE)
}
a<-table.row.end(a)
}
a<-table.end(a)
table.save(a,file='mytable.tab')
a<-table.start()
nr <- nrow(ezout$Mauchly)
nc <- ncol(ezout$Mauchly)
a<-table.row.start(a)
a<-table.element(a,'Mauchlys Test for Sphericity', nc+1,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'effect', 1,TRUE)
a<-table.element(a,'W',1,TRUE)
a<-table.element(a,'p', 1,TRUE)
a<-table.element(a,'p<0.05', 1,TRUE)
a<-table.row.end(a)
for ( i in 1:nr){
a<-table.row.start(a)
a<-table.element(a,ezout$Mauchly$Effect[i], 1, TRUE)
for(j in 2:nc){
if (j != 4)
a<-table.element(a,round(ezout$Mauchly[[j]][i], digits = 3), 1, FALSE)
else
a<-table.element(a,ezout$Mauchly[[j]][i], 1, FALSE)
}
a<-table.row.end(a)
}
a<-table.end(a)
table.save(a,file='mytable1.tab')
a<-table.start()
nr <- nrow(ezout$Spher)
nc <- ncol(ezout$Sphe)
a<-table.row.start(a)
a<-table.element(a,'Sphericity Corrections', nc+1,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'effect', 1,TRUE)
a<-table.element(a,'GGe',1,TRUE)
a<-table.element(a,'p[GG]', 1,TRUE)
a<-table.element(a,'p[GG]<0.05', 1,TRUE)
a<-table.element(a,'HFe', 1,TRUE)
a<-table.element(a,'p[HF]', 1,TRUE)
a<-table.element(a,'p[HF]<0.05', 1,TRUE)
a<-table.row.end(a)
for ( i in 1:nr){
a<-table.row.start(a)
a<-table.element(a,ezout$Spher$Effect[i], 1, TRUE)
for(j in 2:nc){
if ( ! ((j == 4) | (j == 7)) )
a<-table.element(a,round(ezout$Spher[[j]][i], digits=3), 1, FALSE)
else
a<-table.element(a,ezout$Spher[[j]][i], 1, FALSE)
}
a<-table.row.end(a)
}
a<-table.end(a)
table.save(a,file='mytable2.tab')
ezP.between<-ezPlot(data = xdf, dv = .(mean_rt), between = .(group), wid = .(sid), do_lines=FALSE, x_lab='group', y_lab='RT' , x=.(group))
bitmap(file = 'between.cairo')
print(ezP.between)
dev.off()
ezstats_between<-ezStats(data = xdf, dv = .(mean_rt), between =.(group), wid = .(sid))
a<-table.start()
nr <- nrow(ezstats_between)
nc <- ncol(ezstats_between)
a<-table.row.start(a)
a<-table.element(a,'Between Effects Comparisons', nc+1,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
for(i in 1:nc){
a<-table.element(a, names(ezstats_between)[i], 1,TRUE)
}
a<-table.row.end(a)
for ( i in 1:nr){
a<-table.row.start(a)
a<-table.element(a,ezstats_between[[1]][i], 1, TRUE)
for(j in 2:nc){
a<-table.element(a,ezstats_between[[j]][i], 1, FALSE)
}
a<-table.row.end(a)
}
a<-table.end(a)
table.save(a,file='mytable3.tab')
ezP.within<-ezPlot(data = xdf, dv = .(mean_rt), within = .(cue, flanker), wid = .(sid), do_lines=TRUE, x_lab='flanker', y_lab='RT' , x=.(flanker), split=.(cue), split_lab = 'cue')
bitmap(file = 'within.cairo')
print(ezP.within)
dev.off()
ezstats_within <- ezStats(data = xdf, dv = .(mean_rt), within = .(cue, flanker), wid = .(sid))
a<-table.start()
nr <- nrow(ezstats_within)
nc <- ncol(ezstats_within)
a<-table.row.start(a)
a<-table.element(a,'Within Effects Comparisons', nc+1,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
for(i in 1:nc){
a<-table.element(a, names(ezstats_within)[i], 1,TRUE)
}
a<-table.row.end(a)
for ( i in 1:nr){
a<-table.row.start(a)
a<-table.element(a,ezstats_within[[1]][i], 1, TRUE)
for(j in 2:nc){
a<-table.element(a, ezstats_within[[j]][i], 1, FALSE)
}
a<-table.row.end(a)
}
a<-table.end(a)
table.save(a,file='mytable4.tab')
|