shapepca: Principal components analysis for shape

In shapes: Statistical Shape Analysis

Principal components analysis for shape

Description

Provides graphical summaries of principal components for shape.

Usage

shapepca(proc, pcno = c(1, 2, 3), type = "r", mag = 1, joinline = c(1, 1), 
           project=c(1,2),scores3d=FALSE,color=2,axes3=FALSE,rglopen=TRUE,zslice=0)

Arguments

proc	List given by the output from procGPA()
pcno	A vector of the PCs to be plotted
type	Options for the types of plot for the $m=2$ planar case: "r" : rows along PCs evaluated at c = -3,0,3 sd's along PC, "v" : vectors drawn from mean to +3 sd's along PC, "s" : plots along c= -3, -2, -1, 0, 1, 2, 3 superimposed, "m" : movie backward and forwards from -3 to +3 sd's along PC, "g" : TPS grid from mean to +3 sd's along PC.
mag	Magnification of the effect of the PC (scalar multiple of sd's)
joinline	A vector stating which landmarks are joined up by lines, e.g. joinline=c(1:n,1) will start at landmark 1, join to 2, ..., join to n, then re-join to landmark 1.
project	The default orthogonal projections if in higher than 2 dimensions
scores3d	Logical. If TRUE then a 3D scatterplot of the first 3 raw PC scores with labels in ‘pcno’ is given, instead of the default plot of the mean and PC vectors.
color	Color of the spheres used in plotting. Default color = 2 (red). If a vector is given then the points are colored in that order.
axes3	Logical. If TRUE then the axes are plotted in a 3D plot.
rglopen	Logical. If TRUE open a new RGL window, otherwise plot in current window.
zslice	For 3D case, type = "g": the z co-ordinate(s) for the grid slice(s)

Details

The mean and PCs are plotted.

Value

No value is returned

Author(s)

Ian Dryden

References

Dryden, I.L. and Mardia, K.V. (2016). Statistical Shape Analysis, with Applications in R (Second Edition). Wiley, Chichester. Chapter 7.

See Also

procGPA

Examples

#2d example
data(gorf.dat)
data(gorm.dat)

gorf<-procGPA(gorf.dat)
gorm<-procGPA(gorm.dat)
shapepca(gorf,type="r",mag=3)
shapepca(gorf,type="v",mag=3)
shapepca(gorm,type="r",mag=3)
shapepca(gorm,type="v",mag=3)

#3D example
#data(macm.dat)
#out<-procGPA(macm.dat)
#movie
#shapepca(out,pcno=1)

shapes documentation built on Feb. 16, 2023, 8:16 p.m.