[MSNoise] (no subject)

Thu Mar 27 19:51:12 UTC 2014

Alberto,

Le 26/03/2014 19:26, Oscar Alberto Castro Artola a écrit :
> Hello Thomas,
>
> all my data have 1 sps so I changed the sampling_rate value to 1.
Aahaaaa, that can explain the problems... 10 seconds window gives 10 
samples for the FFT, padded to 16 (nearest power of 2), then, you have 
16 points between 0.0 and 0.5 Hz in the FFT:
In [2]: numpy.linspace(0,0.5,16)
Out[2]:
array([ 0.        ,  0.03333333,  0.06666667,  0.1       , 0.13333333,
         0.16666667,  0.2       ,  0.23333333,  0.26666667, 0.3       ,
         0.33333333,  0.36666667,  0.4       ,  0.43333333, 0.46666667,
         0.5       ])
With your band selection (0.03, 0.2), that means you'll have 5 (maybe 6 
if upper bound is inclusive) for the linear regression. In theory, that 
should be enough. There might be some limitation in the statsmodels (the 
package to compute the WLS). I'll investigate that asap !

  You will
find attached the Configuration table, I hope this time is OK, fi not, can
I send it to your personal email?

Yeah, send it to me directly, it seems the attach doesn't pass the 
firewall (will check that later)

Best regards from Brussels

