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(88 + ,95 + ,95 + ,81 + ,84 + ,98 + ,92 + ,87 + ,90 + ,97 + ,89 + ,95 + ,91 + ,91 + ,88 + ,81 + ,97 + ,84 + ,87 + ,87 + ,88 + ,92 + ,92 + ,89 + ,92 + ,95 + ,104 + ,84 + ,85 + ,91 + ,79 + ,79 + ,84 + ,81 + ,80 + ,91 + ,78 + ,69 + ,80 + ,73 + ,80 + ,79 + ,84 + ,75 + ,73 + ,91 + ,86 + ,84 + ,84 + ,77 + ,85 + ,92 + ,85 + ,91 + ,86 + ,94 + ,78 + ,89 + ,87 + ,82 + ,86 + ,86 + ,83 + ,90 + ,82 + ,82 + ,71 + ,86 + ,72 + ,73 + ,107 + ,90 + ,91 + ,96 + ,99 + ,97 + ,93 + ,90 + ,88 + ,90 + ,78 + ,80 + ,77 + ,73 + ,78 + ,91 + ,90 + ,88 + ,86 + ,100 + ,74 + ,77 + ,88 + ,67 + ,73 + ,95 + ,86 + ,88 + ,89 + ,95 + ,97 + ,89 + ,102 + ,94 + ,93 + ,112 + ,88 + ,103 + ,88 + ,91 + ,80 + ,86 + ,84 + ,77 + ,81 + ,98 + ,97 + ,87 + ,89 + ,98 + ,102 + ,81 + ,88 + ,104 + ,88 + ,100 + ,88 + ,93 + ,89 + ,90 + ,97 + ,82 + ,101 + ,107 + ,85 + ,92 + ,82 + ,88 + ,96 + ,79 + ,114 + ,100 + ,113 + ,111 + ,107 + ,88 + ,91 + ,95 + ,84 + ,82 + ,96 + ,111 + ,99 + ,83 + ,98 + ,91 + ,102 + ,85 + ,82 + ,84 + ,114 + ,105 + ,105 + ,105 + ,99 + ,91 + ,92 + ,105 + ,100 + ,89 + ,114 + ,100 + ,96 + ,111 + ,105 + ,89 + ,82 + ,88 + ,90 + ,94 + ,103 + ,91 + ,93 + ,87 + ,92 + ,83 + ,71 + ,93 + ,90 + ,81 + ,74 + ,71 + ,69 + ,67 + ,NA + ,88 + ,89 + ,83 + ,84 + ,86 + ,88 + ,82 + ,83 + ,84 + ,90 + ,98 + ,105 + ,93 + ,96 + ,101 + ,74 + ,71 + ,69 + ,67 + ,73 + ,86 + ,91 + ,86 + ,90 + ,88 + ,93 + ,105 + ,90 + ,83 + ,90 + ,85 + ,82 + ,87 + ,78 + ,84 + ,91 + ,86 + ,92 + ,100 + ,89 + ,75 + ,72 + ,75 + ,67 + ,73 + ,97 + ,98 + ,89 + ,94 + ,98 + ,88 + ,80 + ,88 + ,89 + ,82 + ,79 + ,84 + ,87 + ,89 + ,77 + ,106 + ,109 + ,113 + ,111 + ,107 + ,102 + ,80 + ,102 + ,88 + ,84 + ,98 + ,112 + ,101 + ,94 + ,93 + ,98 + ,86 + ,87 + ,111 + ,96 + ,114 + ,112 + ,105 + ,111 + ,107 + ,95 + ,90 + ,101 + ,105 + ,95 + ,98 + ,96 + ,79 + ,97 + ,89 + ,82 + ,84 + ,86 + ,83 + ,82 + ,97 + ,94 + ,95 + ,105 + ,92 + ,94 + ,112 + ,104 + ,69 + ,100 + ,84 + ,81 + ,79 + ,80 + ,75 + ,88 + ,105 + ,86 + ,87 + ,94 + ,83 + ,77 + ,81 + ,67 + ,83 + ,98 + ,105 + ,92 + ,92 + ,92 + ,106 + ,112 + ,95 + ,102 + ,102 + ,94 + ,91 + ,97 + ,90 + ,98 + ,91 + ,92 + ,87 + ,98 + ,94 + ,89 + ,72 + ,94 + ,94 + ,83 + ,102 + ,100 + ,105 + ,101 + ,107 + ,103 + ,103 + ,94 + ,96 + ,107 + ,96 + ,89 + ,88 + ,94 + ,91 + ,89 + ,86 + ,105 + ,105 + ,83 + ,112 + ,112 + ,113 + ,105 + ,100 + ,74 + ,71 + ,69 + ,84 + ,77 + ,74 + ,73 + ,69 + ,72 + ,74 + ,102 + ,99 + ,100 + ,96 + ,92 + ,106 + ,85 + ,113 + ,94 + ,100 + ,91 + ,80 + ,88 + ,107 + ,92 + ,94 + ,88 + ,104 + ,103 + ,93 + ,98 + ,95 + ,94 + ,95 + ,89 + ,85 + ,86 + ,91 + ,73 + ,86 + ,108 + ,92 + ,103 + ,110 + ,97 + ,90 + ,96 + ,95 + ,84 + ,82 + ,74 + ,71 + ,69 + ,67 + ,73 + ,81 + ,88 + ,86 + ,73 + ,93 + ,104 + ,98 + ,108 + ,84 + ,97 + ,90 + ,108 + ,83 + ,79 + ,93 + ,95 + ,103 + ,102 + ,91 + ,89 + ,85 + ,81 + ,94 + ,92 + ,86 + ,114 + ,96 + ,113 + ,107 + ,89 + ,114 + ,112 + ,113 + ,111 + ,107 + ,93 + ,95 + ,93 + ,89 + ,90 + ,104 + ,92 + ,105 + ,110 + ,104 + ,92 + ,103 + ,88 + ,98 + ,96 + ,100 + ,95 + ,90 + ,95 + ,93 + ,92 + ,94 + ,92 + ,84 + ,104 + ,77 + ,79 + ,81 + ,73 + ,83 + ,86 + ,91 + ,86 + ,86 + ,92 + ,87 + ,87 + ,85 + ,98 + ,83 + ,91 + ,82 + ,86 + ,90 + ,83 + ,102 + ,97 + ,97 + ,97 + ,96 + ,100 + ,92 + ,103 + ,106 + ,84 + ,113 + ,112 + ,102 + ,106 + ,107) + ,dim=c(5 + ,104) + ,dimnames=list(c('WJ10AFS' + ,'WJ10AVA' + ,'WJ10ARD' + ,'WJ10AMA' + ,'WJ10AKN') + ,1:104)) > y <- array(NA,dim=c(5,104),dimnames=list(c('WJ10AFS','WJ10AVA','WJ10ARD','WJ10AMA','WJ10AKN'),1:104)) > 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 = '20' > par2 = 'FALSE' > par1 = 'grey' > ylab = 'Degree Grade' > xlab = 'Placement Student' > main = 'Reddy-Moores Placement Data' > par5 <- '2' > par4 <- '1' > par3 <- '20' > par2 <- 'FALSE' > 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/fisher/rcomp/tmp/16vhg1382448824.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