Displaying 1 result from an estimated 1 matches for "rhlab".
Did you mean:
hlab
2006 Mar 31
1
mutual information for two time series
...}
}
jointplogp[is.nan(jointplogp)]<-0
#sum of joint probabilities
hjoint<-sum(jointplogp)*-1
mut[i]<-hlh+hrh-hjoint
#classic defn of mutual information
maint<-paste("MI by Kernel Density is ",as.character(mut[i]))
rhlab<-paste("rh lag ",as.character(i))
xlabl<-"lh"
contour(kd$x1,kd$x2,kd$fhat,xlab=xlabl,ylab=rhlab,main=maint)
#dammit not producing realistic numbers
}
return(mut)
}
Thanks
Tim Parkin
--
Checked by AVG Free Edition.
[[alternative HTML version del...