Nothing
#### computes the approximate likelihood for given theta, my, sigma
#### data in via y, P and M
compute.approx.lr <-
function(theta, mu, sigma, y, Py, M){
### precompute some things
n = length(y)
-M*theta - 0.5*M*theta^2 - n*log(sigma) - 0.5*sum((y -theta*Py - mu)^2)/sigma^2
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.