#R-code for Comparing Properties of Estimators for Continuous Uniform Distribution

set.seed(372)
theta=50
sample.size=10
no.of.samples=500
momvec=mlevec=adjmlevec=rep(NA,no.of.samples)

for (i in 1:no.of.samples){
obsy=runif(sample.size,0,theta)
print(obsy)
momvec[i]=2*mean(obsy)
mlevec[i]=max(obsy)
adjmlevec[i]=(((n+1)/n)*max(obsy))
}
print(momvec)
print(mlevec)


par(mfrow=c(2,1))
hist(momvec,breaks=20:75)
abline(v=50,lty=2)
mean(momvec)
var(momvec)
hist(mlevec,breaks=20:75)
abline(v=50,lty=2)
mean(mlevec)
var(mlevec)

par(mfrow=c(1,2))
hist(abs(momvec-theta),breaks=0:50,prob=T,ylim=c(0,.25))
hist(abs(mlevec-theta),breaks=0:50,prob=T,ylim=c(0,.25))

par(mfrow=c(2,1))
hist(adjmlevec,breaks=20:75)
abline(v=50,lty=2)
mean(adjmlevec)
var(adjmlevec)
hist(abs(adjmlevec-theta),breaks=0:50,prob=T,ylim=c(0,.25))