.Rhistory

data(iris)
colMeans(iris)
iris
colMeans(iris[,1:4])
q()
n
library(ape)
mat <- vcv.phylo(compute.brlen(rtree(10)))
mat
q()
n
 library(ape)
tr <- compute.brlen(rtree(100))
write.tree(tr, file="tree.phy")
getwd()
q()
n
data(iris)
iris
?stack
?stack
q()
n
?stack
q()
n
library(ape)
tr = rtree(20)
tr
tr = compute.brlen(tr)
plot(tr)
write.tree(tr, "exampletree.phy")
q()
n
?cov2cor
cov2cor
 Is <- sqrt(1/diag(mat))
mat
q()
n
library(ape)
tree <- rtree(4)
vcv.phylo(tree)
cor(iris)
q()
n
library(xtable)
dat <- read.table("anolis.txt.csv", sep=" ")
dat
colnames(dat) <- rownames(dat)
dat
rownames(dat) <- dat[,1]
dat
dat <- dat[,-1]
dat
colnames(dat) <- rownames(dat)
dat
cor(dat)
 JL = jaw length, HW = head width, PECT = pectoral width, PELV = pelvic width, HUM = humerus, UL = ulna, FEM = femur, and TIB = tibia. 
rownames(dat) <- c("Jaw Length", "Head Width", "Pectoral Width", "Pelvic Width", "Humerus", "Ulna", "Femur", "Tibia")
dat
xtable(dat)
q()
n