x <- c(127269 ,75620 ,88255 ,104300 ,53891 ,23637 ,239665 ,26847 ,82376 ,108063 ,130696 ,106111 ,77157 ,156938 ,111186 ,193991 ,137623 ,129655 ,90067 ,118206 ,128389 ,226536 ,67073 ,89947 ,101451 ,121691 ,147649 ,101001 ,88745 ,145509 ,158310 ,82934 ,144300 ,50968 ,169102 ,129402 ,197808 ,102540 ,76806 ,135028 ,158063 ,97980 ,99862 ,49417 ,79864 ,238236 ,66613 ,133589 ,179955 ,104416 ,88737 ,24019 ,178249 ,64672 ,85913 ,149157 ,136059 ,147863 ,157786 ,97096 ,194551 ,166067 ,107928 ,136136 ,105283 ,74153 ,95612 ,108093 ,118089 ,51000 ,97034 ,52125 ,96331 ,106760 ,82657 ,159858 ,161402 ,90287 ,74853 ,87661 ,85840 ,224167 ,107905 ,82275 ,80953 ,90496 ,67181 ,155920 ,124542 ,84443 ,93028 ,106264 ,84922 ,132246 ,145658 ,106443 ,108558 ,138044 ,136773 ,183992 ,114968 ,126632 ,80456 ,39948 ,108674 ,101041 ,114400 ,153398 ,57427 ,78870 ,127197 ,119023 ,20764 ,99305 ,61675 ,72554 ,168643 ,21054 ,151340 ,23175 ,157461 ,103303 ,48355 ,51536 ,38214 ,86725 ,159468 ,90604 ,183334 ,137626 ,204060 ,114015 ,51227 ,147342 ,210790 ,41227 ,198679 ,81901 ,130930 ,82730 ,50656 ,195126 ,115032 ,173260 ,144766 ,57297 ,112567 ,142346 ,0 ,14688 ,98 ,455 ,0 ,0 ,111354 ,148268 ,0 ,203 ,7199 ,46660 ,17547 ,73567 ,969 ,76925) par4 = 'Unknown' par3 = 'FALSE' par2 = 'grey' par1 = '' xlab = 'Total time spent in RFC ' main = 'Histogram' #'GNU S' R Code compiled by R2WASP v. 1.0.44 () #Author: Dr. Ian E. Holliday #To cite this work: Ian E. Holliday, 2009, YOUR SOFTWARE TITLE (vNUMBER) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_YOURPAGE.wasp/ #Source of accompanying publication: #Technical description: par1 <- as.numeric(par1) if (par3 == 'TRUE') par3 <- TRUE if (par3 == 'FALSE') par3 <- FALSE if (par4 == 'Unknown') par1 <- as.numeric(par1) if (par4 == 'Interval/Ratio') par1 <- as.numeric(par1) if (par4 == '3-point Likert') par1 <- c(1:3 - 0.5, 3.5) if (par4 == '4-point Likert') par1 <- c(1:4 - 0.5, 4.5) if (par4 == '5-point Likert') par1 <- c(1:5 - 0.5, 5.5) if (par4 == '6-point Likert') par1 <- c(1:6 - 0.5, 6.5) if (par4 == '7-point Likert') par1 <- c(1:7 - 0.5, 7.5) if (par4 == '8-point Likert') par1 <- c(1:8 - 0.5, 8.5) if (par4 == '9-point Likert') par1 <- c(1:9 - 0.5, 9.5) if (par4 == '10-point Likert') par1 <- c(1:10 - 0.5, 10.5) postscript(file="/var/wessaorg/rcomp/tmp/1pzek1321026986.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) if(is.numeric(x[1])) { if (is.na(par1)) { myhist<-hist(x,col=par2,main=main,xlab=xlab,right=par3) } else { if (par1 < 0) par1 <- 3 if (par1 > 50) par1 <- 50 myhist<-hist(x,breaks=par1,col=par2,main=main,xlab=xlab,right=par3) } } else { plot(mytab <- table(x),col=par2,main='Frequency Plot',xlab=xlab,ylab='Absolute Frequency') } dev.off() if(is.numeric(x[1])) { myhist n <- length(x) #Note: the /var/wessaorg/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab load(file="/var/wessaorg/rcomp/createtable") a<-table.start() a<-table.row.start(a) a<-table.element(a,hyperlink('http://www.xycoon.com/histogram.htm','Frequency Table (Histogram)',''),6,TRUE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Bins',header=TRUE) a<-table.element(a,'Midpoint',header=TRUE) a<-table.element(a,'Abs. Frequency',header=TRUE) a<-table.element(a,'Rel. Frequency',header=TRUE) a<-table.element(a,'Cumul. Rel. Freq.',header=TRUE) a<-table.element(a,'Density',header=TRUE) a<-table.row.end(a) crf <- 0 if (par3 == FALSE) mybracket <- '[' else mybracket <- ']' mynumrows <- (length(myhist$breaks)-1) for (i in 1:mynumrows) { a<-table.row.start(a) if (i == 1) dum <- paste('[',myhist$breaks[i],sep='') else dum <- paste(mybracket,myhist$breaks[i],sep='') dum <- paste(dum,myhist$breaks[i+1],sep=',') if (i==mynumrows) dum <- paste(dum,']',sep='') else dum <- paste(dum,mybracket,sep='') a<-table.element(a,dum,header=TRUE) a<-table.element(a,myhist$mids[i]) a<-table.element(a,myhist$counts[i]) rf <- myhist$counts[i]/n crf <- crf + rf a<-table.element(a,round(rf,6)) a<-table.element(a,round(crf,6)) a<-table.element(a,round(myhist$density[i],6)) a<-table.row.end(a) } a<-table.end(a) table.save(a,file="/var/wessaorg/rcomp/tmp/2yqe71321026986.tab") } else { mytab reltab <- mytab / sum(mytab) n <- length(mytab) #Note: the /var/wessaorg/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab load(file="/var/wessaorg/rcomp/createtable") a<-table.start() a<-table.row.start(a) a<-table.element(a,'Frequency Table (Categorical Data)',3,TRUE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Category',header=TRUE) a<-table.element(a,'Abs. Frequency',header=TRUE) a<-table.element(a,'Rel. Frequency',header=TRUE) a<-table.row.end(a) for (i in 1:n) { a<-table.row.start(a) a<-table.element(a,labels(mytab)$x[i],header=TRUE) a<-table.element(a,mytab[i]) a<-table.element(a,round(reltab[i],4)) a<-table.row.end(a) } a<-table.end(a) table.save(a,file="/var/wessaorg/rcomp/tmp/3och11321026986.tab") } try(system("convert tmp/1pzek1321026986.ps tmp/1pzek1321026986.png",intern=TRUE))