R version 2.8.0 (2008-10-20) Copyright (C) 2008 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(284.4,212.8,226.9,308.4,262,227.9,236.1,320.4,271.9,232.8,237,313.4,261.4,226.8,249.9,314.3,286.1,226.5,260.4,311.4,294.7,232.6,257.2,339.2,279.1,249.8,269.8,345.7,293.8,254.7,277.5,363.4,313.4,272.8,300.1,369.5,330.8,287.8,305.9,386.1,335.2,288,308.3,402.3,352.8,316.1,324.9,404.8,393,318.9,327,442.3,383.1,331.6,361.4,445.9,386.6,357.2,373.6,466.2,409.6,369.8,378.6,487,419.2,376.7,392.8,506.1,458.4,387.4,426.9,565,464.8,444.5,449.5,556.1,499.6,451.9,434.9,553.8,510,432.9,453.2,547.6,485.8,452.6,456.6,565.7,514.8,464.3,430.9,588.3,503.1,442.6,448,554.5,504.5,427.3,473.1,526.2,547.5,440.2,468.7,574.5,492.6,432.6,479.8,575.7,474.6,405.3,434.6,535.1,452.6,429.5,417.2,551.8,464,416.6,422.9,553.6,458.6,427.6,429.2,534.2,481.7,416,440.2,538.7,473.8,439.9,446.8,597.5,467.2,439.4,447.4,568.5,485.9,442.1,430.5,600,464.5,423.6,437,574,443,410,420,532,432,420,411,512,449,382) > par4 = 'Unknown' > par3 = 'FALSE' > par2 = 'grey' > par1 = '' > 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/1km911229626179.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] 200 250 300 350 400 450 500 550 600 $counts [1] 11 16 18 16 42 24 13 14 $intensities [1] 0.001428571 0.002077922 0.002337662 0.002077922 0.005454545 0.003116883 [7] 0.001688312 0.001818181 $density [1] 0.001428571 0.002077922 0.002337662 0.002077922 0.005454545 0.003116883 [7] 0.001688312 0.001818181 $mids [1] 225 275 325 375 425 475 525 575 $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/2y0ap1229626179.tab") > > system("convert tmp/1km911229626179.ps tmp/1km911229626179.png") > > > proc.time() user system elapsed 0.387 0.163 0.464