similar to: nlminb with constraints failing on some platforms

Hello, R 3.5.2 on ubuntu 18.04. sessionInfo() at the end. Works with me, same results, cannot reproduce the error. f <- function(x) sum( log(diff(x)^2+.01) + (x[1]-1)^2 ) opt <- nlminb(rep(0, 10), f, lower=-1, upper=3) str(opt) xhat <- rep(1, 10) all.equal(opt$par, xhat, tol=0) # good: 5.53 e-7 #[1] "Mean relative difference: 5.534757e-07" all.equal(opt$objective,

> On 1 Feb 2019, at 10:00, Martin Maechler <maechler at stat.math.ethz.ch> wrote: > > ........ >>> sessionInfo() >> R version 3.5.2 (2018-12-20) >> Platform: x86_64-pc-linux-gnu (64-bit) >> Running under: Scientific Linux release 6.4 (Carbon) > >> Matrix products: default >> BLAS/LAPACK:

Prof Nash, Prof Galanos Is it possible to use a generic code stub in front of packages that use optimx to improve optimx use or curtail it according to the requirements? Best Regards Amit +91 7899381263 ________________________________________________________________________ Please request Skype as available 5th Year FPM (Ph.D.) in Finance and Accounting Area Indian Institute

>>>>> Kasper Kristensen via R-devel >>>>> on Mon, 28 Jan 2019 08:56:39 +0000 writes: > I've noticed unstable behavior of nlminb on some Linux > systems. The problem can be reproduced by compiling > R-3.5.2 using gcc-8.2 and running the following snippet: > f <- function(x) sum( log(diff(x)^2+.01) + (x[1]-1)^2 ) > opt

No error on Windows 10, R.3.5.2 patched, Rblas compiled with OpenBLAS 0.20, Rlapack is base. > f <- function(x) sum( log(diff(x)^2+.01) + (x[1]-1)^2 ) > opt <- nlminb(rep(0, 10), f, lower=-1, upper=3) > str(opt) List of 6 $ par : num [1:10] 1 1 1 1 1 ... $ objective : num -41.4 $ convergence: int 0 $ iterations : int 66 $ evaluations: Named int [1:2] 96 830 ..-

This is not about the failure on some platforms, which is an important issue. However, what is below may provide a temporary workaround until the source of the problem is uncovered. FWIW, the problem seems fairly straightforward for most optimizers at my disposal in the R-forge (developmental) version of the optimx package at https://r-forge.r-project.org/projects/optimizer/ I used the code ##

I get the failure message. To be specific: adcomp.git>R CMD BATCH --quiet test_nlminb.R adcomp.git>cat test_nlminb.Rout > f <- function(x) sum( log(diff(x)^2+.01) + (x[1]-1)^2 ) > opt <- nlminb(rep(0, 10), f, lower=-1, upper=3) > xhat <- rep(1, 10) > abs( opt$objective - f(xhat) ) < 1e-4? ## Must be TRUE [1] FALSE My system is described by: adcomp.git>uname

Thank you, Brad (and others), >>>>> Brad Bell on Mon, 4 Feb 2019 07:21:18 -0700 writes: > I get the failure message. To be specific: adcomp.git> R CMD BATCH --quiet test_nlminb.R adcomp.git> cat test_nlminb.Rout >> f <- function(x) sum( log(diff(x)^2+.01) + (x[1]-1)^2 ) >> opt <- nlminb(rep(0, 10), f, lower=-1, upper=3) >>

Dear all, I'm experiencing an unusual installation error for one package. Could anyone suggest how I can best investigate this from here please? I'm sorry this isn't very much to go on. Hopefully someone can point me in the right direction. The problem seems to be my library. It is large (3,418 packages) and I'd like to avoid rebuilding it, if possible. All packages are up to

