x <- c(255 ,280.2 ,299.9 ,339.2 ,374.2 ,393.5 ,389.2 ,381.7 ,375.2 ,369 ,357.4 ,352.1 ,346.5 ,342.9 ,340.3 ,328.3 ,322.9 ,314.3 ,308.9 ,294 ,285.6 ,281.2 ,280.3 ,278.8 ,274.5 ,270.4 ,263.4 ,259.9 ,258 ,262.7 ,284.7 ,311.3 ,322.1 ,327 ,331.3 ,333.3 ,321.4 ,327 ,320 ,314.7 ,316.7 ,314.4 ,321.3 ,318.2 ,307.2 ,301.3 ,287.5 ,277.7 ,274.4 ,258.8 ,253.3 ,251 ,248.4 ,249.5 ,246.1 ,244.5 ,243.6 ,244 ,240.8 ,249.8 ,248 ,259.4 ,260.5 ,260.8 ,261.3 ,259.5 ,256.6 ,257.9 ,256.5 ,254.2 ,253.3 ,253.8 ,255.5 ,257.1 ,257.3 ,253.2 ,252.8 ,252 ,250.7 ,252.2 ,250 ,251 ,253.4 ,251.2 ,255.6 ,261.1 ,258.9 ,259.9 ,261.2 ,264.7 ,267.1 ,266.4 ,267.7 ,268.6 ,267.5 ,268.5 ,268.5 ,270.5 ,270.9 ,270.1 ,269.3 ,269.8 ,270.1 ,264.9 ,263.7 ,264.8 ,263.7 ,255.9 ,276.2 ,360.1 ,380.5 ,373.7 ,369.8 ,366.6 ,359.3 ,345.8 ,326.2 ,324.5 ,328.1 ,327.5 ,324.4 ,316.5 ,310.9 ,301.5 ,291.7 ,290.4 ,287.4 ,277.7 ,281.6 ,288 ,276 ,272.9 ,283 ,283.3 ,276.8 ,284.5 ,282.7 ,281.2 ,287.4 ,283.1 ,284 ,285.5 ,289.2 ,292.5 ,296.4 ,305.2 ,303.9 ,311.5 ,316.3 ,316.7 ,322.5 ,317.1 ,309.8 ,303.8 ,290.3 ,293.7 ,291.7 ,296.5 ,289.1 ,288.5 ,293.8 ,297.7 ,305.4 ,302.7 ,302.5 ,303 ,294.5 ,294.1 ,294.5 ,297.1 ,289.4 ,292.4 ,287.9 ,286.6 ,280.5 ,272.4 ,269.2 ,270.6 ,267.3 ,262.5 ,266.8 ,268.8 ,263.1 ,261.2 ,266 ,262.5 ,265.2 ,261.3 ,253.7 ,249.2 ,239.1 ,236.4 ,235.2 ,245.2 ,246.2 ,247.7 ,251.4 ,253.3 ,254.8 ,250 ,249.3 ,241.5 ,243.3 ,248 ,253 ,252.9 ,251.5 ,251.6 ,253.5 ,259.8 ,334.1 ,448 ,445.8 ,445 ,448.2 ,438.2 ,439.8 ,423.4 ,410.8 ,408.4 ,406.7 ,405.9 ,402.7 ,405.1 ,399.6 ,386.5 ,381.4 ,375.2 ,357.7 ,359 ,355 ,352.7 ,344.4 ,343.8 ,338 ,339 ,333.3 ,334.4 ,328.3 ,330.7 ,330 ,331.6 ,351.2 ,389.4 ,410.9 ,442.8 ,462.8 ,466.9 ,461.7 ,439.2 ,430.3 ,416.1 ,402.5 ,397.3 ,403.3 ,395.9 ,387.8 ,378.6 ,377.1 ,370.4 ,362 ,350.3 ,348.2 ,344.6 ,343.5 ,342.8 ,347.6 ,346.6 ,349.5 ,342.1 ,342 ,342.8 ,339.3 ,348.2 ,333.7 ,334.7 ,354 ,367.7 ,363.3 ,358.4 ,353.1 ,343.1 ,344.6 ,344.4 ,333.9 ,331.7 ,324.3 ,321.2 ,322.4 ,321.7 ,320.5 ,312.8 ,309.7 ,315.6 ,309.7 ,304.6 ,302.5 ,301.5 ,298.8 ,291.3 ,293.6 ,294.6 ,285.9 ,297.6 ,301.1 ,293.8 ,297.7 ,292.9 ,292.1 ,287.2 ,288.2 ,283.8 ,299.9 ,292.4 ,293.3 ,300.8 ,293.7 ,293.1 ,294.4 ,292.1 ,291.9 ,282.5 ,277.9 ,287.5 ,289.2 ,285.6 ,293.2 ,290.8 ,283.1 ,275 ,287.8 ,287.8 ,287.4 ,284 ,277.8 ,277.6 ,304.9 ,294 ,300.9 ,324 ,332.9 ,341.6 ,333.4 ,348.2 ,344.7 ,344.7 ,329.3 ,323.5 ,323.2 ,317.4 ,330.1 ,329.2 ,334.9 ,315.8 ,315.4 ,319.6 ,317.3 ,313.8 ,315.8 ,311.3) par4 = 'Unknown' par3 = 'TRUE' par2 = 'red' par1 = '' xlab = 'x' 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/1ii811321125843.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/23l2f1321125843.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/3xukz1321125843.tab") } try(system("convert tmp/1ii811321125843.ps tmp/1ii811321125843.png",intern=TRUE))