23.6 Classification trees with categorical explanatory variables

epilobium <- read.table("c:\\temp\\epilobium.txt",header=T)
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])