R version 2.6.0 (2007-10-03)
Copyright (C) 2007 The R Foundation for Statistical Computing
ISBN 3-900051-07-0

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.

> y <- c(-7315,6179,-7146,2033,3447,-2844,-178,39,-3386,1915,1602,-3011,9165,-3699,777,3033,-12392,2467,1932,-2101,3876,-1876,-1679,-893,-5942,-1194,2869,1663,5207,1902,-912,-1428,2833,-1176,280,-1575,5831,-2561,2152,-7643,-1304,1348,-2029,5939,-5121,-574,2362,-2874,-2936,2707,-3112,4182,4363,-948,1813,-1356,1373,-460,3739,-7519,5797,4664,-1047,-5647,194,-4177,3149,154,-2136,3999,-3206,3455,-5506,-3550,1790,5278,2999,-2759,1227,-1875,-2315,2307,-2419)
> x <- c(-1.48033E-05
+ ,-0.000559728
+ ,-0.001661464
+ ,-0.000206111
+ ,0.000391534
+ ,-0.001306447
+ ,-0.001030865
+ ,-0.000471111
+ ,-0.000514407
+ ,-0.000349283
+ ,-0.00163664
+ ,0.00059502
+ ,-0.001215959
+ ,-0.001610338
+ ,0.000514799
+ ,-0.001358222
+ ,-0.001649181
+ ,9.56556E-05
+ ,0.000571407
+ ,-0.002569009
+ ,0.000609157
+ ,0.000228
+ ,0.001641234
+ ,0.001548685
+ ,-0.001081926
+ ,0.000501447
+ ,-0.001669946
+ ,-0.001643796
+ ,0.000933258
+ ,0.001225474
+ ,0.000671829
+ ,-0.001156959
+ ,0.001554204
+ ,0.000697463
+ ,0.000298945
+ ,0.000118779
+ ,-0.000454979
+ ,-0.000758256
+ ,-0.001883814
+ ,0.000350195
+ ,0.000809793
+ ,-0.0002472
+ ,-0.000114181
+ ,-0.000733486
+ ,-0.000524683
+ ,0.001066399
+ ,8.9809E-05
+ ,-0.000808946
+ ,-0.000936656
+ ,-0.000496893
+ ,0.001386333
+ ,0.001422328
+ ,-0.000277787
+ ,-0.00029437
+ ,-0.001290847
+ ,-0.000512573
+ ,0.000477472
+ ,-0.000111006
+ ,8.17718E-05
+ ,1.75385E-05
+ ,-0.000274125
+ ,-0.000991596
+ ,-0.001164287
+ ,-0.001917012
+ ,0.000216283
+ ,-0.000619477
+ ,-0.000813221
+ ,-0.000210582
+ ,-0.001277084
+ ,0.000931432
+ ,0.001251757
+ ,0.00067761
+ ,-0.001415316
+ ,-0.001332513
+ ,-0.002224636
+ ,0.000261857
+ ,-0.001180315
+ ,-0.001831976
+ ,-0.000562689
+ ,-0.001427773
+ ,0.000379311
+ ,0.002175279
+ ,0.000570308)
> ylab = 'Inschrijvingen PW'
> xlab = 'Tot. CPI gebruik PV'
> #'GNU S' R Code compiled by R2WASP v. 1.0.44 ()
> #Author: Prof. Dr. P. Wessa
> #To cite this work: Wessa P., (2007), Pearson Correlation (v1.0.2) 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
> #Technical description: Write here your technical program description (don't use hard returns!)
> postscript(file="/var/www/html/rcomp/tmp/182s71199534624.ps",horizontal=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)
> 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')
> cxy <- as.matrix(corxy$estimate)[1,1]
> load(file='/var/www/html/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/www/html/rcomp/tmp/2o5mm1199534624.tab") 
> 
> system("convert tmp/182s71199534624.ps tmp/182s71199534624.png")
> 
> 
> proc.time()
   user  system elapsed 
  0.615   0.165   1.052