R/psip.bt.R

psip.bt <-
function(x,c1,M)
{
    x1 <- (x-M)/c1
    ivec1 <- (x1 < 0)
    ivec2 <- (x1 >  1)
    return(ivec1+(1-ivec1-ivec2)*((1-x1^2)^2+4*x*x1*(1-x1^2)/c1))
}
musto101/wilcox_R documentation built on May 23, 2019, 10:52 a.m.