install.packages("rpart")
library(rpart)
?rpart
data(iris)
iris
av1 <- rpart(Species~.,data=iris)
av1
install.packages("DMwR")
library(DMwR)
prettyTree(av1)
?rpart.control
?rpartXse
av2 <- rpart(Species~.,data=iris,control=rpart.control(minbucket=2))
av2
prettyTree(av2)
par(mfrow=c(2,1))
prettyTree(av1)
prettyTree(av2)
dev.off()
iris
150/5
teste <- sample(150,30)
teste
av3 <- rpart(Species~.,data=iris[-teste,],control=rpart.control(minbucket=2))
prettyTree(av3)
saida <- predict(av3,newdata=iris[teste,-5])
saida
?predict
help(rpart.predict)
help(predict.rpart)
saida2 <- predict(av3,newdata=iris[teste,-5],type="c")
saida2
iris[teste,5]
tabel(saida2,iris[teste,5])
table(saida2,iris[teste,5])
install.packages("MASS")
library(MASS)
x1 <- c(rnorm(100,4,2),rnorm(60,0,4))
x2 <- c(rnorm(100,3,2),rnorm(60,-2,2))
x1
grupo c(rep("+",100),rep("-",60))
grupo <- c(rep("+",100),rep("-",60))
grupo
?plot
plot(x1,x2,col=c(rep("blue",100),rep("red",60)))
l <- lda(grupo~x1+x2)
l
z1 <- c(rnorm(10,4,2),rnorm(10,0,4))
z2 <- c(rnorm(10,3,2),rnorm(10,-2,2))
saida <- predict(l,newdata=data.frame(x1=z1, x2=z2))
saida
saida2 <- predict(l,newdata=data.frame(x1=x1, x2=x2))
saida2
table(grupo,saida2$class)
qq <- qda(grupo~x1+x2)
saida3 <- predict(qq,newdata=data.frame(x1=x1, x2=x2))
table(grupo,saida3$class)