newdt: newdt
In evoshape/deformetrics: R functions for shape analysis

Description Usage Arguments Value Author(s) Examples

This function performs the TPS Parallel Transport according to Varano et al (accepted). Deformation occurring within groups are estimate relatively to per-group local references ("locs" argument") and applied to a Common Reference ("CR").

1
2
3

newdt(array, group, CR = NULL, pole = CR, diffpole = F, locs = NULL,
  center = T, CSinit = F, tolrot = 1, tol = 1e-08, doopa = T,
  domopa = T)

`array`	numeric: an array kxmxn of landmark coordinates
`group`	character vector: a factor with group affiliation
`CR=NULL:`	if NULLL CR is estimated as the consensus of the array via Generalized Procrustes Analysis (GPA). Otherwise the specified CR is used.
`pole=CR;`	by default the Riemannian pole coincides with CR.
`diffpole=F`	logical: if TRUE the pole is not coincident with CR. This is not supposed to be used by common users.
`locs=NULL:`	if NULL the local references ar estimated via separate per-group GPAs. Oterwise the specifiesed locs are used.
`center=T`	logical: if TRUE shapes are centered.
`CSinit=F`	logical: if TRUE caling at unit size is performed.
`tolrot=1`	numeric: tolerance for rotation during GPA.
`tol=1e-8`	numeric: tolerance for eigenvalue decomposition
`doopa=T`	logical: if TRUE locs are aligned with CR via Ordinary Procrustes Analyisis (OPA)
`domopa=T`	logical: if TRUE each shape within each group is aligned with its proper loc via Modified Ordinary Procrustes Analysis introduced in Varano et al (accepted)

out matrix: a matrix containing the transported shapes.

Paolo Piras

## Not run:  
data(case1)
group<-factor(rep(1:5,each=21))
plottraj(procSym(case1,CSinit=T)$PCscores[,1:2],group,asp=1)######## a classic approach does not recover the cycle
dtcase1<-newdt(case1,group,CR=procSym(case1[,,firstsfac(group)],CSinit=F,reflect=F,scale=F,pcAlign=F)$mshape,locs=case1[,,firstsfac(group)],tolrot=10)
plottraj(procSym(dtcase1,CSinit=T)$PCscores[,1:2],group,asp=1)######## transported data recover the cycle

## End(Not run)