# NOT RUN {
library(randomForest)
data(iris)
iris.rf <- randomForest(iris[,-5], iris[,5], forest = TRUE)
# For a false positive rate of 1%
iris.sft <- sft(iris.rf, 0.01)
print(iris.sft)
# To iterate through a range of alpha values
alpha <- c(0.01,0.05, 0.1,0.15,0.2, 0.25)
threshold <- NULL
for(i in seq_along(alpha)){
threshold[i] <- sft(iris.rf, alpha[i])$sft
}
plot(alpha, threshold, type = 'b')
# }
Run the code above in your browser using DataLab