Re: [R-pkg-devel] [External] Re: incomplete gamma function Fortran subroutine

2020-07-21 Thread Ben Bolker
  I didn't find it until I said "before I compose a snarky reply to the 
list I'd better search for "Fortran" through *all* of WRE to make sure 
it's not there somewhere ..."


On 7/21/20 5:01 PM, Wang, Zhu wrote:

Indeed, that is the piece I missed completely! Glad that I didn't try to invent 
wheels.
Thanks,
Zhu

-Original Message-
From: Ben Bolker 
Sent: Tuesday, July 21, 2020 3:58 PM
To: luke-tier...@uiowa.edu; Wang, Zhu 
Cc: John P. Nolan ; r-package-devel@r-project.org
Subject: Re: [External] Re: [R-pkg-devel] incomplete gamma function Fortran 
subroutine

    Interesting.  It looks like the way to make use of this is to combine the information in 
section 6.7.1 ("Distribution functions", including pgamma ...) with section 6.6 
"Calling C from Fortran and vice versa" ?

On 7/21/20 4:33 PM, luke-tier...@uiowa.edu wrote:

Looking at this section of WRE may help:

https://cran.r-project.org/doc/manuals/r-release/R-exts.html#Numerical
-analysis-subroutines


Best,

luke

On Tue, 21 Jul 2020, Wang, Zhu wrote:


Thanks to Ben and John. Perhaps the program should call pgamma in R
and pass the number to Fortran. Calling some Fortran subroutines
older than R can trigger concerns when submitting the package to the
CRAN.

Best,
Zhu

-Original Message-
From: John P. Nolan 
Sent: Tuesday, July 21, 2020 1:43 PM
To: Wang, Zhu ; Ben Bolker ;
r-package-devel@r-project.org
Subject: RE: [R-pkg-devel] incomplete gamma function Fortran
subroutine

As others have said, built-in function pgamma is a (normalized)
version of the incomplete gamma function!   John

-Original Message-
From: R-package-devel  On
Behalf Of Wang, Zhu
Sent: Tuesday, July 21, 2020 2:16 PM
To: Ben Bolker ; r-package-devel@r-project.org
Subject: Re: [R-pkg-devel] incomplete gamma function Fortran
subroutine

External Email: Use caution with links and attachments.

Sorry for not making myself clear: The Fortran subroutine in an R
package needs to call incomplete gamma function.

-Original Message-
From: R-package-devel  On
Behalf Of Ben Bolker
Sent: Tuesday, July 21, 2020 12:54 PM
To: r-package-devel@r-project.org
Subject: Re: [R-pkg-devel] incomplete gamma function Fortran
subroutine

     Is there a reason not to use pgamma(), possibly adjusted by a
(de-)normalizing constant?   (See detailed notes in ?pgamma)

On 7/21/20 1:44 PM, Wang, Zhu wrote:

Hello,

In an R function within a package, I would like to call a Fortran
subroutine to compute lower gamma function. Any advice will be
appreciated.

Thanks!

Zhu Wang

   [[alternative HTML version deleted]]

__
R-package-devel@r-project.org mailing list
https://urldefense.proofpoint.com/v2/url?u=https-3A__stat.ethz.ch_ma
il
man_listinfo_r-2Dpackage-2Ddevel=DwIGaQ=U0G0XJAMhEk_X0GAGzCL7Q
=7
rQvU8hscCTWlvO-F5wI2-2eTiW40XI5qUKda0AnbG0=Y5sTjoEyQhEnvYqP-rmi1Pm
vi
Z_5jj7ur9P8ujvLiBc=w0p6b_yBQ1jDH3amMKQGvmEKYJD-BAid_CphFO37yu0=

__
R-package-devel@r-project.org mailing list
https://urldefense.proofpoint.com/v2/url?u=https-3A__stat.ethz.ch_mai
lman_listinfo_r-2Dpackage-2Ddevel=DwIGaQ=U0G0XJAMhEk_X0GAGzCL7Q
=7rQvU8hscCTWlvO-F5wI2-2eTiW40XI5qUKda0AnbG0=Y5sTjoEyQhEnvYqP-rmi1P
mviZ_5jj7ur9P8ujvLiBc=w0p6b_yBQ1jDH3amMKQGvmEKYJD-BAid_CphFO37yu0
= __
R-package-devel@r-project.org mailing list
https://urldefense.proofpoint.com/v2/url?u=https-3A__stat.ethz.ch_mai
lman_listinfo_r-2Dpackage-2Ddevel=DwIGaQ=U0G0XJAMhEk_X0GAGzCL7Q
=7rQvU8hscCTWlvO-F5wI2-2eTiW40XI5qUKda0AnbG0=Y5sTjoEyQhEnvYqP-rmi1P
mviZ_5jj7ur9P8ujvLiBc=w0p6b_yBQ1jDH3amMKQGvmEKYJD-BAid_CphFO37yu0
= __
R-package-devel@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-package-devel



