data(Laurasiatherian)
## empty object
new("multiphyDat")
## simple conversion with nicely ordered output
## Not run:
# genes <- list(gene1=subset(Laurasiatherian,,1:1600, FALSE),
# gene2=subset(Laurasiatherian,,1601:3179, FALSE))
# x <- new("multiphyDat", genes)
# x
# ## End(Not run)
## trickier conversion with missing sequences / wrong order
genes <- list(gene1=subset(Laurasiatherian,1:40),
gene2=subset(Laurasiatherian,8:47))
x <- new("multiphyDat", genes)
x
Run the code above in your browser using DataLab