similar to: singular gradient error in nls

Why fails nls with "singular gradient" here? I post a minimal example on the bottom and would be very happy if someone could help me. Kind regards, ########### # define some constants smallc <- 0.0001 t <- seq(0,1,0.001) t0 <- 0.5 tau1 <- 0.02 # generate yy(t) yy <- 1/2 * ( 1- tanh((t - t0)/smallc) * exp(-t / tau1) ) + rnorm(length(t))*0.01 # show the curve

We obtained some disturbing results from convolve() (inaccuracies and negative probabilities). We'll try to make the context clear in as few lines as possible... Our function panjer() (code below) basically computes recursively the probability mass function of a compound Poisson distribution. When the Poisson parameter lambda is very large, the starting value of the recursive scheme ---

Hello List I'm moving this over from the bioC list as, although the problem I'm working on is biological, the current bottle neck is my poor understanding of R. I wonder if someone would help me with the following function. cumulMetric <- function(deMirPresGenes, deMirs){ ??? #need to match position of each miR in deMirPresGenes with its FC to form a vector of FC in correct order ?

Hi, I have a long list of lists from which I want to efficiently extract and rbind elements. So I'm using the approach below: f <- function(i){ out <- replicate(5, list(matrix(rnorm(80), nc=20))) names(out) <- letters[1:5] out } set.seed(1) lst <- lapply(1:1.5e6, f) (t0 <- system.time(tmp <- do.call(rbind, lapply(lst, '[[', 'b')))) Is there

Un texte encapsul? et encod? dans un jeu de caract?res inconnu a ?t? nettoy?... Nom : non disponible URL : <https://stat.ethz.ch/pipermail/r-help/attachments/20110517/1812d6c6/attachment.pl>

Hi all, My app is going to count sth based on poll id in the url. The url looks like: http://address.com?id=dwNKiItvcyrQ5sycnIhmJablDfXsc9tshaGIVyNIei7.e7&some_other_parameters In a controller, I''ve created a hash where the keys are the ids: @codes = { ''dwNKiItvcyrQ5sycnIhmJablDfXsc9tshaGIVyNIei7.e7'' => {:partner_id => 0,

Dear R-helpers, I am using confint() within a function, and I want to turn off the message it prints: x <- rnorm(100) y <- x^1.1+rnorm(100) nlsfit <- nls(y ~ g0*x^g1, start=list(g0=1,g1=1)) > confint(nlsfit) Waiting for profiling to be done... 2.5% 97.5% g0 0.4484198 1.143761 g1 1.0380479 2.370057 I cannot find any way to turn off 'Waiting for. .." I tried

Hello, I'm trying to fit experimental data with a model and nls. For some experiments, I have data with x from 0 to 1.2 and the fit is quite good. But it can happen that I have data only the [0,0.8] range (see the example below) and, then, the fit is not correct. I would like to add a constraint, for example : the second derivative must be positive. But I don't know how to add this to

Dear R experts, beginners and everyone else, I'm calculating "cumulative periodogram" using the command "cpgram" [1] from the MASS library. Here is a short example with the "lh" (hormone level) dataset: library(MASS) plot(lh,type="l",ylab="value",xlab="time", main="Hormone Levels (lh)") spectrum(lh,

On Fri, June 12, 2015 4:16 pm, jd1008 wrote: > > > On 06/12/2015 03:05 PM, Valeri Galtsev wrote: >> On Fri, June 12, 2015 3:54 pm, jd1008 wrote: >>> >>> On 06/12/2015 02:32 PM, Valeri Galtsev wrote: >>>> On Fri, June 12, 2015 3:20 pm, John R Pierce wrote: >>>>> On 6/12/2015 1:03 PM, Valeri Galtsev wrote: >>>>>> But the

Do you know this? Dario ------- Messaggio inoltrato ------- Da: stan <stanl-fedorauser at vfemail.net> Reply-to: Community support for Fedora users <users at lists.fedoraproject.org> A: users at lists.fedoraproject.org Oggetto: CIA Outlaw Country attack against CentOS / Rhel (and Fedora?) Is this credible? Data: Thu, 29 Jun 2017 15:51:43 -0700 Wikileaks released a document about an

Hi Christian, Chuck (and lists) It seems that the problem may be the strange behaviour of 'unstack' inside a function. See this thread in the R mailing list: http://tolstoy.newcastle.edu.au/R/help/04/03/1160.html Anyway, I got round the problem by using 'aggregate' instead of converting to a list and then tapply to sum values of metric. Probably more efficient as well. Thanks

An embedded and charset-unspecified text was scrubbed... Name: inte tillg?nglig URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20080512/9533b31f/attachment.pl>

Dear R users, It is a pleasure for me to announce the availability of the new package kwb.hantush (0.2.1)? on CRAN. Its objective is the calculation of groundwater mounding beneath an (stormwater) infiltration basin by solving the Hantush (1967) equation. For checking the correct implementation of the algorithm the R modelling results were cross-checked against alternative models assessed in

Dear R users, It is a pleasure for me to announce the availability of the new package kwb.hantush (0.2.1)? on CRAN. Its objective is the calculation of groundwater mounding beneath an (stormwater) infiltration basin by solving the Hantush (1967) equation. For checking the correct implementation of the algorithm the R modelling results were cross-checked against alternative models assessed in

I am trying to fit a rank-frequency distribution with 3 unknowns (a, b and k) to a set of data. This is my data set: y <- c(37047647,27083970,23944887,22536157,20133224, 20088720,18774883,18415648,17103717,13580739,12350767, 8682289,7496355,7248810,7022120,6396495,6262477,6005496, 5065887,4594147,2853307,2745322,454572,448397,275136,268771) and this is the fit I'm trying to do: nlsfit

I have an inquire about the RNG in R It is known that when we use the " rnorm " function , we pass the arguments : 1- number of variables to be generated 2- mean vector of the normal random errors. 3- standard deviation vector of the normal random errors. my question is the following Is the a way (a function) in R that we could specify the covariance matrix in step 3, instead of the

Dear list, I have a problem with the cumsum function. I have a data frame like the following one variable Year value EC01 2005 5 EC01 2006 10 AAO1 2005 2 AAO1 2006 4 what I would like to obtain is variable Year value cumsum EC01 2005 5 5 EC01

I have data vectors x and y both with 179 observations. I'm trying to fit a nonlinear model with five parameters using nls. The formula is only defined within a range of x-values, it should be zero otherwise, thus my attempted use of ifelse: > df<-data.frame(x,y) >

On Fri, June 12, 2015 3:54 pm, jd1008 wrote: > > > On 06/12/2015 02:32 PM, Valeri Galtsev wrote: >> On Fri, June 12, 2015 3:20 pm, John R Pierce wrote: >>> On 6/12/2015 1:03 PM, Valeri Galtsev wrote: >>>> But the bottom line is the same: in both cases you are executing >>>> somebody's else code on your computer. >>> >>> your