suppressMessages(library("wsrf"))
suppressMessages(library("randomForest"))
# prepare parameters
ds <- iris
target <- "Species"
vars <- names(ds)
if (sum(is.na(ds[vars]))) ds[vars] <- na.roughfix(ds[vars])
ds[target] <- as.factor(ds[[target]])
(form <- as.formula(paste(target, "~ .")))
set.seed(500)
length(train <- sample(nrow(ds), 0.7*nrow(ds)))
length(test <- setdiff(seq_len(nrow(ds)), train))
# build model
model.wsrf <- wsrf(form, data=ds[train, vars], parallel=FALSE)
model.wsrf.nw <- wsrf(form, data=ds[train, vars], weights=FALSE, parallel=FALSE)
model.wsrf.nw.vi <- wsrf(form, data=ds[train, vars], weights=FALSE, importance=TRUE, parallel=FALSE)
model.subset <- subset.wsrf(model.wsrf, 1:200)
model.combine <- combine.wsrf(model.wsrf, model.wsrf.nw)
# evaluate
# Note:
# 32bit system and 64bit system will have different results, however,
# if random seed is fixed, the same results will be presented in the
# same system.
cl <- predict(model.wsrf, newdata=ds[test, vars], type="class")$class
cl.nw <- predict(model.wsrf.nw, newdata=ds[test, vars], type="class")$class
cl.subset <- predict(model.subset, newdata=ds[test, vars], type="class")$class
cl.combine <- predict(model.combine, newdata=ds[test, vars], type="class")$class
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.