R version 3.0.2 (2013-09-25) -- "Frisbee Sailing" Copyright (C) 2013 The R Foundation for Statistical Computing Platform: i686-pc-linux-gnu (32-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(85,90,81,56,85,51,73,60,80,121,81,50,99,74,78,74,30,32,81,92,69,73,91,91,94,89,86,80,69,56,49,46,75,57,69,85,76,36,21,57,53,75,38,59,66,27,77,58,34,60,36,78,82,48,40,64,73,38,55,56,96,85,70,85,11,57,78,0,59,65,46,99,70,47,47,76,73,49,71,63,57,18,71,82,49,63,102,55,53,55,45,55,41,42,55,138,47,54,82,69,59,55,59,51,23,82,58,65,133,102,28,69,23,49,65,58,78,53,50,49,72,61,0,51,75,71,29,10,56,0,1,60,89,97,98,94,83,69,89,94,90,78,88,102,82,110,90,90,80,68,89,94,81,65,102,86,84,89,92,72,82,94,81,80,91,85,86,103,85,93,82,60,81,94,98,79,83,86,75,86,94,85,90,64,105,96,82,95,99,82,102,77,75,77,77,98,87,73,73,85,76,94,73,60,99,106,94,89,79,112,94,82,75,86,79,82,92,94,84,93,60,82,71,89,110,74,98,96,106,103,87,98,72,98,84,77,85,85,99,90,87,99,89,97,106,82,90,92,95,99,86,90,84,70,92,98,110,92,93,87,91,104,88,92,107,85,93,89,101,94,82,88,73,112,93,77,93,93,87,95,93,74,69,80,82,98,112,101,87,90,72,85,97,107,77,89,103,87,112,95,92,89,86,85,92,78) > par1 = '10' > xlab = 'N Blogs' > main = 'Blogging Acivity' > par1 <- '10' > #'GNU S' R Code compiled by R2WASP v. 1.2.291 () > #Author: root > #To cite this work: Ian E. Holliday, 2012, Histogram, QQplot and Density Plot (v1.0.1) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/Ian.Holliday/rwasp_Distributional Plots.wasp/ > #Source of accompanying publication: > # > x <- x[!is.na(x)] > postscript(file="/var/fisher/rcomp/tmp/1fbsl1385984070.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > par1 <- as.numeric(par1) > myhist<-hist(x, breaks=par1, col=2, xlab=xlab, main=main) > dev.off() null device 1 > postscript(file="/var/fisher/rcomp/tmp/2jf1s1385984070.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > qqnorm(x) > qqline(x, col=2) > dev.off() null device 1 > postscript(file="/var/fisher/rcomp/tmp/3mwrg1385984070.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(density(x), xlab=xlab, main=main) > dev.off() null device 1 > > #Note: the /var/fisher/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab > load(file="/var/fisher/rcomp/createtable") > > a<-table.start() > a<-table.row.start(a) > a<-table.element(a,'Descriptive Statistics',3,TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'mean',header=TRUE) > a<-table.element(a,signif(mean(x), 6)) > a<-table.element(a,hyperlink('http://www.xycoon.com/arithmetic_mean.htm','formula','click to see the formula')) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'standard deviation',header=TRUE) > a<-table.element(a,signif(sd(x), 6)) > a<-table.element(a,hyperlink('http://www.xycoon.com/unbiased1.htm','formula','click to see the formula')) > a<-table.row.end(a) > a<-table.end(a) > table.save(a,file="/var/fisher/rcomp/tmp/4vla51385984070.tab") > > try(system("convert tmp/1fbsl1385984070.ps tmp/1fbsl1385984070.png",intern=TRUE)) character(0) > try(system("convert tmp/2jf1s1385984070.ps tmp/2jf1s1385984070.png",intern=TRUE)) character(0) > try(system("convert tmp/3mwrg1385984070.ps tmp/3mwrg1385984070.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 2.207 0.521 2.702