## G.3 PROBABILITY AND STATISTICS

### Refine

#### Document Type

- ZIB-Report (3)

#### Language

- English (3)

#### Has Fulltext

- yes (3)

#### Is part of the Bibliography

- no (3)

#### Keywords

For Kendallâ€™s shape space we determine analytically Jacobi fields and parallel transport, and compute geodesic regression. Using the derived expressions, we can fully leverage the geometry via Riemannian optimization and reduce the computational expense by several orders of magnitude. The methodology is demonstrated by performing a longitudinal statistical analysis of epidemiological shape data.
As application example we have chosen 3D shapes of knee bones, reconstructed from image data of the Osteoarthritis Initiative. Comparing subject groups with incident and developing osteoarthritis versus normal controls, we find clear differences in the temporal development of femur shapes. This paves the way for early prediction of incident knee osteoarthritis, using geometry data only.

In many applications one is interested to compute transition probabilities of a Markov chain.
This can be achieved by using Monte Carlo methods with local or global sampling points.
In this article, we analyze the error by the difference in the $L^2$ norm between the true transition probabilities and the approximation
achieved through a Monte Carlo method.
We give a formula for the error for Markov chains with locally computed sampling points. Further, in the case of reversible Markov chains, we will deduce a formula for the error when sampling points are computed globally.
We will see that in both cases the error itself can be approximated with Monte Carlo methods.
As a consequence of the result, we will derive surprising properties of reversible Markov chains.

In this article we present a new idea for approximating exit rates for diffusion processes living in a craggy landscape. We are especially interested in the exit rates of a process living in a metastable regions. Due to the fact that Monte Carlo simulations perform quite poor and are very computational expensive in this setting we create several similar situations with a smoothed potential. For this we introduce a new parameter $\lambda \in [0,1]$ ($\lambda = 1$ very smoothed potential, $\lambda=0$ original potential) into the potential which controls the influence the smoothing. We then sample the exit rate for different parameters $\lambda$ the exit rate from a given region. Due to the fact that $\lambda$ is connected to the exit rate we can use this dependency to approximate the real exit rate. The method can be seen as something between hyperdynamics and temperature accelerated MC.