Tube-Coordinates and [Graphics:../Images/index_gr_452.gif]-formula

In this part, we first give an expression of the smooth surface in [Graphics:../Images/index_gr_453.gif] which specifies the boundary of the region of interest in the η-space. The tube-coordinate system is then defined as a pair of the coordinate system on the surface and the coordinate along the normal direction. The signed distance from the boundary is slightly modified for generalization, and named as  a modified signed distance characterized  by a coefficient vector [Graphics:../Images/index_gr_454.gif]. The [Graphics:../Images/index_gr_455.gif]-formula is derived as an expression of  the distribution function of the modified signed distance, which is obtained up to [Graphics:../Images/index_gr_456.gif] terms ignoring the error of [Graphics:../Images/index_gr_457.gif].

Startup

This section initializes the Mathematica session.

packages
error messages
distribution functions

Exponential family

This section summarizes the results for the exponential family derived in the previous part.

the expectation of the exponential of a polynomial function of the normal vector

Here we give the log of the expectation of  the exponential of [Graphics:../Images/index_gr_458.gif], where a1, a2, and a3 are of order [Graphics:../Images/index_gr_459.gif], and a4 is [Graphics:../Images/index_gr_460.gif]. [Graphics:../Images/index_gr_461.gif] is a multivariate normal random vector [Graphics:../Images/index_gr_462.gif]of dim-dimensions with mean [Graphics:../Images/index_gr_463.gif], and the identity covariance matrix.  [Graphics:../Images/index_gr_464.gif] is obtained up to [Graphics:../Images/index_gr_465.gif]terms.

define tensors
logeexppoly

the canonical form of the density function

Here we give the canonical form of log f(y;η) for the exponential family of distributions, and store it in logdensityy. The metric phi2eta=[Graphics:../Images/index_gr_493.gif]is also given here.

define tensors
logdensityy
phi2eta

Tube-coordinates

First, the expression of the surface is specified in the Taylor series. Then, the tangent vectors and the normal vector are obtained. The tube-coordinates (u,v) are defined and used instead of the η-parametrization. Here u is dim-1 dimensional vector specifying a point on the surface, and v is the signed distance. The density function [Graphics:../Images/index_gr_517.gif] is obtained from [Graphics:../Images/index_gr_518.gif], where the parameter value is specified as [Graphics:../Images/index_gr_519.gif] without losing the generality.

preliminary

Some functions and tensors are defined here.

indices
simplification functions
define tensors

the coordinates around the smooth surface

The surface is specified by [Graphics:../Images/index_gr_562.gif], [Graphics:../Images/index_gr_563.gif] and [Graphics:../Images/index_gr_564.gif].  They are stored in "foo1" and "foo2" or corresponding "rule1" and "rule2".  The region of interest is specified by [Graphics:../Images/index_gr_565.gif]. The tangent vectors are given by foo3=[Graphics:../Images/index_gr_566.gif] and  foo4=[Graphics:../Images/index_gr_567.gif], or corresponding "rule3" and "rule4".  We also obtain phi2bu= [Graphics:../Images/index_gr_568.gif][Graphics:../Images/index_gr_569.gif] , which is the metric in the tangent space. The elements of the normal vector are denoted as [Graphics:../Images/index_gr_570.gif], a=1,...,dim, which are given in foo15=[Graphics:../Images/index_gr_571.gif] and foo16= [Graphics:../Images/index_gr_572.gif], or in the corresponding "rule15" and "rule16". The reparametrization between η-coordinates and (u,v)-coordinates are specified by [Graphics:../Images/index_gr_573.gif], and given in foo21 = [Graphics:../Images/index_gr_574.gif] and foo22 = [Graphics:../Images/index_gr_575.gif], or in "rule21", "rule22", and "rule22b".

smooth surface
tangent vectors
the normal vector
(u,v)-coordinate system

change of variables

The Jacobian of the change of variables η↔(u,v) is [Graphics:../Images/index_gr_741.gif]. The asymptotic expression of  [Graphics:../Images/index_gr_742.gif] is obtained up to [Graphics:../Images/index_gr_743.gif] term in "logdetJ". The density function [Graphics:../Images/index_gr_744.gif] is obtained from [Graphics:../Images/index_gr_745.gif] as shown in "logdensityuv", where the parameter value is specified as [Graphics:../Images/index_gr_746.gif].

