diff --git a/server.r b/server.r index 722fd34baa02328d8c3e9ed086a42f035f043851..e9bea18d932d636672d40ddb1fca82c5442c6ea9 100644 --- a/server.r +++ b/server.r @@ -47,7 +47,7 @@ set.seed(2018+as.numeric(input$submit_model)) population <- rexp(nrow*ncol,rate=lambda) population <- matrix(population,nrow,ncol) -par(mar=c(0.5, 0.5, 0.5, 0.5)) +#par(mar=c(0.5, 0.5, 0.5, 0.5)) x = cbind(rep(1:ncol, nrow), gl(nrow, ncol)) populationIndex = cbind(rep(1:ncol, nrow), gl(nrow, ncol)) @@ -99,8 +99,7 @@ xbar=c(xbar, as.numeric(mean(BootSample))) } cat(file=stderr(), " ", as.numeric(tx_star_b), " ", "\n") - - + CIboot=c(orig.mean-quantile(tx_star_bMINUStheta_fn,1-0.05/2),orig.mean-quantile(tx_star_bMINUStheta_fn,0.05/2)) CIboot2=c(2*orig.mean-quantile(tx_star_b,1-0.05/2),2*orig.mean-quantile(tx_star_b,0.05/2)) CItboot=c(orig.mean-sqrt(orig.var)*quantile(u_star_b,1-0.05/2)/sqrt(n),orig.mean-sqrt(orig.var)*quantile(u_star_b,0.05/2)/sqrt(n)) @@ -153,7 +152,10 @@ xbar=c() -par(mfrow=c(4,2)) +#par(mfrow=c(4,2)) +L <- matrix(1:6, ncol=2, byrow=TRUE) +layout(L) + mu=input$mu n=input$n @@ -268,16 +270,17 @@ u_star_b=c(u_star_b,sqrt(n)*(mean(population[BootSampleIndexMatrix])-mean(popula mycol2<-c(rep("gray",length(as.numeric(guessInput()))-1),"red") -plot(0,type="n",main="",ylab='',yaxt='n',xaxt='n',bty='n') -plot(0,type="n",main="",ylab='',yaxt='n',xaxt='n',bty='n') +#plot(0,type="n",main="",ylab='',yaxt='n',xaxt='n',bty='n') +#plot(0,type="n",main="",ylab='',yaxt='n',xaxt='n',bty='n') -par(mar=c(2, 3, 0, 1),mgp=c(5,1,0)) -hist(as.numeric(guessInput()),freq = TRUE,xlim=c(0,mu+100),main="",ylab='',xaxt='n',bty='n',cex.axis=2) +#par(mar=c(2, 3, 0, 1),mgp=c(5,1,0)) +par(mar=c(2,5,2,1)) +hist(as.numeric(guessInput()),freq = TRUE,xlim=c(0,mu+100),main="",xaxt='n',bty='n',cex.axis=2, cex.lab=2) points(as.numeric(guessInput()),rep(0,length(as.numeric(guessInput()))),cex.lab=2.2,xlim=c(0,mu+100), xlab='',ylab='',pch=17,col=mycol2,cex=1.4, xaxt='n',bty='n') axis(1, seq(0,mu+100,length=5), cex.axis=2, tck=-.01) -mtext("Frequency",at=-100,side=1,line =-22,cex=2,las=2) +#mtext("Frequency",at=-100,side=1,line =-2,cex=2,las=2) points(50,0,cex.lab=2.2,xlim=c(1,3),ylim=c(0,150),xlab='',ylab='',pch=15,col="blue",cex=1.4, yaxt='n', xaxt='n',bty='n') @@ -338,7 +341,3 @@ values4$u_star_b=c() values5$xbar=c() }) } - - - -