R version 2.13.0 (2011-04-13)
Copyright (C) 2011 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
Platform: i486-pc-linux-gnu (32-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.

> y <- c(270018
+ ,179444
+ ,222373
+ ,218443
+ ,162874
+ ,70849
+ ,498732
+ ,33186
+ ,207822
+ ,213274
+ ,296074
+ ,237633
+ ,164107
+ ,358752
+ ,222781
+ ,369889
+ ,305704
+ ,322896
+ ,176082
+ ,263411
+ ,271965
+ ,425544
+ ,179306
+ ,189897
+ ,219420
+ ,207280
+ ,267198
+ ,270750
+ ,155915
+ ,327474
+ ,278019
+ ,204039
+ ,318563
+ ,97717
+ ,369331
+ ,273950
+ ,422946
+ ,215710
+ ,115469
+ ,336047
+ ,324178
+ ,166266
+ ,195153
+ ,173510
+ ,153778
+ ,455168
+ ,78800
+ ,208051
+ ,334657
+ ,175523
+ ,213060
+ ,24188
+ ,372238
+ ,65029
+ ,101097
+ ,279012
+ ,302218
+ ,323485
+ ,339837
+ ,252529
+ ,370483
+ ,303406
+ ,233632
+ ,264889
+ ,228595
+ ,216027
+ ,187965
+ ,237323
+ ,232765
+ ,175699
+ ,239314
+ ,73566
+ ,242585
+ ,187167
+ ,191920
+ ,359644
+ ,341637
+ ,206059
+ ,201783
+ ,182231
+ ,153613
+ ,447353
+ ,145943
+ ,280343
+ ,80953
+ ,150216
+ ,156923
+ ,365370
+ ,318651
+ ,179797
+ ,251466
+ ,254506
+ ,185890
+ ,263577
+ ,314255
+ ,189252
+ ,222504
+ ,285198
+ ,368878
+ ,397681
+ ,287015
+ ,285330
+ ,186856
+ ,43287
+ ,185468
+ ,219475
+ ,259692
+ ,301614
+ ,121726
+ ,154165
+ ,306952
+ ,297982
+ ,23623
+ ,195817
+ ,61857
+ ,163766
+ ,384053
+ ,21054
+ ,252805
+ ,31961
+ ,311281
+ ,240153
+ ,174892
+ ,152043
+ ,38214
+ ,198094
+ ,353021
+ ,196269
+ ,403932
+ ,316105
+ ,396725
+ ,187992
+ ,102424
+ ,283950
+ ,401260
+ ,137843
+ ,383703
+ ,157429
+ ,236370
+ ,282399
+ ,217478
+ ,366774
+ ,236660
+ ,173260
+ ,323545
+ ,168994
+ ,246745
+ ,301703
+ ,1
+ ,14688
+ ,98
+ ,455
+ ,0
+ ,0
+ ,233143
+ ,368328
+ ,0
+ ,203
+ ,7199
+ ,46660
+ ,17547
+ ,116678
+ ,969
+ ,199726)
> x <- c(140824
+ ,110459
+ ,105079
+ ,112098
+ ,43929
+ ,76173
+ ,187326
+ ,22807
+ ,144408
+ ,66485
+ ,79089
+ ,81625
+ ,68788
+ ,103297
+ ,69446
+ ,114948
+ ,167949
+ ,125081
+ ,125818
+ ,136588
+ ,112431
+ ,103037
+ ,82317
+ ,118906
+ ,83515
+ ,104581
+ ,103129
+ ,83243
+ ,37110
+ ,113344
+ ,139165
+ ,86652
+ ,112302
+ ,69652
+ ,119442
+ ,69867
+ ,101629
+ ,70168
+ ,31081
+ ,103925
+ ,92622
+ ,79011
+ ,93487
+ ,64520
+ ,93473
+ ,114360
+ ,33032
+ ,96125
+ ,151911
+ ,89256
+ ,95671
+ ,5950
+ ,149695
+ ,32551
+ ,31701
+ ,100087
+ ,169707
+ ,150491
+ ,120192
+ ,95893
+ ,151715
+ ,176225
+ ,59900
+ ,104767
+ ,114799
+ ,72128
+ ,143592
+ ,89626
+ ,131072
+ ,126817
+ ,81351
+ ,22618
+ ,88977
+ ,92059
+ ,81897
+ ,108146
+ ,126372
+ ,249771
+ ,71154
+ ,71571
+ ,55918
+ ,160141
+ ,38692
+ ,102812
+ ,56622
+ ,15986
+ ,123534
+ ,108535
+ ,93879
+ ,144551
+ ,56750
+ ,127654
+ ,65594
+ ,59938
+ ,146975
+ ,143372
+ ,168553
+ ,183500
+ ,165986
+ ,184923
+ ,140358
+ ,149959
+ ,57224
+ ,43750
+ ,48029
+ ,104978
+ ,100046
+ ,101047
+ ,197426
+ ,160902
+ ,147172
+ ,109432
+ ,1168
+ ,83248
+ ,25162
+ ,45724
+ ,110529
+ ,855
+ ,101382
+ ,14116
+ ,89506
+ ,135356
+ ,116066
+ ,144244
+ ,8773
+ ,102153
+ ,117440
+ ,104128
+ ,134238
+ ,134047
+ ,279488
+ ,79756
+ ,66089
+ ,102070
+ ,146760
+ ,154771
+ ,165933
+ ,64593
+ ,92280
+ ,67150
+ ,128692
+ ,124089
+ ,125386
+ ,37238
+ ,140015
+ ,150047
+ ,154451
+ ,156349
+ ,0
+ ,6023
+ ,0
+ ,0
+ ,0
+ ,0
+ ,84601
+ ,68946
+ ,0
+ ,0
+ ,1644
+ ,6179
+ ,3926
+ ,52789
+ ,0
+ ,100350)
> ylab = 'Time spent in RFC'
> xlab = 'Total Size'
> postscript(file="/var/wessaorg/rcomp/tmp/1r5411324126583.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/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/2m1kh1324126583.tab") 
> 
> try(system("convert tmp/1r5411324126583.ps tmp/1r5411324126583.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  0.512   0.128   0.640