mcmccode_2D_A: function for MCMC sampling for 2D landmark data (Gaussian...

In debashischatterjee111/Bligom: Bayesian Analysis of Landmark based Geometric Morphometrics

Description

Draws posterior from 8 parameters "c1","c2", "b", "theta", "Sigma_0", "Sigma_1", "Sigma_2", "delta"

Usage

mcmccode_2D_A(myData, choice, Nlandmark, tune, Nsample)

Arguments

myData	3D array containing 2 dimensional landmark
choice	Which to compare with 1st object or mean obj
Nlandmark	How many landmarks are there for each object? (usually nrow(myData,,1))
tune	Tuning value of MCMC sampler
Nsample	Number of MCMC sample desired

Details

"c1","c2" = Location parameter #' "b" = dilation parameter "theta" =Rotation angle "Sigma_k" = Anisotropic error variation parameter (k=0,1,2) Note that, here we are assuming Isotropy of error variance

Value

matrix containing samples from posterior density of parameter

Examples

## Not run: 
require(shapes)
data(apes)
myData = apes$x
ape_ress_10000=mcmccode_2D_A(apes$x,10,nrow(myData[,,1]), tune=1,10)

## End(Not run)

debashischatterjee111/Bligom documentation built on Dec. 31, 2020, 11:15 p.m.