Jacobian
density function f(u,v|v0)

[Graphics:../Images/index_gr_891.gif]-formula

We modify the signed distance v to obtain a modified signed distance specified by [Graphics:../Images/index_gr_892.gif]. The density function of w and its cumulants are obtained up to [Graphics:../Images/index_gr_893.gif] terms. The distribution function of w is calculated by applying the Cornish-Fisher expansion to the cumulants of w. We would also take account of the scaling by the factor tau as well as the local coordinates at the projection in the below.

modified signed distance w

The inverse series of the modified signed distance specifies  [Graphics:../Images/index_gr_894.gif]. We have assumed that cbr[0] and cbr[2] are [Graphics:../Images/index_gr_895.gif] and cbr[1], cbr[3], and all br[a',r] are [Graphics:../Images/index_gr_896.gif]. The other coefficients are [Graphics:../Images/index_gr_897.gif] for r>=4. Then the same order applies to cr[r]'s. The modified signed distance w is characterized by the coefficient vector c=(cr[0],cr[1],cr[2],cr[3]) up to [Graphics:../Images/index_gr_898.gif] terms, since we can ignore the linear term in u as explained later. The change of variable is given in "rulevinuw" for v expressed in terms of u and w. The Jacobian is given in logjvw=log [Graphics:../Images/index_gr_899.gif]. The joint density of (u,w) is obtained by f(u,w|v0)=f(u,v(w,u)|v0)J , and log f(u,w|v0) is stored in "logdensityuw".  We then calculate  [Graphics:../Images/index_gr_900.gif] as an application of "logeexppoly" to logdensityuw, and stored in "logdensityw". In fact, the linear term in [Graphics:../Images/index_gr_901.gif], namely [Graphics:../Images/index_gr_902.gif] does not contribute to the argument for deriving the distribution function of w as seen in logdensityw. By using "logeexppoly" again, we obtain the cumulants of [Graphics:../Images/index_gr_903.gif] as shown in "cumulantw".

define the modified signed distance as a series of v.
density function of w
cumulants of w

distribution function of w

The Cornish-Fisher expansion for the standardized random variable is taken from Johnson and Kotz (1994) as shown in "cfexpx" below. We first obtain the same expansion for nonstandardized variable as shown in "cfexpw", and apply it to the cumulants of w. This gives the distribution function of w, and we obtain [Graphics:../Images/index_gr_1044.gif], where [Graphics:../Images/index_gr_1045.gif] is the quantile function of the standard normal distribution. The same expression, but without MathTensor notation, is also given in "zform". Finally, the scaling by the factor "tau" is applied to these results, and zc-formula is stored in "zformulatau" as well as in "zformtau".

Cornish-Fisher expansion  (p.66 of Johnson and Kotz 1994).
Cornish-Fisher expansion of w
zc-formula using a simplified notation
scaling by the factor "tau".

local coordinates at the projection

We consider a local coordinate Δη around a point η(u0,0) on the surface, where u0 indicates any specified value of u.  This will be used for u0 specifying the projection of y onto the surface. The change of variable η↔Δη  is specified by  [Graphics:../Images/index_gr_1156.gif] for each u0.  The expression for [Graphics:../Images/index_gr_1157.gif] is given in "rule93", and that for [Graphics:../Images/index_gr_1158.gif] is in "rule94". The surface is expressed as [Graphics:../Images/index_gr_1159.gif], where the coefficients are shown in foo101 and foo102. Next,  the expression for [Graphics:../Images/index_gr_1160.gif] is obtained and stored in foo114[ua,ub]. This is equated with [Graphics:../Images/index_gr_1161.gif], and the coefficients  [Graphics:../Images/index_gr_1162.gif], [Graphics:../Images/index_gr_1163.gif], and [Graphics:../Images/index_gr_1164.gif] are obtained in foo121. The inverse of the metric=[Graphics:../Images/index_gr_1165.gif]is in rule131, which is used for [Graphics:../Images/index_gr_1166.gif] . These conversion rules are summarized in "rulesproj". zc-formula evaluated at η(u0,0) is shown in zformulau0, and that for scaling tau is in zformulatauu0.

the expression for the surface in the local coordinates
the expressions for the potential derivatives
Geometric quantities at the projection
zc-formula


Converted by Mathematica      July 21, 2003