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() }
|