1

我正在尝试通过R 包将 python 模块limmbo（https://github.com/HannahVMeyer/limmbo）与 R 一起使用。reticulate我已经成功安装limmbo了 Anaconda2。我现在正在尝试使用该函数limmbo$core$vdbootstrap$LiMMBo$runBootstrapCovarianceEstimation，如下面的代码所示。当我运行下面的代码时，我收到一个关于将 float64 转换为 integer64 的错误。

```{r}
library(reticulate)
import("limmbo") -> limmbo
```

然后我运行python代码：

```{python}
import numpy
from numpy.random import RandomState
from numpy.linalg import cholesky as chol
from limmbo.core.vdsimple import vd_reml
from limmbo.io.input import InputData
random = RandomState(15)
N = 100
S = 1000
P = 3
snps = (random.rand(N, S) < 0.2).astype(float)
kinship = numpy.dot(snps, snps.T) / float(10)
y  = random.randn(N, P)
pheno = numpy.dot(chol(kinship), y)
pheno_ID = [ 'PID{}'.format(x+1) for x in range(P)]
samples = [ 'SID{}'.format(x+1) for x in range(N)]
datainput = InputData()
datainput.addPhenotypes(phenotypes = pheno,
phenotype_ID = pheno_ID, pheno_samples = samples)
datainput.addRelatedness(relatedness = kinship,
relatedness_samples = samples)
```

当我尝试运行 R 函数时出现问题limmbo$core$vdbootstrap$LiMMBo$runBootstrapCovarianceEstimation：

```{r}
(limmbo$core$vdbootstrap$LiMMBo(py$datainput, timing = TRUE, iterations = 100, S = 2) -> foo)
limmbo$core$vdbootstrap$LiMMBo$runBootstrapCovarianceEstimation(foo, cpus = 1, seed = 12345)
```



Error in py_call_impl(callable, dots$args, dots$keywords) : 
  TypeError: Cannot cast array from dtype('float64') to dtype('int64') according to the rule 'safe'

Detailed traceback: 
  File "/Users/frederickboehm/anaconda2/lib/python2.7/site-packages/limmbo/core/vdbootstrap.py", line 96, in runBootstrapCovarianceEstimation
    minCooccurrence=minCooccurrence)
  File "/Users/frederickboehm/anaconda2/lib/python2.7/site-packages/limmbo/core/vdbootstrap.py", line 353, in __generateBootstrapMatrix
    rand_state = np.random.RandomState(seed)
  File "mtrand.pyx", line 644, in mtrand.RandomState.__init__
  File "mtrand.pyx", line 687, in mtrand.RandomState.seed
Calls: <Anonymous> ... eval -> eval -> <Anonymous> -> py_call_impl ->     .Call
Execution halted
4

2 回答 2

2

首先，通过 .import 导入你的 numpy 模块 np <- import("numpy", convert = FALSE)。

然后您可以int64使用reticulate::np_array(datainput, dtype = np$int64).

您可以在本教程中了解有关如何操作和创建数组的更多信息。

希望这可以帮助。

于 2018-03-22T22:08:48.247 回答
1

袁的教程（上面答案中的链接）包含允许我回答问题的建议。这是我修改后的 R 代码，到目前为止，它可以工作：

np <- import("numpy", convert = FALSE)
(limmbo$core$vdbootstrap$LiMMBo(datainput, timing = TRUE, iterations = np_array(10, dtype = "int64"), S = np_array(2, dtype = "int64")) -> foo)
limmbo$core$vdbootstrap$LiMMBo$runBootstrapCovarianceEstimation(foo, cpus = np$int(1), seed = np_array(1232, dtype = "int64"))
于 2018-03-23T01:15:53.133 回答