我可以用 R Cubist 结果制作图表吗？答案

【问题标题】：Can I make a diagram out of R Cubist outcome?我可以用 R Cubist 结果制作图表吗？
【发布时间】：2013-07-07 19:52:37
【问题描述】：

我刚开始使用 R 包 Cubist，它会创建如下输出：

Cubist [Release 2.07 GPL 版] 2013 年 7 月 9 日星期二 19:46:48

目标属性`结果'

从 undefined.data 中读取 260 个案例（9 个属性）

型号：

规则 1：[26 例，平均值 0.3，范围 0 到 8，est err 0.3]
if
    B4 <= 54.96766
    B7 > 39.66716
then
    outcome = 0
规则 2：[48 例，平均值 0.3，范围 0 到 8，est err 0.6]
if
    B1 > 56.99043
    B5 > 74.11118
    B6 > 155.996
then
    outcome = 9.1 - 0.17 B5 + 0.25 B7 + 0.15 NDVI

我可以使用 predict 来使用模型，但我想制作树的图形。我不认为可以查看手册，但我想知道是否有人知道这样做的方法。

【问题讨论】：

只是想澄清一下，我正在寻找这样的图表：stackoverflow.com/questions/17460595/…

标签： r diagram regression decision-tree

【解决方案1】：

rpart 和 partykit 都是用于绘制树形图的有用 R 包。

编辑：在此处查看派对工具包的示例here

【讨论】：

据我所知，rpart 不会在节点处输出回归公式，它只会输出平均值。我还没有研究过partykit，你有它可以制作的图形的链接吗？
This give an example of how to use it。虽然仔细阅读我认为它可能与 rpart 非常相似。在我的原始答案中添加了一个示例

【解决方案2】：

您可以使用dotplot.cubist 获得模型条件和系数的直观视图。举个例子：

library(mlbench)
library(Cubist)
library(gridExtra)
data(BostonHousing)
mod1 <- cubist(x = BostonHousing[, -14], y = BostonHousing$medv)
summary(mod1)
p1 <- dotplot(mod1, what = "splits",main='Conditions')
p2 <- dotplot(mod1, what = "coefs",main='Coefs')
grid.arrange(p1,p2)

【讨论】：

很好的答案，但我应该澄清我在寻找什么样的树形图。现在在问题下方的 cmets 中。