inst/doc/MoranI.R

### R code from vignette source 'MoranI.Rnw'

###################################################
### code chunk number 1: MoranI.Rnw:21-22
###################################################
options(width=60)


###################################################
### code chunk number 2: MoranI.Rnw:119-122
###################################################
body <- c(4.09434, 3.61092, 2.37024, 2.02815, -1.46968)
longevity <- c(4.74493, 3.3322, 3.3673, 2.89037, 2.30259)
names(body) <- names(longevity) <- c("Homo", "Pongo", "Macaca", "Ateles", "Galago")


###################################################
### code chunk number 3: MoranI.Rnw:128-134
###################################################
library(ape)
trnwk <- "((((Homo:0.21,Pongo:0.21):0.28,Macaca:0.49):0.13,Ateles:0.62)"
trnwk[2] <- ":0.38,Galago:1.00);"
tr <- read.tree(text = trnwk)
plot(tr)
axisPhylo()


###################################################
### code chunk number 4: MoranI.Rnw:140-142
###################################################
w <- 1/cophenetic(tr)
w


###################################################
### code chunk number 5: MoranI.Rnw:146-147
###################################################
diag(w) <- 0


###################################################
### code chunk number 6: MoranI.Rnw:151-152
###################################################
Moran.I(body, w)


###################################################
### code chunk number 7: MoranI.Rnw:172-173
###################################################
Moran.I(body, w, alt = "greater")


###################################################
### code chunk number 8: MoranI.Rnw:178-179
###################################################
Moran.I(longevity, w)


###################################################
### code chunk number 9: MoranI.Rnw:240-243
###################################################
data(carnivora)
carnivora$log10SW <- log10(carnivora$SW)
carnivora$log10FW <- log10(carnivora$FW)


###################################################
### code chunk number 10: MoranI.Rnw:247-250
###################################################
fm1.carn <- log10SW ~ Order/SuperFamily/Family/Genus
co1 <- correlogram.formula(fm1.carn, data = carnivora)
plot(co1)


###################################################
### code chunk number 11: MoranI.Rnw:265-268
###################################################
fm2.carn <- log10SW + log10FW ~ Order/SuperFamily/Family/Genus
co2 <- correlogram.formula(fm2.carn, data = carnivora)
print(plot(co2))


###################################################
### code chunk number 12: MoranI.Rnw:276-277
###################################################
plot(co2, FALSE)

Try the ape package in your browser

Any scripts or data that you put into this service are public.

ape documentation built on May 30, 2017, 6:11 a.m.