Novel Bayesian Procrustes Variance-based Inferences in Geometric Morphometrics Novel R package: BPviGM1

by   Debashis Chatterjee, et al.

Compared to abundant classical statistics-based literature, to date, very little Bayesian literature exists on Procrustes shape analysis in Geometric Morphometrics, probably because of being a relatively new branch of statistical research and because of inherent computational difficulty associated with Bayesian analysis. Moreover, we may obtain a plethora of novel inferences from Bayesian Procrustes analysis of shape parameter distributions. In this paper, we propose to regard the posterior of Procrustes shape variance as morphological variability indicators. Here we propose novel Bayesian methodologies for Procrustes shape analysis based on landmark data's isotropic variance assumption and propose a Bayesian statistical test for model validation of new species discovery using morphological variation reflected in the posterior distribution of landmark-variance of objects studied under Geometric Morphometrics. We will consider Gaussian distribution-based and heavy-tailed t distribution-based models for Procrustes analysis. To date, we are not aware of any direct R package for Bayesian Procrustes analysis for landmark-based Geometric Morphometrics. Hence, we introduce a novel, simple R package BPviGM1 ("Bayesian Procrustes Variance-based inferences in Geometric Morphometrics 1"), which essentially contains the R code implementations of the computations for proposed models and methodologies, such as R function for Markov Chain Monte Carlo (MCMC) run for drawing samples from posterior of parameters of concern and R function for the proposed Bayesian test of model validation based on significance morphological variation. As an application, we can quantitatively show that primate male-face may be genetically viable to more shape-variation than the same for females.


page 9

page 14


A Bayesian level set method for an inverse medium scattering problem in acoustics

In this work, we are interested in the determination of the shape of the...

Comparison and Bayesian Estimation of Feature Allocations

Feature allocation models postulate a sampling distribution whose parame...

The Efficiency of Geometric Samplers for Exoplanet Transit Timing Variation Models

Transit timing variations (TTVs) are a valuable tool to determine the ma...

Dealing with Stochastic Volatility in Time Series Using the R Package stochvol

The R package stochvol provides a fully Bayesian implementation of heter...

Efficient Bayesian Sampling Using Normalizing Flows to Assist Markov Chain Monte Carlo Methods

Normalizing flows can generate complex target distributions and thus sho...

A New Wald Test for Hypothesis Testing Based on MCMC outputs

In this paper, a new and convenient χ^2 wald test based on MCMC outputs ...

A new method for faster and more accurate inference of species associations from novel community data

Joint Species Distribution models (jSDMs) explain spatial variation in c...