shapes (version 1.2.5)

shapepca: 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)

Value

No value is returned

Details

The mean and PCs are plotted.

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

# NOT RUN {
#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)
# }