Displaying 1 result from an estimated 1 matches for "moment1".
Did you mean:
moment
2006 Aug 10
0
sn package - skew t - code for analytical expressions for first 4 moments
...is on
a multivariate skew-t distribution.
hope it is one day useful for someone else too.
bye
luc
# computing analytical moments for skew t of azzalini
xi=location;
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 use...