if (exists('par1')) {
mylist <- list('093497914', '091633583', '091545633', '092839155', '092500082', '090501706', '092211812', '093779474', '091042523', '081749120', '091786029', '090770058', '085544789', '104936319', '091559401', '093404783', '091048813', '092904486', '084134428', '091230212', '092184549', '091140287', '091733975', '093161037', '091440154', '094925487', '080764980', '091210373', '093014342', '090711109', '092058147', '093156181', '091315885', '093797124', '091696267', '092745188', '092222690', '091066183', '092344394', '091011699', '093290108', '092035274', '092253188', '091729602', '080921642', '090967625', '091694311', '092236815', '091420374', '091209399', '091395027', '092843222', '105704881', '081670488', '091466626', '090449969', '090931774', '092809217', '091651783', '085381472', '091925813', '092023662', '081971709', '092172580', '091658759', '091044697', '091167516', '092098936', '091565387', '092178744', '091348506', '083052093', '091305831', '091243533', '091610779', '090371132', '092028265', '092502307', '090799026', '090933435', '091461414', '091291523', '092150784', '091642792', '093930880', '082333704', '093341857', '092152973', '090754081', '092374438', '091412531', '105841647', '093345110', '090168679', '091212067', '095547972', '093721794', '083161650', '092523128', '093678258', '092141667', '063786965', '109076388', '109070560', '099112808', '109066046', '060240011', '109081223', '109091598', '109059945', '109014896', '109036393', '060942791', '071193407', '109036441', '109013095', '099108786', '093931658')
if (length(which(mylist == par1)) > 0) {
par1 <- paste(par1,'_',runif(1),'_',date(),sep='')
} else {
load(file='createtable')
a<-table.start()
a<-table.row.start(a)
a<-table.element(a,'Error',1,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'The User ID you entered does not exist! ')
a<-table.row.end(a)
a<-table.end(a)
table.save(a,file='mytable1a.tab')
stop('The User ID you entered does not exist!')
}
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
allsd<- c(sd(nx1), sd(nx2), sd(nx3))
names(allsd)<- c('Sample 1', 'Sample 2', 'Sample 3')
allsd
load(file='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('https://automated.biganalytics.eu/Ian.Holliday/rwasp_mcq20102.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('https://automated.biganalytics.eu/Ian.Holliday/rwasp_mcq20102.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('https://automated.biganalytics.eu/Ian.Holliday/rwasp_mcq20102.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('https://automated.biganalytics.eu/Ian.Holliday/rwasp_mcq20102.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='mytable1.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='mytable2.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='mytable3.tab')
bitmap(file='hist1.png')
hist(nx1)
dev.off()
bitmap(file='hist2.png')
hist(nx2)
dev.off()
bitmap(file='hist3.png')
hist(nx3)
dev.off()
bitmap(file='pairs.png')
pairs(x=cbind(nx1, nx2, nx3))
dev.off()
}
|
