# AFS Kansas City 20-Aug-16

# Appropriately set the working directory before this
# This also ran library(FSA) which also provides alkIndivAge(), Summarize(), hist(), col2rgbt()
source("../scripts/ALK_Construction.R")

ls()
headtail(sp.len)

sp.len.mod <- alkIndivAge(ALK.obs,age~tl,data=sp.len)
headtail(sp.len.mod)

sp.comb <- rbind(sp.age,sp.len.mod)
str(sp.comb)

( agefreq <- xtabs(~age,data=sp.comb) )
prop.table(agefreq)

hist(~age,data=sp.comb,w=1,xlab="Age (yrs)")

( sp.sum <- Summarize(tl~age,data=sp.comb,digits=2) )
plot(tl~age,data=sp.comb,ylab=lblTL,xlab="Age (yrs)",pch=19,col=col2rgbt("black",1/10))
lines(mean~fact2num(age),data=sp.sum,col="blue",lwd=2)


# Script created at 2016-07-25 19:44:10