similar to: Fitting data to a sigmoidal curve

I haven't used r in quite a while but would like to get back into it. I have a problem that I would like to solve with r. I have some multinomial data that looks to follow an asymmetric sigmoidal growth pattern. Solving a multinomial gompertz partial differential equation in r is what I’m after. Would anyone be able to provide me the code and packages to do something like this? Regards,

I am trying to understand my population abundance data and am looking into analyses of change point to try and determine, at approximately what point do populations begin to change (either decline or increasing). Can anyone offer suggestions on ways to go about this? I have looked into bcp and strucchange packages but am not completely convinced that these are appropriate for my data. Here is

I am trying to write a function that fits a sigmoid given a X and Y vector guessing the start parameters. I use nls. What I did (enclosed) seems to work well with many data points but if I want to fit small vectors like : pressure <- c(5,15,9,35,45) gas <- c(1000,2000,3000,4000,5000) it do not work. The help page says that it do no not work on zero residual data. Massimo Cressoni

Dear friends, I have asked last few days about cross-random effects using PQL, but I have not receive any answer because might my question was not clear. My question was about analysing the salamander mating data using PQL. This data contain cross-random effects for (male) and for (female). By opining MASS and lme library. I wrote this code sala.glmm <- glmmPQL(fixed=y~WSf*WSM,

Hi, I am rather new to R, but both myself and another much more experience user cannot figure this out. I have a collegue who has a 800+ nonlinear regressions to run for seed germination (different species, treatments, etc.) over time. I created a looping structure to extract the parameters from each regression; I will then use the parameters themselves for further analysis. I would like to

Hello, all. I have co-opted a number of functions that can be used to plot the hazard/survival functions and associated density distribution for a Gompertz mortality model, given known parameters. The Gompertz hazard model has been shown to fit relatively well to the human adult lifespan. For example, if I wanted to plot the hazard (i.e., mortality) functions: pop1 <- function (t) {

Hi, I have three years of species presence/absence data for a set of ~100 ponds, with a list of associated habitat characteristics (park, vegetation, hydroperiod, etc.). The datasets differ slightly by year because some ponds were dry in some years and not in others. My goal is to look at habitat associations for a couple focal species. For each year, I have a data table that looks like

Dear all! I tried to fit Gompertz growth models to describe cummulative germination rates using nls. I used the following code: germ.model<-nls(percent.germ~a*exp(-b*exp(-k*day)),data=tab,start=list(a=100,b=10,k=0.5)) My problem is that I want that the fitted model goes through the origin, since germination cannot start before the experiment was started, and y-max should be 100. Does anyone

Dear list: I am attempting to use what I thought would be a pretty straightforward practical application of Cox regression. I figure users of the survival package must have come across this problem before, so I would like to ask you how you dealt with it. I have set up an illustrative example and included it at the end of this post. I took a sample of 100 data points from each of two populations

Hello, I want to perform a nonlinear curve fit in order to obtain parameter estimates from experimentally determined data (y in dependence of x), but with an implicit function, thus, a function of which I cannot isolate y on the left-hand side of the equation. As far as I understand, the functions I found up to now (nls, optim) all work only for explicit functions. My data looks like

Hi everyone, I'm trying to perform the following regressions in order to compare linear vs. sigmoidal fit of the relationship between my dependent variable (y) and one explaining parameter (x2), both including the confounding effects of a third variable (x1): quasi-pois-lin <- glm(y ~ x1 + x2, family = quasipoisson(link="identity"), data=fit) quasi-pois-sig <- gnm(y ~ x1 +

I'm trying to plot an extrapolated logistic sigmoid curve using glm(..., family=binomial) as follows, but neither the fitted() points or the predict()ed curve are plotting correctly: > year <- c(2003+(6/12), 2004+(2/12), 2004+(10/12), 2005+(4/12)) > percent <- c(0.31, 0.43, 0.47, 0.50) > plot(year, percent, xlim=c(2003, 2007), ylim=c(0, 1)) > lm <- lm(percent ~ year)

Hi list, I was wondering how I should go about fitting a sigmoid curve to a dataset. More specifically how I estimate parameters a and b in the following equation: 1 / 1+exp(-(x-a)*b) with b the steepness of the sigmoid curve and a the shift of the center of the sigmoid curve relative to the center of your dataframe. The fit is in function of x, the location within the input vector and y, an

Hi all, Sorry if this has been answered already, but I couldn't find it in the archives or general internet. Is it possible to implement the gompertz distribution as survreg.distribution to use with survreg of the survival library? I haven't found anything and recent attempts from my side weren't succefull so far. I know that other packages like 'eha' and

Hello i did a stacked barplot using ggplot and R arranged the bars of the items in different orders. i don?t know why. but i want to have the same order in every stacked bar. I used the code data1 <- read.table("N_O_W_MAI.txt", header=TRUE, dec = ",") attach(data1) Teich1<-factor(Teich,levels=c(5,7,9,11,"G") ,ordered=is.ordered(Teich))

I'm trying to fit a Gompertz sigmoid as follows: x <- c(15, 16, 17, 18, 19) # arbitrary example data here; y <- c(0.1, 1.8, 2.2, 2.6, 2.9) # actual data is similar gm <- nls(y ~ a+b*exp(-exp(-c*(x-d))), start=c(a=?, b=?, c=?, d=?)) I have been unable to properly set the starting value '?'s. All of my guesses yield either a "singular gradient" error if they

List gurus, I'm trying to code a Gompertz growth curve function as part of a larger project and have run across a problem due to my ignorance of environments. Some sample data and the function are as follows: growth <- data.frame(age = c(1.92, 3, 5.83, 3.17, 15.5, 1.17, 5.58, 13.33, 14.29, 5.83, 13.79, 6.33, 13.75, 16.83, 13, 11.67, 0.25, 1.73, 9.46, 5.67), length = c(157, 165, 179,

Hi, is there any package (or source code snippet) that will evaluate the Frechet distance for curves represented as sets of points? Searching around only threw up references to a Frechet distribution. Thanks, ------------------------------------------------------------------- Rajarshi Guha <rxg218 at psu.edu> <http://jijo.cjb.net> GPG Fingerprint: 0CCA 8EE2 2EEB 25E2 AB04 06F7 1BB9

hi , i used the folllwing command to scaffold, G:\my\webblog>ruby script/generate scaffold webblog id:integer title:string body :text created_at:datetime after when i migrate with the follwing command rake db:migrate i got the error as (in G:/my/webblog) == 1 CreateWebblogs: migrating ================================================ -- create_table(:webblogs) rake aborted! Mysql::Error:

Hi I am struggling to fit a non-linear trend using the likfit function in geoR. Specifically I want a sigmoidal function, something like SSfpl in the nls package to fit the trend. But it seems trend.spatial in geoR only works with lm or glm type models. Any ideas how I can specify the model to calculate the kriging parameters using REML, including the parameters of a sigmoidal trend function