# NOT RUN {
data(gorilla) # load the data
riem.female = wrap.landmark(gorilla$female) # wrap as RIEMOBJ
opar <- par(no.readonly=TRUE)
for (i in 1:30){
if (i < 2){
plot(riem.female$data[[i]], cex=0.5,
xlim=c(-1,1)/2, ylim=c(-2,2)/5,
main="30 female gorilla skull preshapes",
xlab="dimension 1", ylab="dimension 2")
lines(riem.female$data[[i]])
} else {
points(riem.female$data[[i]], cex=0.5)
lines(riem.female$data[[i]])
}
}
par(opar)
# }
Run the code above in your browser using DataLab