R version 2.9.0 (2009-04-17)
Copyright (C) 2009 The R Foundation for Statistical Computing
ISBN 3-900051-07-0

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 <- c(14861,14583.3,15305.8,17903.9,16379.4,15420.3,17870.5,15912.8,13866.5,17823.2,17872,17420.4,16704.4,15991.2,16583.6,19123.5,17838.7,17209.4,18586.5,16258.1,15141.6,19202.1,17746.5,19090.1,18040.3,17515.5,17751.8,21072.4,17170,19439.5,19795.4,17574.9,16165.4,19464.6,19932.1,19961.2,17343.4,18924.2,18574.1,21350.6,18594.6,19823.1,20844.4,19640.2,17735.4,19813.6,22160,20664.3,17877.4,20906.5,21164.1,21374.4,22952.3,21343.5,23899.3,22392.9,18274.1,22786.7,22321.5,17842.2,16373.5,16087.1,16555.9,17880.2,16764.5,16049,18288.3,17570.4,15133.4,19334.2,19291.8,20176.7)
> par4 = 'Unknown'
> par3 = 'FALSE'
> par2 = 'grey'
> par1 = '50'
> xlab = 'x'
> main = 'Histogram'
> #'GNU S' R Code compiled by R2WASP v. 1.0.44 ()
> #Author: Prof. Dr. P. Wessa
> #To cite this work: AUTHOR(S), (YEAR), 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: Office for Research, Development, and Education
> #Technical description: Write here your technical program description (don't use hard returns!)
> 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/www/html/rcomp/tmp/1dguz1266861824.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
> 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)
+ }
> dev.off()
null device 
          1 
> myhist
$breaks
 [1] 13800 14000 14200 14400 14600 14800 15000 15200 15400 15600 15800 16000
[13] 16200 16400 16600 16800 17000 17200 17400 17600 17800 18000 18200 18400
[25] 18600 18800 19000 19200 19400 19600 19800 20000 20200 20400 20600 20800
[37] 21000 21200 21400 21600 21800 22000 22200 22400 22600 22800 23000 23200
[49] 23400 23600 23800 24000

$counts
 [1] 1 0 0 1 0 1 2 1 1 0 2 3 3 2 2 0 1 2 4 3 8 1 2 3 0 1 2 3 2 2 4 1 0 0 1 2 2 3
[39] 0 0 0 1 2 0 1 1 0 0 0 0 1

$intensities
 [1] 6.944444e-05 0.000000e+00 0.000000e+00 6.944444e-05 0.000000e+00
 [6] 6.944444e-05 1.388889e-04 6.944444e-05 6.944444e-05 0.000000e+00
[11] 1.388889e-04 2.083333e-04 2.083333e-04 1.388889e-04 1.388889e-04
[16] 0.000000e+00 6.944444e-05 1.388889e-04 2.777778e-04 2.083333e-04
[21] 5.555556e-04 6.944444e-05 1.388889e-04 2.083333e-04 0.000000e+00
[26] 6.944444e-05 1.388889e-04 2.083333e-04 1.388889e-04 1.388889e-04
[31] 2.777778e-04 6.944444e-05 0.000000e+00 0.000000e+00 6.944444e-05
[36] 1.388889e-04 1.388889e-04 2.083333e-04 0.000000e+00 0.000000e+00
[41] 0.000000e+00 6.944444e-05 1.388889e-04 0.000000e+00 6.944444e-05
[46] 6.944444e-05 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00
[51] 6.944443e-05

$density
 [1] 6.944444e-05 0.000000e+00 0.000000e+00 6.944444e-05 0.000000e+00
 [6] 6.944444e-05 1.388889e-04 6.944444e-05 6.944444e-05 0.000000e+00
[11] 1.388889e-04 2.083333e-04 2.083333e-04 1.388889e-04 1.388889e-04
[16] 0.000000e+00 6.944444e-05 1.388889e-04 2.777778e-04 2.083333e-04
[21] 5.555556e-04 6.944444e-05 1.388889e-04 2.083333e-04 0.000000e+00
[26] 6.944444e-05 1.388889e-04 2.083333e-04 1.388889e-04 1.388889e-04
[31] 2.777778e-04 6.944444e-05 0.000000e+00 0.000000e+00 6.944444e-05
[36] 1.388889e-04 1.388889e-04 2.083333e-04 0.000000e+00 0.000000e+00
[41] 0.000000e+00 6.944444e-05 1.388889e-04 0.000000e+00 6.944444e-05
[46] 6.944444e-05 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00
[51] 6.944443e-05

$mids
 [1] 13900 14100 14300 14500 14700 14900 15100 15300 15500 15700 15900 16100
[13] 16300 16500 16700 16900 17100 17300 17500 17700 17900 18100 18300 18500
[25] 18700 18900 19100 19300 19500 19700 19900 20100 20300 20500 20700 20900
[37] 21100 21300 21500 21700 21900 22100 22300 22500 22700 22900 23100 23300
[49] 23500 23700 23900

$xname
[1] "x"

$equidist
[1] TRUE

attr(,"class")
[1] "histogram"
> n <- length(x)
> 
> #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,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/www/html/rcomp/tmp/2p6bo1266861824.tab") 
> 
> try(system("convert tmp/1dguz1266861824.ps tmp/1dguz1266861824.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  0.414   0.190   0.536