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.

> par1 = 'a-101530532'
> par1 <- 'a-101530532'
> #'GNU S' R Code compiled by R2WASP v. 1.2.291 ()
> #Author: aston2
> #To cite this work: I. Holliday, 2012, MC Exam Q1 (v1.3) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_mcq20101.wasp/
> #Source of accompanying publication: 
> #
> library(digest)
> myran <- runif(1)
> if (myran <= 1) {
+ nx1 <- rnorm(100)
+ nx2 <- rnorm(100, sd=2.0, mean=0.0)
+ nx3 <- rnorm(100, sd=2.0, mean=-2.0)
+ answer <- 'nx1'
+ }
> if (myran < 0.66) {
+ nx2 <- rnorm(100)
+ nx3 <- rnorm(100, sd=2.0, mean=0.0)
+ nx1 <- rnorm(100, sd=2.0, mean=-2.0)
+ answer <- 'nx2'
+ }
> if (myran < 0.33) {
+ nx3 <- rnorm(100)
+ nx1 <- rnorm(100, sd=2.0, mean=0.0)
+ nx2 <- rnorm(100, sd=2.0, mean=-2.0)
+ answer <- 'nx3'
+ }
> answer <- digest(answer,ser=F)
> allmeans<- c(mean(nx1),mean(nx2),mean(nx3))
> names(allmeans)<-c('Sample 1', 'Sample 2', 'Sample 3')
> allmeans
   Sample 1    Sample 2    Sample 3 
-0.03616298 -2.21009170  0.18237323 
> allsd<- c(sd(nx1), sd(nx2), sd(nx3))
> names(allsd)<- c('Sample 1', 'Sample 2', 'Sample 3')
> allsd
 Sample 1  Sample 2  Sample 3 
1.9254320 2.1174114 0.8928418 
> 
> #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,'Identify the sample most similar to a standard normal  distribution.',1,TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,hyperlink(paste('http://www.wessa.net/rwasp_May_mock2.wasp?par1=',par1,'&par2=nx1&par3=',answer,'&par0=1',sep=''),'The sample called nx1','',target=''))
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,hyperlink(paste('http://www.wessa.net/rwasp_May_mock2.wasp?par1=',par1,'&par2=nx2&par3=',answer,'&par0=1',sep=''),'The sample called nx2','',target=''))
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,hyperlink(paste('http://www.wessa.net/rwasp_May_mock2.wasp?par1=',par1,'&par2=nx3&par3=',answer,'&par0=1',sep=''),'The sample called nx3','',target=''))
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,hyperlink(paste('http://www.wessa.net/rwasp_May_mock2.wasp?par1=',par1,'&par2=wnga&par3=',answer,'&par0=1',sep=''),'Avoid Guessing Penalty','',target=''))
> a<-table.row.end(a)
> a<-table.end(a)
> table.save(a,file="/var/wessaorg/rcomp/tmp/1c6kn1337529345.tab") 
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a,'Means of the three samples',4,TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'Sample',1,TRUE)
> a<-table.element(a,'nx1',1,TRUE)
> a<-table.element(a,'nx2',1,TRUE)
> a<-table.element(a,'nx3',1,TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'Value',1,TRUE)
> a<-table.element(a,round(mean(nx1),3))
> a<-table.element(a,round(mean(nx2),3))
> a<-table.element(a,round(mean(nx3),3))
> a<-table.row.end(a)
> a<-table.end(a)
> table.save(a,file="/var/wessaorg/rcomp/tmp/2csyb1337529345.tab") 
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a,'Standard deviations of the three samples',4,TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'Sample',1,TRUE)
> a<-table.element(a,'nx1',1,TRUE)
> a<-table.element(a,'nx2',1,TRUE)
> a<-table.element(a,'nx3',1,TRUE)
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'Value',1,TRUE)
> a<-table.element(a,round(sd(nx1),3))
> a<-table.element(a,round(sd(nx2),3))
> a<-table.element(a,round(sd(nx3),3))
> a<-table.row.end(a)
> a<-table.end(a)
> table.save(a,file="/var/wessaorg/rcomp/tmp/3l0n11337529345.tab") 
> postscript(file="/var/wessaorg/rcomp/tmp/4jdgx1337529345.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
> hist(nx1)
> dev.off()
null device 
          1 
> postscript(file="/var/wessaorg/rcomp/tmp/51lkk1337529345.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
> hist(nx2)
> dev.off()
null device 
          1 
> postscript(file="/var/wessaorg/rcomp/tmp/6xqaz1337529345.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
> hist(nx3)
> dev.off()
null device 
          1 
> postscript(file="/var/wessaorg/rcomp/tmp/76y3c1337529345.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
> pairs(x=cbind(nx1, nx2, nx3))
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/4jdgx1337529345.ps tmp/4jdgx1337529345.png",intern=TRUE))
character(0)
> try(system("convert tmp/51lkk1337529345.ps tmp/51lkk1337529345.png",intern=TRUE))
character(0)
> try(system("convert tmp/6xqaz1337529345.ps tmp/6xqaz1337529345.png",intern=TRUE))
character(0)
> try(system("convert tmp/76y3c1337529345.ps tmp/76y3c1337529345.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  1.096   0.270   1.362