R version 3.0.2 (2013-09-25) -- "Frisbee Sailing" Copyright (C) 2013 The R Foundation for Statistical Computing Platform: i686-pc-linux-gnu (32-bit) R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type 'license()' or 'licence()' for distribution details. R is a collaborative project with many contributors. Type 'contributors()' for more information and 'citation()' on how to cite R or R packages in publications. Type 'demo()' for some demos, 'help()' for on-line help, or 'help.start()' for an HTML browser interface to help. Type 'q()' to quit R. > x <- array(list(58 + ,62 + ,64 + ,63 + ,60 + ,58 + ,62 + ,64 + ,63 + ,60 + ,58 + ,62 + ,64 + ,63 + ,60 + ,58 + ,62 + ,64 + ,63 + ,60 + ,58 + ,68 + ,66 + ,64 + ,63 + ,59 + ,70 + ,68 + ,67 + ,63 + ,62 + ,71 + ,69 + ,68 + ,66 + ,63 + ,71 + ,70 + ,69 + ,67 + ,68 + ,72 + ,70 + ,69 + ,67 + ,68 + ,72 + ,71 + ,69 + ,69 + ,70 + ,72 + ,72 + ,70 + ,70 + ,70 + ,72 + ,72 + ,70 + ,70 + ,71 + ,72 + ,73 + ,72 + ,71 + ,71 + ,72 + ,74 + ,72 + ,71 + ,73 + ,72 + ,74 + ,73 + ,73 + ,73 + ,72 + ,75 + ,73 + ,73 + ,74 + ,72 + ,77 + ,75 + ,74 + ,75 + ,72 + ,77 + ,77 + ,75 + ,75 + ,73 + ,77 + ,77 + ,75 + ,75 + ,74 + ,78 + ,77 + ,75 + ,75 + ,74 + ,78 + ,78 + ,77 + ,77 + ,75 + ,78 + ,79 + ,78 + ,77 + ,75 + ,78 + ,79 + ,78 + ,77 + ,75 + ,78 + ,80 + ,79 + ,77 + ,75 + ,78 + ,80 + ,79 + ,78 + ,75 + ,78 + ,81 + ,79 + ,78 + ,75 + ,79 + ,81 + ,79 + ,78 + ,75 + ,79 + ,81 + ,79 + ,80 + ,77 + ,79 + ,82 + ,80 + ,80 + ,77 + ,79 + ,82 + ,81 + ,80 + ,78 + ,80 + ,82 + ,81 + ,80 + ,78 + ,81 + ,82 + ,81 + ,80 + ,78 + ,81 + ,82 + ,82 + ,80 + ,78 + ,81 + ,83 + ,82 + ,80 + ,79 + ,82 + ,83 + ,82 + ,80 + ,79 + ,82 + ,83 + ,82 + ,80 + ,79 + ,82 + ,83 + ,82 + ,81 + ,80 + ,82 + ,83 + ,82 + ,81 + ,80 + ,83 + ,83 + ,82 + ,81 + ,81 + ,83 + ,83 + ,82 + ,82 + ,81 + ,83 + ,83 + ,82 + ,82 + ,81 + ,83 + ,84 + ,83 + ,82 + ,82 + ,84 + ,84 + ,83 + ,82 + ,82 + ,84 + ,84 + ,83 + ,82 + ,82 + ,84 + ,84 + ,84 + ,82 + ,82 + ,84 + ,85 + ,84 + ,83 + ,82 + ,84 + ,86 + ,84 + ,83 + ,82 + ,84 + ,86 + ,84 + ,84 + ,82 + ,84 + ,86 + ,85 + ,84 + ,82 + ,85 + ,86 + ,85 + ,84 + ,82 + ,86 + ,87 + ,85 + ,84 + ,82 + ,86 + ,87 + ,85 + ,84 + ,84 + ,86 + ,87 + ,85 + ,84 + ,85 + ,86 + ,88 + ,85 + ,84 + ,85 + ,86 + ,88 + ,85 + ,84 + ,85 + ,86 + ,88 + ,85 + ,85 + ,86 + ,87 + ,88 + ,85 + ,85 + ,87 + ,87 + ,88 + ,85 + ,86 + ,87 + ,87 + ,88 + ,85 + ,86 + ,87 + ,87 + ,88 + ,86 + ,86 + ,87 + ,87 + ,88 + ,86 + ,86 + ,88 + ,87 + ,88 + ,86 + ,86 + ,88 + ,87 + ,89 + ,86 + ,86 + ,88 + ,88 + ,89 + ,87 + ,86 + ,88 + ,88 + ,89 + ,87 + ,86 + ,88 + ,88 + ,89 + ,87 + ,86 + ,88 + ,88 + ,89 + ,87 + ,86 + ,89 + ,88 + ,89 + ,87 + ,86 + ,89 + ,89 + ,89 + ,87 + ,86 + ,89 + ,89 + ,89 + ,87 + ,86 + ,90 + ,89 + ,90 + ,88 + ,87 + ,90 + ,90 + ,90 + ,88 + ,87 + ,90 + ,90 + ,90 + ,88 + ,87 + ,90 + ,90 + ,91 + ,88 + ,87 + ,90 + ,90 + ,91 + ,88 + ,87 + ,91 + ,90 + ,91 + ,89 + ,87 + ,91 + ,90 + ,91 + ,89 + ,88 + ,91 + ,91 + ,91 + ,89 + ,88 + ,91 + ,91 + ,91 + ,89 + ,88 + ,91 + ,92 + ,92 + ,90 + ,88 + ,91 + ,92 + ,92 + ,90 + ,89 + ,91 + ,92 + ,92 + ,90 + ,89 + ,91 + ,92 + ,92 + ,91 + ,89 + ,91 + ,93 + ,92 + ,91 + ,89 + ,92 + ,94 + ,92 + ,92 + ,89 + ,92 + ,95 + ,93 + ,92 + ,89 + ,92 + ,95 + ,93 + ,92 + ,89 + ,92 + ,95 + ,93 + ,92 + ,89 + ,92 + ,95 + ,93 + ,92 + ,90 + ,92 + ,95 + ,94 + ,92 + ,90 + ,92 + ,95 + ,94 + ,92 + ,90 + ,93 + ,95 + ,94 + ,92 + ,90 + ,93 + ,96 + ,94 + ,93 + ,90 + ,93 + ,96 + ,94 + ,93 + ,90 + ,93 + ,96 + ,95 + ,93 + ,90 + ,93 + ,96 + ,95 + ,93 + ,90 + ,93 + ,96 + ,95 + ,93 + ,90 + ,95 + ,97 + ,95 + ,94 + ,91 + ,95 + ,97 + ,96 + ,94 + ,91 + ,95 + ,97 + ,96 + ,94 + ,91 + ,95 + ,98 + ,96 + ,95 + ,91 + ,95 + ,98 + ,96 + ,95 + ,91 + ,95 + ,98 + ,96 + ,95 + ,92 + ,95 + ,98 + ,97 + ,95 + ,92 + ,95 + ,99 + ,97 + ,95 + ,92 + ,95 + ,99 + ,97 + ,95 + ,92 + ,95 + ,99 + ,97 + ,95 + ,92 + ,95 + ,99 + ,97 + ,95 + ,92 + ,95 + ,100 + ,97 + ,97 + ,92 + ,96 + ,100 + ,97 + ,97 + ,94 + ,97 + ,100 + ,98 + ,97 + ,94 + ,97 + ,101 + ,98 + ,97 + ,94 + ,97 + ,102 + ,98 + ,97 + ,95 + ,97 + ,102 + ,98 + ,97 + ,95 + ,97 + ,102 + ,98 + ,97 + ,95 + ,97 + ,102 + ,98 + ,97 + ,96 + ,97 + ,103 + ,98 + ,97 + ,96 + ,97 + ,103 + ,98 + ,98 + ,97 + ,97 + ,103 + ,99 + ,99 + ,97 + ,97 + ,104 + ,99 + ,99 + ,97 + ,98 + ,104 + ,99 + ,99 + ,98 + ,98 + ,105 + ,99 + ,99 + ,98 + ,99 + ,105 + ,100 + ,100 + ,98 + ,100 + ,105 + ,100 + ,100 + ,98 + ,100 + ,106 + ,100 + ,100 + ,98 + ,100 + ,106 + ,100 + ,100 + ,98 + ,102 + ,106 + ,100 + ,101 + ,98 + ,102 + ,106 + ,101 + ,102 + ,98 + ,102 + ,108 + ,101 + ,103 + ,98 + ,102 + ,108 + ,101 + ,104 + ,98 + ,102 + ,108 + ,101 + ,104 + ,98 + ,102 + ,108 + ,102 + ,104 + ,98 + ,102 + ,109 + ,102 + ,104 + ,99 + ,102 + ,109 + ,102 + ,105 + ,99 + ,102 + ,109 + ,102 + ,105 + ,99 + ,104 + ,109 + ,102 + ,105 + ,99 + ,104 + ,109 + ,102 + ,106 + ,102 + ,104 + ,109 + ,102 + ,107 + ,102 + ,104 + ,109 + ,103 + ,108 + ,102 + ,105 + ,109 + ,103 + ,109 + ,102 + ,105 + ,110 + ,103 + ,109 + ,102 + ,105 + ,110 + ,103 + ,109 + ,102 + ,105 + ,111 + ,104 + ,109 + ,104 + ,105 + ,111 + ,104 + ,109 + ,104 + ,105 + ,112 + ,105 + ,110 + ,104 + ,106 + ,112 + ,106 + ,111 + ,105 + ,106 + ,112 + ,108 + ,112 + ,105 + ,108 + ,114 + ,108 + ,112 + ,106 + ,109 + ,114 + ,108 + ,113 + ,107 + ,109 + ,114 + ,108 + ,114 + ,110 + ,110 + ,116 + ,108 + ,114 + ,112 + ,110 + ,116 + ,108 + ,114 + ,112 + ,110 + ,117 + ,109 + ,114 + ,112 + ,110 + ,118 + ,110 + ,117 + ,113 + ,111 + ,120 + ,110 + ,118 + ,113 + ,113 + ,121 + ,113 + ,118 + ,113 + ,113 + ,121 + ,113 + ,118 + ,114 + ,113 + ,122 + ,114 + ,118 + ,117 + ,114 + ,123 + ,119 + ,121 + ,117 + ,114 + ,123 + ,119 + ,121 + ,117 + ,114 + ,123 + ,119 + ,121 + ,117 + ,114 + ,123 + ,119 + ,121) + ,dim=c(5 + ,162) + ,dimnames=list(c('WJ7ARD' + ,'WJ7AMA' + ,'WJ7AWL' + ,'WJ7AKN' + ,'WJ7ASK') + ,1:162)) > y <- array(NA,dim=c(5,162),dimnames=list(c('WJ7ARD','WJ7AMA','WJ7AWL','WJ7AKN','WJ7ASK'),1:162)) > for (i in 1:dim(x)[1]) + { + for (j in 1:dim(x)[2]) + { + y[i,j] <- as.numeric(x[i,j]) + } + } > par5 = '2' > par4 = '1' > par3 = '5' > par2 = 'TRUE' > par1 = 'grey' > ylab = 'Test Scores' > xlab = 'CARE Variable Names' > main = 'CARE Age 7 Data' > par5 <- '2' > par4 <- '1' > par3 <- '5' > par2 <- 'TRUE' > par1 <- 'grey' > #'GNU S' R Code compiled by R2WASP v. 1.2.291 () > #Author: aston2 > #To cite this work: Ian E. Holliday, 2012, Histogram and QQplot (v1.0.6) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/Ian.Holliday/rwasp_Reddy-Moores%20Data%20Boxplot.wasp/ > #Source of accompanying publication: > # > postscript(file="/var/wessaorg/rcomp/tmp/1nih51382291855.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > par2<- as.logical(par2) # notches > par3<-as.numeric(par3) # percentage trim > if(par3>45){par3<-45;warning('trim limited to 45%')} > if(par3<0){par3<-0;warning('negative trim makes no sense. Trim is zero.')} > par4 <- as.numeric(par4) #factor column > par5 <- as.numeric(par5) # response column > x <- t(x) > x1<-as.numeric(x[,par5]) # response > f1<-as.character(x[,par4]) # factor > x2<-x1[f1=='no'] > f2 <- f1[f1=='no'] > lotrm<-as.integer(length(x2)*par3/100) > hitrm<-as.integer(length(x2)*(100-par3)/100) > srt<-order(x2,f2) > trmx1<-x2[srt[lotrm:hitrm]] > trmf1<-f2[srt[lotrm:hitrm]] > x3<-x1[f1=='yes'] > f3 <- f1[f1=='yes'] > lotrm<-as.integer(length(x3)*par3/100) > hitrm<-as.integer(length(x3)*(100-par3)/100) > srt<-order(x3,f3) > trmx2<-x3[srt[lotrm:hitrm]] > trmf2<-f3[srt[lotrm:hitrm]] > xtrm<-c(trmx1,trmx2) > ftrm<-c(trmf1,trmf2) > xtrm[1:6] [1] NA NA NA NA NA NA > ftrm[1:6] [1] NA NA NA NA NA NA > r<-boxplot(xtrm~as.factor(ftrm), col=par1, notch=par2, main=main, xlab=xlab, ylab=ylab) Error in boxplot.default(split(mf[[response]], mf[-response]), ...) : invalid first argument Calls: boxplot -> boxplot.formula -> boxplot -> boxplot.default In addition: Warning message: In is.na(rows) : is.na() applied to non-(list or vector) of type 'NULL' Execution halted