__
R-package-devel@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-package-devel


Re: [R-pkg-devel] [External] Re: incomplete gamma function Fortran subroutine

2020-07-21 Thread Wang, Zhu
Indeed, that is the piece I missed completely! Glad that I didn't try to invent 
wheels.
Thanks,
Zhu

-Original Message-
From: Ben Bolker  
Sent: Tuesday, July 21, 2020 3:58 PM
To: luke-tier...@uiowa.edu; Wang, Zhu 
Cc: John P. Nolan ; r-package-devel@r-project.org
Subject: Re: [External] Re: [R-pkg-devel] incomplete gamma function Fortran 
subroutine

   Interesting.  It looks like the way to make use of this is to combine the 
information in section 6.7.1 ("Distribution functions", including pgamma ...) 
with section 6.6 "Calling C from Fortran and vice versa" ?

On 7/21/20 4:33 PM, luke-tier...@uiowa.edu wrote:
> Looking at this section of WRE may help:
>
> https://cran.r-project.org/doc/manuals/r-release/R-exts.html#Numerical
> -analysis-subroutines
>
>
> Best,
>
> luke
>
> On Tue, 21 Jul 2020, Wang, Zhu wrote:
>
>> Thanks to Ben and John. Perhaps the program should call pgamma in R 
>> and pass the number to Fortran. Calling some Fortran subroutines 
>> older than R can trigger concerns when submitting the package to the 
>> CRAN.
>>
>> Best,
>> Zhu
>>
>> -Original Message-
>> From: John P. Nolan 
>> Sent: Tuesday, July 21, 2020 1:43 PM
>> To: Wang, Zhu ; Ben Bolker ; 
>> r-package-devel@r-project.org
>> Subject: RE: [R-pkg-devel] incomplete gamma function Fortran 
>> subroutine
>>
>> As others have said, built-in function pgamma is a (normalized) 
>> version of the incomplete gamma function!   John
>>
>> -Original Message-
>> From: R-package-devel  On 
>> Behalf Of Wang, Zhu
>> Sent: Tuesday, July 21, 2020 2:16 PM
>> To: Ben Bolker ; r-package-devel@r-project.org
>> Subject: Re: [R-pkg-devel] incomplete gamma function Fortran 
>> subroutine
>>
>> External Email: Use caution with links and attachments.
>>
>> Sorry for not making myself clear: The Fortran subroutine in an R 
>> package needs to call incomplete gamma function.
>>
>> -Original Message-
>> From: R-package-devel  On 
>> Behalf Of Ben Bolker
>> Sent: Tuesday, July 21, 2020 12:54 PM
>> To: r-package-devel@r-project.org
>> Subject: Re: [R-pkg-devel] incomplete gamma function Fortran 
>> subroutine
>>
>>     Is there a reason not to use pgamma(), possibly adjusted by a 
>> (de-)normalizing constant?   (See detailed notes in ?pgamma)
>>
>> On 7/21/20 1:44 PM, Wang, Zhu wrote:
>>> Hello,
>>>
>>> In an R function within a package, I would like to call a Fortran 
>>> subroutine to compute lower gamma function. Any advice will be 
>>> appreciated.
>>>
>>> Thanks!
>>>
>>> Zhu Wang
>>>
>>>   [[alternative HTML version deleted]]
>>>
>>> __
>>> R-package-devel@r-project.org mailing list 
>>> https://urldefense.proofpoint.com/v2/url?u=https-3A__stat.ethz.ch_ma
>>> il
>>> man_listinfo_r-2Dpackage-2Ddevel=DwIGaQ=U0G0XJAMhEk_X0GAGzCL7Q
>>> =7 
>>> rQvU8hscCTWlvO-F5wI2-2eTiW40XI5qUKda0AnbG0=Y5sTjoEyQhEnvYqP-rmi1Pm
>>> vi 
>>> Z_5jj7ur9P8ujvLiBc=w0p6b_yBQ1jDH3amMKQGvmEKYJD-BAid_CphFO37yu0=
>>
>> __
>> R-package-devel@r-project.org mailing list 
>> https://urldefense.proofpoint.com/v2/url?u=https-3A__stat.ethz.ch_mai
>> lman_listinfo_r-2Dpackage-2Ddevel=DwIGaQ=U0G0XJAMhEk_X0GAGzCL7Q
>> =7rQvU8hscCTWlvO-F5wI2-2eTiW40XI5qUKda0AnbG0=Y5sTjoEyQhEnvYqP-rmi1P
>> mviZ_5jj7ur9P8ujvLiBc=w0p6b_yBQ1jDH3amMKQGvmEKYJD-BAid_CphFO37yu0
>> = __
>> R-package-devel@r-project.org mailing list 
>> https://urldefense.proofpoint.com/v2/url?u=https-3A__stat.ethz.ch_mai
>> lman_listinfo_r-2Dpackage-2Ddevel=DwIGaQ=U0G0XJAMhEk_X0GAGzCL7Q
>> =7rQvU8hscCTWlvO-F5wI2-2eTiW40XI5qUKda0AnbG0=Y5sTjoEyQhEnvYqP-rmi1P
>> mviZ_5jj7ur9P8ujvLiBc=w0p6b_yBQ1jDH3amMKQGvmEKYJD-BAid_CphFO37yu0
>> = __
>> R-package-devel@r-project.org mailing list 
>> https://stat.ethz.ch/mailman/listinfo/r-package-devel
>>
>
__
R-package-devel@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-package-devel


