Problem using printf
by Jan Tille
Hello to all and a (maybe late) Happy New Year.
I have a problem saving the contents of matrix into a file .
I used the tiny script:
outfile "matrix_test.txt" --write
printf "%#24.10g", test
outfile --close
While the matrix in Gretl is 497 x 1317 the printed matrix contains the first 133 rows only.
I have alternatively tried a right click on the matrix and use the "save as csv" option (with decimal separator set to ","). Pasting this into Excel, gives me 406 rows.
A right click on the matrix properties shows 497 rows and 1317 columns and when I open the matrix, Gretl also shows that all rows a filled with data (or nan).
In the most probable case, I have just made an error but and I would be glad, if someone has a hint.
Thanks in advance and kind regards
Jan
7 years, 2 months
Different Behavior Script and GUI
by juergen.malitte＠t-online.de
Hello,
I would make a note respective the different behavior of the command
"setobs n 1:1 --stacked-cross-section" by comparison with the appropriate
functionality in the GUI. I wouldn't say it's a mistake, but it evokes some
confusion.
Suppose you have a dataset "Umsatz_Panel.txt", which contains the sales of
i =4 units over t = 2 periods, for example:
1980
2110
2200
2312
2455
2380
4200
4560
Opening the file in the GUI you choose the data-structure "panel" and then
"gestapelte Querschnitte" (stacked cross-sectional). The result will be the
following:
1:1 1980
1:2 2110
1:3 2200
1:4 2312
2:1 2455
2:2 2380
2:3 4200
2:4 4560
That looks good.and the representation of the identifier has the format
"t:i" as expected (because the first staple of cross-sectional data belongs
to the period 1 etc.).
If you open the file "Umsatz_Panel.txt" in a script with the commands:
open "Umsatz_Panel.txt"
setobs 4 1:1 --stacked-cross-section
then you get the result
1:1 1980
1:2 2455
2:1 2110
2:2 2380
3:1 2200
3:2 4200
4:1 2312
4:2 4560
This is not the same representation as in the GUI, Gretl has changed the
format "t:i" into "i:t" !! I know, that "i:t" is the preferred and
standard representation for panel data in Gretl. If you store the above
representation "t:i" (received over the GUI) in a gdt-file and load this
again, you will then get the same result, will say the representation
"i:t".
First I was a little bit confused concerning the different behavior in the
GUI and in the script. But i think it should work identically !!
In respect of the command "setobs n 1:1 --stacked-time-series" the behavior
is the same as in the GUI !
Best wishes
Jürgen Malitte
7 years, 2 months
A question on squares()
by oleg_komashko＠ukr.net
Dear all, The example to explain the situation with comments and a question attached Oleh
7 years, 2 months
Regression coefficients
by Marc O'Callaghan
Dear all,
I am back with another, perhaps very nooby, question.
Namely:
I have successfully extracted 3rd-degree regressions from my temperature data. My question is now, how can I find out the exact coefficient of those regressions' x^3 ? More exactly, (how) can I make a new series out of the regression coefficients for each year ?
Hoping I am being clear and thanks in advance,
Best to all,
Marc O'C.
7 years, 2 months
gretl crashes during quantile regression
by Jan Tille
Dear list members,
I am trying to estimate several quantile regressions, each with a constant and a single regressor, on a rolling basis (497 steps). The quantile vector is set to q={0.1, 0.25, 0.5, 0.75, 0.90) yielding 5 quantile regressions per step and individual (1317). Using a simple weighting scheme, where the weights sum to 1, the estimated conditional quantiles ("y-hats") are aggregated so that this yield one estimate per step.
Overall, I have 16 regressors and at the end, I intend to obtain one matrix for each regressor. At the end of the day, this should amount to 497 (steps) * 1317(indiviuals) * 16(regressors) * 5(quantiles) = 5,236,392,0 quantile regressions.
Gretl works quite well for some time, but breaks down (i.e. does not react any more) somewhere during the process.
I suppose this occurs at the end, prior to returning the script output, which is quite lengthy (I do use --quiet option after quantreg command).
I tried to eliminate the 16th regressors, after which Gretl broke down the first time and redid the estimation using the other 15 regressors. This time, I saved the matrices using mwrite()-function. After the 15th round, Gretl broke down again. I opend the matrix number 15 in excel and saw that all rows and columns were filled with either data or, if the y-series was to short, with "nan".
Is there an option to tell Gretl to suppress the entire script output or save it as .txt-file instead?
I am not sure, maybe it is an issue with my data, because system performance (RAM 8GB) is not supposed to be the limiting factor, is it?
Another possibility is that there is a mistake in my script. But you already were so kind, to hint, that, last time, I forget an endloop. This time I also counted my loop - endloop combinations (and they match).
Thanks in advance.
Best regard,
Jan
7 years, 2 months
names in boxplot graph is incorrect (panel data)
by cociuba mihai
After setting the name to be used in graphs the boxplot shows the name of
the series.
<hansl>
clear
open nls_panel.gdt
#no description , no graph-name
boxplot lwage --output="display"
gnuplot lwage --time-series --with-lines --output="display"
#setting the graph-name
setinfo lwage --graph-name="lwage name to display"
#the name of the series is displayed correctly as per graph-name
gnuplot lwage --time-series --with-lines --output="display"
##the name of the series is displayed incorrectly as series name
boxplot lwage --output="display"
Mihai
7 years, 2 months
error when using the summary command (core dumped)
by cociuba mihai
Hello,
I've attached the dataset, the script that I use in order to import the
data and also the log, Gretl is build from sources (git, date of build
29.12.2015) and running on a Ubuntu 14.04.
Steps to reproduce the error:
-import the dataset, run
#######
<hansl>
clear
set skip_missing off
open "bazaextinsabankscopeFINAL.xls" --sheet=1 --coloffset=5 --quiet
matrix M={dataset}
#print M
summary
# hello gretl, this is a panel! unde N este periodicitatea
scalar C=cols(M)
scalar R=rows(M)
scalar N=R/9
print C R N
setobs N 1.1 --stacked-cross-section
series year = time + 2004
summary
#no error#
< end hansl>
###########
But when using the Summary option from the GUI menu gretl exits with the
error "core dumped"
The same dataset in R is working without any error (tested only the same
summary command).
Also when trying to run the script from
http://lists.wfu.edu/pipermail/gretl-users/2014-January/009553.html in
order to set the description to the variables I receive the following error:
#########
? string line
? scalar i = 1
Generated scalar i = 1
? loop while getline(D, line) -q
> setinfo i --description="@line"
> i++
> endloop
Command has insufficient arguments
>> setinfo i --description="@line"
##########
Any ideas of what seems to be problem?
Mihai
7 years, 2 months
list memory question
by oleg_komashko＠ukr.net
Dear all, A very simple question: how much memory take lists? Oleh
7 years, 2 months
OpenBLAS warning
by oleg_komashko＠ukr.net
Dear all, The warning after current GIT:
gretl is built using OpenMP, but is linked against
OpenBLAS parallelized via pthreads. This combination
of threading mechanisms is not recommended. Ideally,
OpenBLAS should also use OpenMP.
Configure options used (as always) ./configure --enable-openmp --enable-build-doc --enable-build-addons --with-odbc
Oleh
7 years, 2 months
Happy New Year!
by oleg_komashko＠ukr.net
Dear all, Happy New Year!
P.S. In my opinion it's absolutely senseless, but that's it's advantage: it could unite people through all over the World!
7 years, 2 months
