3

我想使用 pmml 库导出 Caret 随机森林模型，以便可以将其用于 Java 中的预测。这是我得到的错误的再现。

data(iris)
require(caret)
require(pmml)
rfGrid2 <- expand.grid(.mtry = c(1,2))
fitControl2 <- trainControl(
  method = "repeatedcv",
  number = NUMBER_OF_CV, 
  repeats = REPEATES)

model.Test <- train(Species ~ .,
  data = iris,
  method ="rf",
  trControl = fitControl2,
  ntree = NUMBER_OF_TREES,
  importance = TRUE,  
  tuneGrid = rfGrid2)

print(model.Test)
pmml(model.Test)

Error in UseMethod("pmml") : 
  no applicable method for 'pmml' applied to an object of class "c('train', 'train.formula')"

我在谷歌上搜索了一段时间，发现实际上很少有关于导出到 PMML 的信息，一般来说，pmml 库的随机森林位于：

methods(pmml)
 [1] pmml.ada          pmml.coxph        pmml.cv.glmnet    pmml.glm          pmml.hclust       pmml.itemsets     pmml.kmeans      
 [8] pmml.ksvm         pmml.lm           pmml.multinom     pmml.naiveBayes   pmml.nnet         pmml.randomForest pmml.rfsrc       
[15] pmml.rpart        pmml.rules        pmml.svm

它使用直接随机森林模型工作，但不是插入符号训练的模型。

library(randomForest)
iris.rf <- randomForest(Species ~ ., data=iris, ntree=20)
# Convert to pmml
pmml(iris.rf)
# this works!!!
str(iris.rf)

List of 19
 $ call           : language randomForest(formula = Species ~ ., data = iris, ntree = 20)
 $ type           : chr "classification"
 $ predicted      : Factor w/ 3 levels "setosa","versicolor",..: 1 1 1 1 1 1 1 1 1 1 ...
...

str(model.Test)
List of 22
 $ method      : chr "rf"
 $ modelInfo   :List of 14
  ..$ label     : chr "Random Forest"
  ..$ library   : chr "randomForest"
  ..$ loop      : NULL
  ..$ type      : chr [1:2] "Classification" "Regression"
...
4

2 回答 2

5

您不能使用或类型调用该pmml方法（即，这是您的对象的类型）。traintrain.formulamodel.Test

该方法的插入符号文档train说您可以访问最佳模型作为finalModel字段。pmml然后，您可以在该对象上调用该方法。

rf = model.Test$finalModel
pmml(rf)

不幸的是，Caret 使用“矩阵接口”（即通过设置x和y字段）指定 RF 模型，而不是使用更常见的“公式接口”（即通过设置formula字段）。AFAIK，“pmml”包不支持导出此类 RF 模型。

所以，看起来你最好的选择是使用两级方法。首先，使用 Caret 包为您的数据集找到最合适的 RF 参数化。其次，使用具有此参数化的“公式接口”手动训练最终的 RF 模型。

于 2014-12-11T20:08:04.807 回答
2

您可以使用该r2pmml软件包来完成这项工作：

library("caret")
library("r2pmml")

data(iris)

train.rf = train(Species ~ ., data = iris, method = "rf")
print(train.rf)
r2pmml(train.rf, "/tmp/train-rf.pmml")
于 2015-02-20T09:17:14.047 回答