[Rd] Silent failure with NA results in fligner.test()

2020-12-19 Thread Karolis K
Hello, In certain cases fligner.test() returns NaN statistic and NA p-value. The issue happens when, after centering with the median, all absolute values become constant, which ten leads to identical ranks. Below are a few examples: # 2 groups, 2 values each # issue is caused by residual values

Re: [Rd] From .Fortran to .Call?

2020-12-19 Thread Ivan Krylov
On Sat, 19 Dec 2020 17:04:59 + "Koenker, Roger W" wrote: > There are comments in various places, including R-extensions §5.4 > suggesting that .Fortran is (nearly) deprecated and hinting that use > of .Call is more efficient and now preferred for packages. My understanding of §5.4 and 5.5 is

Re: [Rd] From .Fortran to .Call?

2020-12-19 Thread Bill Dunlap
To make C prototypes for routines defined in a *.f file you can use gfortran's -fc-prototypes-external flag. You still have to convert 'name_' to 'F77_NAME(name). bill@Bill-T490:~/packages/quantreg/src$ gfortran -fc-prototypes-external -fsyntax-only boot.f ... [elided defs of complex types] ... /

[Rd] From .Fortran to .Call?

2020-12-19 Thread Koenker, Roger W
There are comments in various places, including R-extensions §5.4 suggesting that .Fortran is (nearly) deprecated and hinting that use of .Call is more efficient and now preferred for packages. I understand and greatly appreciate its use with dyn.load() and R CMD SHLIB for development workflow.