Mann-Kendall test
by Sven Schreiber
Hi,
I've just read about the non-parametric Mann-Kendall test for a trend
(hence my recent postings about Kendall's tau and ties in a series).
It looked pretty easy, so I'm attaching an attempt at a basic
implementation. Full disclosure: I've used the description from
https://cran.r-project.org/web/packages/trend/vignettes/trend.pdf (only
the basic version in section 1.1).
All you have to do is call "MannKendall(x)", where x is your series of
interest.
I'd be happy about feedback regarding cross-checks and correctness (for
example in the above document in equ. (4) the case S>0 wrongly appears
twice, and I had to make an educated guess where to switch the
signs...), whether it should become a function package and so on.
thanks,
sven
8 years, 9 months
another ternary woe
by Sven Schreiber
Hi,
very short minimal script here:
<hansl>
matrix heyho = ("a" == "b") ? garbl[3,] : I(2)
</hansl>
There is no object "garbl", but that's by design in the case that the
condition isn't true. Gretl says "[" would be invalid here.
thanks,
sven
8 years, 10 months
args in lags() vs. mlag()
by Sven Schreiber
Hi,
I find it confusing in routine work that the lags() function wants to
have the lag number as the first arg, while the mlag() function wants to
have the input matrix first and afterwards the lag number.
Don't know how, maybe this could be aligned somehow in the future.
cheers,
sven
9 years
Variable names in summary statistics (and similar)
by Schaff, Frederik
Dear Gretl-Devs,
when "printing" summary statistics, correlation matrixes, etc., the usual limit of 31 characters for variables is not held. Instead, variables are truncated to 15 characters... In some cases making them indistinguishable. If it was possible to also use the limit of 31 chars here, that would be really nice.
Kind Regards
Frederik Schaff
------------
Dr. Frederik Schaff
FernUniversität in Hagen
Department of Business and Economics
Chair of Economic Theory (Prof. A. Endres)
Universitätsstraße 11 (TGZ)
58097 Hagen
Voice: +49 (0) 2331 987-4454
E-mail: Frederik.Schaff(a)FernUni-Hagen.de<mailto:Frederik.Schaff@FernUni-Hagen.de>
WWW: e.feu.de/schaff
9 years
minor report on formatted text import
by Sven Schreiber
Hi,
concerning the heuristics gretl uses when parsing formatted text data: I
just had a case with tab-delimited data where gretl gave up, and all I
had to do is to replace all tabs with semicolons.
The other format information: daily dates were given as dd.mm.yyyy, the
decimal separator was a comma. The date column had the label "Date".
I'm not saying this is a bug and I know heuristics can never be perfect.
Just wanted to let you know that gretl was almost getting there, there
was only one step left.
cheers,
sven
9 years
RFC: gretl calendar
by Allin Cottrell
Further follow-up to Ignacio's comment in
http://lists.wfu.edu/pipermail/gretl-devel/2017-February/007439.html
(As a reminder, Ignacio remarked that gretl was considering 1700 as a
leap year, which it wasn't on the Gregorian calendar, and also that
our weekday() function was producing some bad values for eighteenth
century dates.)
In git and snapshots I've now revamped gretl's calendar so that it's
consistently Gregorian by default. At the same time I've added some
apparatus for handling Julian dates, in case that should be needed for
dealing with historical data. Internally, we're now using GLib's
GDate API (with some additional algorithms for handling Julian dates)
in place of the old code based on the unix "cal" program.
The current setup is discussed in a new chapter 16, "Calendar dates",
in the User's Guide. I've partially updated the help for the relevant
calendrical functions in line with the new chapter, but not yet fully
since I'd like to wait and get people's reactions.
So, comments welcome.
Allin
9 years
gretl website question
by oleg_komashko＠ukr.net
Dear all,
Is Ukrainian translation official?
No mention about it at the site.
Oleh
9 years
function package warnings
by Allin Cottrell
For package authors: here's a listing of deprecation warnings produced
by packages currently on the server. If you can find time to update
your package(s) that would be much appreciated!
Note that the warning for using '=' as an equality test is very
explicit so it should be easy to find the location in the code.
* Missing "include" in sample script:
DHF_test.gfn
plot_scatter_label.gfn
* Use of "catch" on calls to user-defined functions
(generates a warning, really should be an error by now):
coint2rec.gfn
GHegy.gfn
* Use of '=' as test for equality (count of warnings
from running sample script):
ADFGLSTest.gfn 2
a_eff.gfn 34425
armax.gfn 9
bartlett.gfn 10
BMA.gfn 18000
BreitungCandelonTest.gfn 2
BSMestim.gfn 359
buys_ballot.gfn 1
Canova_Hansen.gfn 2
carr.gfn 49
cffilt.gfn 1
coint2rec.gfn 7
DiebMar.gfn 11
DIF_panelcoint.gfn 1003
dynMultipl.gfn 1
genr_dates.gfn 7
getQuandl.gfn 45
GJR-garchm.gfn 69
hetero.gfn 2
lagreg.gfn 12109
levene.gfn 2
LLTestim.gfn 427
margeff_el.gfn 13872
MGARCH.gfn 5
mwu.gfn 1
Robinson.gfn 1
robust_RESET.gfn 1
RVNR.gfn 2
SB.gfn 200
SETAR.gfn 50
sfa_mod.gfn 9
ShapleyCont.gfn 1
sparse.gfn 102
StrucBreak.gfn 1
summary_xy.gfn 191
tobithetm.gfn 1
tramolin.gfn 1
var_chow.gfn 2
VARrec.gfn 1
waldTest.gfn 26
Wilcoxon.gfn 1
Allin
9 years
Is SF git down?
by Riccardo (Jack) Lucchetti
or is it just me?
-------------------------------------------------------
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
-------------------------------------------------------
9 years
Re: [Gretl-devel] Gretl-devel Digest, Vol 121, Issue 20
by Ignacio Diaz-Emparanza
OK. I understand. Good to know that gretl is doing well.
It is interesting that in this url the calendar for Spain in 1700 is
different from the US calendar.
I assume the "cal" code maybe is doing some adjustments also for some
days in 1756 (and also some other years in the XVIII century). In
particular the gretl function 'weekday' produces negative values for
some days in may 1756 (they should be values only from 0(sunday) to
6(saturday)).
<script>
nulldata 1500
setobs 7 1755-01-01 --time-series
series d = weekday($obsmajor,$obsminor,$obsmicro)
</script>
El 16/02/17 a las 18:00, gretl-devel-request(a)lists.wfu.edu escribió:
> Gretl follows the unix "cal" program on this. In lib/src/calendar.c:
>
> /* leap year -- account for gregorian reformation in 1752 */
> #define leap_year(yr) \
> ((yr) <= 1752 ? !((yr) % 4) : \
> (!((yr) % 4) && ((yr) % 100)) || !((yr) % 400))
>
> See also
>
> https://www.timeanddate.com/calendar/index.html?year=1700
--
Ignacio Díaz-Emparanza
Departamento de Economía Aplicada III (Econometría y Estadística)
Universidad del País Vasco - Euskalherriko Unibertsitatea, UPV/EHU
Tfno: (+34) 94 601 3732
http://www.ehu.eus/ea3
9 years
