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 |
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
1 2 3 4 5 6 7 8 | ## 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)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.