R version 2.8.0 (2008-10-20)
Copyright (C) 2008 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(727,817,918,786,803,756,725,523,538,587,505,521,498,550,637,622,668,669,670,499,539,593,429,622,533,655,835,686,706,869,777,739,637,597,629,940,444,496,801,659,767,876,601,697,745,655,572,628,650,677,900,780,896,1092,823,735,770,915,645,566,707,785,762,712,714,823,609,620,619,638,483,535,617,698,804,824,878,1019,974,773,734,827,804,721,659,732,839,994,828,1039,1072,803,1035,922,834,1739,359,513,699,741,793,877,750,752,675,682,583,632,606,645,980,847,941,1066,936,880,808,741,780,675,782,795,873,727,998,768,714,782,578,664,560,516,752,597,716,691,752,718,737,621,472,719,497,536,653,605,637,743,719,653,675,590,527,534,463,542,568,501,678,774,665,742,715,638,656,606,498,587,677,547,871,731,752,862,619,700,667,667,650,547,637,655,703,886,896,831,741,833,750,779,655,739,845,795,1021,726,1045,915,852,772,729,755,691,729,702,702,894,765,753,876,781,776,606,775,663,649,821,771,635,1070,693,779)
> x <- c(2560
+ ,2491
+ ,2380
+ ,2291
+ ,2079
+ ,1929
+ ,1851
+ ,1607
+ ,1661
+ ,2259
+ ,1668
+ ,2011
+ ,1944
+ ,1958
+ ,1844
+ ,1868
+ ,1701
+ ,2338
+ ,2018
+ ,1302
+ ,2168
+ ,2139
+ ,1560
+ ,2093
+ ,1973
+ ,2090
+ ,2811
+ ,1984
+ ,1849
+ ,2433
+ ,2071
+ ,1855
+ ,1756
+ ,1898
+ ,1770
+ ,1969
+ ,1769
+ ,2139
+ ,3013
+ ,2061
+ ,2132
+ ,2973
+ ,2081
+ ,2257
+ ,2075
+ ,2084
+ ,1747
+ ,2092
+ ,1919
+ ,2551
+ ,2643
+ ,2153
+ ,2496
+ ,2645
+ ,2035
+ ,2294
+ ,2205
+ ,2044
+ ,1762
+ ,1897
+ ,1821
+ ,1905
+ ,2111
+ ,1643
+ ,1956
+ ,1977
+ ,1685
+ ,1393
+ ,1574
+ ,1793
+ ,1562
+ ,1510
+ ,1675
+ ,1965
+ ,2173
+ ,2395
+ ,2197
+ ,2257
+ ,2885
+ ,1594
+ ,1950
+ ,1772
+ ,1280
+ ,1724
+ ,1473
+ ,1461
+ ,1576
+ ,1900
+ ,1618
+ ,2303
+ ,1994
+ ,1575
+ ,1893
+ ,1788
+ ,1817
+ ,3233
+ ,727
+ ,1121
+ ,1665
+ ,1401
+ ,1415
+ ,2058
+ ,1544
+ ,1379
+ ,1402
+ ,1313
+ ,1296
+ ,1398
+ ,1288
+ ,1563
+ ,1972
+ ,1496
+ ,1481
+ ,1819
+ ,1479
+ ,1635
+ ,1511
+ ,1547
+ ,1388
+ ,1958
+ ,1390
+ ,1597
+ ,1842
+ ,1396
+ ,1671
+ ,1385
+ ,1632
+ ,1313
+ ,1300
+ ,1431
+ ,1398
+ ,1198
+ ,1292
+ ,1434
+ ,1660
+ ,1837
+ ,1455
+ ,1315
+ ,1642
+ ,1069
+ ,1209
+ ,1586
+ ,1122
+ ,1063
+ ,1125
+ ,1414
+ ,1347
+ ,1403
+ ,1299
+ ,1547
+ ,1515
+ ,1247
+ ,1639
+ ,1296
+ ,1063
+ ,1282
+ ,1365
+ ,1268
+ ,1532
+ ,1455
+ ,1393
+ ,1515
+ ,1510
+ ,1225
+ ,1577
+ ,1417
+ ,1224
+ ,1693
+ ,1633
+ ,1639
+ ,1914
+ ,1586
+ ,1552
+ ,2081
+ ,1500
+ ,1437
+ ,1470
+ ,1849
+ ,1387
+ ,1592
+ ,1589
+ ,1798
+ ,1935
+ ,1887
+ ,2027
+ ,2080
+ ,1556
+ ,1682
+ ,1785
+ ,1869
+ ,1781
+ ,2082
+ ,2570
+ ,1862
+ ,1936
+ ,1504
+ ,1765
+ ,1607
+ ,1577
+ ,1493
+ ,1615
+ ,1700
+ ,1335
+ ,1523
+ ,1623
+ ,1540
+ ,1637
+ ,1524
+ ,1419
+ ,1821
+ ,1593
+ ,1357
+ ,1263
+ ,1750
+ ,1405
+ ,1393
+ ,1639
+ ,1679
+ ,1551
+ ,1744
+ ,1429
+ ,1784)
> ylab = 'Wallonië'
> xlab = 'Vlaanderen'
> #'GNU S' R Code compiled by R2WASP v. 1.0.44 ()
> #Author: Prof. Dr. P. Wessa
> #To cite this work: AUTHOR(S), (YEAR), YOUR SOFTWARE TITLE (vNUMBER) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_YOURPAGE.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/freestat/rcomp/tmp/1y41k1229180914.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]
> 
> #Note: the /var/www/html/freestat/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab
> load(file="/var/www/html/freestat/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/freestat/rcomp/tmp/2q7lb1229180915.tab") 
> 
> system("convert tmp/1y41k1229180914.ps tmp/1y41k1229180914.png")
> 
> 
> proc.time()
   user  system elapsed 
  0.581   0.241   0.630