[MSNoise] s07

[Thomas Lecocq]

Kasper,

I think it is related to using resample('D', how='mean') while, in 
theory, the input data is already of Daily frequency. The goal of this 
call is to have a continuous record. So, adding ""fill_method='pad' " 
should prevent NaN appearing:

py1_wmean = py1_wmean.resample('D', how='mean', fill_method='pad')
py1_wstd = py1_wstd.resample('D', how='mean', fill_method='pad'
)



Thomas



Le 05/12/2014 02:56, Kasper van Wijk a écrit :
> Dear MSnoisers:
>
> In step s07, I have a problem with the resampling in pandas. The result is
> that I lose the errorbars on my dv/v's.
>
> py1_wmean, py1_wstd = get_wavgwstd(allbut)
> py1_wmean = py1_wmean.resample('D', how='mean')
> py1_wstd = py1_wstd.resample('D', how='mean').fillna(0.0)
> data = detrend(py1_wmean)
> print 'data= ', data
>
> [kasper at localhost MSNoise-master]$ python s07plot_dtt.py
> data=  [ nan  nan  nan  nan  nan  nan  nan  nan  nan  nan  nan  nan  nan
>   nan  nan
>    nan  nan  nan  nan  nan  nan  nan  nan  nan  nan  nan  nan  nan  nan  nan
>    nan  nan  nan  nan  nan  nan  nan  nan  nan  nan  nan  nan  nan  nan  nan
>    nan  nan  nan  nan  nan  nan  nan  nan  nan  nan  nan  nan  nan  nan  nan
>    nan  nan  nan  nan  nan  nan  nan  nan  nan  nan  nan  nan  nan  nan  nan
>    nan  nan  nan  nan  nan  nan  nan  nan  nan  nan  nan  nan  nan  nan  nan
>    nan  nan  nan  nan  nan  nan  nan  nan  nan  nan  nan  nan  nan  nan  nan
>    nan  nan  nan  nan  nan  nan  nan  nan  nan  nan  nan  nan  nan  nan  nan
>    nan  nan  nan  nan  nan  nan  nan  nan  nan  nan  nan  nan  nan  nan  nan
>    nan  nan  nan  nan  nan  nan  nan  nan  nan  nan  nan  nan  nan  nan  nan
>    nan  nan  nan  nan  nan  nan  nan  nan  nan  nan  nan  nan  nan  nan  nan
>    nan  nan  nan  nan  nan  nan  nan  nan  nan  nan  nan  nan  nan  nan  nan
>    nan  nan  nan  nan  nan  nan  nan  nan  nan  nan  nan  nan  nan  nan  nan
>    nan  nan  nan  nan  nan  nan  nan  nan  nan  nan  nan  nan  nan  nan  nan
>    nan  nan  nan  nan]
>
> If I don't resample, life is good:
>
>   py1_wmean, py1_wstd = get_wavgwstd(allbut)
>   #py1_wmean = py1_wmean.resample('D', how='mean')
>   #py1_wstd = py1_wstd.resample('D', how='mean').fillna(0.0)
>   data = detrend(py1_wmean)
>   print 'data= ', data
>
> [kasper at localhost MSNoise-master]$ python s07plot_dtt.py
> data=  [  1.05332051e-02   1.24279614e-02  -4.84688864e-01   1.11580621e-02
>     8.28039210e-03   1.04757327e-02   1.07993328e-02  -5.96110125e-04
>     7.23556870e-03   9.96330072e-03   9.54712307e-03   1.10273139e-02
>     1.06786539e-02   1.11155446e-02   9.51223128e-03   1.05238797e-02
>     1.00996973e-02   9.37587307e-03   7.18620810e-03   8.97052425e-03
>     9.08720458e-03   5.78305838e-03   8.33841882e-03   9.67250220e-03
>     8.61533043e-03   9.87386678e-03   8.22291954e-03   1.00831938e-02
>     9.36887461e-03   6.17043697e-03   8.22567034e-03   8.65890679e-03
>     8.22813898e-03   7.76614092e-03   7.84201987e-03   6.62494706e-03
>     8.02592632e-03   7.57269605e-03   8.68610949e-03   8.13714761e-03
>     7.87601472e-03   7.59459241e-03   7.61133449e-03   7.23192283e-03
>
> Cheers,
> kasper
> _______________________________________________
> MSNoise mailing list
> MSNoise at mailman-as.oma.be
> http://mailman-as.oma.be/mailman/listinfo/msnoise