| Stickleback data {paleoTS} | R Documentation |
Data of dorsal spine data from a fossil stickleback lineage.
data(dorsal.spines)
a paleoTS object for the time-series of dorsal spine counts
These data are already in the form of a paleoTS object, with sample means (mm), variances (vv)
sample sizes (nn), and ages in Kyr (tt).
Bell, M.A., M.P. Travis and D.M. Blouw 2006. Inferring natural selection in a fossil threespine stickleback. Paleobiology 32:562-577.
data(dorsal.spines)
# get subset of samples from invading lineage (tt>=4.5 Kyr), only those with nn>=5
ok<- dorsal.spines$tt >= 4.5 & dorsal.spines$nn >=5
ds2<- sub.paleoTS(dorsal.spines, ok=ok)
# convert time scale to generations (500 gen per Kyr)
ds2$tt<- ds2$tt*(1000/2)
plot(ds2, pool=TRUE)
# compare URW (drift) and OU (adaptive) models
m.urw<- opt.joint.URW(ds2, pool=TRUE)
m.ou<- opt.joint.OU(ds2, pool=TRUE)
aw<- akaike.wts(c(m.urw$AICc, m.ou$AICc))
cat ("Akaike wts: ", round(aw, 3),"\t", "[URW, OU]", "\n")