R devel - Apr 2006 - contribution offer: df() for non-centrality != 0

Hi R-devels,

I noticed that for the scheduled R 2.3.0,
there is still no function df() for
non-centrality != 0.

Easter is quickly approaching, so as my little
gift, I would like to offer you some code to fill
this gap:

After excluding some "unusual" cases,
it amounts to differentiating the corresponding
function pnf() --- which in turn leads to the
function dnbeta() already available in the
stats package for some time.

You may find this code in

http://www.uni-bayreuth.de/departments/math/org/mathe7/R-devel/dnf.tar.gz

This archive includes the C-File    dnf.c      
as well as some adapted versions     Fdist_new.Rd   
and     distn_new.R;
in the latter two files I indicated all modifications
w.r.t.  Fdist.Rd    and    distn.R     by
a corresponding %%-tag.

You may modify these files as you like in order to
conform your design principles

--- and hopefully we will soon have df() for ncp, too...

Best,
Peter

Thank you, Peter,

but note that this is clearly too late for 2.3.0 :
If you look at the release schedule at
   http://developer.R-project.org/ 
you see that  April 10 was "Feature Freeze".

If it looks good -- and it does after a quick glimpse --
it should be possible to be integrated for 2.3.1 though -- but
that is not at all scheduled yet. 

With regards, 
Martin

>>>>> "PetRd" == Peter Ruckdeschel
<Peter.Ruckdeschel at uni-bayreuth.de>
>>>>>     on Thu, 13 Apr 2006 15:15:11 +0200 writes:
    PetRd> Hi R-devels, I noticed that for the scheduled R
    PetRd> 2.3.0, there is still no function df() for
    PetRd> non-centrality != 0.

    PetRd> Easter is quickly approaching, so as my little gift,
    PetRd> I would like to offer you some code to fill this gap:

    PetRd> After excluding some "unusual" cases, it amounts to
    PetRd> differentiating the corresponding function pnf() ---
    PetRd> which in turn leads to the function dnbeta() already
    PetRd> available in the stats package for some time.

    PetRd> You may find this code in

    PetRd>
http://www.uni-bayreuth.de/departments/math/org/mathe7/R-devel/dnf.tar.gz

    PetRd> This archive includes the C-File dnf.c as well as
    PetRd> some adapted versions Fdist_new.Rd and distn_new.R;
    PetRd> in the latter two files I indicated all modifications
    PetRd> w.r.t.  Fdist.Rd and distn.R by a corresponding
    PetRd> %%-tag.

    PetRd> You may modify these files as you like in order to
    PetRd> conform your design principles

    PetRd> --- and hopefully we will soon have df() for ncp,
    PetRd> too...

Hi Martin,

Martin Maechler wrote:
> but note that this is clearly too late for 2.3.0 :
> If you look at the release schedule at
>    http://developer.R-project.org/ 
> you see that  April 10 was "Feature Freeze".
>   
Yes, don't worry, I have been aware of this;

I did not mean to submit the code for 2.3.0, but just used
this version as a reference  to be
"up-to-date".
> If it looks good -- and it does after a quick glimpse --
> it should be possible to be integrated for 2.3.1 though 
This would be fine !
Thank you

Peter