2017-08-09 7 views
-1

Lorsque j'essaie d'appeler le paquetage BayesTree R de python avec des données simples, j'obtiens l'erreur "NA/NaN/Inf dans un appel de fonction étranger" même si toutes les données sont des nombres réels positifs.rpy2 conversion NumPy automatique NA/NaN/Inf dans une fonction étrangère appel

Source Code

import numpy as np 
# R interface for python 
import rpy2 
# For importing R packages 
from rpy2.robjects.packages import importr 
# Activate conversion from numpy to R 
import rpy2.robjects.numpy2ri 
rpy2.robjects.numpy2ri.activate() 

train_x_py = np.array([[0.0, 0.0], 
         [0.0, 1.0], 
         [1.0, 1.0]]) 

# Any 3-length float vector fails for training y 
train_y_py = np.array([1.0,2.0,3.0]) 

test_x_py = np.array([[0.2, 0.0], 
         [0.2, 0.2], 
         [1.0, 0.2]]) 

# Create R versions of the training and testing data 
train_x = rpy2.robjects.r.matrix(train_x_py, nrow=3, ncol=2) 
train_y = rpy2.robjects.vectors.FloatVector(train_y_py) 
test_x = rpy2.robjects.r.matrix(test_x_py, nrow=3, ncol=2) 

print(train_x) 
print(train_y) 
print(test_x) 

BayesTree = importr('BayesTree') 
response = BayesTree.bart(train_x, train_y, test_x, 
          verbose=False, ntree=100) 

# The 7th return value is the estimated response 
response = response[7] 

print(response) 

Code de sortie/Erreur

 [,1] [,2] 
[1,] 0 0 
[2,] 0 1 
[3,] 1 1 

[1] 1 2 3 

    [,1] [,2] 
[1,] 0.2 0.0 
[2,] 0.2 0.2 
[3,] 1.0 0.2 

Traceback (most recent call last): 
    File "broken_rpy2.py", line 32, in <module> 
    verbose=False, ntree=100) 
    File "/Library/Frameworks/Python.framework/Versions/3.6/lib/python3.6/site-packages/rpy2/robjects/functions.py", line 178, in __call__ 
    return super(SignatureTranslatedFunction, self).__call__(*args, **kwargs) 
    File "/Library/Frameworks/Python.framework/Versions/3.6/lib/python3.6/site-packages/rpy2/robjects/functions.py", line 106, in __call__ 
    res = super(Function, self).__call__(*new_args, **new_kwargs) 
rpy2.rinterface.RRuntimeError: Error in (function (x.train, y.train, x.test = matrix(0, 0, 0), sigest = NA, : 
    NA/NaN/Inf in foreign function call (arg 7) 

L'erreur sur la ligne 32 à laquelle il fait référence est:

response = BayesTree.bart(train_x, train_y, test_x, 
          verbose=False, ntree=100) 

Configuration du système

Système d'exploitation: Mac OS X Sierra 10.12.6

Python Version: Python 3.6.1

R Version: R 3.4.1

Packages Python : pip 9.0.1, rpy2 2.8.6, numpy 1.13.0

Question

Est-ce ma propre erreur de l'utilisateur, ou est-ce un bug dans rpy2?

+0

Correction, j'ai testé ce problème sans utiliser numpy. Tant que les matrices R produisent la même sortie de code, l'erreur se produit. –

+0

En éditant le titre dans une minute, j'ai reproduit exactement le même problème dans R. Cela se passe à l'intérieur du paquet BayesTree. Il n'y a pas de correctif sur le côté Python/Numpy/rpy2 des choses. –

+1

Au lieu de modifier le titre, postez simplement une réponse (ou ce que vous considérez être la solution) et acceptez-la. :) – MSeifert

Répondre

0

Ceci est un problème dans le paquet R "BayesTree". Vous pouvez reproduire le problème dans R directement avec le code suivant (en supposant que vous avez installé le paquet BayesTree).

train_x = matrix(c(0,0,1,0,1,1),nrow=3,ncol=2) 
train_y = as.vector(c(1,2,3)) 
test_x = matrix(c(.2,.2,1.,.0,.2,.2),nrow=3,ncol=2) 
result = bart(train_x,train_y,test_x,verbose=FALSE,ntree=100)