Juergen Hedderich
2001-Jan-09 08:44 UTC
[R] quantile function for noncentral f-distribution
hello R-friends, I'm looking for a quantile function for the noncentral f-distribution in the area of equivalence hypotheses testing. Can somebody help me? Many thanks ----------------------------------------------------------------- Dipl. Inform. J. Hedderich Institut f?r Medizinische Informatik Phone : 0431 / 5973182 und Statistik im Klinikum an der CAU Fax : 0431 / 5973193 Brunswiker Str. 10 24105 Kiel E-mail:hedderich at medinfo.uni-kiel.de ----------------------------------------------------------------- -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the subject !) To: r-help-request at stat.math.ethz.ch _._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._
Prof Brian Ripley
2001-Jan-09 09:00 UTC
[R] quantile function for noncentral f-distribution
On Tue, 9 Jan 2001, Juergen Hedderich wrote:> > hello R-friends, > I'm looking for a quantile function for the noncentral > f-distribution in the area of equivalence hypotheses testing. > Can somebody help me?Since pf has a non-centrality parameter, you can easily invert it. myqf <- function(p, df1, df2, ncp) { uniroot(function(x) pf(x, df1, df2, ncp) - p, , 0, 100)$root } adjust 100 to suit (one can surely come up with an useful computed upper bound if ncp might be large), and add a tol arg to uniroot if you care.> > Many thanks > > ----------------------------------------------------------------- > Dipl. Inform. J. Hedderich > Institut für Medizinische Informatik Phone : 0431 / > 5973182 > und Statistik im Klinikum an der CAU Fax : 0431 / 5973193 > > Brunswiker Str. > 10 > 24105 Kiel > E-mail:hedderich at medinfo.uni-kiel.de > ----------------------------------------------------------------- > -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- > r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html > Send "info", "help", or "[un]subscribe" > (in the "body", not the subject !) To: r-help-request at stat.math.ethz.ch > _._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._ >-- Brian D. Ripley, ripley at stats.ox.ac.uk Professor of Applied Statistics, http://www.stats.ox.ac.uk/~ripley/ University of Oxford, Tel: +44 1865 272861 (self) 1 South Parks Road, +44 1865 272860 (secr) Oxford OX1 3TG, UK Fax: +44 1865 272595 -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the subject !) To: r-help-request at stat.math.ethz.ch _._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._
Hi R-users! I'm building R package using R 1.2.0 for windows. To check the completeness of the program I used the following command: from c:/rw1020/rw1020 cd src/library Rcmd check meanscore (# `meanscore' is the name of a subdirectory contains the package.) However, I got the following error message (see the last line) * using log directory C:/rw1020/rw1020/src/meanscore.Rcheck * checking for file `library/meanscore/DESCRIPTION' ... OK * checking package directory ... OK * checking DESCRIPTION Package entry ... OK * checking DESCRIPTION Version entry ... OK * checking DESCRIPTION Title entry ... OK * checking DESCRIPTION License entry ... OK * checking DESCRIPTION Author entry ... OK * checking Rd files ... OK * checking for undocumented objects ... OK * creating meanscore-Ex.R ... ERROR But when I'm using the primitive command `make pkgcheck' from c:/rw1020/rw1020/src/gnuwin32, the meanscore-Ex.R file is successfully created and the output is correct. I've already checked the AUTOEXEC file and all the required paths are there. I'm using Perl 5.6.0 build 623 on windows 98 OS. The set of rtools I'm using is the one from http://www.stats.ox.ac.uk/pub/Rtools/tools.zip (downloaded 09/01/01). And I don't have C/Fortran codes in the package. Is this something to do with the `check' perlscript at R_HOME/bin or do I need to edit the Mkrules file? Any suggestion will be very appreciated, Agus ------------------------ Agus Salim PhD student Dept. of Statistics University College Cork, Ireland Tel: 353-21-4902197 ------------------------ -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the subject !) To: r-help-request at stat.math.ethz.ch _._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._