R version 2.15.2 (2012-10-26) -- "Trick or Treat"
Copyright (C) 2012 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
Platform: i686-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(2.465174278
+ ,0.931192069
+ ,0.711452013
+ ,0.377183262
+ ,2.236740543
+ ,0.713040346
+ ,5.398554225
+ ,0.836420224
+ ,1.066455767
+ ,0.539451529
+ ,1.148668486
+ ,0.887746415
+ ,0.982298258
+ ,1.833295857
+ ,1.213173841
+ ,0.875413916
+ ,1.886406807
+ ,2.501723932
+ ,20.68464671
+ ,1.273047089
+ ,1.351081227
+ ,0.510491792
+ ,1.02685726
+ ,0.9435998
+ ,0.83565382
+ ,2.323293959
+ ,1.215575543
+ ,0.682603434
+ ,2.298259358
+ ,0.495497906
+ ,1.449464423
+ ,1.230981011
+ ,1.946545359
+ ,41.09568458
+ ,1.166003988
+ ,1.039045527
+ ,1.727527078
+ ,0.57014764
+ ,6.258638225
+ ,1.456811622
+ ,1.20360678
+ ,11.6070412
+ ,8.86957821
+ ,1.606921259
+ ,8.450290249
+ ,3.500918366
+ ,4.358933039
+ ,5.803955039
+ ,3.271119553
+ ,13.79015521
+ ,19.1986642
+ ,12.66054399
+ ,6.604966516
+ ,2.283512797
+ ,5.588444625
+ ,14.48342128
+ ,6.047639561
+ ,2.265745018
+ ,3.73875744
+ ,3.655522621
+ ,10.73115923
+ ,51.9968201
+ ,11.32211324
+ ,1.422567613
+ ,22.78075733
+ ,1.753938078
+ ,3.028556706
+ ,36.79690672
+ ,11.70628687
+ ,2.146786529
+ ,35.86205879
+ ,24.76992116
+ ,1.141761387
+ ,4.227539945
+ ,1.74873164
+ ,2.923244097
+ ,1.625318562
+ ,1.498778844
+ ,5.821176903
+ ,18.54626322
+ ,1.115054801
+ ,59.47912753
+ ,8.348444665
+ ,5.885859465
+ ,51.86285367
+ ,17.2267059
+ ,34.4208771
+ ,1.012470134
+ ,12.87098797
+ ,1.819797467
+ ,18.99327441
+ ,2.695604249
+ ,0.998056641
+ ,7.267825576
+ ,5.985410164
+ ,15.11979791
+ ,1.360090159
+ ,45.9577151
+ ,2.8079531
+ ,7.947593866
+ ,8.071039311
+ ,6.454250195
+ ,61.25096573
+ ,2.038822398
+ ,3.437673071
+ ,21.02243546
+ ,6.492329285
+ ,3.099361502
+ ,3.723876516
+ ,13.05988721
+ ,7.201885679
+ ,13.03962519
+ ,14.71813705
+ ,12.22706716
+ ,8.012906737
+ ,19.7240993
+ ,3.353902438
+ ,3.862711115
+ ,1.763453356
+ ,3.389291708
+ ,24.21336095
+ ,8.712897965
+ ,15.67216479
+ ,2.81760235
+ ,20.10899999
+ ,32.99783963
+ ,15.10295621
+ ,38.04642666
+ ,15.82539576
+ ,24.50336329
+ ,36.28979145
+ ,17.95651738
+ ,3.17242222
+ ,20.22346238
+ ,17.07301372
+ ,2.386059461
+ ,5.0429669
+ ,22.51255609
+ ,12.1945271
+ ,25.47503117
+ ,8.001388659
+ ,10.85369035
+ ,8.825689919
+ ,8.124307221
+ ,35.97971684
+ ,14.91412617
+ ,29.25935373
+ ,3.820067354
+ ,21.7287718
+ ,17.76647335
+ ,10.14087895)
> x <- c(15.3,19.2,31.3,0.5,24.8,45.7,88.8,22.1,3.9,30.4,8.6,31.6,12,30.8,15,30.3,17.2,32,19.9,6.2,31.2,0.6,11.9,10.4,36.4,1.1,22.6,25.5,7.3,13.5,11.9,22.5,16.3,1.6,8.5,35,8.3,24.8,9.2,5.5,1.3,8.9,0.7,21.7,7.8,10.7,19,10.2,12.5,14,21.7,15,31.8,0.7,7.3,4.9,21,6.2,14.7,23.1,1.3,10.7,9.9,58.2,8,24.1,7.3,12.7,18.8,3.5,5.9,5.1,0.1,3.8,17.3,8.6,6.5,7.6,15.6,1.3,4.9,8.1,13.4,6.6,5.8,12.6,10,3.5,6.4,11.5,17.3,12.5,3.5,1.5,7.3,17,5,13.1,0.5,1.9,0.5,4.7,6.2,34.2,4.2,14.6,7.1,1.5,0.4,7.7,11.6,1.2,1.6,1.5,0.9,2.8,4.3,54.8,1,2.8,5.5,2,0.4,0.5,6.4,2.7,1.4,8.9,1.1,1.4,1.4,1.1,0.8,1.6,0.6,4.6,0.6,1.2,1.8,0.4,1.1,1,5.1,0.8,0.7,0.6,1,0.7,1.2,0.5,1.6)
> ylab = 'Intentional homicide per 100,000'
> xlab = 'Guns per 100 residents'
> #'GNU S' R Code compiled by R2WASP v. 1.2.291 ()
> #Author: root
> #To cite this work: Wessa P., (2012), Pearson Correlation (v1.0.6) 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/fisher/rcomp/tmp/1arda1355256421.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', na.rm = T)
> cxy <- as.matrix(corxy$estimate)[1,1]
> 
> #Note: the /var/fisher/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab
> load(file="/var/fisher/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/fisher/rcomp/tmp/2fclk1355256421.tab") 
> 
> try(system("convert tmp/1arda1355256421.ps tmp/1arda1355256421.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  0.852   0.200   1.032