R 語言和本福德定律：程序執行中的警告

Question

誰能幫我？ 我的代碼有一些警告。 我該如何解決這些問題？ 輸入是從 Matlab 中提取的簡單向量 1*196608。

警告：

1: In empirical.distribution$dist.freq - benford.dist.freq :
  longer object length is not a multiple of shorter object length

2: In ((empirical.distribution$dist.freq - benford.dist.freq)^2)/benford.dist.freq :
  longer object length is not a multiple of shorter object length

3: In empirical.distribution$dist - benford.dist :
  longer object length is not a multiple of shorter object length

4: In as.data.table.list(x, keep.rownames = keep.rownames,  ... :
  Item 1 has 9 rows but longest item has 10; recycled with remainder

我的代碼：

library(benford.analysis)
library(R.matlab)
f=readMat('vectorimage.mat')
s=f$vectorimage
y=as.numeric(s);#convert to numeric type
d=benford(y,number.of.digits = 1)
print(d)
plot(d)

Answer 1

根據您的數據，我認為您的負值和正值都在 y 中變化，這可能是正態分布的，因此您可以修改 benford，如下所示：

benford(y, number.of.digits = 3, sign = "both", discrete=TRUE, round=3)