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
7 months
BACE Estimation difficulties
by Brian Revell
I thought i would try this option estimating a simple model (always the
best way to gain experience and confidence with the routine).
Independent variable, 1 indep variable . No lag on Y. . Zero selected for
out of sample forecast Constant -can be omitted. Best 4 models selected
-again as the default.
This was the output:-
Dependent viariable was removed from the X list.* It is not there -there is
only a single indep variable. *
Error message from BACE_GUI():
Average model size should be greater then 0 and lower than K !!! *Average
model size specified is the default 1*
*Some guidance wlecomed please*
*Brian*
8 months
Re: How to apply the HAC option in Gretl?
by Cottrell, Allin
On Mon, Nov 7, 2022 at 12:29 PM Torbjørn Lorentzen
<torbjorn.lorentzen(a)outlook.com> wrote:
>
> I'm a new user of Gretl. I'm using ols in the estimation of a time series model. T want to us ethe HAC-option to correct for autocorrelation and heteroscedasticity. I tried the robust option, but to what I can see it only corrects for heteroscedasticity? Great if anyone could help me with the Gretl script for HAC and how to put it together with the ols command. All suggestions are appreciated 😊
Hi Tobben. If you go to the menu item /Tools/Preferences/General
you'll see an "HCCME" tab. That lets you set your preferences for
robust covariance matrix estimation. Note that the HAC option is only
available when your dataset is recognized as time-series.
In scripting usage the relevant option is --robust. What exactly that
does depends on the aforementioned setting. But if the data are time
series gretl will use the HAC estimator by default. In that case you
should see a line like the following in the estimation output:
HAC standard errors, bandwidth 4 (Bartlett kernel)
Allin Cottrell
8 months, 2 weeks
Some suggestions on the "Number of cases ’correctly predicted'" from a logit/probit outputs
by Fred Engst
Hi Jack, Allin, Sven and all others on the gretl team,
As you know, "Number of cases ’correctly predicted’" in a logit/probit model can be miss-leading even in a 50/50 split case.
What we should be comparing is not zero ‘correctly predicted’, but rather random assignments based on sample mean.
If a sample is 50/50 split, a random assignment would get 50% "correctly predicted", in theory. If our model's 'correctly predicted' is 70%, we are only 20 percentage points higher than a model based on random assignment, representing an improvement over the random assignment model by only 40%.
Thus, I would like to propose an alternative output from gretl, i.e. the “Extra number of cases 'correctly predicted' over random assignment” (or something like that), call this dot_R-square perhaps.
Dot_R-saure = (Y_hat_model - Y_hat_random) / (1-Y_hat_random)
where Y_hat_model = sum(Y_hat_model_i=Y_i)/N
Y_hat_random = Y_hat^2 + (1-Y_hat)^2
Y_hat is the sample mean
Y_hat_model_i = Pro(Y_i = 1) >Y_hat
Pro(Y_i = 1) >Y_hat = 1, if Pro(Y_i = 1) >Y_hat is true
Unlike the McFadden R-squared, the interpretation of this is fairly straight forward, i.e. the percent that our model is better off than a model based on random assignment.
Best,
Fred
8 months, 2 weeks
deseas function arima option
by Javier Sansa
I have been trying to use the different options of the deseas function in
gretl. All the options I used seem to work as intended except the arima
option. To be more specific I have tried to build a deseasonalizing
procedure in which a certain arima model is imposed. According to the gretl
reference manual this should be done by building a bundle and using "arima"
and a 6x1 matrix which defines the model, like in the following example:
matrix start ={0,2,1,0,1,1}
bundle b=defbundle(verbose=2,seats=1,outliers=7, arima=start)
The bundle definition works, i.e a bundle b is built, but function deseas
fails because it does not recognize the arima key
I attach a file with data and a short script should something be unclear.
Please, let me know if I failed to follow the manual or there is something
strange going on with the program.
And since we are at the end of 2022, a happy new year!
Javier Sansa
8 months, 2 weeks
count panel data - poisson model with dummies?
by lars.ahnland＠outlook.com
Hi,
I have count data with what appears to be poisson distribution of the dependent variable in panel data. There is no automatic method for doing this in Gretl as I have understood it, so I am using a count, poisson, model with cross-sectional dummies. Gretl deletes the last dummy so it becomes a fixed effects model, so intuitively, it looks as if the model should be ok. Do you agree? One problem is that Gretl fails to include an overdispersion (Chi-square) test, so I am not sure the poisson model is valid. Can you help me with this too?
Best regards,
Lars Ahnöand
8 months, 2 weeks
Function package DHF_test is going to be retired
by Sven Schreiber
Hello everybody,
this is to let you know that the contributed function package "DHF_test"
is going to be removed from the public package server soon. If you rely
on that package for whatever reason, you can keep a local copy on your
system(s). (And also perhaps write us about the fact that you're using
it, which would be interesting as well.)
The package DHF_test implements an early test for seasonal unit roots,
its name stands for Dickey-Hasza-Fuller. But for seasonal unit root
testing we currently recommend to use the "GHegy" package instead, which
also has been available for many years on the gretl package server, and
which refers to the HEGY publication (Hylleberg/Engle/Granger/Yoo, 1990).
DHF_test package author Marcin obviously knows about this
reorganization, as he is also a member of the gretl team.
If you have any remaining questions about this (and perhaps also about
seasonal unit root testing in gretl), please go ahead and ask them.
cheers
sven
8 months, 2 weeks
Gretl no longer reads panel data
by lars.ahnland＠outlook.com
Hi,
Suddenly, Gretl no longer reads panel data. It stopped working today, with a file that always formed before. The error started appearing when I added a variable, but it was still there when I deleted the variable. Gretl says it discovered the first column to contain a string variable and suggests to remove it. But if I do it will no longer be panel data. Anyone else experiencing this? How do I fix it?
Best,
Lars
9 months
Gretl no longer reads panel data
by lars.ahnland＠outlook.com
Hi,
Suddenly, Gretl no longer reads panel data. It stopped working today, with a file that always formed before. The error started appearing when I added a variable, but it was still there when I deleted the variable. Gretl says it discovered the first column to contain a string variable and suggests to remove it. But if I do it will no longer be panel data. Anyone else experiencing this? How do I fix it?
Best,
Lars
9 months
chatGPT on gretl
by Artur T.
Dear all,
chatGPT has written a poem for our favorite econometrics software which
I would like to share with you ;-)
Gretl, the econometrics software,
Makes it easy to analyze data galore.
Time series or panel, it's all the same,
With Gretl, the possibilities are tame.
It offers many commands and functions,
To estimate and test models with precision.
With Gretl, we become econometrics kings,
Always ready to learn new things.
So we install Gretl today,
And soon we'll be econometrics at play.
Gretl, the software that takes us there,
In the world of data and econometrics, without a care.
Artur
9 months, 1 week
