similar to: avoiding termination of nls given convergence failure

Hello. I am trying to fit a nonlinear mixed model involving 3 parameters. I have successfully made a self-starting function. getInitial() correctly outputs the initial estimates. I can also use the nlsList with this function to get the separate nonlinear fits by group. However, I get an error message when using the nlme function. Here is the relevent code:

Hello, I have come across a result that I cannot explain, and am hoping that someone else can provide an answer. A student fitted a mixed model using the lme function: out<- lme(fixed=Y~A+B+A:B, random=~1|Site). Y is a continuous variable while A and B are factors. The data set is balanced with the same number of observations in each combination of A and B. There are two hierarchical

I am getting an error message in a call to nlme and cannot understand what is happening. I explain the steps below in the hope that someone can explain the error and how to correct it. STEP 1: Data set: name: marouane.data. This is a data frame whose first few lines are as follows: > marouane.data[1:13,] species plant leaf irradiance photosynthesis chlorophyll 1

Hello. I am trying to understand the syntax of the nonlinear least squares function (nls) when the function definition is made outside of the call. Here is the context. 1. If I specify the following command, it works fine: > fit2<-nls( + A~Am*(1-exp(-alpha*(I-LCP))),data=dat1, + start=list(Am=3.6,alpha=0.01,LCP=20)) 2. Now, I want to be able to specify the function definition

Hello. I am using the nonlinear least squares function (nls). The function that I am trying to fit seems to be very sensitive to the starting values and, if these are not chosen properly, the nls function stops and gives an error message: Error in numericDeriv(form[[3]], names(ind), env) : Missing value or an Infinity produced when evaluating the model In addition: Warning

Hello. Following pages 342-347 of Pinheiro & Bates, I am trying to write a self-starting nonlinear function (a non-rectagular hyperbola) to be used in nonlinear least squares regression (and eventually for a mixed model). When I use the getInitial function for my self-starting function I get the following error message: > getInitial(photo~NRhyperbola(Irr,theta,Am,alpha,Rd),dat) Error

Hello, I am trying to do a nonlinear model using the "nls" command in R software. The data I am using is as follows: A<-c(7.132000,8.668667,9.880667,8.168000,10.863333,10.381333,11.059333,7.589333,4.716667,4.268667,7.265333,10.309333,8.456667,13.359333,8.624000,13.571333,12.523333,4.084667 ,NaN,NaN)

Hello. I know how to make a bar plot in which a numeric y variable is plotted against some grouping variable X (say, groups A, B, C) when this grouping variable is subdivided into each of two subgroups; so the bars would be: (group A subgroup 1) beside (group A subgroup 2), then (group B subgroup 1) beside (group B subgroup 2), and so on. This is done using the beside=TRUE argument in the

I'm wondering if I'm seeing a regression of 5616 with 3.4.8. I'm trying to set-up pptpd with winbind, which I'm doing on two machines on Debian lenny, and I'm trying on Debian Squeeze now. The Windows client is saying "Error 778: It was not possible to verify the identity of the server." The logs say that everything is ok, and that the client is hanging up the

Hello. I am working on an analysis involving the nonlinear mixed model function (nlme) in R. The data consist of measures of carbon fixation by leaves as a function of light intensity and the parametric function (standard in this area because it has a biological interpretation) is a non-rectangular hyperbola. I cannot get the nonlinear mixed model (nlme) function to converge cleanly. I am

Hi R users, I have a question about adding uncertainty bars to stacked bar plots. DF: year A B C Amin Amax Bmin Bmax Cmin Cmax 2009 40 45 15 30 61 23 56 14 17 2010 36 41 23 26 54 22 51 22 24 I use the code below: DF.refm = melt(subset(DF[,c(1:4)]),id.vars='year',variable_name='Legend') fig1 =

On 18/06/17 12:10, lily li wrote: > Hi R users, > > I have a question about adding uncertainty bars to stacked bar plots. > > DF: > year A B C Amin Amax Bmin Bmax Cmin Cmax > 2009 40 45 15 30 61 23 56 14 17 > 2010 36 41 23 26 54 22 51 22 24 > > I use the code below: > >

Hello everyone! I have a data frame of 136 variables with 270 observations. I would like to create a new data frame such that each element of that data frame contains the maximum value of the 6 prior lags of the initial data frame. So for example, if my original data frame, A, were A1=c(7.72, 7.94, 7.56, 7.54, 0.93, 0.59, 7.21, 8.00, 7.19, 7.57)A2=c(4.27, 3.70, 3.80, 3.67, 3.83, 3.95, 4.02, 2.06,

My goal is to make VPN access to our Samba PDC (FreeBSD 7.0) so that users can access there home shares from Windows clients. I have read the instructions at http://samba.org/ftp/unpacked/lorikeet/pppd/final-report.pdf, but I can't make it work. Don't know if is due to my lack of skills or has something to do with the Popop functionality in FreeBSD. Following the instructions, i made a

I am trying to set up a secure PPP tunnel between an OpenSSH client and server, and am having problems establishing the tunnel. ----------------------------------------------------------------------------- Server information: Stock Redhat 6.1 machine running a 2.2.12 kernel OpenSSH version 2.2.0p1 (downloaded as Redhat RPMs, revision 2) OpenSSL version 0.9.5a (downloaded as Redhat RPMs, revision

After upgrading to R 2.6.0, I'm having trouble running lmer: model <- lmer(primed ~ log(dist.time)*role + 1|target.utt, data=data.utts) Error in UseMethod("as.logical") : no applicable method for "as.logical" So I thought I'd upgrade lme4 to the latest version, but unfortunately the compilation fails - perhaps there's a missing #include: R CMD INSTALL

Hello. Douglas Bates has explained in a previous posting to R why he does not output residual degrees of freedom, F values and probabilities in the mixed model (lmer) function: because the usual degrees of freedom (obs - fixed df -1) are not exact and are really only upper bounds. I am interpreting what he said but I am not a professional statistician, so I might be getting this wrong... Does

Hi, I am having problems establishing an adsl connection in dom0 (but ONLY in dom0, without xen everything is ok). Following is my setup/situation: 1. Installed clean Fedora 7 (with virtualization option selected) 2. Booted with the non-xen kernel & setup adsl. Adsl connection works without any problems. 3. Rebooted with xen kernel into dom0. 4. adsl connection can no longer be activated

Dear R users, I want to do a fitted.contour plot of selected columns of a dataframe M with M$AM and M$Irradiance as x and y axes respectively. The level of the contour shall be determined by M$PR. Some words on my data first. Dataframe M looks like: head(M$Irradiance) [1] 293 350 412 419 477 509 head(M$AM) [1] 2.407 2.161 1.964 1.805 1.673 1.563 head(M$PR) [1] 70.102 72.600 75.097 80.167

Hi all, I read the setup at http://www.blackh0le.net/articles/vpn-dun-howto.html to setup my VPN. However, I'm having a problem which I think is proxy-ARP not working. I like to ask you to see if you know what's going on. When I ping 10.77.1.1 from windows XP machine the packets get to the 10.77.1.1 machine, but they don't have a return path to get back. When I do ping the windows