Nothing
SPU_E <- function(gamma, n1, n2, cov){
if(gamma %% 2 == 1){
e <- 0
}
if(gamma %% 2 == 0){
ga.half <- gamma/2
P <- 0
for(d in seq(0, ga.half, 1)){
P <- P + 1/(factorial(d)*factorial(ga.half - d)*
n1^d*n2^(ga.half - d))
}
e <- factorial(gamma)/2^ga.half*P*sum(diag(cov)^ga.half)
}
return(e)
}
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.