[Rd] gfortran 9 quantreg bug
bhh @end|ng |rom x@4@||@n|
Sun Aug 4 16:26:40 CEST 2019
I have run
gfortran -c -fsyntax-only -fimplicit-none -Wall -pedantic rqbr.f
in the src folder of quantreg.
There are many warnings about defined but not used labels.
Also two errors such as "Symbol ‘in’ at (1) has no IMPLICIT type".
And warnings such as: Warning: "Possible change of value in conversion from REAL(8) to INTEGER(4) at ..."
No offense intended but this fortran code is awful. I wouldn't want to debug this before an extensive cleanup by
getting rid of as many numerical labels as possible, indenting and doing something about the warnings "Possible change of value ...".
This is going to be very difficult.
> On 4 Aug 2019, at 08:48, Koenker, Roger W <rkoenker using illinois.edu> wrote:
> I’d like to solicit some advice on a debugging problem I have in the quantreg package.
> Kurt and Brian have reported to me that on Debian machines with gfortran 9
> f = summary(rq(foodexp ~ income, data = engel, tau = 1:4/5))
> fails because summary() produces bogus estimates of the coefficient bounds.
> This example has been around in my R package from the earliest days of R, and
> before that in various incarnations of S. The culprit is apparently rqbr.f which is
> even more ancient, but must have something that gfortran 9 doesn’t approve of.
> I note that in R-devel there have been some other issues with gfortran 9, but these seem
> unrelated to my problem. Not having access to a machine with an R/gfortran9
> configuration, I can’t apply my rudimentary debugging methods. I’ve considered
> trying to build gfortran on my mac air and then building R from source, but before
> going down this road, I wondered whether others had other suggestions, or
> advice about my proposed route. As far as I can see there are not yet
> binaries for gfortran 9 for osx.
> Roger Koenker
> r.koenker using ucl.ac.uk<mailto:r.koenker using ucl.ac.uk>
> Department of Economics, UCL
> London WC1H 0AX.
> [[alternative HTML version deleted]]
> R-devel using r-project.org mailing list
More information about the R-devel