Re: [R-pkg-devel] [External] Re: incomplete gamma function Fortran subroutine

2020-07-21 Thread Ben Bolker
  Interesting.  It looks like the way to make use of this is to combine 
the information in section 6.7.1 ("Distribution functions", including 
pgamma ...) with section 6.6 "Calling C from Fortran and vice versa" ?


On 7/21/20 4:33 PM, luke-tier...@uiowa.edu wrote:

Looking at this section of WRE may help:

https://cran.r-project.org/doc/manuals/r-release/R-exts.html#Numerical-analysis-subroutines 



Best,

luke

On Tue, 21 Jul 2020, Wang, Zhu wrote:

Thanks to Ben and John. Perhaps the program should call pgamma in R 
and pass the number to Fortran. Calling some Fortran subroutines 
older than R can trigger concerns when submitting the package to the 
CRAN.


Best,
Zhu

-Original Message-
From: John P. Nolan 
Sent: Tuesday, July 21, 2020 1:43 PM
To: Wang, Zhu ; Ben Bolker ; 
r-package-devel@r-project.org

Subject: RE: [R-pkg-devel] incomplete gamma function Fortran subroutine

As others have said, built-in function pgamma is a (normalized) 
version of the incomplete gamma function!   John


-Original Message-
From: R-package-devel  On 
Behalf Of Wang, Zhu

Sent: Tuesday, July 21, 2020 2:16 PM
To: Ben Bolker ; r-package-devel@r-project.org
Subject: Re: [R-pkg-devel] incomplete gamma function Fortran subroutine

External Email: Use caution with links and attachments.

Sorry for not making myself clear: The Fortran subroutine in an R 
package needs to call incomplete gamma function.


-Original Message-
From: R-package-devel  On 
Behalf Of Ben Bolker

Sent: Tuesday, July 21, 2020 12:54 PM
To: r-package-devel@r-project.org
Subject: Re: [R-pkg-devel] incomplete gamma function Fortran subroutine

    Is there a reason not to use pgamma(), possibly adjusted by a 
(de-)normalizing constant?   (See detailed notes in ?pgamma)


On 7/21/20 1:44 PM, Wang, Zhu wrote:

Hello,

In an R function within a package, I would like to call a Fortran 
subroutine to compute lower gamma function. Any advice will be 
appreciated.


Thanks!

Zhu Wang

  [[alternative HTML version deleted]]

__
R-package-devel@r-project.org mailing list
https://urldefense.proofpoint.com/v2/url?u=https-3A__stat.ethz.ch_mail
man_listinfo_r-2Dpackage-2Ddevel=DwIGaQ=U0G0XJAMhEk_X0GAGzCL7Q=7
rQvU8hscCTWlvO-F5wI2-2eTiW40XI5qUKda0AnbG0=Y5sTjoEyQhEnvYqP-rmi1Pmvi
Z_5jj7ur9P8ujvLiBc=w0p6b_yBQ1jDH3amMKQGvmEKYJD-BAid_CphFO37yu0=


__
R-package-devel@r-project.org mailing list 
https://urldefense.proofpoint.com/v2/url?u=https-3A__stat.ethz.ch_mailman_listinfo_r-2Dpackage-2Ddevel=DwIGaQ=U0G0XJAMhEk_X0GAGzCL7Q=7rQvU8hscCTWlvO-F5wI2-2eTiW40XI5qUKda0AnbG0=Y5sTjoEyQhEnvYqP-rmi1PmviZ_5jj7ur9P8ujvLiBc=w0p6b_yBQ1jDH3amMKQGvmEKYJD-BAid_CphFO37yu0=

