使用rpy2复现R的density()函数却得到不同输出
复现R中density()函数的问题
我尝试复现R中的density()函数,此前已参考过Stack Overflow上的相关帖子。我编写了Python代码(借助rpy2调用R的stats包),但运行结果与R原生代码的结果不一致;同时使用sklearn.neighbors.KernelDensity按相同参数配置也无法复现一致结果。
R原生代码
set.seed(1234) x <- rnorm(10) dx <- density(x, adjust=1)
Python代码(基于rpy2)
from rpy2 import robjects from rpy2.robjects.packages import importr from rpy2.robjects import vectors import numpy as np X = [-1.2070657, 0.2774292, 1.0844412, -2.3456977, 0.4291247, 0.5060559, -0.5747400, -0.5466319, -0.5644520, -0.8900378] print(" my X: ", X) print("\n") stats = importr("stats") column = vectors.IntVector(X) output = stats.density(column, adjust=1) x = np.array(output[0]) y = np.array(output[1]) print("dx$x from R : ", x[:10] ) print("\n") print("dx$y from R : ", y[:10])
内容的提问来源于stack exchange,提问作者Paga
相关产品推荐
相关产品推荐

