পূর্ববর্তী উত্তরে যেমন এটি ইতিমধ্যে উল্লিখিত হয়েছে, রিগ্রেশন / রিগ্রেশন গাছের জন্য এলোমেলো বনগুলি ডেটা পয়েন্টগুলির জন্য প্রশিক্ষণের ডেটা রেঞ্জের পরিধি ছাড়িয়ে ডেটা পয়েন্টগুলির জন্য প্রত্যাশিত পূর্বাভাস দেয় না কারণ তারা এক্সট্রোপোলেট (ভাল) করতে পারে না। একটি রিগ্রেশন ট্রি নোডের একটি শ্রেণিবিন্যাস নিয়ে গঠিত, যেখানে প্রতিটি নোড একটি গুণ মানের উপর সঞ্চালনের জন্য একটি পরীক্ষা নির্দিষ্ট করে এবং প্রতিটি পাত (টার্মিনাল) নোড একটি পূর্বাভাসিত আউটপুট গণনা করার জন্য একটি বিধি নির্দিষ্ট করে। আপনার ক্ষেত্রে পরীক্ষার পর্যবেক্ষণ গাছের মধ্য দিয়ে গাছের পাতাগুলিতে প্রবাহিত হয়, উদাহরণস্বরূপ, "যদি x> 335, তবে y = 15", যা এলোমেলো বন দ্বারা গড় হয়।
এখানে একটি আর স্ক্রিপ্ট রয়েছে যা এলোমেলো বন এবং লিনিয়ার রিগ্রেশন উভয়ই পরিস্থিতিটির চিত্রিত করছে। এলোমেলো বনের ক্ষেত্রে, ডেটা পয়েন্টগুলি পরীক্ষা করার জন্য ভবিষ্যদ্বাণীগুলি ধ্রুবক হয় যা হয় সর্বনিম্ন প্রশিক্ষণের ডেটা এক্স-মানের নীচে বা সর্বোচ্চ প্রশিক্ষণের ডেটা এক্স-মানের উপরে।
library(datasets)
library(randomForest)
library(ggplot2)
library(ggthemes)
# Import mtcars (Motor Trend Car Road Tests) dataset
data(mtcars)
# Define training data
train_data = data.frame(
x = mtcars$hp, # Gross horsepower
y = mtcars$qsec) # 1/4 mile time
# Train random forest model for regression
random_forest <- randomForest(x = matrix(train_data$x),
y = matrix(train_data$y), ntree = 20)
# Train linear regression model using ordinary least squares (OLS) estimator
linear_regr <- lm(y ~ x, train_data)
# Create testing data
test_data = data.frame(x = seq(0, 400))
# Predict targets for testing data points
test_data$y_predicted_rf <- predict(random_forest, matrix(test_data$x))
test_data$y_predicted_linreg <- predict(linear_regr, test_data)
# Visualize
ggplot2::ggplot() +
# Training data points
ggplot2::geom_point(data = train_data, size = 2,
ggplot2::aes(x = x, y = y, color = "Training data")) +
# Random forest predictions
ggplot2::geom_line(data = test_data, size = 2, alpha = 0.7,
ggplot2::aes(x = x, y = y_predicted_rf,
color = "Predicted with random forest")) +
# Linear regression predictions
ggplot2::geom_line(data = test_data, size = 2, alpha = 0.7,
ggplot2::aes(x = x, y = y_predicted_linreg,
color = "Predicted with linear regression")) +
# Hide legend title, change legend location and add axis labels
ggplot2::theme(legend.title = element_blank(),
legend.position = "bottom") + labs(y = "1/4 mile time",
x = "Gross horsepower") +
ggthemes::scale_colour_colorblind()