reiser.2000.REML <- function(x,y,sigmaX,sigmaY){
muX <- mean(x)
muY <- mean(y)
delta.hat <- pnorm((muX-muY)/(sqrt(sigmaX^2+sigmaY^2)))
delta.hat
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.