Description: |
|
Auxiliary routine for rICfil:
calculates for dimension p>(=)2
diag(E[ YY' u min(b/|aIhY|,u) ])
and diag(E[ YY' min(b/|aIhY|,u)^2 ])
for u square root of a Chi^2_p-variable,
and Y~ufo(S_2) indep of u
by using a polar representation of
Lambda:= I^{1/2} Y u, u = | I^{-1/2} Lambda |,
Y=I^{-1/2} Lambda /u
The integrals are evaluated stepwise, first conditioning on Y
and calculated "analytically" using Ewinn, Ew2inn
and then the outer integration is done by MC-Integration
along the directions Y, parametrized by a
sin-cos-representation.
|