[Top][All Lists]

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]


From: Przemek Klosowski
Subject: hurst.m
Date: Wed, 28 Sep 2005 16:26:08 -0400 (EDT)

I ran into this problem with 'hurst':

hurst is the user-defined function from the file

 -- Function File:  hurst (X)
     Estimate the Hurst parameter of sample X via the rescaled range
     statistic.  If X is a matrix, the parameter is estimated for every
     single column.

but there's a problem with the matrix: I get this for "hurst(rand(5))"

    error: operator -: nonconformant arguments (op1 is 5x5, op2 is 1x5)

I think it can be fixed by 

*** hurst.m     2005-09-28 16:12:59.000000000 -0400
--- /usr/share/octave/2.1.71/m/signal/hurst.m   2005-08-03 15:54:46.000000000 
*** 42,48 ****
    [xr, xc] = size (x);

    s = std (x);
!   w = cumsum (x - repmat(mean (x),xr,1));
    RS = (max(w) - min(w)) ./ s;
    H = log (RS) / log (xr);

--- 42,48 ----
    [xr, xc] = size (x);

    s = std (x);
!   w = cumsum (x - mean (x));
    RS = (max(w) - min(w)) ./ s;
    H = log (RS) / log (xr);

because mean() sums column-wise, and I think the intent of hurst()
is to process the results column by column. Please have a look at this
and forward the patch to John if it's OK.

reply via email to

[Prev in Thread] Current Thread [Next in Thread]