similar to: (no subject)

While developing updates to optimization tools for R (I was behind 3 of the 5 codes in optim(), but Brian Ripley did the implementation), I've been seeing this kind of error when the objective function cannot be computed so returns NA. Examples: attempts to divide by 0 or sqrt(-ve) or log(0). A crude work-around is to detect bad situations when evaluating the function and return a very

Dear R users I have a data frame that looks like this: m j X1 X2 1 1 1 2 1 ... 1 J 2 1 2 2 2 ... 2 J . . . M 1 M 2 M ... M J The data frame essentially looks like a panel: m and j are indicators where m is 1:M and j is 1:J (and M > J); X1 and X2 are variables. I

Dear R community, I am using the R package `interval` in order to perform some modelling tests of the NPMLE convergence in the case of censoring. So all I am doing is drawing a sample from exponential distribution, making it a censored sample and computing the NPMLE of its distribution function. But when run on Linux Calculate 10.4 the program keeps crashing and reporting a segmentation fault

> version _ platform i686-pc-linux-gnu arch i686 os linux-gnu system i686, linux-gnu status major 1 minor 7.1 year 2003 month 06 day 16 language R Bug: integrate(f,lower,upper,extra_args) where f <- function(x,extra_args) { body } integrate doesn't pass the extra arguments when calling f. As a first check of this finding I integrated dnorm from

When I run gnls I get the error: Error in nls(y ~ cbind(1, 1/(1 + exp((xmid - x)/exp(lscal)))), data = xy, : step factor 0.000488281 reduced below 'minFactor' of 0.000976563 My first thought was to decrease minFactor but gnlsControl does not contain minFactor nor nlsMinFactor (see below). It does however contain nlsMaxIter and nlsTol which I assume are the analogs of

Rembert.DeBlander@econ.kuleuven.ac.be wrote: >>version > PLEASE: a) Use a sensible subject line b) Please post only post bug reports if you are sure it is a bug (and here it is not a bug!) c) Check the recent R version (1.9.1) if posting bugs, not an outdated one! d) Only named arguments are passed through the ... argument: integrate(dnorm, lower=-Inf, upper=3, mean=-1, sd=7) Uwe

Hi, I'm having a problem converting a Matlab program into R. The R code works almost all the time, but about 4% of the time R's optim function gets stuck on a local minimum whereas matlab's fminsearch function does not (or at least fminsearch finds a better minimum than optim). My understanding is that both functions default to Nelder-Mead optimization, but what's different about

Hi everyone, I am looking to use R as a MATLAB replacement for linear algebra. I've done a fairly good job for finding replacements for most of the functions I'm interested in, I John Fox wrote a program for implementing the reduced row echelon form of a matrix (by doing the Gauss-Jordan elimination). I modified it a bit: rref <- function(A,

Dear Michael, One key is adjustment of nls optimizer tolerance. I notice it has to be higher than usual, but, I recovered your noisy "known" parameter values with an error of K1 (-7%) and k1 (-6%): #### Miller problem with Dalgaard modifications ## Linares 1/22/2004 ## Solution 1 nls(noisy ~ lsoda(xstart, time, one.compartment.model, c(K1=K1, k2=k2))[,2], data=C1.lsoda,

Hi, I got the following problem. I am doing a maximum likelihood estimation for a Kalman Filter. For this purpose, I have to invert an error matrix Ffast of dimension "no. parameters X no.parameters". The usualy optim methods often find only local minima, so I decided to make the optimization using the SANN algorithm, which is very slow already. However, this becomes a real problem

Hello Trying to develop variograms and kriged surfaces from a point file. Here is what I've done so far. library(gstat) # also loads library(sp) library(lattice) soilpts$x <- soilpts$UTM_X soilpts$y <- soilpts$UTM_Y soil.dat <- subset(soilpts, select=c(x, y, Area, BulkDensity, LOI, TP, TN, TC, Total_Mg)) dim(soil.dat) [1] 1292 7 coordinates(soil.dat) <- ~ x+y

Hi All, Can someone tell me if solve function shown below for my version of R is proper or not? I am using R 2.7.2 .Wherever i have used this function ,i got results which were different from the expected results as computed using SPSS. Description of this function says:-- Solve a System of EquationsDescription This generic function solves the equation a %*% x = b for x, where b can be either

Hi Rick Thanks to Dieter Menne I did manage to solve the problem of imposing bounds on the parameter space duirng an nlsList fit. He suggested using optim to optimize the parameters prior to each fit. This worked well for me as I had a customized selfStart function that then optimized the parameters for each individual separately. if you rewrote your selfStart function as: powrDpltInit <-

I have been told by the CRAN administrators that the following code generated an error on 64-bit Fedora Linux (gcc, clang) and on Solaris machines (sparc, x86), but runs well on all other systems): > fn <- function(x, y) ifelse(x^2 + y^2 <= 1, 1 - x^2 - y^2, 0) > tol <- 1.5e-8 > fy <- function(x) integrate(function(y) fn(x, y), 0, 1,

To Whom It May Concern: I have an R package and want to put this package be part of R-project and available to anyone who is interested in. The R package is created for my paper, titled "Acceptance Sampling Plans from Truncated Life Tests Based on the Birnbaum-Saunders Distribution for Percentiles". The paper has been accepted by Communications in Statistics: Simulation and

The only observation I can make is that the change to round() was made in r77727 whereas your R-devel appears to be r77715 (so would not exhibit the fixed behaviour). My guess is that there was a perpetual installation failure after r77715 but that the test folder was still retrieved and used. On Sat, 8 Feb 2020 at 19:27, Berwin A Turlach <berwin.turlach at gmail.com> wrote: > >

Dear list, I have a problem with nested functions and I don't manage to get it solved. I know I should be looking in environments, and I have tried a lot, but it keeps on erroring. An easy version of the problem is as follows: innerfunction<-function() { print(paste(a, " from inner function")) print(paste(b, " from inner function")) setwd(wd) }

On Sat, Feb 8, 2020 at 9:27 AM Berwin A Turlach <berwin.turlach at gmail.com> wrote: > > G'day all, > > I have daily scripts running to install the patched version of the > current R version and the development version of R on my linux box > (Ubuntu 18.04.4 LTS). > > The last development version that was successfully compiled and > installed was "R Under

Solaris 2.6, R Version 0.63.0 (November 14, 1998) Version 0.62.4 (October 24, 1998) > seq(0.15, 0.70, 0.05) [1] 0.15 0.20 0.25 0.30 0.35 0.40 0.45 0.50 0.55 0.60 0.65 The much reviled (on r-help) 0.62.3 got such simple cases right (even if 0.62.1 did not). Jim Lindsey insisted that such bugs are reported here, so I am. And yes, this is a real example and leads to an

Greetings, I would like to calculate tolerancelimits for a series of predicted values from nonlinear regression models. I've been using the tolerance package but the self-starting functions are not in the derivitive tables. When I spellout the functions and supply starting values I repeatedly get an error message regardless of the starting values I use (see output below). Does anyone