boxcount               package:RTisean               R Documentation

_R_e_n_y_i _e_n_t_r_o_p_y _e_s_t_i_m_a_t_e

_D_e_s_c_r_i_p_t_i_o_n:

     Estimates the Renyi entropy using a partition of the phase space.

_U_s_a_g_e:

     boxcount(series, l, x = 0, c, d = 1, M, Q = 2, R, r, scale = 20)

_A_r_g_u_m_e_n_t_s:

  series: a vector or a matrix. 

       l: number of data to use. 

       x: number of lines to be ignored. 

       c: column to be read.  

       d: delay for the delay vectors. 

       M: number of components, maximal embedding dimension. 

       Q: order of the entropy. 

       R: maximal length scale. 

       r: minimal length scale. 

   scale: number of epsilon values. 

_D_e_t_a_i_l_s:

     This function also can handle multivariate data, so that the phase
     space  is built of the components of the time series plus a
     temporal embedding, if desired.

_V_a_l_u_e:

     A list containing as many lists as the number of components, each
     containing as many matrices as the number of dimensions. Each
     matrix contains: the value of epsilon in the first column, the
     'Q'th order entropy (H_Q(dimension,epsilon)) in the second column
     and the 'Q'th order differential entropy (
     H_Q(dimension,epsilon)-H_Q(dimension-1,epsilon)) in the third
     column.

_R_e_f_e_r_e_n_c_e_s:

     http://www.mpipks-dresden.mpg.de/~tisean/

_S_e_e _A_l_s_o:

     'd2', 'c1'

_E_x_a_m_p_l_e_s:

     ## Not run: 

     dat <- henon(10000)
     boxout<- boxcount(dat,d=2)
     plot(boxout[[1]][,1],boxout[[1]][,2],ylim=c(0,8), t="l",
     xlab="Epsilon",ylab="Entropy",main="Renyi Entropy of Henon Data")
     lines(boxout[[2]][,1],boxout[[2]][,2],t="l",col=2)
     legend(1.5,7, c("Embedding dimension 1", "Embedding dimension 2"),fill=c(1,2), bty="n")

     ## End(Not run)