__
R-package-devel@r-project.org mailing list 
https://urldefense.proofpoint.com/v2/url?u=https-3A__stat.ethz.ch_mailman_listinfo_r-2Dpackage-2Ddevel=DwIGaQ=U0G0XJAMhEk_X0GAGzCL7Q=7rQvU8hscCTWlvO-F5wI2-2eTiW40XI5qUKda0AnbG0=Y5sTjoEyQhEnvYqP-rmi1PmviZ_5jj7ur9P8ujvLiBc=w0p6b_yBQ1jDH3amMKQGvmEKYJD-BAid_CphFO37yu0=

__
R-package-devel@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-package-devel





__
R-package-devel@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-package-devel


Re: [R-pkg-devel] [External] Re: incomplete gamma function Fortran subroutine

2020-07-21 Thread luke-tierney

Looking at this section of WRE may help:

https://cran.r-project.org/doc/manuals/r-release/R-exts.html#Numerical-analysis-subroutines

Best,

luke

On Tue, 21 Jul 2020, Wang, Zhu wrote:


Thanks to Ben and John. Perhaps the program should call pgamma in R and pass 
the number to Fortran. Calling some Fortran subroutines older than R can 
trigger concerns when submitting the package to the CRAN.

Best,
Zhu

-Original Message-
From: John P. Nolan 
Sent: Tuesday, July 21, 2020 1:43 PM
To: Wang, Zhu ; Ben Bolker ; 
r-package-devel@r-project.org
Subject: RE: [R-pkg-devel] incomplete gamma function Fortran subroutine

As others have said, built-in function pgamma is a (normalized) version of the 
incomplete gamma function!   John

-Original Message-
From: R-package-devel  On Behalf Of 
Wang, Zhu
Sent: Tuesday, July 21, 2020 2:16 PM
To: Ben Bolker ; r-package-devel@r-project.org
Subject: Re: [R-pkg-devel] incomplete gamma function Fortran subroutine

External Email: Use caution with links and attachments.

Sorry for not making myself clear: The Fortran subroutine in an R package needs 
to call incomplete gamma function.

-Original Message-
From: R-package-devel  On Behalf Of Ben 
Bolker
Sent: Tuesday, July 21, 2020 12:54 PM
To: r-package-devel@r-project.org
Subject: Re: [R-pkg-devel] incomplete gamma function Fortran subroutine

Is there a reason not to use pgamma(), possibly adjusted by a 
(de-)normalizing constant?   (See detailed notes in ?pgamma)

On 7/21/20 1:44 PM, Wang, Zhu wrote:

Hello,

In an R function within a package, I would like to call a Fortran subroutine to 
compute lower gamma function. Any advice will be appreciated.

Thanks!

Zhu Wang

  [[alternative HTML version deleted]]

__
R-package-devel@r-project.org mailing list
https://urldefense.proofpoint.com/v2/url?u=https-3A__stat.ethz.ch_mail
man_listinfo_r-2Dpackage-2Ddevel=DwIGaQ=U0G0XJAMhEk_X0GAGzCL7Q=7
rQvU8hscCTWlvO-F5wI2-2eTiW40XI5qUKda0AnbG0=Y5sTjoEyQhEnvYqP-rmi1Pmvi
Z_5jj7ur9P8ujvLiBc=w0p6b_yBQ1jDH3amMKQGvmEKYJD-BAid_CphFO37yu0=


__
R-package-devel@r-project.org mailing list 
https://urldefense.proofpoint.com/v2/url?u=https-3A__stat.ethz.ch_mailman_listinfo_r-2Dpackage-2Ddevel=DwIGaQ=U0G0XJAMhEk_X0GAGzCL7Q=7rQvU8hscCTWlvO-F5wI2-2eTiW40XI5qUKda0AnbG0=Y5sTjoEyQhEnvYqP-rmi1PmviZ_5jj7ur9P8ujvLiBc=w0p6b_yBQ1jDH3amMKQGvmEKYJD-BAid_CphFO37yu0=
__
R-package-devel@r-project.org mailing list 
https://urldefense.proofpoint.com/v2/url?u=https-3A__stat.ethz.ch_mailman_listinfo_r-2Dpackage-2Ddevel=DwIGaQ=U0G0XJAMhEk_X0GAGzCL7Q=7rQvU8hscCTWlvO-F5wI2-2eTiW40XI5qUKda0AnbG0=Y5sTjoEyQhEnvYqP-rmi1PmviZ_5jj7ur9P8ujvLiBc=w0p6b_yBQ1jDH3amMKQGvmEKYJD-BAid_CphFO37yu0=
__
R-package-devel@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-package-devel



--
Luke Tierney
Ralph E. Wareham Professor of Mathematical Sciences
University of Iowa  Phone: 319-335-3386
Department of Statistics andFax:   319-335-3017
   Actuarial Science
241 Schaeffer Hall  email:   luke-tier...@uiowa.edu
Iowa City, IA 52242 WWW:  http://www.stat.uiowa.edu

__
R-package-devel@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-package-devel