R version 2.9.0 (2009-04-17)
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(66.67,66.33,64.33,64,63.33,61.33,64.67,63,60.67,63.67,60.67,61.67,62.33,60.33,59.67,60.33,59.33,58.67,58.67,59.33,57.33,59.33,56,53.67,58.67,49.33,71.33,70.33,69,66,66,63.33,65.33,64.33,64,61.67,63.67,64.67,61.67,62,61.33,63.67,61.33,62.33,59.67,59.33,61.67,58.67,58,56.67,59.67,58,57,57.67,58.67,55.33,56,55.67,53.33,53.67,51,47,4.33,70,68.67,67.67,66,65.67,65.67,63.67,63.67,64,62,62,61.67,61.67,63.33,61,62.33,60.33,60.33,60.67,57.67,58.33,58,57.33,56.67,58,55.33,55.67,54.67,56.33,55,55,54.67,54.33,49,48.33,49.67,43.67,6.33,3,66.67,67.33,65.33,66,65.67,66.67,65.67,65,64.67,66.67,63.67,63.33,63.67,63.33,63.67,63,61.67,61.33,60.67,60,61.67,61.33,58.67,60.33,59.67,59.33,59.67,61,61,60,60,58.67,58.33,58,56.33,54.67,55.33,54,52.67,44,65.67,65,66.33,64,62.33,61.33,63,63.67,62,61.33,64.67,62.67,64,61,60.67,59.67,60.33,56.67,56.67,54.33,51,51,47,68,65,64,64,64,62,61,60,60,62,60,59,61,60,60,58,58,60,58,59,56,54,51,47)
> x <- c(70.8
+ ,69.6
+ ,69.87
+ ,67.47
+ ,67.6
+ ,67.13
+ ,66.27
+ ,66.73
+ ,68.07
+ ,67.8
+ ,64.8
+ ,64.6
+ ,64.2
+ ,64.2
+ ,63.67
+ ,61
+ ,59.67
+ ,59.67
+ ,59.8
+ ,60.73
+ ,59.4
+ ,58.07
+ ,57.47
+ ,70.73
+ ,72.87
+ ,66
+ ,66.07
+ ,66
+ ,66.27
+ ,64
+ ,63.67
+ ,63.73
+ ,63.33
+ ,63.53
+ ,63.53
+ ,62.87
+ ,59.53
+ ,62.8
+ ,60.8
+ ,59.8
+ ,56.67
+ ,57.67
+ ,58.4
+ ,55.47
+ ,56.2
+ ,71.53
+ ,68.67
+ ,65.67
+ ,66.73
+ ,67.33
+ ,66.73
+ ,66.87
+ ,65.8
+ ,64.73
+ ,65.47
+ ,63.6
+ ,64.07
+ ,64.67
+ ,63.73
+ ,62.53
+ ,61.93
+ ,62.67
+ ,62.8
+ ,61.33
+ ,62.6
+ ,59.13
+ ,61.27
+ ,59.47
+ ,57.87
+ ,59.73
+ ,61.4
+ ,58.8
+ ,58.33
+ ,57.47
+ ,57.13
+ ,55
+ ,51.53
+ ,72.73
+ ,73
+ ,70.8
+ ,70.07
+ ,71.67
+ ,71.07
+ ,70.67
+ ,70.73
+ ,70.73
+ ,68.6
+ ,69.6
+ ,66.47
+ ,67.07
+ ,68.67
+ ,66.93
+ ,65.93
+ ,68.87
+ ,66.53
+ ,65.8
+ ,66.6
+ ,66
+ ,65
+ ,66.8
+ ,65.6
+ ,66
+ ,65.67
+ ,64.67
+ ,65.07
+ ,64.67
+ ,65.07
+ ,65.2
+ ,64.87
+ ,63.47
+ ,62.6
+ ,64.07
+ ,63.73
+ ,64.67
+ ,61.6
+ ,61.6
+ ,60.47
+ ,61.27
+ ,63
+ ,61.47
+ ,60.87
+ ,61.67
+ ,62.87
+ ,62.4
+ ,59.73
+ ,60.13
+ ,58.8
+ ,59.6
+ ,58.93
+ ,60.13
+ ,58.2
+ ,58.27
+ ,58.27
+ ,55.07
+ ,53.87
+ ,52.33
+ ,47.2
+ ,37.93
+ ,72.73
+ ,70.07
+ ,70.67
+ ,72.07
+ ,68.8
+ ,68.8
+ ,67.47
+ ,66.73
+ ,66.53
+ ,66
+ ,67.6
+ ,66
+ ,66
+ ,66.53
+ ,65.8
+ ,64.27
+ ,64.67
+ ,64.6
+ ,64.13
+ ,65.47
+ ,62.93
+ ,63.53
+ ,62.13
+ ,63.87
+ ,64.67
+ ,63.33
+ ,63.13
+ ,62.8
+ ,62.4
+ ,62.4
+ ,62.6
+ ,61.47
+ ,62.2
+ ,63
+ ,61.8
+ ,59.73
+ ,60.33
+ ,60.13
+ ,59.53
+ ,59
+ ,55.93
+ ,41.87
+ ,36.33
+ ,71.67
+ ,71.47
+ ,70.47
+ ,69.53
+ ,70.73
+ ,69.93
+ ,68.73
+ ,67.53)
> par2 = 'less'
> par1 = 'pearson'
> ylab = 'No Placement'
> xlab = 'Placement'
> #'GNU S' R Code compiled by R2WASP v. 1.0.44 ()
> #Author: Dr. Ian E. Holliday
> #To cite this work: Ian E. Holliday, 2009, 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: 
> #Technical description: 
> postscript(file="/var/www/html/rcomp/tmp/1mguo1289306747.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
> plot(x,y, xlab=xlab, ylab=ylab)
> abline(coef=(lm(y ~ x))$coeff, col='red')
> ctp<-cor.test(x, y, use='pair', method='pearson')
> legend('bottomright', c('cor =', as.character(signif(ctp$estimate, digits=4))))
> dev.off()
null device 
          1 
> (sw1<-shapiro.test(x))

	Shapiro-Wilk normality test

data:  x 
W = 0.8979, p-value = 4.242e-10

> (sw2<-shapiro.test(y))

	Shapiro-Wilk normality test

data:  y 
W = 0.6249, p-value < 2.2e-16

> (cxy<-cor.test(x, y, method=par1))

	Pearson's product-moment correlation

data:  x and y 
t = -1.0387, df = 187, p-value = 0.3003
alternative hypothesis: true correlation is not equal to 0 
95 percent confidence interval:
 -0.2161336  0.0677226 
sample estimates:
        cor 
-0.07573978 

> 
> #Note: the /var/www/html/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab
> load(file="/var/www/html/rcomp/createtable")
> 
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a,paste('Shapiro Test (', xlab,')',sep=''),3,TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'Statistic W',header=TRUE)
> a<-table.element(a,'Probability',header=TRUE)
> a<-table.element(a,'Link',header=TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,signif(sw1$statistic, digits=4))
> a<-table.element(a,signif(sw1$p.value, digits=4))
> a<-table.element(a,hyperlink('http://en.wikipedia.org/wiki/Shapiro%E2%80%93Wilk_test','Shapiro-Wilk Test','Shapiro-Wilk Test'))
> a<-table.row.end(a)
> a<-table.end(a)
> table.save(a,file="/var/www/html/rcomp/tmp/27zbb1289306747.tab") 
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a, paste('Shapiro Test (', ylab,')',sep=''),3,TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'Statistic W',header=TRUE)
> a<-table.element(a,'probability',header=TRUE)
> a<-table.element(a,'Link',header=TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,signif(sw2$statistic, digits=4))
> a<-table.element(a,signif(sw2$p.value, digits=4))
> a<-table.element(a,hyperlink('http://en.wikipedia.org/wiki/Shapiro%E2%80%93Wilk_test','Shapiro-Wilk Test','Shapiro-Wilk Test'))
> a<-table.row.end(a)
> a<-table.end(a)
> table.save(a,file="/var/www/html/rcomp/tmp/30qsw1289306747.tab") 
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a,'Correlation Test X-Y',6,TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'Statistic',header=TRUE)
> a<-table.element(a,'Value',header=TRUE)
> a<-table.element(a,'probability',header=TRUE)
> a<-table.element(a,'Method',header=TRUE)
> a<-table.element(a,'alternative',header=TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,names(cxy$estimate))
> a<-table.element(a,signif(cxy$estimate, digits=4))
> a<-table.element(a,signif(cxy$p.value, digits=4))
> a<-table.element(a,cxy$method)
> a<-table.element(a,cxy$alternative)
> a<-table.row.end(a)
> a<-table.end(a)
> table.save(a,file="/var/www/html/rcomp/tmp/4389k1289306747.tab") 
> 
> try(system("convert tmp/1mguo1289306747.ps tmp/1mguo1289306747.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  0.414   0.172   1.050