R version 3.0.2 (2013-09-25) -- "Frisbee Sailing"
Copyright (C) 2013 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.

> 
> RC.capture <- function (expression, collapse = NULL) {
+     resultConn <- textConnection('RC.resultText', open = 'w', local=TRUE)
+     sink(resultConn)
+     on.exit(function() {
+         sink()
+         close(resultConn)
+     })
+     expression
+     on.exit(NULL)
+     sink()
+     close(resultConn)
+     return(paste(c(RC.resultText, ''), collapse = collapse, sep = ''))
+ }
> RC.texteval <- function (sourceText, collapse = NULL, echo = TRUE) {
+     sourceConn <- textConnection(sourceText, open = 'r')
+     on.exit(close(sourceConn))
+     result <- RC.capture(source(file = sourceConn, local = FALSE, echo = echo, print.eval = TRUE), collapse = collapse)
+     on.exit(NULL)
+     close(sourceConn)
+     res <- ''
+     for(i in 1:length(result)) {
+       if (result[i]!='') res <- paste(res,result[i],'
+ ',sep='')
+     }
+     return(res)
+ }
> y <- c(1,2,2,2,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,6,6,6,6,6,6,6,6,6,6,6,6,6,6,6,6,6,6,6,6,6,6,6,6,6,6,6,6,6,6,6,6,6,6,6,6,6,6,6,6,6,6,6,6,6,6,6,6,7,7,7,7,7,7,7,7,7,7,7,7,7,7,7,7,7,7,7,7,7,7,7,7,7,7,7,7,7,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,8,9,9,9,9,9,9,9,9,9,9,9,9,9,9,9,9,9,9,9,10,10,10,10,10,11,11,11,12,12,12,12,13,14,17,23)
> x <- c(3.55
+ ,3.52
+ ,3.53
+ ,3.55
+ ,3.43
+ ,3.46
+ ,3.47
+ ,3.48
+ ,3.48
+ ,3.49
+ ,3.5
+ ,3.5
+ ,3.51
+ ,3.54
+ ,3.55
+ ,3.55
+ ,3.56
+ ,3.56
+ ,3.57
+ ,3.57
+ ,3.59
+ ,3.63
+ ,3.37
+ ,3.39
+ ,3.43
+ ,3.45
+ ,3.45
+ ,3.45
+ ,3.48
+ ,3.48
+ ,3.48
+ ,3.48
+ ,3.49
+ ,3.49
+ ,3.5
+ ,3.5
+ ,3.51
+ ,3.51
+ ,3.51
+ ,3.52
+ ,3.52
+ ,3.52
+ ,3.52
+ ,3.52
+ ,3.52
+ ,3.52
+ ,3.53
+ ,3.53
+ ,3.54
+ ,3.54
+ ,3.54
+ ,3.55
+ ,3.55
+ ,3.56
+ ,3.56
+ ,3.56
+ ,3.57
+ ,3.57
+ ,3.57
+ ,3.57
+ ,3.58
+ ,3.58
+ ,3.58
+ ,3.59
+ ,3.6
+ ,3.6
+ ,3.42
+ ,3.43
+ ,3.45
+ ,3.45
+ ,3.47
+ ,3.48
+ ,3.48
+ ,3.48
+ ,3.48
+ ,3.48
+ ,3.49
+ ,3.49
+ ,3.49
+ ,3.49
+ ,3.49
+ ,3.49
+ ,3.49
+ ,3.49
+ ,3.5
+ ,3.5
+ ,3.5
+ ,3.5
+ ,3.5
+ ,3.5
+ ,3.5
+ ,3.51
+ ,3.51
+ ,3.51
+ ,3.51
+ ,3.52
+ ,3.52
+ ,3.53
+ ,3.53
+ ,3.53
+ ,3.53
+ ,3.54
+ ,3.54
+ ,3.56
+ ,3.56
+ ,3.56
+ ,3.56
+ ,3.57
+ ,3.57
+ ,3.57
+ ,3.58
+ ,3.58
+ ,3.58
+ ,3.6
+ ,3.69
+ ,3.33
+ ,3.34
+ ,3.39
+ ,3.4
+ ,3.44
+ ,3.45
+ ,3.47
+ ,3.47
+ ,3.47
+ ,3.48
+ ,3.48
+ ,3.48
+ ,3.48
+ ,3.49
+ ,3.49
+ ,3.49
+ ,3.49
+ ,3.49
+ ,3.49
+ ,3.5
+ ,3.5
+ ,3.5
+ ,3.5
+ ,3.5
+ ,3.5
+ ,3.51
+ ,3.51
+ ,3.51
+ ,3.51
+ ,3.51
+ ,3.51
+ ,3.52
+ ,3.52
+ ,3.52
+ ,3.53
+ ,3.54
+ ,3.54
+ ,3.54
+ ,3.55
+ ,3.55
+ ,3.55
+ ,3.55
+ ,3.56
+ ,3.58
+ ,3.58
+ ,3.58
+ ,3.59
+ ,3.62
+ ,3.29
+ ,3.4
+ ,3.41
+ ,3.43
+ ,3.44
+ ,3.44
+ ,3.44
+ ,3.45
+ ,3.45
+ ,3.46
+ ,3.46
+ ,3.46
+ ,3.47
+ ,3.48
+ ,3.48
+ ,3.48
+ ,3.49
+ ,3.5
+ ,3.5
+ ,3.51
+ ,3.51
+ ,3.52
+ ,3.52
+ ,3.53
+ ,3.54
+ ,3.56
+ ,3.56
+ ,3.57
+ ,3.58
+ ,3.32
+ ,3.4
+ ,3.41
+ ,3.42
+ ,3.42
+ ,3.42
+ ,3.43
+ ,3.43
+ ,3.48
+ ,3.5
+ ,3.5
+ ,3.51
+ ,3.51
+ ,3.53
+ ,3.53
+ ,3.53
+ ,3.53
+ ,3.54
+ ,3.54
+ ,3.54
+ ,3.39
+ ,3.45
+ ,3.45
+ ,3.47
+ ,3.47
+ ,3.47
+ ,3.48
+ ,3.48
+ ,3.48
+ ,3.48
+ ,3.49
+ ,3.49
+ ,3.51
+ ,3.52
+ ,3.53
+ ,3.53
+ ,3.54
+ ,3.54
+ ,3.56
+ ,3.43
+ ,3.45
+ ,3.47
+ ,3.47
+ ,3.54
+ ,3.43
+ ,3.5
+ ,3.55
+ ,3.37
+ ,3.48
+ ,3.5
+ ,3.5
+ ,3.51
+ ,3.46
+ ,3.49
+ ,3.5)
> ylab = 'y'
> xlab = 'x'
> main = 'Scatterplot and histograms'
> #'GNU S' R Code compiled by R2WASP v. 1.2.327 ()
> #Author: root
> #To cite this work: Wessa P., (2014), Pearson Correlation (v1.0.9) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_correlation.wasp/
> #Source of accompanying publication: Office for Research, Development, and Education
> #
> x <- x[!is.na(y)]
> y <- y[!is.na(y)]
> y <- y[!is.na(x)]
> x <- x[!is.na(x)]
> postscript(file="/var/wessaorg/rcomp/tmp/133m51430538290.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
> histx <- hist(x, plot=FALSE)
> histy <- hist(y, plot=FALSE)
> maxcounts <- max(c(histx$counts, histx$counts))
> xrange <- c(min(x),max(x))
> yrange <- c(min(y),max(y))
> nf <- layout(matrix(c(2,0,1,3),2,2,byrow=TRUE), c(3,1), c(1,3), TRUE)
> par(mar=c(4,4,1,1))
> plot(x, y, xlim=xrange, ylim=yrange, xlab=xlab, ylab=ylab, sub=main)
> par(mar=c(0,4,1,1))
> barplot(histx$counts, axes=FALSE, ylim=c(0, maxcounts), space=0)
> par(mar=c(4,0,1,1))
> barplot(histy$counts, axes=FALSE, xlim=c(0, maxcounts), space=0, horiz=TRUE)
> dev.off()
null device 
          1 
> lx = length(x)
> makebiased = (lx-1)/lx
> varx = var(x)*makebiased
> vary = var(y)*makebiased
> corxy <- cor.test(x,y,method='pearson', na.rm = T)
> cxy <- as.matrix(corxy$estimate)[1,1]
> 
> #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,'Pearson Product Moment Correlation - Ungrouped Data',3,TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'Statistic',1,TRUE)
> a<-table.element(a,'Variable X',1,TRUE)
> a<-table.element(a,'Variable Y',1,TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,hyperlink('http://www.xycoon.com/arithmetic_mean.htm','Mean',''),header=TRUE)
> a<-table.element(a,mean(x))
> a<-table.element(a,mean(y))
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,hyperlink('http://www.xycoon.com/biased.htm','Biased Variance',''),header=TRUE)
> a<-table.element(a,varx)
> a<-table.element(a,vary)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,hyperlink('http://www.xycoon.com/biased1.htm','Biased Standard Deviation',''),header=TRUE)
> a<-table.element(a,sqrt(varx))
> a<-table.element(a,sqrt(vary))
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,hyperlink('http://www.xycoon.com/covariance.htm','Covariance',''),header=TRUE)
> a<-table.element(a,cov(x,y),2)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,hyperlink('http://www.xycoon.com/pearson_correlation.htm','Correlation',''),header=TRUE)
> a<-table.element(a,cxy,2)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,hyperlink('http://www.xycoon.com/coeff_of_determination.htm','Determination',''),header=TRUE)
> a<-table.element(a,cxy*cxy,2)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,hyperlink('http://www.xycoon.com/ttest_statistic.htm','T-Test',''),header=TRUE)
> a<-table.element(a,as.matrix(corxy$statistic)[1,1],2)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'p-value (2 sided)',header=TRUE)
> a<-table.element(a,(p2 <- as.matrix(corxy$p.value)[1,1]),2)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'p-value (1 sided)',header=TRUE)
> a<-table.element(a,p2/2,2)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'Degrees of Freedom',header=TRUE)
> a<-table.element(a,lx-2,2)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'Number of Observations',header=TRUE)
> a<-table.element(a,lx,2)
> a<-table.row.end(a)
> a<-table.end(a)
> table.save(a,file="/var/wessaorg/rcomp/tmp/23uqx1430538290.tab") 
> library(moments)
> library(nortest)
> jarque.x <- jarque.test(x)
> jarque.y <- jarque.test(y)
> if(lx>7) {
+ ad.x <- ad.test(x)
+ ad.y <- ad.test(y)
+ }
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a,'Normality Tests',1,TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,paste('<pre>',RC.texteval('jarque.x'),'</pre>',sep=''))
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,paste('<pre>',RC.texteval('jarque.y'),'</pre>',sep=''))
> a<-table.row.end(a)
> if(lx>7) {
+ a<-table.row.start(a)
+ a<-table.element(a,paste('<pre>',RC.texteval('ad.x'),'</pre>',sep=''))
+ a<-table.row.end(a)
+ a<-table.row.start(a)
+ a<-table.element(a,paste('<pre>',RC.texteval('ad.y'),'</pre>',sep=''))
+ a<-table.row.end(a)
+ }
> a<-table.end(a)
> table.save(a,file="/var/wessaorg/rcomp/tmp/3ptfp1430538290.tab") 
> library(car)
> postscript(file="/var/wessaorg/rcomp/tmp/4hcob1430538290.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
> qq.plot(x,main='QQplot of variable x')
Warning message:
'qq.plot' is deprecated.
Use 'qqPlot' instead.
See help("Deprecated") and help("car-deprecated"). 
> dev.off()
null device 
          1 
> postscript(file="/var/wessaorg/rcomp/tmp/561q71430538290.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
> qq.plot(y,main='QQplot of variable y')
Warning message:
'qq.plot' is deprecated.
Use 'qqPlot' instead.
See help("Deprecated") and help("car-deprecated"). 
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/133m51430538290.ps tmp/133m51430538290.png",intern=TRUE))
character(0)
> try(system("convert tmp/4hcob1430538290.ps tmp/4hcob1430538290.png",intern=TRUE))
character(0)
> try(system("convert tmp/561q71430538290.ps tmp/561q71430538290.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  1.355   0.192   1.559