您的訂閱是我製作影片的動力
訂閱點這裡~
若內容有誤,還請留言指正,謝謝您的指教
#影片程式碼
data(iris)
iris <- subset(iris,Species!="setosa")
iris$Species <- gsub("virginica",1,iris$Species)
iris$Species <- gsub("versicolor",0,iris$Species)
iris$Species <- as.numeric(iris$Species)
library(rminer)
model<-glm(Species~Sepal.Length+Sepal.Width+Petal.Length+Petal.Width,
data=iris, family=binomial(link="logit"))
glm <- Importance(model,iris,measure="AAD")
glm$imp
par(mfrow=c(1,4))
boxplot(iris$Sepal.Length)$out
boxplot(iris$Sepal.Width)$out
boxplot(iris$Petal.Length)$out
boxplot(iris$Petal.Width)$out
#因為離群值不多,所以先不轉
n <- nrow(iris)
set.seed(1117)
subiris <- sample(seq_len(n), size = round(0.7 * n))
traindata <- iris[subiris,]
testdata <- iris[ - subiris,]
model<-glm(Species~Sepal.Length+Sepal.Width+Petal.Length+Petal.Width,
family=binomial(link="logit"), data=traindata)