Преобразование модели rpart в PMML (с использованием пакета 'pmml')

Question

При попытке преобразовать мою rpart модель в pmml я получаю следующую ошибку:

Fehler in if (ff$nsurrogate[parent_ii] > 0) { :   
    Fehlender Wert, wo TRUE/FALSE nötig ist
    (Missing value where TRUE / FALSE is needed)

Эта ошибка может быть воспроизведена с помощью кода ниже:

library(rpart)
library(pmml)
df <- structure(list(a = structure(c(15L, 1L, 13L, 8L, 11L, 25L, 6L, 
                                     24L, 27L, 9L, 2L, 18L, 28L, 14L, 5L, 17L, 20L, 21L, 16L, 7L, 
                                     22L, 19L, 23L, 26L, 3L, 10L, 12L, 4L), .Label = c("013", "018", 
                                                                                       "063", "073", "122", "173", "212", "216", "296", "355", "410", 
                                                                                       "415", "423", "428", "453", "481", "534", "586", "678", "701", 
                                                                                       "735", "746", "778", "812", "818", "855", "864", "998"), class = "factor"), 
                     y = c(1.029993, 0.95987, 0.95987, 0.95987, 0.95987, 0.95987, 
                           0.95987, 0.969903, 0.95987, 0.860644, 0.95987, 0.969903, 
                           0.900669, 0.95987, 0.95987, 0.95987, 1.12018, 0.95987, 0.95987, 
                           0.95987, 0.95987, 0.880656, 0.95987, 0.939858, 0.95987, 0.939858, 
                           0.95987, 0.95987)), row.names = c(NA, -28L), class = "data.frame")

model <- rpart(y ~ a, df, control = rpart.control(minsplit = 1, minbucket = 2, cp=-1))
pmml.rpart(model)

sessionInfo () R версия 3.5.1 (2018-07-02) Платформа: x86_64-redhat-linux-gnu (64-разрядная версия) Работает под: Red Hat Enterprise Linux

Продукты матрицы:BLAS / LAPACK по умолчанию: /usr/lib64/R/lib/libRblas.so

языковой стандарт: [1] LC_CTYPE = de_DE.UTF-8 LC_NUMERIC = C LC_TIME = de_DE.UTF-8 LC_COLLATE = de_DE.UTF-8 LC_MONETARY = de_DE.UTF-8
[6] LC_MESSAGES = de_DE.UTF-8 LC_PAPER = de_DE.UTF-8 LC_NAME = C LC_ADDRESS = C LC_TELEPHONE = C
[11] LC_MEASUREMENT = de_DE.F8 LC_IDENTIFICATION = C

прикрепленные базовые пакеты: [1] stats graphics grDevices использует наборы методов базы данных

другие присоединенные пакеты: [1] pmml_1.5.7 XML_3.98-1.16 rpart_4.1-13

загружено через пространство имен (и не подключено): [1] compiler_3.5.1 magrittr_1.5 tools_3.5.1 yaml_2.2.0 stringi_1.2.4 stringr_1.3.1

Answer 1

Попробуйте вместо этого использовать пакет r2pmml: https://github.com/jpmml/r2pmml

Преобразование успешно выполняется с приведенным выше кодом как есть, и сгенерированный файл модели PMML меньше, чище и достоверно корректен:

library("r2pmml")
r2pmml(model, "model.pmml")

Answer 2

В настоящее время df$a является фактором, который на самом деле не имеет смысла, учитывая, что количество строк равно числу уровней фактора.Исправление этого с

df$a <- as.numeric(as.character(df$a))

также позволяет запускать

pmml.rpart(model)