pgls.print {mmodely}

Print the results of a PGLS model fit

Description

Print the results of a PGLS model fit

Usage

pgls.print(pgls, all.vars=names(pgls$data$data)[-1], 
           model.no=NA, mtx.out=NA, write=TRUE, print=FALSE)

Arguments

pgls

a fit PGLS model

all.vars

the names of all the variables to be reported

model.no

the model number (can be the order that models were run

mtx.out

should a matrix of the tabular summary results be returned

write

should the matrix of summary results be written to disk?

print

should the matrix of summary results be printed to screen?

Value

A matrix of summary results of a fit PGLS model

Examples

data.path <- system.file("extdata","primate-example.data.csv", package="mmodely")
data <- read.csv(data.path, row.names=1)
data$gn_sp <- rownames(data)

tree.path <- system.file("extdata","primate-springer.2012.tre", package="mmodely")

#5. RAxML phylogram based on the 61199 bp concatenation of 69 nuclear and ten mitochondrial genes.
phyl <- ape::read.tree(tree.path)[[5]]

phyl <- trim.phylo(phylo=phyl, gs.vect=data$gn_sp)

comp <- comp.data(phylo=phyl, df=data)

a.PGLS <- caper::pgls(formula('OC~mass.Kg + DPL.km'),  data=comp)

pgls.print(a.PGLS, all.vars=names(a.PGLS$data$data)[-1],
           model.no=NA, mtx.out='', write=FALSE, print=FALSE)