Stickleback data from Bell et al. (2006)

Description

Dorsal spine, pterygiophore, and pelvic score data from a fossil stickleback lineage.

Usage

1
2
3

Format

paleoTS objects for time-series of stickleback data.

Details

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).

Source

Bell, M.A., M.P. Travis and D.M. Blouw 2006. Inferring natural selection in a fossil threespine stickleback. Paleobiology 32:562-577.

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
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), reset starting time to t=0
ds2$tt<- ds2$tt*(1000/2)	
ds2$tt<- ds2$tt - ds2$tt[1]

# compare URW (drift) and OU (adaptive) models
m.urw<- opt.joint.URW(ds2, pool=TRUE)
m.ou<- opt.joint.OU(ds2, pool=TRUE)
compareModels(m.urw, m.ou)

plot(ds2, modelFit=m.ou, pool=TRUE)

Want to suggest features or report bugs for rdrr.io? Use the GitHub issue tracker.