R version 2.10.1 (2009-12-14)
Copyright (C) 2009 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(0.0004619415
+ ,0.0004675645
+ ,0.0004537342
+ ,0.0004593573
+ ,0.0004406834
+ ,0.0004500941
+ ,0.0004685962
+ ,0.0004496081
+ ,0.0004688521
+ ,0.0004590911
+ ,0.0004587804
+ ,0.0004592529
+ ,0.0005336746
+ ,0.0005338313
+ ,0.0005517349
+ ,0.0005520146
+ ,0.0005523253
+ ,0.0006278578
+ ,0.000599586
+ ,0.0006089967
+ ,0.0006643019
+ ,0.000718216
+ ,0.0006441434
+ ,0.0006253128
+ ,0.0006257724
+ ,0.0006254927
+ ,0.0006258068
+ ,0.0005502048
+ ,0.00057832
+ ,0.0005782507
+ ,0.0005234128
+ ,0.0005252564
+ ,0.000515605
+ ,0.0005342788
+ ,0.0005158461
+ ,0.0005344337
+ ,0.0005154457
+ ,0.0005341892
+ ,0.0005904195
+ ,0.0005811171
+ ,0.0006371907
+ ,0.0005808047
+ ,0.0005902459
+ ,0.0006458549
+ ,0.0006303712
+ ,0.0006025294
+ ,0.0005839626
+ ,0.0005558475
+ ,0.0004809432
+ ,0.0004753202
+ ,0.0004003465
+ ,0.0003723871
+ ,0.0003905962
+ ,0.0002975526
+ ,0.0003130363
+ ,0.0003128797
+ ,0.0003314465
+ ,0.0003258235
+ ,0.0003257539
+ ,0.0003313769
+ ,0.0003163822
+ ,0.0003257539
+ ,0.0002924804
+ ,0.0002793078
+ ,0.0004292553
+ ,0.0004292553
+ ,0.0004292553
+ ,0.0004292553
+ ,0.0004292553
+ ,0.0004292553
+ ,0.0004292553
+ ,0.0004292553
+ ,0.0004442501
+ ,0.0004611192
+ ,0.0003017999
+ ,0.0003112106
+ ,0.0003111411
+ ,0.0003111411
+ ,0.0003017694
+ ,0.0003111411
+ ,0.0003261359
+ ,0.0003167641
+ ,0.0003017694
+ ,0.0003167641
+ ,0.0003261359
+ ,0.0003167251
+ ,0.0003197241
+ ,0.0003321938
+ ,0.0003321938
+ ,0.0003321938
+ ,0.0003171991
+ ,0.0003265708
+ ,0.0003321938
+ ,0.0003322634
+ ,0.0003321941
+ ,0.0003322331
+ ,0.0003292341
+ ,0.0003505026
+ ,0.0003598743
+ ,0.0003598743
+ ,0.0003654973
+ ,0.0003655669
+ ,0.0003748693
+ ,0.000365428
+ ,0.0004404711
+ ,0.0004404711
+ ,0.0004405406
+ ,0.0004217279
+ ,0.0004123562
+ ,0.0004873299
+ ,0.0004873299
+ ,0.0004873299
+ ,0.0005059869
+ ,0.0005153586
+ ,0.0004404544
+ ,0.0004684138
+ ,0.0004683443
+ ,0.0004590418
+ ,0.0005433873
+ ,0.0004590418
+ ,0.0005433873
+ ,0.0005339461
+ ,0.0005153584
+ ,0.0005342585
+ ,0.0005341196
+ ,0.0005342754
+ ,0.0005343449
+ ,0.0005623043
+ ,0.0004873306
+ ,0.0005154457
+ ,0.000440472
+ ,0.0004687438
+ ,0.0004686745
+ ,0.0004497744
+ ,0.000440472
+ ,0.0004217285
+ ,0.000440559
+ ,0.0004219021
+ ,0.0004219021
+ ,0.0003937869
+ ,0.0003937869
+ ,0.0003936303
+ ,0.0003936996
+ ,0.0003936996
+ ,0.0003936996
+ ,0.0004686733
+ ,0.0004404015
+ ,0.0004404708
+ ,0.0004404708
+ ,0.000468586
+ ,0.000468586
+ ,0.0004404708
+ ,0.0005154446
+ ,0.000506112
+ ,0.0005154144
+ ,0.0004404406
+ ,0.0004687124
+ ,0.0004686431
+ ,0.0004592714
+ ,0.000440528
+ ,0.0004594281
+ ,0.0004687305
+ ,0.0003937567
+ ,0.0003937176
+ ,0.0004693889
+ ,0.0004600172
+ ,0.0004410479
+ ,0.0004411172
+ ,0.0004411172
+ ,0.0005160909
+ ,0.0004878191
+ ,0.0004878884
+ ,0.0004878884
+ ,0.0004878884
+ ,0.0004122864
+ ,0.0004217277
+ ,0.0004124252
+ ,0.0004123559
+ ,0.0004404711
+ ,0.0003655669
+ ,0.0003747315
+ ,0.0003654774
+ ,0.0003654081
+ ,0.0003747798
+ ,0.0003654081
+ ,0.0003653386
+ ,0.0003934537
+ ,0.0003935926
+ ,0.0003747798
+ ,0.0003653386
+ ,0.0003655457
+ ,0.0003747998
+ ,0.0003654973
+ ,0.0003842408
+ ,0.0003842408
+ ,0.0003748691
+ ,0.0003561256
+ ,0.0003560561
+ ,0.0003467536
+ ,0.0003561949
+ ,0.0003561256
+ ,0.0003467539
+ ,0.0003561256
+ ,0.0003373822
+ ,0.0003373822
+ ,0.0003468234
+ ,0.0003467542
+ ,0.0003373824
+ ,0.0003467542
+ ,0.0003466846
+ ,0.0003373822
+ ,0.0003467539
+ ,0.0003317591
+ ,0.0003317591
+ ,0.0003317591
+ ,0.0003316896
+ ,0.0003317589
+ ,0.0003317589
+ ,0.0003317589
+ ,0.0003318284
+ ,0.0003356926
+ ,0.0003304397
+ ,0.0003454344
+ ,0.0003360627
+ ,0.0003454344
+ ,0.0003454344
+ ,0.0003361018
+ ,0.0003454735
+ ,0.0003454042
+ ,0.000335963
+ ,0.0003264758
+ ,0.0003317982
+ ,0.0003316599
+ ,0.0003260369
+ ,0.0003167043
+ ,0.0003017095
+ ,0.0003111117
+ ,0.000296092
+ ,0.0002961613
+ ,0.0002905632
+ ,0.0003085424
+ ,0.0003010955
+ ,0.0002920467
+ ,0.0003016361)
> x <- c(22854
+ ,23026
+ ,22927
+ ,22838
+ ,22686
+ ,22799
+ ,22820
+ ,22835
+ ,22854
+ ,22650
+ ,22657
+ ,22744
+ ,22664
+ ,22688
+ ,22447
+ ,22044
+ ,21796
+ ,21400
+ ,21470
+ ,21569
+ ,21711
+ ,21850
+ ,21943
+ ,21959
+ ,22008
+ ,22044
+ ,22024
+ ,22047
+ ,22067
+ ,22106
+ ,22005
+ ,21909
+ ,21793
+ ,21775
+ ,21830
+ ,21814
+ ,21802
+ ,21851
+ ,21812
+ ,21780
+ ,21771
+ ,21757
+ ,21759
+ ,21638
+ ,21510
+ ,21411
+ ,21385
+ ,21490
+ ,19065
+ ,19919
+ ,20919
+ ,20912
+ ,20913
+ ,20810
+ ,20808
+ ,20746
+ ,20601
+ ,20383
+ ,20091
+ ,20099
+ ,19896
+ ,19793
+ ,19421
+ ,19187
+ ,18969
+ ,18708
+ ,18566
+ ,18200
+ ,18200
+ ,17013
+ ,16998
+ ,16894
+ ,16937
+ ,16952
+ ,16882
+ ,16864
+ ,16944
+ ,16987
+ ,17060
+ ,17143
+ ,17218
+ ,17280
+ ,17278
+ ,17310
+ ,17389
+ ,17502
+ ,17609
+ ,17816
+ ,17982
+ ,17996
+ ,18247
+ ,18489
+ ,18510
+ ,18670
+ ,18823
+ ,18967
+ ,18946
+ ,19275
+ ,19546
+ ,19514
+ ,19836
+ ,19731
+ ,20165
+ ,20448
+ ,20624
+ ,20969
+ ,21135
+ ,21316
+ ,21642
+ ,21743
+ ,22021
+ ,22144
+ ,22416
+ ,22454
+ ,22584
+ ,22676
+ ,22878
+ ,22905
+ ,23049
+ ,23288
+ ,23442
+ ,23570
+ ,23791
+ ,23769
+ ,23921
+ ,23353
+ ,24162
+ ,24407
+ ,24552
+ ,24782
+ ,24783
+ ,25071
+ ,24496
+ ,25156
+ ,25228
+ ,25384
+ ,25740
+ ,25860
+ ,23791
+ ,24306
+ ,25278
+ ,24727
+ ,25568
+ ,25672
+ ,25405
+ ,25404
+ ,25858
+ ,26160
+ ,26088
+ ,26192
+ ,26177
+ ,26286
+ ,25606
+ ,26227
+ ,26143
+ ,26084
+ ,26119
+ ,26038
+ ,26038
+ ,26031
+ ,26010
+ ,25946
+ ,25088
+ ,25670
+ ,25425
+ ,25338
+ ,25157
+ ,25029
+ ,24876
+ ,24851
+ ,24747
+ ,24523
+ ,24301
+ ,24182
+ ,24085
+ ,23961
+ ,23699
+ ,23850
+ ,23669
+ ,23340
+ ,23394
+ ,23227
+ ,23263
+ ,23253
+ ,23239
+ ,23206
+ ,23131
+ ,23083
+ ,23088
+ ,22953
+ ,22955
+ ,22845
+ ,22834
+ ,22848
+ ,22784
+ ,22769
+ ,22741
+ ,22625
+ ,20344
+ ,20832
+ ,21083
+ ,21242
+ ,21287
+ ,21184
+ ,21044
+ ,20872
+ ,20819
+ ,20698
+ ,20424
+ ,20249
+ ,20151
+ ,19898
+ ,19652
+ ,19516
+ ,19416
+ ,19293
+ ,19015
+ ,18971
+ ,18902
+ ,18778
+ ,18622
+ ,18574
+ ,18480
+ ,18337
+ ,18279
+ ,18327
+ ,18343
+ ,18405
+ ,18464
+ ,18051
+ ,18591
+ ,18774
+ ,18859
+ ,18915
+ ,18463
+ ,18668
+ ,18836
+ ,19251
+ ,16875
+ ,17379
+ ,17859
+ ,18118
+ ,18373
+ ,18668)
> ylab = 'progress(%)/second'
> xlab = 'num CPUs'
> #'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/yougetitorg/rcomp/tmp/1vf801303056220.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)
> 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/yougetitorg/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab
> load(file="/var/yougetitorg/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/yougetitorg/rcomp/tmp/205r41303056220.tab") 
> 
> try(system("convert tmp/1vf801303056220.ps tmp/1vf801303056220.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  0.570   0.270   0.653