Displaying 1 result from an estimated 1 matches for "jointplogp".
2006 Mar 31
1
mutual information for two time series
...formation at. I hope to use it to set up optimal bivariate
embeddings but also to measure the strength of the relationship between the
time series
mutual2<-function(datin,lmax) {
library(KernSmooth)
lmax<-lmax+1
lhplogp<-numeric(100)
rhplogp<-numeric(100)
mut<-numeric(lmax)
jointplogp<-matrix(0,nrow=100,ncol=100)
lh<-as.vector(datin[,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),trunc...