Gretl-Hansl "IDE" for Sublime editor
by Artur Tarassow
Dear all,
some weeks ago I've started to switch to the sublime editor
(https://www.sublimetext.com/) for writing gretl code.
Even though I like the simplicity of the gretl editor and its features
such as syntax highlighting and auto-itendation, it lacks some features
of modern IDEs such as "goto-anything", "goto-definition", custom
keybindings, fancy themes, git-implementation, snippets etc. which make
life much easier when working on larger projects. Don't get me wrong,
the gretl editor is great but was _never_ supposed to become a proper
software-development IDE but rather has another focus which is totally fine.
So I started to write the "Hansl-Gretl-Language" package for sublime
which includes the following features:
- 3 gretl build-systems (client mode, batch mode, and REPL mode) for
executing hansl code by means of sublime (plots are also working!)
- syntax-highlighting
- completion of gretl commands, accessors and keywords
- some snippet examples for speeding up coding
The project still has the following (known) issues:
- no auto-itendation (still have to figure out how this works)
- issues with some corner-cases which are not syntax-highlighted (regex
can become so hard!)
The package can be downloaded through sublimes package control system,
and can be found here:
https://packagecontrol.io/packages/Hansl-Gretl-Language
If somebody wants to participate on this project, check out the code on
my github repository:
https://github.com/atecon/Hansl-Gretl-Language
For those interesting in the sublime editor, check out "OdatNurd"'s
brilliant tutorials on youtube:
https://www.youtube.com/user/nurdz
Enjoy the package,
Artur
4 months
Gretlwiki announcement
by DeFeroci
Dear all,
I would like to announce that the gretl wiki is back up and running at:
https://gretlwiki.econ.univpm.it/index.php/Main_Page
For security reasons (and previous experience) only registered users can
modify the wiki and self-registration is no longer possible. So, if you
want to be added, please write to me privately or reply to this email.
Information on how to be registered and using the wiki can be found here:
https://gretlwiki.econ.univpm.it/index.php/Wiki
At the moment the wiki administrators are me and Cristián Arturo Ducoing
Ruiz.
If you have any suggestions or requests, you are welcome to contact us.
Happy holidays.
Federico
5 months, 3 weeks
Simulated maximum likelihood using Halton numbers
by Alecos Papadopoulos
I have a sample of size n=1789. I need to run simulated maximum
likelihood on it, generating series of pairs of independent (w,u) per
observation to compute the value of the likelihood as an simuluated average.
Instead of using the RNG of gretl, I was thinking of using Halton
numbers, in sequences of length, say, 100. It would appear that this
requires to construct (once), 1789 (x) 2 = 3578 sequences of Halton
numbers of length 100.
Given my current understanding of the "halton" function in gretl, It
appears that I can only construct 40 distinct sequences maximum.
Naively, one could think "then, construct 40 sequences each of length
9000" and then chop them in pieces of length 100 to get the number of
required distinct sequences (or generate 1 sequence with length 357800
for that matter)... but how would one preserve the "well-balanced
spacing" in (0,1) that is the advantage of Halton numbers? (again this
"idea" is just to make small talk -it may be totally silly).
Is there any way to generate what I need inside gretl?
--
Alecos Papadopoulos PhD
Affiliate Researcher
Dpt of Economics, Athens University of Economics and Business
Foundation for Economic and Industrial Research (IOBE)
web: alecospapadopoulos.wordpress.com/
ORCID:0000-0003-2441-4550
5 months, 4 weeks
Fwd: Norm of gradient in mle command
by Riccardo (Jack) Lucchetti
On 27/10/2023 01:16, Alecos Papadopoulos wrote:
> In an mle command with numerical derivatives I get in the script output
> (actual example)
>
> <<
>
> Gradients: 2.4007 6.7772 -5.3557 -0.89559 (norm
> 1.96e+000)
>
> >>
>
> Can somebody please advise which norm is this, preferably giving also
> its computation formula?
Call b is the parameter vector ad g the gradient vector, both with n
elements. Then the norm is computed as
gradnorm = sqrt(abs(b'g)/n)
I don't remember where this formula came from, Allin and I worked on
this piece of code long ago.
If you're curious to see the rest of the code, see
https://sourceforge.net/p/gretl/git/ci/master/tree/lib/src/gretl_bfgs.c
from line 1147 onwards.
-------------------------------------------------------
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
-------------------------------------------------------
--
-------------------------------------------------------
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
-------------------------------------------------------
6 months
Fwd: Inconsistency in info related to Gauss-Legendre
by Riccardo (Jack) Lucchetti
On 27/10/2023 01:36, Alecos Papadopoulos wrote:
> In the help text for function "quadtable" we read
>
> <<When the Gauss–Legendre type is selected, the optional arguments a and
> b can be used to control the lower and upper limits of integration, the
> default values being –1 and 1. >>
>
> In Gretl guide ch. 37.4 page 384, in a table summarizing numerical
> integration we read
>
> <<Gauss–Legendre by default, a = 0 and b = 1>>
>
> Which one is the actual default value, a=-1, or a=0?
The guide is wrong, I'll correct it in a few minutes: the default is -1.
The pertinent code is at
https://sourceforge.net/p/gretl/git/ci/master/tree/lib/src/matrix_extra.c
line 2866.
Thanks!!!
-------------------------------------------------------
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
-------------------------------------------------------
--
-------------------------------------------------------
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
-------------------------------------------------------
6 months
Inconsistency in info related to Gauss-Legendre
by Alecos Papadopoulos
In the help text for function "quadtable" we read
<<When the Gauss–Legendre type is selected, the optional arguments a and
b can be used to control the lower and upper limits of integration, the
default values being –1 and 1. >>
In Gretl guide ch. 37.4 page 384, in a table summarizing numerical
integration we read
<<Gauss–Legendre by default, a = 0 and b = 1>>
Which one is the actual default value, a=-1, or a=0?
Thanks.
--
Alecos Papadopoulos PhD
Affiliate Researcher
Dpt of Economics, Athens University of Economics and Business
Foundation for Economic and Industrial Research (IOBE)
web: alecospapadopoulos.wordpress.com/
ORCID:0000-0003-2441-4550
6 months
Norm of gradient in mle command
by Alecos Papadopoulos
In an mle command with numerical derivatives I get in the script output
(actual example)
<<
Gradients: 2.4007 6.7772 -5.3557 -0.89559 (norm
1.96e+000)
>>
Can somebody please advise which norm is this, preferably giving also
its computation formula?
Thanks.
--
Alecos Papadopoulos PhD
Affiliate Researcher
Dpt of Economics, Athens University of Economics and Business
Foundation for Economic and Industrial Research (IOBE)
web: alecospapadopoulos.wordpress.com/
ORCID:0000-0003-2441-4550
6 months
Gretl Crash on Windows
by Federico Fiorani
Hi everybody,
the following script crashes on Windows (gretl 2023c-git build date
2023-10-07) when I try to convert to a bundle.
<code>
bundle bun = null
bun.header = "text/html,application/xhtml+xml,application/xml"
bun.URL = "http://dataservices.imf.org/REST/SDMX_JSON.svc/Dataflow"
curl(&bun)
a = bread(bun.output)
</code>
Does anyone have the same problem?
Best
Federico
6 months, 1 week
