n<-20
p1<-0.5
p2<-0.6
k <- seq(0, n, by = 1)
par(mfrow=c(3,1)) 
plot (k, dbinom(k, n, p1), type='l',col='red', ylab="density",
      main = "binomial likelihood; red: p=0.5; blue: p=0.6; black: ratio")
lines(k, (dbinom(k, n, p2)), col='blue')
plot(k, dbinom(k, n, p1)/dbinom(k, n, p2), ylab="ratio", main="likelihood
ratio for 0.5 vs 0.6") 
k <- seq(10, n, by = 1)
plot(k, dbinom(k, n, p1)/dbinom(k, n, p2), ylab="ratio", main="likelihood
ratio for 0.5 vs 0.6")