R version 3.2.3 (2015-12-10) -- "Wooden Christmas-Tree" Copyright (C) 2015 The R Foundation for Statistical Computing Platform: x86_64-pc-linux-gnu (64-bit) 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(7612 + ,7381 + ,6978 + ,6819 + ,6688 + ,6454 + ,6679 + ,6921 + ,7807 + ,7898 + ,7832 + ,7384 + ,7620 + ,7281 + ,6929 + ,6587 + ,6071 + ,5928 + ,5964 + ,6374 + ,7160 + ,7213 + ,6890 + ,6525 + ,6739 + ,6580 + ,6391 + ,6254 + ,6114 + ,5978 + ,6315 + ,6427 + ,7132 + ,7292 + ,7708 + ,7525 + ,7450 + ,7526 + ,7263 + ,7070 + ,6893 + ,6781 + ,7188 + ,7015 + ,8273 + ,8470 + ,8230 + ,8137 + ,8122 + ,8367 + ,8141 + ,7750 + ,7504 + ,7330 + ,7608 + ,7647 + ,8942 + ,8865 + ,8320 + ,8207 + ,8105 + ,8290 + ,8162 + ,8051 + ,7699 + ,7440 + ,7656 + ,7549 + ,9086 + ,8942 + ,8764 + ,8500 + ,8239 + ,8443 + ,8349 + ,8288 + ,7970 + ,7496 + ,7745 + ,7543 + ,9036 + ,9075 + ,8859 + ,8605 + ,8419 + ,8495 + ,8284 + ,7582 + ,7691 + ,7046 + ,7442 + ,7596 + ,8597 + ,8436 + ,7881 + ,7477 + ,7508 + ,7361 + ,7299 + ,6914 + ,6768 + ,6746 + ,7052 + ,7139 + ,7714 + ,7750 + ,7622 + ,7424 + ,7444 + ,7208 + ,7128 + ,7022 + ,6688 + ,6199 + ,6400 + ,6474 + ,7182 + ,7330 + ,7410 + ,7442 + ,7753 + ,7762 + ,7814 + ,7838 + ,7298 + ,7155 + ,7076 + ,7450 + ,8216 + ,8246 + ,8335 + ,8171 + ,8485 + ,8435 + ,8369 + ,8210 + ,7888 + ,8061 + ,8139 + ,7837 + ,8943 + ,8523 + ,8104 + ,7969 + ,7921 + ,7930 + ,7706 + ,7552 + ,7379 + ,6946 + ,7128 + ,7393 + ,8092 + ,8004 + ,7903 + ,7710 + ,7867 + ,7860 + ,7723 + ,7477 + ,7126 + ,7161 + ,7162 + ,7406 + ,7944 + ,8084 + ,8088 + ,7972 + ,8184 + ,7914 + ,7845 + ,7610 + ,7278 + ,6883 + ,7123 + ,7182 + ,7912 + ,7893 + ,7671 + ,7403 + ,7663 + ,7589 + ,7450 + ,7069 + ,6670 + ,6285 + ,6506 + ,6539 + ,7291 + ,7391 + ,7126 + ,6752 + ,6835 + ,6664 + ,6562 + ,6174 + ,5741 + ,5398 + ,5203 + ,5673 + ,6379 + ,6418 + ,6272 + ,6059) > par4 = 'Unknown' > par3 = 'FALSE' > par2 = 'grey' > par1 = '' > xlab = 'Hoeveelheid' > main = 'NWWz, -25jr, Man, Brussel' > par4 <- 'Unknown' > par3 <- 'FALSE' > par2 <- 'grey' > par1 <- '' > #'GNU S' R Code compiled by R2WASP v. 1.2.327 (Wed, 26 Aug 2015 15:05:17 +0100) > #Author: root > #To cite this work: Wessa P., (2015), Histogram (v1.0.15) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_histogram.wasp/ > #Source of accompanying publication: > # > 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/1hvak1455976427.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 { + barplot(mytab <- sort(table(x),T),col=par2,main='Frequency Plot',xlab=xlab,ylab='Absolute Frequency') + } > dev.off() null device 1 > 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/2ik5b1455976427.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/3mfe01455976427.tab") + } > > try(system("convert tmp/1hvak1455976427.ps tmp/1hvak1455976427.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 0.545 0.132 0.668