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

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

"run2D" MCMC sampling function for running MCMC sampler

Description

2D landmark data(Gaussian likelihood with Isotropic Error Variance) Accepts new parameter 5*1 point with probability alpha

Usage

run2D(t, tune, myData, choice, Nsample)

Arguments

t	An array containing 5*1 parameters
tune	array containing 5*1 tuning value
myData	3D array containing 2 dimensional landmark
choice	Which to compare with 1st object or mean obj
Nsample	number of MCMC sample to generated

Details

"c1","c2" = Location parameter "b" = dilation parameter "theta" =Rotation angle "Sigma" = Isotropic error variation parameter

Value

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

Examples

## Not run: 
require(shapes)
data(apes)
myData = apes$x
r=run2D(rnorm(5, mean=1, sd=1),rep(0.1,5),
 myData, 10, 5000)
head(r); dim(r)

## End(Not run)

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