Re: [Rd] unstable corner of parameter space for qbeta?

2020-03-26 Thread Martin Maechler
> Ben Bolker > on Wed, 25 Mar 2020 21:09:16 -0400 writes: > I've discovered an infelicity (I guess) in qbeta(): it's not a bug, > since there's a clear warning about lack of convergence of the numerical > algorithm ("full precision may not have been achieved"). I can work

Re: [Rd] unstable corner of parameter space for qbeta?

2020-03-26 Thread peter dalgaard
It's a pretty extreme case, try e.g. curve(pbeta(x, shape1, shape2), n=10001), and (probably -- I can't be bothered to work out the relation between beta shapes and F df parameters this morning...) outside what is normally encountered in statistical analyses. Notice though, that you have lower=F

Re: [Rd] unstable corner of parameter space for qbeta?

2020-03-26 Thread J C Nash
Given that a number of us are housebound, it might be a good time to try to improve the approximation. It's not an area where I have much expertise, but in looking at the qbeta.c code I see a lot of root-finding, where I do have some background. However, I'm very reluctant to work alone on this, an

Re: [Rd] unstable corner of parameter space for qbeta?

2020-03-26 Thread Martin Maechler
> J C Nash > on Thu, 26 Mar 2020 09:29:53 -0400 writes: > Given that a number of us are housebound, it might be a good time to try to > improve the approximation. It's not an area where I have much expertise, but in > looking at the qbeta.c code I see a lot of root-findi

Re: [Rd] unstable corner of parameter space for qbeta?

2020-03-26 Thread J C Nash
Despite the need to focus on pbeta, I'm still willing to put in some effort. But I find it really helps to have 2-3 others involved, since the questions back and forth keep matters moving forward. Volunteers? Thanks to Martin for detailed comments. JN On 2020-03-26 10:34 a.m., Martin Maechler w

Re: [Rd] unstable corner of parameter space for qbeta?

2020-03-26 Thread Ravi Varadhan
This is also strange: qbeta <- function (p, shape1, shape2, ncp = 0, lower.tail = TRUE, log.p = FALSE) { if (missing(ncp)) .Call(C_qbeta, p, shape1, shape2, lower.tail, log.p) else .Call(C_qnbeta, p, shape1, shape2, ncp, lower.tail, log.p) } Since the default value is

[Rd] Rebuilding and re-checking of downstream dependencies on CRAN Mac build machines

2020-03-26 Thread Winston Chang
I have two questions about the CRAN machines that build binary packages for Mac. When a new version of a package is released, (A) Do the downstream dependencies get re-checked? (B) Do the downstream dependencies get re-built? I have heard (but do not know for sure) that the answer to (A) is no

Re: [Rd] unstable corner of parameter space for qbeta?

2020-03-26 Thread Martin Maechler
> Ravi Varadhan > on Thu, 26 Mar 2020 18:33:43 + writes: > This is also strange: > qbeta <- function (p, shape1, shape2, ncp = 0, lower.tail = TRUE, log.p = FALSE) > { > if (missing(ncp)) > .Call(C_qbeta, p, shape1, shape2, lower.tail, log.p) > else .Call(

Re: [Rd] Rebuilding and re-checking of downstream dependencies on CRAN Mac build machines

2020-03-26 Thread Simon Urbanek
Winston, the Mac CRAN build builds a package only if either is true: 1) the package has not passed checks 2) there is a new version of the package since last successful build+check The old build machine doesn't have the capacity to do full re-builds (it would take over 24h - currently the nightl

Re: [Rd] Rebuilding and re-checking of downstream dependencies on CRAN Mac build machines

2020-03-26 Thread Hadley Wickham
If I do install.packages("dplyr", type = "source"), I see: Installing package into ‘/Users/hadley/R’ (as ‘lib’ is unspecified) trying URL 'https://cran.rstudio.com/src/contrib/dplyr_0.8.5.tar.gz' Content type 'application/x-gzip' length 1378766 bytes (1.3 MB) ==

Re: [Rd] Rebuilding and re-checking of downstream dependencies on CRAN Mac build machines

2020-03-26 Thread Winston Chang
Simon, The link I provided to the httpuv issue has detailed information about the error: https://github.com/rstudio/httpuv/issues/260 The error occurs when I run the following (although note, depending on which mirror you use, the recent CRAN server issues may result in problems downloading the

Re: [Rd] unstable corner of parameter space for qbeta?

2020-03-26 Thread Ben Bolker
On 2020-03-26 4:02 a.m., Martin Maechler wrote: >> Ben Bolker >> on Wed, 25 Mar 2020 21:09:16 -0400 writes: > > > I've discovered an infelicity (I guess) in qbeta(): it's not a bug, > > since there's a clear warning about lack of convergence of the numerical > > algorit

[Rd] Expressions from boxplot() passed to bxp()

2020-03-26 Thread Marius Hofert
Hi, Is this expected behavior (R-3.6.0)? dat <- cbind(x = 1:10, y = 10:1) ylab <- substitute(X[t], list(t = 2)) plot(dat, ylab = ylab) # works (correctly displays ylab) boxplot(dat, ylab = ylab) # fails boxplot(dat, ylab = as.expression(ylab)) # works Thanks & cheers, M