if (FALSE) {
data(plethodon)
Y.gpa <- gpagen(plethodon$land) #GPA-alignment
A <- Y.gpa$coords
A[[1]] <- NA # make a missing value, just for example
mshape(Y.gpa$coords) # mean (consensus) configuration
# mshape(A, na.action = 1) # will return an error
mshape(A, na.action = 2) # returns NA in spot of missing value
mshape(A, na.action = 3) # finds mean values from all possible values
}
Run the code above in your browser using DataLab