Nothing
M21p <-
function(RR, alpha.f=c(0.5,0.5), alpha.m=c(0.5,0.5)) {
x.m <- RR[1:3]
x.f <- RR[4:6]
a.f <- alpha.f[1] + 2*x.f[1] + x.f[2]
b.f <- alpha.f[2] + 2*x.f[3] + x.f[2]
a.m <- alpha.m[1] + 2*x.m[1] + x.m[2]
b.m <- alpha.m[2] + 2*x.m[3] + x.m[2]
dens <- lfactorial(sum(x.f)) + lfactorial(sum(x.m)) - sum(lfactorial(x.f)) - sum(lfactorial(x.m)) + (x.f[2]+x.m[2])*log(2) +
lgamma(sum(alpha.f)) + lgamma(sum(alpha.m)) - sum(lgamma(alpha.f)) - sum(lgamma(alpha.m)) +
lgamma(a.f) + lgamma(b.f) - lgamma(a.f+b.f) + lgamma(a.m) + lgamma(b.m) - lgamma(a.m+b.m)
dens <- exp(dens)
return(dens)
}
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.