'print' mixture with series
by Sven Schreiber
Hi, not sure if the following should work:
<hansl>
open denmark
series se
string s
print se s # error
</hansl>
The cause seems to be mixing a series with another type in the print
command. (Even though the error msg somewhat misleadingly complains
about the string type.)
thanks
sven
1 day, 12 hours
please unsubscrime me
by Dr.Leonardo- Romero
Dr. Leonardo Suárez Romero.
*Consultor en Finanzas. *
Cel: 3333687873
Este correo y cualquiera de sus archivos adjuntos pueden contener
información privilegiada y confidencial. Si Ud. no es el destinatario
previsto, por favor notifique inmediatamente al remitente; borre este
correo y destruya cualquier copia en su poder. Cualquier uso o divulgación
de esta información por una persona distinta al destinatario previsto queda
prohibida y puede considerarse ilegal. Protección de datos personales.- La
persona física con actividad empresarial, Leonardo Suárez Romero con número
de cédula profesional G2010509 expedida por la Comisión Nacional de Seguros
y Fianzas, con domicilio ubicado en Av. De Los Pinos 480, Colonia Rinconada
de Santa Rita, C.P. 45120 en Zapopan Jalisco, le informa que los datos
personales de personas físicas o morales que reciba en calidad de encargado
para el tratamiento de los mismos, serán únicamente utilizados para
promover los servicios de intermediación técnica y financiera ante las
aseguradoras del mercado asegurador Mexicano, así como la administración de
pólizas, gestión administrativa, servicios de cobranza , ayuda en la
determinación de siniestros amparados por la intermediación de nuestros
contratos, tal como lo establece la ley. Para mayor información acerca del
tratamiento y de sus derechos, usted puede solicitar en el domicilio
indicado, en el correo electrónico leonardosuarezromero(a)gmail.com y en la
página de Internet, el aviso de privacidad completo o comunicarse al
teléfono (33) 33 68 78 73.
1 day, 13 hours
Optional arguments (in user-defined hansl functions) and nested function calls
by Sven Schreiber
Hi,
looking at some new hansl code I'm reminded again that we have a
fundamental problem for (mildly) advanced programming, I'm afraid. An
example is probably the best way to show what I mean:
<hansl>
function void level2 (string s2)
eval exists(s) ? s : "nothing"
end function
function void level1 (string s1[null])
if !exists(s1)
string s1 = null # this second null-setting is necessary
endif
level2(s1)
end function
level1()
</hansl>
What's the problem here: In the level1 function the arg s1 gets a
default 'null' value, but this kind of null value cannot be passed on to
another function. (Gretl complains if you comment out the if-block.)
Note that this is true even though the level2 function properly handles
the possibility that its arg s2 might be 'null' (== might not exist) !
Instead a hansl programmer is forced to check for the null value of s1,
and then --a little absurdly-- the programmer must set it to 'null' yet
again. This feels quite stupid, to be honest. (No offense intended, of
course!) And it also wastes lines of code, making everything harder to read.
Wouldn't it be possible that for every function argument --I know,
except for the list type with its special role-- that has a null value
from the outer scope, this null value is automatically "copied" to the
function's inner scope? Basically that gretl automatically executes the
equivalent of the if-block above? I think this would make it much easier
/ more elegant to write real-world hansl code with nested function calls.
thanks
sven
1 day, 16 hours
glitch with switching decimal separator
by Sven Schreiber
Hi,
with a snapshot from March 11th, it seems that when (running in German
which natively would have a decimal comma) I switch off the comma in the
settings, it works, but when I try to switch it back _on_, then the
decimal point stays. Only a gretl restart corrects the issue. (Actually
I seem to remember that a restart was always required in the past, but
I'm not getting that alert message here.)
I'm talking about small scripts in the script editor which print out
decimal numbers.
thanks
sven
4 days, 4 hours
search in packages' descriptive labels
by Sven Schreiber
Hi,
it would be useful if - in the package listing window - using the search
text box would not only search the names of the packages (as it seems to
be doing), but also in the descriptive strings. We may have much less
packages than R, but we have enough that doing the eyeball scan is not a
good alternative.
OK, I acknowledge that one could also search on the web page
(https://gretl.sourceforge.net/cgi-bin/gretldata.cgi?opt=SHOW_FUNCS),
but by that rationale we wouldn't need any built-in search facility.
thanks
sven
4 days, 4 hours
LIML glitch
by Riccardo (Jack) Lucchetti
Hi all,
I just tried running tsls with the --liml option on a just-identified
equation. As is well-known, this should return the same results as
ordinary 2sls. For some reason, I'm getting an error: I suspect this is
linked to the fact that there's no constant term in the model. I don't
have time to debug this myself today, but maybe somebody else can.
Here's s simple script exposing the problem:
<hansl>
nulldata 100
series w = uniform()
series x = normal()
series y = normal()
tsls y const x ; const w
tsls y const x ; const w --liml # this works
tsls y x ; w
tsls y x ; w --liml # this doesn't
</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
-------------------------------------------------------
1 week, 2 days
ERROR: aggregate() w. user-defined functions
by Artur T.
Hi all,
I obtain an error with the aggregate() function when calling a
user-defined function for aggregation.
The following example works fine, if one groups the data by a single
variable but it fails with a second group-by variable:
<hansl>
set verbose off
clear
open grunfeld -q
function scalar quantile_range (series y)
/* Compute the interquartile range. */
scalar ret = quantile(y, 0.75) - quantile(y, 0.25)
print ret
return ret
end function
# Works
list groupby = firmname
list L = invest value
matrix agg = aggregate(L, groupby, quantile_range)
print agg
# Fails
list groupby = firmname ticker # we add "ticker"
list L = invest value
matrix agg = aggregate(L, groupby, quantile_range)
print agg
</hansl>
The error I get with current git-version using Ubuntu, is:
<error>
Data error
*** error in function quantile_range, line 3
> scalar ret = quantile(y, 0.75) - quantile(y, 0.25)
</error
Best,
Artur
1 week, 4 days
Win7 64bit, gretl 32bit and MS-MPI 10
by Marcin Błażejowski
Hi,
I need to run some test on Windows, so I have configuration like in
subject: Win7 Pro 64bit + gretl 2023b-git 32bit + MS-MPI 10. When
calling simple script:
<hansl>
mpi
eval $mpirank
end mpi
</hansl>
I'm getting error message that "there is no exit from MPI_Abort8@
procedure". In gretl output I have: "C:\Program Files\Microsoft
MPI\Bin\mpiexec.exe" /np 4 "C:\Program Files (x86)\gretl\gretlmpi"
"C:\Users\Admin2\AppData\Roaming\gretl\gretltmp-mpi.inp": exit code
3221225785
Is this related with the fact that gretl is 32 bit while Windows 64 bit?
Marcin
--
Marcin Błażejowski
2 weeks
alias for "scatters" in time series?
by Sven Schreiber
Hi,
the "scatters" command has the nice feature that in a time-series
context it makes it possible to have multiple plots of series against
time. But I find the word "scatters" completely unintuitive, arguably
misleading, and I'm surprised every time I read it in somebody's code.
Can we have a time-series-ish alias for that use case?
thanks
sven
2 weeks
lrvar bug (?) in panels
by Sven Schreiber
Hi,
I suspect that the lrvar function (long-run variance estimation) is
slightly wrong in the panel case. Consider the following example:
<hansl>
function matrix lrvarFE(series x, int b "bandwidth")
N = max($unit)
out = 0
loop i = 1..N
smpl i i --unit
errorif(abs(mean(x)) > 1e-10, "whoops, found non-zero within mean")
out += lrvar(x, b)
endloop
return out / N
end function
# -- test case
open grunfeld
# do an explicit FE demeaning first, for easier comparison
panel kstock const
series kres = $uhat
eval lrvarFE(kres, 5) # 2.0608e+005
# comparison
eval lrvar(kres, 5) # 177215.62
# the following wrong on purpose:
eval lrvar({kres}, 5) # same result as the previous line
</hansl>
Explanation: The last calculation (which is wrong) converts the kres
series directly into a matrix. Given gretl's panel storage format of
stacked time-series, this means that the values for all units are simply
stacked in the resulting single column. Then the Bartlett kernel of the
lrvar() function partly connects the different panel units, which is
spurious and distorts the result.
Given that gretl's result in the pen-ultimate line is identical to that,
I guess that's what happens.
Instead, the output of my lrvarFE function is what I believe to be the
correct one.
The similar lrcovar() function is not directly affected, because it
takes a matrix instead of (a list of) series as input, and so it's the
caller's responsibility to handle the data correctly. (Although it would
be nice if it accepted a list, but that's a different issue.)
But I could imagine that the filter functions could suffer from a
similar problem (bkfilt, bwfilt, hpfilt, perhaps pergm?, kdensity?); but
haven't checked that.
Does all that sound right?
Sven
2 weeks, 4 days
