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(109,88,87,78,72,85,73,77,90,78,96,84,87,105,89,78,96,79,90,96,95,79,86,99,111,105,123,102,82,86,117,94,110,86,109,118,64,78,86,98,64,84,87,88,106,78,79,95,98,101,104,91,78,123,109,77,77,102,88,81,83,72,95,109,91,87,112,95,103,108,116,71,70,105,112,86,109,112,97,122,97,64,84,110,70,100,102,108,123,78,109,106,95,108,74,86,82,88,96,114,111,88,123,89,114,121,100,83,92,108,87,84,95,90,82,84,96,109,84,83,93,120,83,103,103,102,90,69,114,64,87,116,121,90,92,81,84,106,75,79,80,74,66,81,98,99,68,100,97,99,104,82,92,99,109,106,89,87,98,95,92,90) > par1 = '15' > xlab = 'Test Score' > main = 'Written Language, Year 7' > par1 <- '15' > #'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/1adcr1382376387.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/2ecmf1382376387.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/3e9w01382376387.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/4dfe41382376387.tab") > > try(system("convert tmp/1adcr1382376387.ps tmp/1adcr1382376387.png",intern=TRUE)) character(0) > try(system("convert tmp/2ecmf1382376387.ps tmp/2ecmf1382376387.png",intern=TRUE)) character(0) > try(system("convert tmp/3e9w01382376387.ps tmp/3e9w01382376387.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.629 0.394 2.005