p1 <- c(10,0.1, 0.1)
p2 <- c(10,0.4, 0.1)
data <- rdavies(100,p1)
likelihood(p1,data)
likelihood(p2,data) #should be smaller.
neg.log.likelihood(p1,rstupid(100)) #should be large negative.
Run the code above in your browser using DataLab