attach(epilobium)
epilobium
model <- tree(species ~ .,epilobium,mindev=1e-6,minsize=2)
plot(model)
text(model,cex=0.7)
taxonomy <- read.table("c:\\temp\\taxonomy.txt",header=T)
attach(taxonomy)
names(taxonomy)
model1 <- tree(Taxon~., taxonomy)
plot(model1)
text(model1)
summary(model1)
print(model1)
model2 <- tree(Taxon~Sepal+Leaf,taxonomy);
partition.tree(model2)
label <- ifelse(Taxon=="I", "a",
ifelse(Taxon=="II","b",ifelse(Taxon=="III","c","d")))
text(Sepal,Leaf,label,col=1+as.numeric(factor(label)))
hump <- function(x,y){
library(tree)
model <- tree(y~x)
xs <- grep("[0-9]",model[[1]][[5]])
xv <- as.numeric(substring(model[[1]][[5]][xs],2,10))
xv <- xv[1:(length(xv)/2)]
xv <- c(min(x),sort(xv),max(x))
yv <- model[[1]][[4]][model[[1]][[1]]=="<leaf>"]
plot(x,y,col="red",pch=16,
xlab=deparse(substitute(x)),ylab=deparse(substitute(y)))
i <- 1
j <- 2
k <- 1
b <- 2*length(yv)+1
for (a in 1:b){
lines(c(xv[i],xv[j]),c(yv[k],yv[i]))
if (a %% 2 == 0){
j <- j+1
k <- k+1}
else{
i <- i+1
}}}
library(lattice)
attach(ethanol)
names(ethanol)
hump(E,NOx)
hump(E[E<1.007],NOx[E<1.007])
hump(E[E<1.006],NOx[E<1.006])