如何計算 FNN package 中互信息的 p 值

Question

I have two columns data where I want to calculate the correlation from mutinfo() in FNN package in R ( https://www.rdocumentation.org/packages/FNN/versions/1.1.3.1/topics/mutinfo ) and then calculate the p 值。 在 R 中使用mutinfo()后如何計算 p 值？ 提前感謝您的幫助。

Answer 1

您可以使用maigesPack包中的bootstrapMI function：

此 function 采用數值矩陣（或兩個向量）並計算自舉（通過排列）p 值以測試互信息值是否等於零。 如果第一個參數是矩陣，則在矩陣的所有行對之間計算 p 值。

這是一個可重現的示例：

#BiocManager::install("maigesPack")
library(maigesPack)
library(FNN)
set.seed(1) # for reproducibility 
x <- runif(50, 0, 1)
y <- rbeta(50, 1, 2)
mutinfo(x, y, k = 2)  
#> [1] 0.04406588
bootstrapMI(x, y, bRep=100, ret="p-value")
#> [1] 0.31

^{由代表 package (v2.0.1) 於 2022 年 7 月 28 日創建}