shapes (version 1.0-2)

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

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 backw

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.

Value

• No value is returned

Details

For $m=3$ the mean and PCs are plotted with orthogonal projections. The display is in the type="m" format only at the moment.

References

Dryden, I.L. and Mardia, K.V. (1998) Statistical Shape Analysis. Wiley, Chichester.

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)