这是我不久前提出的上一个问题的后续,该问题最近得到了回答。
我已经建立了几个gbm
模型dismo::gbm.step
,这些模型依赖gbm
于 R 包中的拟合函数gbm
,以及 R 包中的交叉验证工具splines
。
作为我分析的一部分,我想使用 R 中可用的一些图形工具(例如透视图)来可视化数据中的成对交互。gbm
和dismo
包都具有检测和建模数据交互的功能。
dismo
Elith et 中解释了其中的实现。al ( 2008 ) 并返回一个统计数据,该统计数据表明模型预测偏离预测变量的线性组合,同时将所有其他预测变量保持在其平均值。
中的实现gbm
使用 Friedman 的 H 统计量 ( Friedman & Popescue, 2005 ),并返回不同的度量标准,并且也没有按其平均值设置其他变量。
建模和绘制的交互dismo::gbm.interactions
非常棒,并且提供了非常丰富的信息。但是,我也想使用gbm::interact.gbm
,部分是为了发表强度,也是为了比较两种方法的结果。
如果我尝试gbm::interact.gbm
在gbm.object
created with中运行dismo
,则会返回错误...</p>
"Error in is.factor(data[, x$var.names[j]]) :
argument "data" is missing, with no default"
我知道dismo::gmb.step
添加了作者认为对gbm model
.
我也明白我的问题的答案在源代码的某个地方。
我的问题是...
是否可以修改在
gbm
中创建的对象dismo
以在中使用gbm::gbm.interact
?如果是这样,这将通过......一个。修改
gbm
创建的对象dismo::gbm.step
?湾。修改源代码
gbm::interact.gbm
?C。做别的事?
我将通过源代码尝试自己解决这个问题,如果我在任何人回答之前提出解决方案,我将回答我自己的问题。