######################
##  Harjoitukset 5  ##
######################

#########
## 5.6 ##
#########

n<-10
x1<-rbinom(100,n,0.1)
x2<-rbinom(100,n,0.5)

lf<-function(theta,n,x){x*log(theta)+(n-x)*log(1-theta)}
max1<-lf(x1/n,n,x1)
max2<-lf(x2/n,n,x2)

max1[is.na(max1)]<-0
max2[is.na(max2)]<-0
# max1<-replace(max1,which(is.na(max1)==T),0)
# max2<-replace(max2,which(is.na(max2)==T),0)

D1<--2*(lf(0.1,n,x1)-max1)
D2<--2*(lf(0.5,n,x1)-max2)

par(mfrow=c(1,2))
hist(D1)
hist(D2)