Prun2D: "Prun2D" MCMC sampling function for running MCMC sampler

In debashischatterjee111/BPviGM1: Novel Bayesian Procrustes Variance-based Inferences in Geometric Morphometrics

"Prun2D" MCMC sampling function for running MCMC sampler

Description

2D landmark data in Pre Shape Space (Gaussian likelihood with Isotropic Error Variance) Accepts new parameter 5*1 point with probability alpha "Sigma" = Isotropic error variation parameter

Usage

Prun2D(t, tune, myData, Nsample)

Arguments

t	An array containing 1 parameters
tune	array containing 1 tuning value
myData	3D array containing 2 dimensional landmark
Nsample	number of MCMC sample to generated

Value

Returns a large matrix containing the actual samples from posterior of 1 parameter point

Examples

## Not run: 
require(shapes)
data(apes)
myData = apes$x
for(i in 1:dim(myData)[3])
{myData[,,i]=
Morpho::rotonto(myData[,,1],myData[,,i])$yrot}
r=Prun2D(runif(1, 0.1,1),rep(0.01,1),
 myData, 5000)
head(r); dim(r)

## End(Not run)

debashischatterjee111/BPviGM1 documentation built on April 8, 2023, 7:28 p.m.