Если матрица, аргумент x
для cubist()
, кажется, требует атрибута dimnames
.
Настройка в R:
library(Cubist)
dt = mtcars
Z = dt[, 4]
X = dt[, 6]
X1 = dt[, 7]
Теперь сравните это (воспроизводит вашу ошибку):
> T = cbind(dt[, 6], dt[, 7])
> str(T)
num [1:32, 1:2] 2.62 2.88 2.32 3.21 3.44 ...
> cubist(x=T, y=Z, committees=10)
cubist code called exit with value 1
Error in strsplit(tmp, "\"")[[1]] : subscript out of bounds
против
> T = cbind(X, X1)
> str(T)
num [1:32, 1:2] 2.62 2.88 2.32 3.21 3.44 ...
- attr(*, "dimnames")=List of 2
..$ : NULL
..$ : chr [1:2] "X" "X1"
> cubist(x=T, y=Z, committees=10)
Call:
cubist.default(x = T, y = Z, committees = 10)
Number of samples: 32
Number of predictors: 2
Number of committees: 10
Number of rules per committee: 1, 1, 1, 1, 1, 1, 1, 1, 1, 1
Естьнесколько способов гарантировать, что dimnames присоединяются через rpy2
.Один из простых способов с вашим кодом - просто явно назвать переменные:
In [15]: T = r['cbind'](X=X,X1=X1)
In [16]: print(r['str'](T))
num [1:32, 1:2] 2.62 2.88 2.32 3.21 3.44 ...
- attr(*, "dimnames")=List of 2
..$ : NULL
..$ : chr [1:2] "X" "X1"
<rpy2.rinterface.NULLType object at 0x7f0d7c0f5608> [RTYPES.NILSXP]
In [17]: print(r['cubist'](x=T,y=Z,committees=10))
Call:
cubist.default(x = structure(c(2.62, 2.875, 2.32, 3.215, 3.44, 3.46,
205, 215, 230, 66, 52, 65, 97, 150, 150, 245, 175, 66, 91, 113, 264, 175,
335, 109), committees = 10L)
Number of samples: 32
Number of predictors: 2
Number of committees: 10
Number of rules per committee: 1, 1, 1, 1, 1, 1, 1, 1, 1, 1