[Gretl-devel] Re: x-range for scatters using daily data
On Mon, 2 Mar 2020, Allin Cottrell wrote:
On Mon, 2 Mar 2020, Artur Tarassow wrote:
> Am 02.03.20 um 10:06 schrieb Riccardo (Jack) Lucchetti:
>> On Sun, 1 Mar 2020, Artur Tarassow wrote:
>>
>>>> Maybe we could find a way to show dates, but there's not much space
>>>> available in "scatters". At least now we don't over-extend
the x-axis.
>>>
>>> Maybe its beneficial to rotate the x-labels to save some space. But it
>>> has no priority I guess.
>>
>> Something like this, perhaps?
>>
>> <hansl>
[nice example script]
>> </hansl>
>
> I don't get any xtics here :-/ But apart from this: This is a very nice
> example of exploiting "scatters" for enhancing multiplots.
Hmm, why aren't you getting any xtics? They're shown here OK. Anyway, Jack's
script might be a good approach to getting dates in place for "scatters",
maybe better than trying to use gnuplot's "timefmt".
I think I found it. You should get tyhe xtics only using current git. With
2019d, the x scale is in "years", with current git is just an integer
sequence. This version of the script should give you the xtics with 2019d
as well:
<hansl>
open IT_companies.gdt --frompkg=MGARCH
scatters 2 3 4 5 --output="(a)dotdir/x.gp"
s = readfile("(a)dotdir/x.gp")
# generate daily xtics
ntics = 6
v = round(seq(1, $nobs, $nobs/ntics))
compat = $version < 20200
xtics = "set xtics ("
loop i = 1 .. ntics --quiet
vi = v[i]
if compat
# rough!!!
x = $obsmajor[vi] + ($obsminor[vi] - 1)/12 + ($obsmicro[vi] - 1)/365
else
x = vi
endif
xtics ~= sprintf("\"%s\" %8.4f", obslabel(vi), x)
if i<ntics
xtics ~= ", "
endif
endloop
xtics ~= ") rotate by -45"
# modify gnuplot script and plot
s = regsub(s, "set xtics (.*)", xtics)
outfile "(a)dotdir/y.gp" --quiet
print s
end outfile
gnuplot --input="(a)dotdir/y.gp" --output=display
</hansl>
-------------------------------------------------------
Riccardo (Jack) Lucchetti
Dipartimento di Scienze Economiche e Sociali (DiSES)
Università Politecnica delle Marche
(formerly known as Università di Ancona)
r.lucchetti(a)univpm.it
http://www2.econ.univpm.it/servizi/hpp/lucchetti
-------------------------------------------------------