R/j.maad.R

`j.maad` <-
function (x) 
{

DNAME=deparse(substitute(x))

## Strip NAs
x<-na.omit(x)


### Robust Standard Deviation J
J<-sqrt(pi/2)*mean(abs(x-median(x))) 

# return(J)

paste("MAAD estimated J =", J, "for data", DNAME)

}

Try the lawstat package in your browser

Any scripts or data that you put into this service are public.

lawstat documentation built on Nov. 23, 2017, 5:05 p.m.