gamma1.to.lambda(gamma1)abs(gamma1)<0.5*(4-pi)*(2 (pi-2))^1.5<="" code="">,
which is about 0.99527.
If some values of gamma1 are not in the feasible region, a warning
message is issued, and NAs are returned.
See the reference below for the expression of the index of skewnnes
of a skew-normal distribution.
0.5*(4-pi)*(2>dsn