9:49 a.m.
Dear Gretl Masters (Allin and Jack ;)),
I did a review of hansl-primer-a4.pdf (learning hansl), and found very few
defects (one or two):
Chapter 7. Control flow, page 27
scalar x = 1
loop foreach f sqrt exp ln
scalar y = $f(1)
print y
endloop
(Unused "scalar x = 1" )
Chapter 10. Series and Lists, page 41
Another option for appending terms to, or dropping terms from, an existing
list is to use += or +=, respectively, as in
xlist += cpi
zlist -= cpi
(you meant -=)
----
That's it.
Thank you for the excellent book.
Best Regards,
Hélio
10:35 a.m.
On Tue, 29 Apr 2014, Hélio Guilherme wrote:
Dear Gretl Masters (Allin and Jack ;)),
I did a review of hansl-primer-a4.pdf (learning hansl), and found very few
defects (one or two):
Chapter 7. Control flow, page 27
scalar x = 1
loop foreach f sqrt exp ln
scalar y = $f(1)
print y
endloop
(Unused "scalar x = 1" )
Chapter 10. Series and Lists, page 41
Another option for appending terms to, or dropping terms from, an existing
list is to use += or +=, respectively, as in
xlist += cpi
zlist -= cpi
(you meant -=)
Thanks, Hélio, those points are now fixed in CVS.
Allin
10:03 a.m.
New subject: Gretl crash
the last line in this script makes gretl to crash: (gretl from today CVS
on Ubuntu Linux 12.04)
<hansl>
nulldata 8
series vari= { 1; 2; 3; 4; 5; 6; 7; 8}
series w = { 1; 1; 1; 1; 2; 2; 2; 2 }
summary vari --weight=w
</hansl>
--
Firma Arista
*Ignacio Díaz-Emparanza*
Zuzendaria/Director
ignacio.diaz-emparanza(a)ehu.es <mailto:ignacio.diaz.emparanza@ehu.es>
94 6013732
*EKONOMIA APLIKATUA III SAILA (EKONOMETRIA ETA ESTATISTIKA)/
DEPARTAMENTO DE ECONOMÍA APLICADA III (ECONOMETRÍA Y ESTADÍSTICA)
UPV/EHU*
Avda. Lehendakari Aguirre, 83 | 48015 BILBAO
*T.: +34 946013740* | *F.: +34 946013754*
*www.ea3.ehu.es* <http://www.ea3.ehu.es>
ERNE! Baliteke mezu honen zatiren bat edo mezu osoa legez babestuta
egotea. Mezuak badu bere hartzailea. Okerreko helbidera heldu bada
(helbidea gaizki idatzi, transmisioak huts egin) eman abisu igorleari,
korreo honi erantzunda. Kontuz! Mezua ez bada zuretzat, ez erabili, ez
zabaldu beste inori, ez kopiatu eta ez baliatu.
¡ATENCIÓN! Este mensaje contiene información privilegiada o confidencial
a la que sólo tiene derecho a acceder el destinatario. Si usted lo
recibe por error le agradeceríamos que no hiciera uso de la información
y que se pusiese en contacto con el remitente.
E-mail hau inprimatu baino lehen egiaztatu inprimatzeko beharra.
Antes de imprimir este e-mail piense bien si es necesario hacerlo.
11:32 a.m.
New subject: Gretl crash
On Wed, 30 Apr 2014, Ignacio Diaz-Emparanza wrote:
the last line in this script makes gretl to crash: (gretl from today
CVS on
Ubuntu Linux 12.04)
<hansl>
nulldata 8
series vari= { 1; 2; 3; 4; 5; 6; 7; 8}
series w = { 1; 1; 1; 1; 2; 2; 2; 2 }
summary vari --weight=w
</hansl>
Thanks, Ignacio: this is now fixed in CVS.
The crash was related to the fact that, with so few datapoints, the 5% and
95% percentiles are effectively undefined, and this was unhandled (my
fault). Fun fact: this is one of the few fixes that was physically coded
jointly by Allin and myself, cause we happen to be in the same room at
the moment! :)
-------------------------------------------------------
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
-------------------------------------------------------
11:45 a.m.
New subject: Gretl crash
We can say that the next release of gretl is so important that Jack and
Allin get in the same room to work on it! ;)
I just wanted to say that this small script is a good candidate for a
Regression Test (software speaking). I have confirmed the crash and the fix.
It would be very useful if we all collect a set of scripts to do Regression
Tests whenever a new build is done.
Hélio
On Wed, Apr 30, 2014 at 12:32 PM, Riccardo (Jack) Lucchetti <
r.lucchetti(a)univpm.it> wrote:
On Wed, 30 Apr 2014, Ignacio Diaz-Emparanza wrote:
the last line in this script makes gretl to crash: (gretl from today CVS
> on Ubuntu Linux 12.04)
>
> <hansl>
> nulldata 8
> series vari= { 1; 2; 3; 4; 5; 6; 7; 8}
> series w = { 1; 1; 1; 1; 2; 2; 2; 2 }
> summary vari --weight=w
> </hansl>
>
Thanks, Ignacio: this is now fixed in CVS.
The crash was related to the fact that, with so few datapoints, the 5% and
95% percentiles are effectively undefined, and this was unhandled (my
fault). Fun fact: this is one of the few fixes that was physically coded
jointly by Allin and myself, cause we happen to be in the same room at the
moment! :)
-------------------------------------------------------
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
-------------------------------------------------------
_______________________________________________
Gretl-devel mailing list
Gretl-devel(a)lists.wfu.edu
http://lists.wfu.edu/mailman/listinfo/gretl-devel
1:05 p.m.
New subject: Gretl crash
Am 30.04.2014 13:45, schrieb Hélio Guilherme:
We can say that the next release of gretl is so important that Jack
and
Allin get in the same room to work on it! ;)
I just wanted to say that this small script is a good candidate for a
Regression Test (software speaking). I have confirmed the crash and the fix.
It would be very useful if we all collect a set of scripts to do
Regression Tests whenever a new build is done.
Yes building up some unit tests (and / or regression tests) over time
would be very good. Allin has a bunch of stuff he runs for himself I
think, but over the years I think there were various occasions where
this did not catch the problem IIRC.
-sven
1:05 p.m.
New subject: Gretl crash
On 30/04/14 13:32, Riccardo (Jack) Lucchetti wrote:
On Wed, 30 Apr 2014, Ignacio Diaz-Emparanza wrote:
> the last line in this script makes gretl to crash: (gretl from today
> CVS on Ubuntu Linux 12.04)
>
> <hansl>
> nulldata 8
> series vari= { 1; 2; 3; 4; 5; 6; 7; 8}
> series w = { 1; 1; 1; 1; 2; 2; 2; 2 }
> summary vari --weight=w
> </hansl>
Thanks, Ignacio: this is now fixed in CVS.
The crash was related to the fact that, with so few datapoints, the 5%
and 95% percentiles are effectively undefined, and this was unhandled
(my fault). Fun fact: this is one of the few fixes that was physically
coded jointly by Allin and myself, cause we happen to be in the same
room at the moment! :)
Thanks. By the way I prepared a function for doing the same thing by
means of using the --matrix option of 'sumary' and I see there are
several results different: median, 5% and 95% percentiles and
Interquartile range ¿is this correct?
(and note: the option --weight is not explained in the reference for the
suumary command)
<hansl>
function void summaryw(series y, series w)
string cname=argname(y)
matrix mvari = { y }
matrix mw = { w }
matrix A = {}
loop i=1..rows(mvari)
matrix A = A | mshape(mvari[i,1], mw[i,1],1)
endloop
colnames(A, "@cname")
A
summary --matrix=A
end function
/*The table below shows the times taken by a group of walkers to
complete a 15-mile walk. Their times have been recorded to the nearest hour.
Time (hours) 3 4 5 6 7 8
Frequency 2 5 12 11 4 3
*/
nulldata 6
series hours= { 3; 4; 5; 6; 7; 8}
series Freq = { 2; 5; 12; 11; 4; 3 }
summaryw(hours, Freq)
summary hours --weight=Freq
</hansl>
--
Firma Arista
*Ignacio Díaz-Emparanza*
Zuzendaria/Director
ignacio.diaz-emparanza(a)ehu.es <mailto:ignacio.diaz.emparanza@ehu.es>
94 6013732
*EKONOMIA APLIKATUA III SAILA (EKONOMETRIA ETA ESTATISTIKA)/
DEPARTAMENTO DE ECONOMÍA APLICADA III (ECONOMETRÍA Y ESTADÍSTICA)
UPV/EHU*
Avda. Lehendakari Aguirre, 83 | 48015 BILBAO
*T.: +34 946013740* | *F.: +34 946013754*
*www.ea3.ehu.es* <http://www.ea3.ehu.es>
ERNE! Baliteke mezu honen zatiren bat edo mezu osoa legez babestuta
egotea. Mezuak badu bere hartzailea. Okerreko helbidera heldu bada
(helbidea gaizki idatzi, transmisioak huts egin) eman abisu igorleari,
korreo honi erantzunda. Kontuz! Mezua ez bada zuretzat, ez erabili, ez
zabaldu beste inori, ez kopiatu eta ez baliatu.
¡ATENCIÓN! Este mensaje contiene información privilegiada o confidencial
a la que sólo tiene derecho a acceder el destinatario. Si usted lo
recibe por error le agradeceríamos que no hiciera uso de la información
y que se pusiese en contacto con el remitente.
E-mail hau inprimatu baino lehen egiaztatu inprimatzeko beharra.
Antes de imprimir este e-mail piense bien si es necesario hacerlo.
1:49 p.m.
New subject: Gretl crash
On Wed, 30 Apr 2014, Ignacio Diaz-Emparanza wrote:
On 30/04/14 13:32, Riccardo (Jack) Lucchetti wrote:
> On Wed, 30 Apr 2014, Ignacio Diaz-Emparanza wrote:
>
>> the last line in this script makes gretl to crash: (gretl from today CVS
>> on Ubuntu Linux 12.04)
>>
>> <hansl>
>> nulldata 8
>> series vari= { 1; 2; 3; 4; 5; 6; 7; 8}
>> series w = { 1; 1; 1; 1; 2; 2; 2; 2 }
>> summary vari --weight=w
>> </hansl>
>
> Thanks, Ignacio: this is now fixed in CVS.
>
> The crash was related to the fact that, with so few datapoints, the 5% and
> 95% percentiles are effectively undefined, and this was unhandled (my
> fault). Fun fact: this is one of the few fixes that was physically coded
> jointly by Allin and myself, cause we happen to be in the same room at the
> moment! :)
Thanks. By the way I prepared a function for doing the same thing by means of
using the --matrix option of 'sumary' and I see there are several results
different: median, 5% and 95% percentiles and Interquartile range ¿is this
correct?
I see your point, but the problem has no easy solution, because to "solve"
the problem (whatever that may mean), we need an acceptable definiton of
an order statistc (especially so, in the case of weighted data).
In gretl, we've always used an "interpolation" approach to the calculation
of quantiles, which is NOT the traditional definition, but I believe there
are very good reasons to stand by it. For example, according to the
traditional definition of a quantile, the median of (-x) needs not be the
same as -median(x). This is something I personally consider utterly
ridiculous. From this a number of other choices follow.
For example, what is the median if x for a population in which x_i =
min(x) for 60% (say) of observations? Well, in my eyes that's open to
interpretation. The current approach of the summary command with the
--weights option is to return NA, but I'm open to be convinced that other
conventions are preferable.
-------------------------------------------------------
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
-------------------------------------------------------
3889
days inactive
3890
days old
7 comments
5 participants
participants (5)
-
Allin Cottrell -
Hélio Guilherme -
Ignacio Diaz-Emparanza -
Riccardo (Jack) Lucchetti -
Sven Schreiber