similar to: Compare parameter estimates of a nlsList object

I have made many different attempts to run a R program in Batch mode. I run a DOS command like: C:\Progra~1\R\rw1061\bin\rterm.exe -q --no-restore --no-save myprogram.R out.txt What happens: 1. DOS Shell starts up well 2. R starts 3. I get the messages: "ARGUMENT 'myprogram.R' __ignored__" "ARGUMENT 'out.txt' __ignored__" 4. R Shell is ready to

Dear All, Using:- R : Copyright 2005, The R Foundation for Statistical Computing Version 2.1.1 (2005-06-20), ISBN 3-900051-07-0 and Jose Pinheiro, Douglas Bates, Saikat DebRoy and Deepayan Sarkar (2005). nlme: Linear and nonlinear mixed effects models. R package version 3.1-65. on a WINDOWS 2000 machine I am trying to run the scripts from the Mixed Models book and am running into some

Dear All, I have some code that works in S-Plus for writing saving a graphics file to disk :- graphsheet(type = "auto", format = "WMF", file = "G:\\north0l.wmf", pages = "auto", print.background = F, orientation="landscape", color.style="color") plot(x,y) dev.off() This works fine in S-Plus. I have tried playing

Dear listers, I am trying to fit a model using nlsList() using alternately a SSfol() selfstart function or its developped equivalent formulae. This preliminary trial works well mydata<-groupedData(Conc~Tps|Organ,data=mydata) mymod1<-nls(Conc~SSfol(Dose,Tps,lKe,lKa,lCl),data=mydata) as well as a developped form: mymod2<-nls(Conc~Dose * exp(lKe+lKa-lCl) *

Hello, nls-users, is it possible to pass additional parameters to the model function that are known and groupwise constant with nlsList? I could not find something like a "keep this fixed" option in the documentation and the code (my fault...?) The current workaround is to break the problem down into groups and use globals to pass the constant parameters, but it is ugly code and

Hi All I'm analysing growth rates using a gompertz (logistic) curve and am attempting to fit parameters for all of my study birds using nlsList. I've been looking for an option in nlsList to set min and max limits for parameter values during estimation, although I have failed to find any under 'control'. Other users seem to have reported that using control in the port algorithm

Hi All. I'm trying to run nlsList an getting an error that makes no sense to me. I have accuracy and reaction time data over many trials for each person (id) When I use nlsList code that is virtually identical to the example in the doc file I get the following error. I've tried everything I could think of and can't get around it. Any ideas what I'm doing wrong? **************

Dear R-helpers, I'm using R 2.5.0 under Windows and am trying to use nlsList from nlme 3.1-80 with the selfstart function for the four parametric logistic function. My first test went well, but now I'm trying to do some more sophisticated things and it does not work anymore. I simulate my data from a five parametric logistic function like this:

Dear R users, Does anybody have a function to calculate logLik or AIC for nlsList objects? After receiving error messages, another user helped me ascertain that this function is not currently written into R. Many thanks Lindsay

DeaR I am trying to use a dynamically create formula with nlsList and nlme, but I cannot get the environment of the string-generated formal to work similarly to the manually entered one. Any idea? Dieter #----- library(nlme) # Pinheiro/Bates p 280 fm1Indom.lis = nlsList(conc~SSbiexp(time,A1,lrc1,A2,lrc2), data=Indometh) nlme(fm1Indom.lis,random=pdDiag(A1+lrc1+A2~1)) # works... # Simulating

I hope that the mailing list is the correct forum for the question below. I have trouble calling functions nlsList and nlme from another function. Any help would be greatly appreciated. Jens Praestgaard Human Genome Sciences Rockville MD. I have a data set v with two components, v$mixeddat and v$init. They are listed below: > v$mixeddat conc result rep sample z 11 20.00000000

Dear all, is there a way to extract the individual variance-covariance matrix of the parameter estimates of the i-th individual using nlsList ? Thanks a lot Peter -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the

Hi All. I'd like to send some control arguments to the nls function when performing a nlsList analysis. I'm fitting a power model to some grouped data and would like to impose lower bounds on the estimates using the "port" algorithm. Obtaining the lower bound constraint works fine with a direct call to nls for a single level of the grouping variable. ?However, the bounds

Howdee, I am able to fit a 4-parameter logistic growth curve to a dataset which comprise many individuals (using R v. 2.3.1). Yet, if I want to obtain the parameters for each individual (i.e., for each 'id') using nlsList, then I obtain an Error message which I have trouble interpreting. Any advice as to how I can solve this problem? Thanks for your time, Marc > reg <-nls(mass ~

Hi list, In the object returned by summary.nlsList, what's the difference between "coefficients" and "parameters"? The have the same "Estimate", different se (therefore t value), but same p values. R.2.8.0 on winxp with nlme_3.1-89 Thanks, ...Tao +++++++++++++++++++++++++++++++++++++++++++++++++ > library(nlme) > fm1 <- nlsList(uptake ~

Hi I am using nlsList from the nlme package to estimate group specific coefficient sets to get a first impression of the structure of coefficients in relation to potential predictor variables. In some cases I got coefficients in other there are conversion problems due to bad start values I guess, resulting in typical error messages as known from nls: Error in nls(formula = formula, data = data,

Dear R users I am using GLMMpql to analyse some nested negative binomial response data. How do I summarise the significance of my main effects? For example, in a standard linear mixed model (lme), I would use anova.lme to obtain an F statistic and P value for each of my main effects: how do I achieve a similar goal using GLMMpql? many thanks Sarah Richardson Sarah Richardson Plant

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 <-

Hello all. I was wondering if the logLike.nls() and logLike.nlme() functions are still being used. Neither function seems to be available in the most recent release of R (1.9.1). The following is contained in the help file for logLik(): "classes which already have methods for this function include: 'glm', 'lm', 'nls' and 'gls', 'lme' and others in

Dear all, I an running R on a Windows 2000 machine (1.5Gb RAM) and am trying to load the lme4 package, however, whenever I attempt to load the library "lme4" I get the following error message I have installed lme4 and Matrix locally from the zip file with no problems. R version 2.6.2 (2008-02-08) Copyright (C) 2008 The R Foundation for Statistical Computing ISBN 3-900051-07-0