x <- array(list(1 ,1 ,6 ,36 ,88 ,3 ,1 ,8 ,56 ,94 ,3 ,1 ,8 ,48 ,90 ,3 ,1 ,7 ,32 ,73 ,1 ,1 ,5 ,44 ,68 ,2 ,1 ,7 ,39 ,80 ,3 ,1 ,8 ,34 ,86 ,3 ,1 ,9 ,41 ,86 ,2 ,1 ,9 ,50 ,91 ,1 ,1 ,3 ,39 ,79 ,1 ,1 ,9 ,62 ,96 ,3 ,1 ,7 ,52 ,92 ,3 ,1 ,9 ,37 ,72 ,3 ,1 ,8 ,50 ,96 ,1 ,1 ,6 ,41 ,70 ,3 ,1 ,7 ,55 ,86 ,3 ,1 ,8 ,41 ,87 ,3 ,1 ,9 ,56 ,88 ,2 ,1 ,7 ,39 ,79 ,2 ,1 ,6 ,52 ,90 ,1 ,1 ,8 ,46 ,95 ,1 ,1 ,7 ,44 ,85 ,3 ,1 ,8 ,41 ,90 ,1 ,1 ,9 ,50 ,115 ,1 ,1 ,9 ,50 ,84 ,3 ,1 ,7 ,44 ,79 ,1 ,1 ,4 ,52 ,94 ,2 ,1 ,7 ,54 ,97 ,3 ,1 ,7 ,44 ,86 ,2 ,1 ,9 ,52 ,111 ,2 ,1 ,7 ,37 ,87 ,3 ,1 ,9 ,52 ,98 ,1 ,1 ,10 ,50 ,87 ,3 ,1 ,5 ,36 ,68 ,3 ,1 ,6 ,50 ,88 ,1 ,1 ,9 ,52 ,82 ,2 ,1 ,9 ,55 ,111 ,1 ,1 ,8 ,31 ,75 ,1 ,1 ,6 ,36 ,94 ,2 ,1 ,6 ,49 ,95 ,1 ,1 ,5 ,42 ,80 ,1 ,1 ,8 ,37 ,95 ,3 ,1 ,8 ,41 ,68 ,2 ,1 ,5 ,30 ,94 ,2 ,1 ,6 ,52 ,88 ,1 ,1 ,9 ,30 ,84 ,1 ,1 ,4 ,44 ,101 ,1 ,1 ,8 ,66 ,98 ,2 ,1 ,9 ,48 ,78 ,3 ,1 ,7 ,43 ,109 ,3 ,1 ,7 ,57 ,102 ,3 ,1 ,6 ,46 ,81 ,3 ,2 ,9 ,54 ,97 ,2 ,2 ,9 ,48 ,75 ,2 ,2 ,8 ,48 ,97 ,3 ,2 ,6 ,62 ,101 ,3 ,2 ,10 ,58 ,101 ,2 ,2 ,8 ,58 ,95 ,2 ,2 ,7 ,62 ,95 ,3 ,2 ,8 ,46 ,95 ,2 ,2 ,3 ,34 ,90 ,3 ,2 ,8 ,66 ,107 ,3 ,2 ,10 ,52 ,92 ,2 ,2 ,7 ,55 ,86 ,2 ,2 ,5 ,55 ,70 ,3 ,2 ,10 ,57 ,95 ,3 ,2 ,5 ,56 ,96 ,2 ,2 ,8 ,55 ,91 ,1 ,2 ,9 ,56 ,87 ,2 ,2 ,6 ,54 ,92 ,2 ,2 ,9 ,55 ,97 ,2 ,2 ,8 ,46 ,102 ,2 ,2 ,5 ,52 ,91 ,3 ,2 ,8 ,32 ,68 ,2 ,2 ,3 ,44 ,88 ,3 ,2 ,7 ,46 ,97 ,1 ,2 ,8 ,59 ,90 ,3 ,2 ,10 ,46 ,101 ,3 ,2 ,9 ,46 ,94 ,1 ,2 ,10 ,54 ,101 ,1 ,2 ,9 ,66 ,109 ,2 ,2 ,8 ,56 ,100 ,1 ,2 ,8 ,59 ,103 ,3 ,2 ,8 ,57 ,94 ,1 ,2 ,9 ,52 ,97 ,3 ,2 ,4 ,48 ,85 ,2 ,2 ,6 ,44 ,75 ,3 ,2 ,7 ,41 ,77 ,1 ,2 ,4 ,50 ,87 ,1 ,2 ,9 ,48 ,78 ,2 ,2 ,7 ,48 ,108 ,3 ,2 ,8 ,59 ,97 ,3 ,2 ,8 ,46 ,106 ,1 ,2 ,7 ,54 ,107 ,2 ,2 ,7 ,55 ,95 ,2 ,2 ,9 ,54 ,107 ,3 ,2 ,8 ,59 ,115 ,3 ,2 ,8 ,44 ,101 ,1 ,2 ,9 ,54 ,85 ,3 ,2 ,9 ,52 ,90 ,3 ,2 ,10 ,66 ,115 ,3 ,2 ,7 ,44 ,95 ,2 ,2 ,8 ,57 ,97 ,3 ,2 ,5 ,39 ,112 ,3 ,2 ,9 ,60 ,97 ,3 ,2 ,8 ,45 ,77 ,3 ,2 ,7 ,41 ,90 ,1 ,2 ,8 ,50 ,94 ,3 ,2 ,8 ,39 ,103 ,1 ,2 ,7 ,43 ,77 ,3 ,2 ,6 ,48 ,98 ,3 ,2 ,7 ,37 ,90 ,3 ,2 ,7 ,58 ,111 ,2 ,2 ,6 ,46 ,77 ,3 ,2 ,6 ,43 ,88 ,2 ,2 ,7 ,44 ,75 ,3 ,2 ,9 ,34 ,92 ,2 ,2 ,6 ,30 ,78 ,2 ,2 ,10 ,50 ,106 ,2 ,2 ,4 ,39 ,80 ,3 ,2 ,8 ,37 ,87 ,3 ,2 ,7 ,55 ,92 ,3 ,2 ,5 ,39 ,86 ,3 ,2 ,9 ,36 ,85 ,2 ,2 ,8 ,43 ,90 ,1 ,2 ,9 ,50 ,101 ,3 ,2 ,8 ,55 ,94 ,3 ,2 ,8 ,43 ,86 ,3 ,3 ,9 ,60 ,86 ,2 ,3 ,8 ,48 ,90 ,1 ,3 ,9 ,30 ,75 ,3 ,3 ,7 ,43 ,86 ,3 ,3 ,6 ,39 ,91 ,3 ,3 ,8 ,52 ,97 ,2 ,3 ,6 ,39 ,91 ,3 ,3 ,5 ,39 ,70 ,2 ,3 ,3 ,56 ,98 ,3 ,3 ,6 ,59 ,96 ,3 ,3 ,8 ,46 ,95 ,2 ,3 ,7 ,57 ,100 ,3 ,3 ,8 ,50 ,95 ,3 ,3 ,6 ,54 ,97 ,3 ,3 ,9 ,50 ,97 ,3 ,3 ,9 ,60 ,92 ,1 ,3 ,10 ,59 ,115 ,3 ,3 ,7 ,41 ,88 ,2 ,3 ,5 ,48 ,87 ,3 ,3 ,8 ,59 ,100 ,3 ,3 ,9 ,60 ,98 ,1 ,3 ,8 ,56 ,102 ,3 ,3 ,4 ,51 ,96) ,dim=c(5 ,151) ,dimnames=list(c('MOMAGE' ,'MVRBIQ0' ,'MWARM30' ,'MC30VRB' ,'WISCRY7V') ,1:151)) y <- array(NA,dim=c(5,151),dimnames=list(c('MOMAGE','MVRBIQ0','MWARM30','MC30VRB','WISCRY7V'),1:151)) for (i in 1:dim(x)[1]) { for (j in 1:dim(x)[2]) { y[i,j] <- as.numeric(x[i,j]) } } par3 = 'TRUE' par2 = '2' par1 = '4' ylab = 'Y Variable Name' xlab = 'X Variable Name' main = 'Title Goes Here' #'GNU S' R Code compiled by R2WASP v. 1.0.44 () #Author: Dr. Ian E. Holliday #To cite this work: Ian E. Holliday, 2009, STARS Bullying Study (v1.0.1) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/Ian.Holliday/rwasp_STARS_Bullying_Study_alt.wasp/ #Source of accompanying publication: #Technical description: cat1 <- as.numeric(par1) # cat2<- as.numeric(par2) # intercept<-as.logical(par3) x <- t(x) x1<-as.numeric(x[,cat1]) f1<-as.character(x[,cat2]) xdf<-data.frame(x1,f1) (V1<-dimnames(y)[[1]][cat1]) (V2<-dimnames(y)[[1]][cat2]) names(xdf)<-c('Response', 'Treatment') if(intercept == FALSE) (lmxdf<-lm(Response ~ Treatment - 1, data = xdf) ) else (lmxdf<-lm(Response ~ Treatment, data = xdf) ) (aov.xdf<-aov(lmxdf) ) (anova.xdf<-anova(lmxdf) ) #Note: the /var/www/html/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab load(file="/var/www/html/rcomp/createtable") a<-table.start() a<-table.row.start(a) a<-table.element(a,'ANOVA Model', length(lmxdf$coefficients)+1,TRUE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a, paste(V1, ' ~ ', V2), length(lmxdf$coefficients)+1,TRUE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a, 'means',,TRUE) for(i in 1:length(lmxdf$coefficients)){ a<-table.element(a, round(lmxdf$coefficients[i], digits=3),,FALSE) } a<-table.row.end(a) a<-table.end(a) table.save(a,file="/var/www/html/rcomp/tmp/1hbiw1290686234.tab") a<-table.start() a<-table.row.start(a) a<-table.element(a,'ANOVA Statistics', 5+1,TRUE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a, ' ',,TRUE) a<-table.element(a, 'Df',,FALSE) a<-table.element(a, 'Sum Sq',,FALSE) a<-table.element(a, 'Mean Sq',,FALSE) a<-table.element(a, 'F value',,FALSE) a<-table.element(a, 'Pr(>F)',,FALSE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a, V2,,TRUE) a<-table.element(a, anova.xdf$Df[1],,FALSE) a<-table.element(a, round(anova.xdf$'Sum Sq'[1], digits=3),,FALSE) a<-table.element(a, round(anova.xdf$'Mean Sq'[1], digits=3),,FALSE) a<-table.element(a, round(anova.xdf$'F value'[1], digits=3),,FALSE) a<-table.element(a, round(anova.xdf$'Pr(>F)'[1], digits=3),,FALSE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a, 'Residuals',,TRUE) a<-table.element(a, anova.xdf$Df[2],,FALSE) a<-table.element(a, round(anova.xdf$'Sum Sq'[2], digits=3),,FALSE) a<-table.element(a, round(anova.xdf$'Mean Sq'[2], digits=3),,FALSE) a<-table.element(a, ' ',,FALSE) a<-table.element(a, ' ',,FALSE) a<-table.row.end(a) a<-table.end(a) table.save(a,file="/var/www/html/rcomp/tmp/2elyn1290686234.tab") postscript(file="/var/www/html/rcomp/tmp/3elyn1290686234.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) boxplot(Response ~ Treatment, data=xdf, xlab=V2, ylab=V1) dev.off() if(intercept==TRUE){ thsd<-TukeyHSD(aov.xdf) postscript(file="/var/www/html/rcomp/tmp/4elyn1290686234.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) plot(thsd) dev.off() a<-table.start() a<-table.row.start(a) a<-table.element(a,'Tukey Honest Significant Difference Comparisons', 5,TRUE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a, ' ', 1, TRUE) for(i in 1:4){ a<-table.element(a,colnames(thsd[[1]])[i], 1, TRUE) } a<-table.row.end(a) for(i in 1:length(rownames(thsd[[1]]))){ a<-table.row.start(a) a<-table.element(a,rownames(thsd[[1]])[i], 1, TRUE) for(j in 1:4){ a<-table.element(a,round(thsd[[1]][i,j], digits=3), 1, FALSE) } a<-table.row.end(a) } a<-table.end(a) table.save(a,file="/var/www/html/rcomp/tmp/5hlft1290686234.tab") } if(intercept==FALSE){ a<-table.start() a<-table.row.start(a) a<-table.element(a,'TukeyHSD Message', 1,TRUE) a<-table.row.end(a) a<-table.start() a<-table.row.start(a) a<-table.element(a,'Must Include Intercept to use Tukey Test ', 1, FALSE) a<-table.row.end(a) a<-table.end(a) table.save(a,file="/var/www/html/rcomp/tmp/6auwe1290686234.tab") } library(car) lt.lmxdf<-levene.test(lmxdf) a<-table.start() a<-table.row.start(a) a<-table.element(a,'Levenes Test for Homogeneity of Variance', 4,TRUE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,' ', 1, TRUE) for (i in 1:3){ a<-table.element(a,names(lt.lmxdf)[i], 1, FALSE) } a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Group', 1, TRUE) for (i in 1:3){ a<-table.element(a,round(lt.lmxdf[[i]][1], digits=3), 1, FALSE) } a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,' ', 1, TRUE) a<-table.element(a,lt.lmxdf[[1]][2], 1, FALSE) a<-table.element(a,' ', 1, FALSE) a<-table.element(a,' ', 1, FALSE) a<-table.row.end(a) a<-table.end(a) table.save(a,file="/var/www/html/rcomp/tmp/7dvc21290686234.tab") try(system("convert tmp/3elyn1290686234.ps tmp/3elyn1290686234.png",intern=TRUE)) try(system("convert tmp/4elyn1290686234.ps tmp/4elyn1290686234.png",intern=TRUE))