Re: [R-pkg-devel] [External] Re: incomplete gamma function Fortran subroutine
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
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
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
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