# Load librarys and mean shape function
library(shapes)
library(car)
library(heplots)

raw_landmarks<-read.csv("fundulus_gm.csv")
fish_config<-array(raw_landmarks[,2],dim=c(17,2,100))

# view raw data
fish_config
# view dimensions
dim(fish_config)

# Create a vector that lists the species for the samples
spvector<-c(1:100)
spvector[1:50]<-"F.olivaceus"
spvector[51:100]<-"F.notatus"
spvector<-as.factor(spvector)
spvector