Я пытаюсь передать Numpy массив пакета GAMLSS в R., не созвучные массивы при переходе Numpy массива R через rpy2
import numpy as np
import rpy2.robjects as robjects
from rpy2.robjects import numpy2ri
numpy2ri.activate()
r = robjects.r
r.library("gamlss")
r.library("gamlss.mx")
L = r['data.frame'](np.array(np.random.normal(size=1000),
dtype=([('x', np.float), ('y', np.float), ('z', np.float)])))
r.gamlssMX(robjects.Formula('z~1'), data=L)
Запуск этого возвращает
Error in y0 - f0 : non-conformable arrays
Но я могу передать кадр данных в функцию линейной модели R.
lm = r.lm(robjects.Formula('x~y'), data=L)
print r.summary(lm.rx())
Я получил груз кода, который считывает двоичный файл в Python, но хотел бы использовать R пакет, следовательно, потребность в rpy2.
- EDIT -
В качестве примера в R:
x <- data.frame(z=c(rnorm(1000), rnorm(1000, mean=4)))
gamlssMX(z~1, K=1, data=x)