Hi,
I try to find a function that calculates the mutual information between
two arrays and i found mutualInfo() from bioDist package.
when i run the example:
x <- matrix(rnorm(100),nrow=5)
mutualInfo(x, nbin=3)
Then just to try out, i have assigned all the rest rows to the same
values as the first row.
x[2,]<-x[1,]
x[3,]<-x[1,]
x[4,]<-x[1,]
x[5,]<-x[1,]
> x
[,1] [,2] [,3] [,4] [,5] [,6]
[,7] [,8] [,9] [,10]
[1,] -0.7542452 0.8033849 1.513750 1.143666 -1.815935 -0.3781354
0.02026738 -0.8870625 0.3354863 -1.775121
[2,] -0.7542452 0.8033849 1.513750 1.143666 -1.815935 -0.3781354
0.02026738 -0.8870625 0.3354863 -1.775121
[3,] -0.7542452 0.8033849 1.513750 1.143666 -1.815935 -0.3781354
0.02026738 -0.8870625 0.3354863 -1.775121
[4,] -0.7542452 0.8033849 1.513750 1.143666 -1.815935 -0.3781354
0.02026738 -0.8870625 0.3354863 -1.775121
[5,] -0.7542452 0.8033849 1.513750 1.143666 -1.815935 -0.3781354
0.02026738 -0.8870625 0.3354863 -1.775121
[,11] [,12] [,13] [,14] [,15] [,16]
[,17] [,18] [,19] [,20]
[1,] -1.264028 0.9292861 1.592756 0.5829724 -0.004288089 0.03062974
0.3517299 0.4434983 -1.622165 -0.2684567
[2,] -1.264028 0.9292861 1.592756 0.5829724 -0.004288089 0.03062974
0.3517299 0.4434983 -1.622165 -0.2684567
[3,] -1.264028 0.9292861 1.592756 0.5829724 -0.004288089 0.03062974
0.3517299 0.4434983 -1.622165 -0.2684567
[4,] -1.264028 0.9292861 1.592756 0.5829724 -0.004288089 0.03062974
0.3517299 0.4434983 -1.622165 -0.2684567
[5,] -1.264028 0.9292861 1.592756 0.5829724 -0.004288089 0.03062974
0.3517299 0.4434983 -1.622165 -0.2684567
Then i expect the output of mutualInfo to be all zeros. However,
> mutualInfo(x, nbin=3, diag=T)
1 2 3 4 5
1 0.0000
2 1.0889 0.0000
3 1.0889 1.0889 0.0000
4 1.0889 1.0889 1.0889 0.0000
5 1.0889 1.0889 1.0889 1.0889 0.0000
Has anyone come across problem like this when using this function? Is
there any recommendations for other functions calculating pairwise
mutual information between two arrays?
Thanks a lot,
Best,
Runxuan
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
