if (exists('par1')) { mylist <- list('aston2', 'a-101401634', 'a-100457812', 'a-102409703', 'a-100998212', 'a-101135359', 'a-101334088', 'a-101566014', 'a-101726641', 'a-105491925', 'a-105092722', 'a-100184457', 'a-92029778', 'a-102356832', 'a-100811922', 'a-104926413', 'a-101530532', 'a-103588056', 'a-109101705', 'a-100760723', 'a-101980911', 'a-91586304', 'a-101891606', 'a-101556956', 'a-101849951', 'a-92481068', 'a-100970919', 'a-101583952', 'a-102004287', 'a-101102247', 'a-102250482', 'a-119159116', 'a-101642385', 'a-102646168', 'a-100976184', 'a-101632847', 'a-101466161', 'a-100800100', 'a-102020117', 'a-102187409', 'a-102079665', 'a-109104511', 'a-102236497', 'a-101906988', 'a-100959040', 'a-101224516', 'a-102461468', 'a-102046634', 'a-103541992', 'a-101687906', 'a-105143781', 'a-101308045', 'a-101337768', 'a-109111645', 'a-119140055', 'a-31053516', 'a-102727504', 'a-102226834', 'a-102947412', 'a-101037695', 'a-101695543', 'a-100931378', 'a-109053895', 'a-94645367', 'a-100615744', 'a-100181087', 'a-109110659', 'a-101128720', 'a-102478101', 'a-103750826', 'a-93248615', 'a-101021155', 'a-93902728', 'a-105738646', 'a-101597524', 'a-100433218', 'a-100984455', 'a-100771660', 'a-101098261', 'a-101638137', 'a-101251615', 'a-104354801', 'a-101017424', 'a-103914037', 'a-101181943', 'a-100984341', 'a-102848832', 'a-100920512', 'a-102557521', 'a-101037190', 'a-101303224', 'a-90792337', 'a-101203533', 'a-71006521', 'a-102120864', 'a-101695314', 'a-102737846', 'a-103080273', 'a-102217957', 'a-119040124', 'a-119191161', 'a-103512844', 'a-119205011', 'a-103558071', 'a-102198623', 'a-101217688', 'a-103093505', 'a-100217755', 'a-103059350', 'a-101580205', 'a-101366953', 'a-102240896', 'a-103082093', 'a-101498760', 'a-93236461', 'a-109111221', 'a-102189388', 'a-119197026', 'a-102804436', 'a-105761172', 'a-103247717', 'a-119122958', 'a-100979473', 'a-103818551', 'a-103165798', 'a-83161650', 'a-102189735', 'a-91676465', 'a-119217689', 'Amanda') 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/rwasp_mcq20122.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/rwasp_mcq20122.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/rwasp_mcq20122.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/rwasp_mcq20122.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() }
|