joint.like <-
function(theta, ys, m0,ma,null = FALSE,const=1){
#find the joint likelihood
if(null) {
ma<-m0
}
l=likelihoodpercond(theta, mu=ma, ys=ys,const)
prod(l)
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.