how to get the time periodicity of panel data (in a script)
by Sven Schreiber
Hi,
for a given panel dataset with embedded information about the time
dimension structure I want to obtain the periodicity in the time
direction; for example the value 4 if the data are quarterly (for
several cross units, of course, because we're in a panel). I know that
the $pd accessor doesn't work here because in a panel context it gives
me the total length of the time dimension (as documented).
I'm wondering whether there is a better (shorter) way than the following
hansl function and test cases:
<hansl>
function scalar pd_paneltime(void)
catch pd = nelem(seasonals())
if $error
return 1 # assume the case of non-seasonal data
# (but could also be simply no panel-time info)
else
return pd
endif
end function
# artificial example
nulldata 20
setobs 5 1:1 --stacked-time-series
setobs 4 1900:1 --panel-time # quarterly
eval pd_paneltime() # gives 4
# real example
open grunfeld # annual
eval pd_paneltime() # gives 1
</hansl>
thanks for any hints
Sven
1 year, 7 months
ketvals 0.9 released
by Riccardo (Jack) Lucchetti
OK, time for some plugging :)
Francesco Valentini and I have just released version 0.9 of the "ketvals"
package, that provides OLS and IV estimation of linear model with
time-varying parameters using a non-parametric technique developed by
Giraitis, Kapetnaios and others.
Download, enjoy, and send us bug reports!
-------------------------------------------------------
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 year, 8 months
Panel Regression Model: Fixed Effect, Random Effects, or OLS Regression?
by 3J LEMA
Dear Experts,
1. If there is heterogeneity in the panel model/data, is it true that the
appropriate model is either the Fixed Effect or Random Effects model?
2. I want to choose the appropriate panel regression model when there are
autocorrelations, heterogeneity, multicollinearity, and non-normality. Any
guiding non-highly mathematical materials that you can suggest?
Thank you.
1 year, 8 months
Strange behavior of the ksmmoth function
by Paolo Chirico
Dears,
for a simulation study I want to repeat the estimate of a state space
model /N/ times. Repeating the estimate is not a problem, but if I call
the ksmooth function after each estimate (MLE block), gretl allows me it
on the first iteration, but not the following ones.
It appears that the function cannot find the bundle in subsequent iterations
This happens if the repetitions occur in a single script run.
Obviously it doesn't happen if every rep has its own run.
It's really strange!
Is it a bag or am I missing something?
Thanks,
Paolo
--
Paolo Chirico
Università del Piemonte Orientale
Dip. di Giurisprudenza e Scienze Politiche,
Economiche e Sociali (DIGSPES)
Alessandria, Italia
--
Questa email è stata esaminata alla ricerca di virus da AVG.
http://www.avg.com
1 year, 8 months
Is calling R script thread-safe in MPI context?
by Marcin Błażejowski
Hi,
suppose I have a R code which load gretldata via '--send-data' flag and
stores the results RES via 'gretl.export(RES, base)' and 'base' is
$mpirank based string. Is this thread-safe?
Marcin
--
Marcin Błażejowski
1 year, 8 months
On the use of Fixed Effect and Random Effects models for Panel Data
by 3J LEMA
I have read a blog stating that if there is heterogeneity in the panel
data the appropriate model is either the Fixed Effect or Random Effects
model, not the OLS Regression. Do you agree with the blogger?
I'm interested to know to determine the appropriate panel regression model
when there are autocorrelations, heterogeneity, multicollinearity, and
non-normality. Appreciate if you suggest a textbook that deals with the
OLS, FE, and RE models for panel data.
Thank you.
Lema
1 year, 8 months
PairPlot package v0.8
by Artur T.
Hi gretl users,
there has been an update of the PairPlot package which may be of
interest for some of you. The recent version is 0.8, and it includes a
couple of bug fixes. The changelog is:
Version 0.8 (January 2022):
- Handle eventual missing values in passed series.
- Bugfix of gnuplot's "offset" option: More than at maximum 4
parameters were passed in case of more than 4 features to plot.
- Bugfix: Works now with more than 8 features.
The following commands will install the recent version and show the help
text:
<hansl>
pkg install PairPlot
include PairPlot.gfn
help PairPlot
</hansl>
Have fun!
Artur
1 year, 8 months
A'la selifc() for list
by Marcin Błażejowski
Hi,
do we have a function similar to selifc() but defined for lists?
Marcin
--
Marcin Błażejowski
1 year, 8 months
Fwd: Graphing Impulse response for AR(1) univariate model
by Cottrell, Allin
The OP's message below was sent to gretl-users-owner, a backdoor way
of trying to get the attention of gretl people ;-)
I'm forwarding it to the gretl-users list, where it should have been
posted, and where might be of interest to others.
Here's Jack Lucchetti's reply to the original email:
<jack>
On Tue, 18 Jan 2022, Matt Marcoux wrote:
> Hi all,
>
> Looking over the documentation for GRETL, and the IRF function will only
> generate an irf after a var/vecm model is estimated.
> However, I am looking to graph the irf for an AR(1) model, Here is my code:
>
> #open data file
>
> open ARMAexamples.csv
>
> #univariate time series AR(1)
>
> arma 1 0; X --nc
>
> #matrix ir1 = irf(0,1)
>
> In stata, we are able to do something similar by using this code. Anyone know
> how one might go about graphing an IRF for an AR(1) model in gretl?
>
> I am looking to create a graph that looks something like this:
>
> Thanks for taking the time to read this!
You'll get what you want from the ARMA_IRF function package.
</jack>
I'll just add: if the term "function package" doesn't resonate with
anyone, please take a look at the "Packages" item under the Help menu
in gretl's main window.
--
Allin Cottrell
Department of Economics
Wake Forest University
1 year, 8 months
