Pfratio2D: "Pfratio2D" function for MCMC posterior sampling

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

"Pfratio2D" function for MCMC posterior sampling

Description

2D landmark data in Pre-shape space (Gaussian likelihood with Isotropic Error Variance) as the name suggest, it evaluates fratio for two parameter vector t1, t2

Usage

Pfratio2D(t1, t2, myData)

Arguments

t1	An array containing 1 parameter sigma
t2	An array containing another 1 parameter sigma
myData	3D array containing 2 dimensional landmark

Value

fratio value

Examples

## Not run: 
require(shapes);
require(Morpho);
data(apes)
myData = apes$x
for(i in 1:dim(myData)[3])
{myData[,,i]=Morpho::rotonto(myData[,,1],myData[,,i])$yrot
}
for(k in 1:10)
{tempp2=runif(1, 0.1, 2);
tempp1=Ppurturb2D(tempp2,rep(0.01,1))
p=Pfratio2D (tempp1, tempp2,  myData)
print(p)
}

## End(Not run)

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