R/dnormt3.R

"dnormt3" <-
function(x, mean=0, sd=1){

if (sd==0)	{
	f <- rep(0, length(x))
	f[x==0] <- Inf
	}
else	{
	f <- sqrt(2)*normt3ip(mu=(x-mean)*sqrt(2)/sd, sigma=1)/sd
	sv <- abs(x - mean)*sqrt(2)/sd >= 37
	f[sv] <- 0
	}
f
}

Try the NORMT3 package in your browser

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

NORMT3 documentation built on April 25, 2022, 5:06 p.m.