source('/home/pw/wessanet/cretab') y <- c(20,16,19.8,18.4,17.1,15.5,14.7,17.1,15.4,16.2,15,17.2,16,17,14.4) x <- c(88.6,71.6,93.3,84.3,80.6,75.2,69.7,82,69.4,83.3,79.6,82.6,80.6,83.5,76.3) par1 = '0' ylab = 'y' xlab = 'x' par1 <- '0' #'GNU S' R Code compiled by R2WASP v. 1.2.327 (Thu, 22 Sep 2016 14:29:32 +0200) #Author: root #To cite this work: Wessa P. (2016) QQ Plots and Histograms (v1.0.5) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_qqplot.wasp/ #Source of accompanying publication: # library(car) if(par1 == '') par1 <- 0 par1 <- as.numeric(par1) if (par1 < 0) par1 <- 0 postscript(file="/home/pw/wessanet/rcomp/tmp/1onb81477482583.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) if(par1 > 0) { myhist<-hist(x, breaks=par1, col=2,main='Histogram (series X)') } else { myhist <- hist(x, col=2) } dev.off() postscript(file="/home/pw/wessanet/rcomp/tmp/23jyv1477482583.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) if(par1 > 0) { myhist<-hist(y, breaks=par1, col=2,main='Histogram (series Y)') } else { myhist <- hist(y, col=2) } dev.off() postscript(file="/home/pw/wessanet/rcomp/tmp/3agfx1477482583.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) qqPlot(x,main='Normal Q-Q Plot (series X)',ylab='Sample Quantiles',xlab='Normal Quantiles') dev.off() postscript(file="/home/pw/wessanet/rcomp/tmp/49jr21477482583.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) qqPlot(y,main='Normal Q-Q Plot (series Y)',ylab='Sample Quantiles',xlab='Normal Quantiles') dev.off() postscript(file="/home/pw/wessanet/rcomp/tmp/5n0z91477482583.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) qqplot(x,y,main='Q-Q Plot (series X vs Y)',ylab='Sample Quantiles of series Y',xlab='Sample Quantiles of series X') qqline(x,y) grid() dev.off() a<-table.start() a<-table.row.start(a) a<-table.element(a,'Descriptive Statistics',2,TRUE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'mean of X',header=TRUE) a<-table.element(a,mean(x)) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'standard deviation of X',header=TRUE) a<-table.element(a,sd(x)) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'mean of Y',header=TRUE) a<-table.element(a,mean(y)) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'standard deviation of Y',header=TRUE) a<-table.element(a,sd(y)) a<-table.row.end(a) a<-table.end(a) table.save(a,file="/home/pw/wessanet/rcomp/tmp/6a7b71477482583.tab") try(system("convert /home/pw/wessanet/rcomp/tmp/1onb81477482583.ps /home/pw/wessanet/rcomp/tmp/1onb81477482583.png",intern=TRUE)) try(system("convert /home/pw/wessanet/rcomp/tmp/23jyv1477482583.ps /home/pw/wessanet/rcomp/tmp/23jyv1477482583.png",intern=TRUE)) try(system("convert /home/pw/wessanet/rcomp/tmp/3agfx1477482583.ps /home/pw/wessanet/rcomp/tmp/3agfx1477482583.png",intern=TRUE)) try(system("convert /home/pw/wessanet/rcomp/tmp/49jr21477482583.ps /home/pw/wessanet/rcomp/tmp/49jr21477482583.png",intern=TRUE)) try(system("convert /home/pw/wessanet/rcomp/tmp/5n0z91477482583.ps /home/pw/wessanet/rcomp/tmp/5n0z91477482583.png",intern=TRUE))