Hi, Is there a faster way to do this? It takes forever, even on a moderately sized dataset. n <- dim(dsn)[1] dsn2 <- dsn[order(-dsn$xhat),] dsn2[1, "cumx"] <- dsn2[1, "xhat"] for (i in 2:n) { dsn2[i, "cumx"] <- dsn2[i - 1, "cumx"] + dsn2[i, "xhat"] } [[alternative HTML version deleted]]

Hello, I experienced a compiler error when I tried to compile the latest version of R i.e. R4.0.2 making iosupport.d from iosupport.c making lapack.d from lapack.c making list.d from list.c making localecharset.d from localecharset.c grep.c(74): catastrophic error: cannot open source file "pcre2.h" # include<pcre2.h> (The pcre2.h header file is actually present!) I used the

Hi all, I am noticing some strange behaviour so I am bringing to this list. In the past when I have submitted bugs to bugzilla, I have come here first for confirmation/advice. Hopefully this is appropriate. Upgrading from R 3.5.2 to R 3.5.3 seems to have elicited some strange behaviour on my Windows machine. R seems to have trouble spawning a new process on my machine. You can noticing with all

On Wed, Sep 9, 2015 at 11:26 PM, Nathan Kurz <nate at verse.com> wrote: > > As a short and simple approach, I just compiled the current R release > on Ubuntu with ICC and MKL using just this: > > $ tar -xzf R-3.2.2.tar.gz > $ cd R-3.2.2 > $ CC=icc CXX=icpc AR=xiar LD=xild CFLAGS="-g -O3 -xHost" CXXFLAGS="-g > -O3 -xHost" ./configure

On Tue, Sep 29, 2015 at 8:11 PM, Nathan Kurz <nate at verse.com> wrote: > Hi Arnaud -- > > I'm glad it's working for you. I'm not sure I understand your final > answer. Are you saying that the version I posted worked for you as > given, or that you had to remove some of the other options? I say it works perfectly when using the single dynamic library (lmkl_rt):

Current R release (3.5.2) and devel return a 0 exit status on error, while prior versions returned a non-zero exit status. On Linux and MacOs, the following line returns TRUE for R-3.5.2 and R-devel, and FALSE for R-3.5.1 and R-3.5.0: system2("Rscript", c("-e", shQuote("stop('foo')"))) == 0 I didn't find this in the NEWS, so I believe this is a bug.

So I am trying to use the mle command (from stats4 package) to estimate a number of parameters using data but it keeps throwing up this error message: Error in solve.default(oout$hessian) : Lapack routine dgesv: system is exactly singular: U[1,1] = 0 This error sometimes indicates that the list of starting values is too far from optimum but this is unlikely since I picked values close to where

> On 6 Feb 2019, at 10:58, Martin Maechler <maechler at stat.math.ethz.ch> wrote: > ..... > --------------------------------------------------------------------------- > > I summarize what has been reported till: > > Failure in these cases > ======== > 1. Kasper K ("Scientific Linux", self compiled R, using Intel's MKL > for BLAS/LAPACK)

Reporting a possible inconsistency or bug in handling stringsAsFactors in as.data.frame.table() Here is a simple test > options()$stringsAsFactors [1] TRUE > x<-c("a","b","c","a","b") > d<-as.data.frame(table(x)) > d x Freq 1 a 2 2 b 2 3 c 1 > class(d$x) [1] "factor" >

I'm trying to fit a model to solve a biological problem. There are multiple independent variables, and also there are multiple responses. Each response is a function of all the independent variables, plus a set of parameters. All the responses depend on the same variables and parameters - just the form of the function changes to define each seperate response. Any ideas how I can fit

Hola! I am on debian testing, all updates up to date. I am running R via emacs and ess, so I canot know if this is an R or emacs/ess problem. What happens is that I want to read an vignette, and calls something like vignette("mle2", package="bbmle") the pdf file opens (for me in Foxit reader), no problem, but then emacs starts to spew a lot of spam, interfering with the