Nothing
RefSigmaG <- function(sigma,Beta,X,y,delta,tol=0.0001,maxit=100,nitmon)
{
# Fixed point algorithm for scale
nit <- 1
repeat{
sigmao <- sigma
sigma <- ( RefAve2G(sigmao,Beta,X,y,delta)*sigmao^2/0.5 )^0.5; d <- sigma-sigmao
if (nit==maxit | abs(d)<tol) break
if(nitmon) cat(nit,sigma,"\n")
nit <- nit+1}
list(sigma=sigma,nit=nit)}
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.