Displaying 1 result from an estimated 1 matches for "kdmat".
Did you mean:
dmat
2006 Mar 31
1
mutual information for two time series
...,1])
rh<-as.vector(datin[,2])
rhemb<-embed(rh,lmax)
lhc<-lh[1:length(rhemb[,1])]
for (i in 1:lmax) {
rhc<-rhemb[,i]
kd<-bkde2D(cbind(lhc,rhc),bandwidth=c(dpik(lhc),dpik(rhc)),gridsize=c(100,10
0),truncate=T)
#2d kernel density estimate in 2 dimensions)
kdmat<-as.matrix(kd$fhat)
for (j in 1:100){
p<-sum(kdmat[j,])
lhplogp[j]<-p*log(p)
}
lhplogp[is.nan(lhplogp)]<-0
hlh<- -1*sum(lhplogp)
#entropy of left hemsisphere
for (j in 1:1...