Re: [Gretl-devel] INFO: $obsmicro refers to day instead of week for weekly series
Am 18.03.2019 um 20:53 schrieb Artur Tarassow:
Just as a follow up: Python's Pandas library also keeps the
structure
<year:month:day> for weekly data. Also <pandas.DatetimeIndex.freq>
involves a series with 52 unique discrete values exactly corresponding
to what gretl's seasonals() function returns (even though you get a list
of series here). So I guess things are fine as they are in gretl ;-)
Well, seasonals() is helpful but the large collection of dummies is
certainly no week indicator. Maybe something like my workaround function
could be added to the extra.zip add-on. (If the seasonals are created in
a function they wouldn't be visible and not be flooding the active dataset.)
To be explicit, this is what I mean - Artur, perhaps you could test it
also with a dataset that doesn't start at week 1?
<hansl>
function series obsweek(void)
if $pd != 52
funcerr "Not a weekly dataset"
endif
list s = seasonals() # gives S1...S52
series obsweek = NA
loop i=1..$pd -q
obsweek = (S$i == 1) ? i : obsweek
endloop
setinfo obsweek --description="Week indicator"
setinfo obsweek --discrete # doesn't work !?
return obsweek
end function
</hansl>
cheers
sven