minor path printing bug
by Sven Schreiber
maybe the following is related to the recent workdir changes;
each time I re-run a script with the line
open de_6dim77q1.gdt
at the top, I get one more slash in the echoed path, like so:
Lese Datendatei
<...>/permaconf///////////////////de_6dim77q1.gdt
However, so far this does not cause any real problems (that's why no bug
tracker item)
cheers,
sven
16 years, 2 months
for loops too restrictive?
by Sven Schreiber
Hi,
why is the following script invalid ("end index must be greater than
start index")
<script1>
loop for i=1..1
print "check"
end loop
</script1>
when the next equivalent one works fine?
<script2>
lagorder = 1
loop for i=1..lagorder
print "check"
end loop
</script2>
thanks,
sven
16 years, 2 months
Improvement of logit regression modelling in gretl
by andreas.rosenblad＠ltv.se
A few suggestions for improving the handling of logit regression models in
gretl:
1. Multinomial models, when the dependent variable is discrete but not
binary and does not have a natural ordering does not have a native
implementation in gretl yet. Please implement it. For example, include a
radio button in the GUI where the user can choose between using a
multinomial logit or an ordered logit (the current implementation) when the
dependent variable is discrete but not binary.
2. Currently, it is necessary that the lowest value for the dependent
variable is 0 and this is used as the reference category. Change this so
that any integer values can be used for the dependent variable and any of
these can be used as reference category.
3. Make it easy for the user to choose reference category by including
radio buttons for First category, Last category or Custom value, as in
e.g. SPSS. See the two attached screenshots.
Best regards
Andreas Rosenblad
(See attached file: logit suggestion.JPG)
(See attached file: SPSS logit.jpg)
16 years, 2 months
change to quantile()
by Riccardo (Jack) Lucchetti
Folks,
as a preliminary step towards fixing #1855731, I re-wrote the quantile()
function in a much more efficient way. Example script:
<script>
nulldata 1000000
set seed 34789023
x = round(normal()*100)/100
# insert a few missing values here and there
# as a robustness check
m = uniform()<0.05
x = m ? NA : x
set stopwatch
loop for (p=0.05; p<0.99; p+=0.05)
q = quantile(x, p)
printf "q = %5.2f\n", q
end loop
et = $stopwatch
</script>
results are the same as with the old function, but speed is about 10 times
higher. I ran quite a few tests, and it seems that results are correct;
however, you can never be too cautious, so if you could test the
quantile() implementation you find in CVS and report any bugs you find,
that would be very helpful.
Thanks.
Riccardo (Jack) Lucchetti
Dipartimento di Economia
Università Politecnica delle Marche
r.lucchetti(a)univpm.it
http://www.econ.univpm.it/lucchetti
16 years, 2 months