>   Thank you very much! Greetings from
> Mexico!
>
> Oscar
>
>
>
> On Wed, Mar 26, 2014 at 12:22 AM, Thomas Lecocq <thlecocq at gmail.com> wrote:
>
>> Hello Oscar,
>>
>>
>>
>> 2014-03-25 15:31 GMT+01:00 Oscar Alberto Castro Artola <
>> oscar.cas.art at gmail.com>:
>>
>>> Hello Thomas,
>>>
>>> first of all, thank you for your time!
>> My pleasure, I'm really happy to see people are using MSNoise.
>>
>>
>>> The filters that I am using are the
>>> following:
>>>
>>> ref low mwcs_low high mwcs_high rms_threshold mwcs_wlen mwcs_step used
>>> 1 0.03 0.03 0.2 0.2 3 10 5 1
>>>
>> That seems to be OK, there should be enough points for the linear
>> regression... Well, depending on the sampling_rate of the CCFs ? Did you
>> leave the default 20Hz value ?
>>
>>> Attached you will find a copy of my Configuration table.
>> This attachment did not make it..
>>
>> Cheers
>>
>> Thomas
>>
>>
>>
>>> I really
>>> appreciate your help! Greetings from Mexico!
>>>
>>
>>> Oscar
>>>
>>>
>>> On Tue, Mar 25, 2014 at 4:40 AM, Thomas Lecocq <thomas.lecocq at oma.be>
>>> wrote:
>>>
>>>> Hello Oscar,
>>>>
>>>> for Esteban, it seems the problems comes from a too-narrow frequency
>> band
>>>> in the MWCS, which doesn't provide enough (any) points for computing
>> the
>>>> phase shift (weighted linear regression).
>>>>
>>>> What are the filters you use ?
>>>>
>>>> Best regards,
>>>>
>>>> Thomas
>>>>
>>>> Le 21/03/2014 23:30, Thomas Lecocq a écrit :
>>>>
>>>>   Hi Oscar,
>>>>> It's the same problem as Esteban's . I'll have a look this weekend !
>>>>>
>>>>> BY then, could you tell me the parameters used (in the filter
>>>>> configuration) ?
>>>>>
>>>>> Cheers,
>>>>>
>>>>> Thom
>>>>>
>>>>> Le 21/03/2014 20:14, Oscar Alberto Castro Artola a écrit :
>>>>>
>>>>>> Hi, I was trying tu run *s05compute_mwcs.py* but I get this error
>>>>>> message:
>>>>>>
>>>>>> *Traceback (most recent call last):*
>>>>>> *  File "s05compute_mwcs.py", line 118, in <module>*
>>>>>> *    cur, ref, f.mwcs_low, f.mwcs_high, goal_sampling_rate, -maxlag,
>>>>>> f.mwcs_wlen, f.mwcs_step)*
>>>>>> *  File "/home/oscar/MSNoise-1.2.3/MWCS.py", line 145, in mwcs*
>>>>>> *    res = sm.regression.linear_model.WLS(phi, v, w**2).fit()*
>>>>>> *  File
>>>>>> "/home/oscar/anaconda/lib/python2.7/site-packages/
>>>>>> statsmodels/regression/linear_model.py",
>>>>>> line 381, in __init__*
>>>>>> *    weights=weights, hasconst=hasconst)*
>>>>>> *  File
>>>>>> "/home/oscar/anaconda/lib/python2.7/site-packages/
>>>>>> statsmodels/regression/linear_model.py",
>>>>>> line 79, in __init__*
>>>>>> *    super(RegressionModel, self).__init__(endog, exog, **kwargs)*
>>>>>> *  File
>>>>>>
>> "/home/oscar/anaconda/lib/python2.7/site-packages/statsmodels/base/model.py",
>>>>>> line 136, in __init__*
>>>>>> *    super(LikelihoodModel, self).__init__(endog, exog, **kwargs)*
>>>>>> *  File
>>>>>>
>> "/home/oscar/anaconda/lib/python2.7/site-packages/statsmodels/base/model.py",
>>>>>> line 52, in __init__*
>>>>>> *    self.data = handle_data(endog, exog, missing, hasconst,
>> **kwargs)*
>>>>>> *  File
>>>>>>
>> "/home/oscar/anaconda/lib/python2.7/site-packages/statsmodels/base/data.py",
>>>>>> line 397, in handle_data*
>>>>>> *    return klass(endog, exog=exog, missing=missing,
>> hasconst=hasconst,
>>>>>> **kwargs)*
>>>>>> *  File
>>>>>>
>> "/home/oscar/anaconda/lib/python2.7/site-packages/statsmodels/base/data.py",
>>>>>> line 78, in __init__*
>>>>>> *    self._check_integrity()*
>>>>>> *  File
>>>>>>
>> "/home/oscar/anaconda/lib/python2.7/site-packages/statsmodels/base/data.py",
>>>>>> line 246, in _check_integrity*
>>>>>> *    if len(self.exog) != len(self.endog):*
>>>>>> *TypeError: len() of unsized object*
>>>>>>
>>>>>> Thank you in advance!
>>>>>>
>>>>>> Oscar
>>>>>> _______________________________________________
>>>>>> MSNoise mailing list
>>>>>> MSNoise at mailman-as.oma.be
>>>>>> http://mailman-as.oma.be/mailman/listinfo/msnoise
>>>>>>
>>>>> _______________________________________________
>>>>> MSNoise mailing list
>>>>> MSNoise at mailman-as.oma.be
>>>>> http://mailman-as.oma.be/mailman/listinfo/msnoise
>>>>>
>>>> _______________________________________________
>>>> MSNoise mailing list
>>>> MSNoise at mailman-as.oma.be
>>>> http://mailman-as.oma.be/mailman/listinfo/msnoise
>>>>
>>>>
>>> _______________________________________________
>>> MSNoise mailing list
>>> MSNoise at mailman-as.oma.be
>>> http://mailman-as.oma.be/mailman/listinfo/msnoise
>>>
>>>
>> _______________________________________________
>> MSNoise mailing list
>> MSNoise at mailman-as.oma.be
>> http://mailman-as.oma.be/mailman/listinfo/msnoise
>>
>>
>>
>>
>> _______________________________________________
>> MSNoise mailing list
>> MSNoise at mailman-as.oma.be
>> http://mailman-as.oma.be/mailman/listinfo/msnoise