shapepca function

Principal components analysis for shape

Provides graphical summaries of principal components for shape.

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.

Returns

No value is returned

References

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

Author(s)

Ian Dryden

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)