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.

> par4 = '12'
> par3 = '889201a3dd5cfc1ceae68e330da06215'
> par2 = 'p13'
> par1 = 'a-101530532'
> par4 <- '12'
> par3 <- '889201a3dd5cfc1ceae68e330da06215'
> par2 <- 'p13'
> par1 <- 'a-101530532'
> #'GNU S' R Code compiled by R2WASP v. 1.2.291 ()
> #Author: aston2
> #To cite this work: I. Holliday, 2012, MC Mock Q9 (v1.3.1) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/Ian.Holliday/rwasp_May_mock5.wasp/
> #Source of accompanying publication: 
> #
> library(digest)
> par4 <- as.numeric(par4)
> if(exists('par2')) {
+ if (par2 != 'wnga') {
+ if (digest(par2,ser=F) == par3) {
+ par4 = par4 + 4
+ lastanswer <- 'correct'
+ }
+ if (digest(par2,ser=F) != par3) {
+ par4 = par4 - 2
+ lastanswer <- 'Incorrect'
+ }
+ }
+ if (par2 == 'wnga') {
+ lastanswer <- 'Skipped Question'
+ }
+ }
> myran <- runif(1)
> x1<-rnorm(100)
> x2<-rnorm(100)
> x3<-rnorm(100)
> x4<-rnorm(100)
> if (myran <= 1) {
+ y1<-x1+0.75*x1*x1*x1+rnorm(100)
+ y2<-x2+1+rnorm(100)
+ y3<-x3+0.75*x3*x3+rnorm(100)
+ y4<-x4-0.75*x4*x4+rnorm(100)
+ answer <- 'y2'
+ }
> if (myran <= 0.75) {
+ y2<-x2+0.75*x2*x2*x2+rnorm(100)
+ y3<-x3+1+rnorm(100)
+ y4<-x4+0.75*x4*x4+rnorm(100)
+ y1<-x1-0.75*x1*x1+rnorm(100)
+ answer <- 'y3'
+ }
> if (myran <= 0.5) {
+ y3<-x3+0.75*x3*x3*x3+rnorm(100)
+ y4<-x4+1+rnorm(100)
+ y1<-x1+0.75*x1*x1+rnorm(100)
+ y2<-x2-0.75*x2*x2+rnorm(100)
+ answer <- 'y4'
+ }
> if (myran <= 0.25) {
+ y4<-x4+0.75*x4*x4*x4+rnorm(100)
+ y1<-x1+1+rnorm(100)
+ y2<-x2+0.75*x2*x2+rnorm(100)
+ y3<-x3-0.75*x3*x3+rnorm(100)
+ answer <- 'y1'
+ }
> answer <- digest(answer,ser=F)
> postscript(file="/var/wessaorg/rcomp/tmp/1v74g1337522301.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
> layout(matrix(c(1,2,3,4), 2,2))
> qqnorm(x1, main='X1')
> qqline(x1)
> qqnorm(x2, main='X2')
> qqline(x2)
> qqnorm(x3, main='X3')
> qqline(x3)
> qqnorm(x4, main='X4')
> qqline(x4)
> dev.off()
null device 
          1 
> postscript(file="/var/wessaorg/rcomp/tmp/27udp1337522301.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
> layout(matrix(c(1,2,3,4), 2,2))
> plot(x1,y1, main='y1 vs. x1')
> abline(reg=lm(y1~x1), col='red')
> mtext('Linear Regression')
> plot(x2,y2,  main='y2 vs. x2')
> abline(reg=lm(y2~x2), col='red')
> mtext('Linear Regression')
> plot(x3,y3,  main='y3 vs. x3')
> abline(reg=lm(y3~x3), col='red')
> mtext('Linear Regression')
> plot(x4,y4,  main='y4 vs. x4')
> abline(reg=lm(y4~x4), col='red')
> mtext('Linear Regression')
> dev.off()
null device 
          1 
> postscript(file="/var/wessaorg/rcomp/tmp/36kzw1337522301.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
> layout(matrix(c(1,2,3,4), 2,2))
> plot(x1, resid(lm(y1~x1)), main='X1, y1')
> abline(0,0, col='red')
> mtext('Regression Residuals')
> plot(x2, resid(lm(y2~x2)), main='X2, Y2')
> abline(0,0, col='red')
> mtext('Regression Residuals')
> plot(x3, resid(lm(y3~x3)), main='X3, Y3')
> abline(0,0, col='red')
> mtext('Regression Residuals')
> plot(x4, resid(lm(y4~x4)), main='X4, Y4')
> abline(0,0, col='red')
> mtext('Regression Residuals')
> dev.off()
null device 
          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,'Your last Question',1,TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,lastanswer,1,TRUE)
> a<-table.row.end(a)
> a<-table.end(a)
> table.save(a,file="/var/wessaorg/rcomp/tmp/4h7kz1337522301.tab") 
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a,'Based on the evidence below which pair of samples is best represented by a simple linear regression?',1,TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,hyperlink(paste('http://www.wessa.net/rwasp_May_mock6.wasp?par1=',par1,'&par2=y1&par3=',answer,'&par0=1&par4=',par4,sep=''),'X1, Y1','',target=''))
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,hyperlink(paste('http://www.wessa.net/rwasp_May_mock6.wasp?par1=',par1,'&par2=y2&par3=',answer,'&par0=1&par4=',par4,sep=''),'X2, Y2','',target=''))
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,hyperlink(paste('http://www.wessa.net/rwasp_May_mock6.wasp?par1=',par1,'&par2=y3&par3=',answer,'&par0=1&par4=',par4,sep=''),'X3, Y3','',target=''))
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,hyperlink(paste('http://www.wessa.net/rwasp_May_mock6.wasp?par1=',par1,'&par2=y4&par3=',answer,'&par0=1&par4=',par4,sep=''),'X4, Y4','',target=''))
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,hyperlink(paste('http://www.wessa.net/rwasp_May_mock6.wasp?par1=',par1,'&par2=wnga&par3=',answer,'&par0=1&par4=',par4,sep=''),'Avoid Guessing Penalty','',target=''))
> a<-table.row.end(a)
> a<-table.end(a)
> table.save(a,file="/var/wessaorg/rcomp/tmp/5ouo71337522301.tab") 
> 
> try(system("convert tmp/1v74g1337522301.ps tmp/1v74g1337522301.png",intern=TRUE))
character(0)
> try(system("convert tmp/27udp1337522301.ps tmp/27udp1337522301.png",intern=TRUE))
character(0)
> try(system("convert tmp/36kzw1337522301.ps tmp/36kzw1337522301.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  1.192   0.207   1.395