Displaying 1 result from an estimated 1 matches for "moment4".
Did you mean:
moment
2006 Aug 10
0
sn package - skew t - code for analytical expressions for first 4 moments
...n;
alpha=shape;
omega=sd(X);
delta=alpha/sqrt(1+alpha^2);
mu=delta*sqrt(df/pi)*gamma(0.5*(df-1))/gamma(0.5*df);
# the 4 first moments
moment1=xi+omega*mu;
moment2=xi^2 + 2*omega*mu*xi + omega^2 * df/(df-2);
moment3=xi^3 + 3*omega*mu*xi^2 + 3*omega^2*df/(df-2)*xi +
omega^3*mu*(3-delta^2)*df/(df-3);
moment4=xi^4 + 4*omega*mu*xi^3 + 6*omega^2*df/(df-2)*xi^2 +
4*omega^3*mu*(3-delta^2)*df/(df-3)*xi+omega^4*3*df^2/((df-2)*(df-4));
# the 4 useful measures
mean=moment1;
var=moment2;
skew=moment3/var^(3/2);
kurt=moment4/var^2 - 3;