477
477
For the transform, Stan uses the first part of an isometric log ratio
478
478
transform; see [ @egozcue +etal:2003] for the basic definitions and
479
479
Chapter 3 of [ @filzmoser +etal:2018] for the pivot coordinate version
480
- used here. Stan uses the isometric log ratio transform because it
481
- induces a geometry with zero correlation among the dimensions, making
482
- it easier for HMC to explore than simpler alternatives such as setting
483
- the final element to the negative sum of the first elements; see, e.g.,
484
- [ @seyboldt :2024] .
485
-
486
-
487
-
480
+ used here. Stan uses the isometric log ratio transform because it
481
+ results in equal variances of the the constrained sum to zero
482
+ vector see, e.g.,[ @seyboldt :2024] . Simpler alternatives, such as setting the
483
+ final element to the negative sum of the first elements, do not result in
484
+ equal variances. The $N - 1$ unconstrained parameters are independent, however,
485
+ the sum-to-zero constraint induces a negative correlation across the
486
+ constrained vector values.
488
487
489
488
### Zero sum transform {-}
490
489
@@ -499,7 +498,7 @@ S_N = 0
499
498
$$
500
499
and
501
500
$$
502
- y_N = -x_{N + 1} \cdot \frac{\ sqrt{N \cdot (N + 1) }{N}}.
501
+ y_N = -x_{N + 1} \cdot \sqrt{1 + \frac{1 }{N}}.
503
502
$$
504
503
The for each $n$ from $N - 1$ down to $1$, let
505
504
$$
0 commit comments