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(196.9 + ,192.1 + ,201.8 + ,186.9 + ,218.0 + ,214.4 + ,227.5 + ,204.1 + ,225.8 + ,223.7 + ,244.7 + ,243.9 + ,257.3 + ,234.5 + ,251.4 + ,243.8 + ,247.4 + ,245.3 + ,262.5 + ,270.0 + ,259.9 + ,262.2 + ,244.9 + ,249.3 + ,268.2 + ,231.2 + ,264.3 + ,252.7 + ,275.5 + ,261.5 + ,275.5 + ,272.3 + ,268.6 + ,270.4 + ,267.7 + ,275.0 + ,272.6 + ,248.6 + ,279.4 + ,270.5 + ,292.8 + ,297.8 + ,296.8 + ,290.9 + ,282.8 + ,312.8 + ,303.2 + ,301.4 + ,289.8 + ,279.6 + ,302.2 + ,299.1 + ,319.7 + ,310.9 + ,315.2 + ,338.5 + ,315.6 + ,321.2 + ,318.5 + ,342.7 + ,261.4 + ,287.0 + ,331.5 + ,326.9 + ,338.6 + ,337.0 + ,358.4 + ,344.5 + ,345.7 + ,344.1 + ,317.4 + ,354.5 + ,345.2 + ,314.1 + ,352.5 + ,361.2 + ,365.9 + ,332.5 + ,364.0 + ,359.1 + ,345.6 + ,366.9 + ,370.2 + ,359.9 + ,366.6 + ,336.3 + ,368.5 + ,374.2 + ,384.3 + ,358.9 + ,407.7 + ,433.3 + ,404.7 + ,392.7 + ,409.7 + ,416.5 + ,414.3 + ,404.3 + ,421.4 + ,372.6 + ,404.7 + ,420.2 + ,438.4 + ,449.1 + ,445.8 + ,413.8 + ,420.5 + ,442.3 + ,438.9 + ,394.5 + ,416.8 + ,402.9 + ,424.5 + ,432.3 + ,484.1 + ,492.7 + ,496.3 + ,471.9 + ,491.2 + ,512.9 + ,482.4 + ,407.9 + ,448.5 + ,431.1 + ,498.8 + ,497.1 + ,517.1 + ,487.7 + ,512.5 + ,550.1 + ,532.5 + ,524.1 + ,515.7 + ,461.0 + ,529.3 + ,467.4 + ,559.8 + ,536.5 + ,531.9 + ,546.5 + ,547.4 + ,536.1 + ,482.8 + ,551.0 + ,532.9 + ,484.1 + ,554.8 + ,537.0 + ,558.0 + ,511.4 + ,502.9 + ,558.6 + ,545.1 + ,574.3 + ,542.2 + ,600.0 + ,588.6 + ,524.4 + ,618.5 + ,580.9 + ,557.2 + ,571.2 + ,597.5 + ,601.7 + ,558.9 + ,600.9 + ,601.0 + ,615.7 + ,578.1 + ,495.9 + ,526.8 + ,522.1 + ,605.1 + ,574.4 + ,609.7 + ,580.7 + ,565.1 + ,590.7 + ,571.5 + ,601.3 + ,567.3 + ,467.9 + ,588.9 + ,579.4 + ,502.6 + ,568.7 + ,616.0 + ,586.2 + ,575.5 + ,599.9 + ,568.2 + ,516.0 + ,493.4 + ,496.8 + ,529.9 + ,491.7 + ,543.2 + ,490.8 + ,554.7 + ,625.7 + ,605.0 + ,645.2 + ,645.2 + ,611.8 + ,600.3 + ,549.8 + ,635.5 + ,617.7 + ,643.5 + ,485.7 + ,689.5 + ,692.0 + ,677.3 + ,704.7 + ,668.6 + ,717.8 + ,689.8 + ,640.4 + ,675.2 + ,528.1 + ,538.0 + ,527.2 + ,655.6 + ,650.6 + ,623.7 + ,748.4 + ,727.4 + ,750.5 + ,678.9 + ,659.5 + ,691.9 + ,639.8 + ,663.8 + ,572.9 + ,592.5 + ,734.8 + ,696.1 + ,589.2 + ,662.9 + ,661.2 + ,672.1 + ,583.7 + ,705.5 + ,631.0 + ,733.3 + ,674.9 + ,695.5 + ,634.1 + ,630.6 + ,635.2 + ,554.1 + ,623.9 + ,679.3 + ,565.6 + ,564.1 + ,637.2 + ,650.8 + ,602.7 + ,587.5 + ,619.2 + ,616.5 + ,637.9 + ,557.9 + ,594.0 + ,668.7 + ,603.3 + ,674.5 + ,573.4 + ,706.0 + ,693.7 + ,627.5 + ,550.7 + ,592.3 + ,660.2 + ,597.3 + ,641.0 + ,663.6 + ,595.9 + ,638.4 + ,665.4 + ,671.4 + ,637.0 + ,685.7 + ,705.8 + ,704.8 + ,734.4 + ,674.2 + ,748.6 + ,763.4 + ,658.0 + ,627.5 + ,528.9 + ,488.3 + ,575.5 + ,735.6 + ,685.3 + ,613.6 + ,629.5 + ,634.7 + ,652.6 + ,728.3 + ,634.3 + ,690.7 + ,676.3 + ,675.4 + ,595.6 + ,712.4 + ,735.8 + ,544.4 + ,567.0 + ,510.0 + ,564.0 + ,630.7 + ,496.7 + ,660.9 + ,601.2 + ,655.2 + ,591.6 + ,606.1 + ,560.7 + ,368.3 + ,371.6 + ,413.9 + ,413.9 + ,389.0 + ,399.2 + ,429.8 + ,395.6 + ,472.0 + ,486.0 + ,525.0 + ,396.0 + ,511.0 + ,525.0 + ,492.0 + ,517.0 + ,525.0 + ,474.0 + ,539.0 + ,468.0 + ,543.0 + ,532.0 + ,565.0 + ,535.0 + ,534.0 + ,546.0 + ,494.0 + ,552.0 + ,511.0 + ,451.0 + ,537.0 + ,494.0 + ,549.0 + ,544.0 + ,598.0 + ,583.0 + ,582.0 + ,589.0 + ,578.0 + ,561.0 + ,592.0 + ,504.0 + ,545.0 + ,547.0 + ,585.0 + ,562.0 + ,520.0 + ,581.0 + ,590.0 + ,562.0 + ,548.0 + ,567.0 + ,542.0 + ,473.0 + ,531.0 + ,462.0 + ,479.0 + ,533.0 + ,552.0 + ,547.0 + ,562.0 + ,524.0 + ,479.0 + ,445.0 + ,406.0 + ,475.0 + ,589.0 + ,495.0 + ,484.0 + ,536.0 + ,555.0 + ,565.0 + ,564.0 + ,573.0 + ,569.0 + ,588.0 + ,546.0 + ,508.0 + ,560.0 + ,558.0 + ,516.0 + ,549.0 + ,595.0 + ,586.0 + ,597.0 + ,592.0 + ,538.0 + ,590.0 + ,576.0 + ,451.0 + ,538.0 + ,555.0 + ,532.0 + ,530.0 + ,553.0 + ,626.0 + ,601.0 + ,573.0 + ,569.0 + ,562.0 + ,468.0 + ,483.0 + ,460.0 + ,411.0 + ,458.0 + ,455.0 + ,600.0 + ,605.0 + ,545.0 + ,549.0 + ,415.0 + ,568.0 + ,577.0 + ,517.0 + ,558.0 + ,518.0 + ,489.0 + ,502.0 + ,569.0 + ,540.0 + ,550.0 + ,557.0 + ,542.0 + ,542.0 + ,582.0 + ,525.0 + ,584.0 + ,562.0 + ,639.0 + ,613.0 + ,604.0 + ,613.0 + ,625.0 + ,654.0 + ,638.0) > par4 = 'Unknown' > par3 = 'FALSE' > par2 = 'grey' > par1 = '6' > xlab = 'hoeveelheid staal in 1000 ton' > 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/1fsgg1267024371.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] 100 200 300 400 500 600 700 800 $counts [1] 3 47 49 75 170 93 18 $intensities [1] 6.593407e-05 1.032967e-03 1.076923e-03 1.648352e-03 3.736264e-03 [6] 2.043956e-03 3.956043e-04 $density [1] 6.593407e-05 1.032967e-03 1.076923e-03 1.648352e-03 3.736264e-03 [6] 2.043956e-03 3.956043e-04 $mids [1] 150 250 350 450 550 650 750 $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/2qpzm1267024371.tab") > > try(system("convert tmp/1fsgg1267024371.ps tmp/1fsgg1267024371.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 0.412 0.157 0.512