if (requireNamespace("dbarts", quietly = TRUE)) {
# Load the dbarts package to access the bart function
library(dbarts)
# Get Data
df <- na.omit(airquality)
# Create Simple dbarts Model For Regression:
set.seed(1701)
dbartModel <- bart(df[2:6],
df[, 1],
ntree = 5,
keeptrees = TRUE,
nskip = 10,
ndpost = 10
)
# Tree Data
trees_data <- extractTreeData(model = dbartModel, data = df)
# Cretae Porximity Matrix
bmProx <- proximityMatrix(
trees = trees_data,
reorder = TRUE,
normalize = TRUE,
iter = 1
)
# MDS plot
mdsBart(
trees = trees_data, data = df, target = bmProx,
plotType = "interactive", level = 0.25, response = "Ozone"
)
}
Run the code above in your browser using DataLab