Anthropometry (version 1.1)

shapes3dMod: 3D shapes plot

Description

This function is a slight modification of the original shapes3d function of the shapes R package so that the resulting plot has customized title and axes. Specifically, the changing lines regarding the original function are those related to its argument axes3 when it is fixed to TRUE.

Usage

shapes3dMod(x,loop=0,type="p",color=2,joinline=c(1:1),axes3=FALSE,rglopen=TRUE,main=main)

Arguments

x

See shapes3d.

loop

See shapes3d.

type

See shapes3d.

color

See shapes3d.

joinline

See shapes3d.

axes3

See shapes3d.

rglopen

See shapes3d.

main

Allows us to give the plot a title if axes3=TRUE.

Value

• A device with the desired plot.

References

Dryden, I. L., (2012). shapes package. R Foundation for Statistical Computing, Vienna, Austria. Contributed package.

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

See Also

shapes3d

Examples

#CLUSTERING WOMEN ACCORDING TO THEIR SHAPE:
landmarks1 <- na.exclude(landmarks)
dim(landmarks1) 
#[1] 574 198 
(num.points <- (dim(landmarks1)[2]) / 3) 
#[1] 66
landmarks2 <- landmarks1[1:50,] #In the interests of simplicity of the computation involved.
(n <- dim(landmarks2)[1]) 
#[1] 100         
    
dg <- array(0,dim = c(num.points,3,n))
for(k in 1:n){            
 for(l in 1:3){            
  dg[,l,k] <-  as.matrix(as.vector(landmarks2[k,][seq(l,
                     dim(landmarks2)[2]+(l-1),by=3)]),ncol=1,byrow=T)
 }
}
 
Nclusters <- 3 ; Nsteps <- 5 ; niter <- 5 ; stopCr <- 0.0001
resLL <- LloydShapes(dg,Nclusters,Nsteps,niter,stopCr,FALSE,TRUE)
copt <- resLL$copt
shapes3dMod(copt[,,1], loop = 0, type = "p", color = 2, joinline = c(1:1), 
            axes3 = TRUE, rglopen = TRUE, main = "Mean shape cluster 1")