similar to: Support for compartment models / matrix exponential method

I try to fit subsequentially different nonlinear models to data. Up to now I tried a first model and if I could not fit this one proceeded to another one. Now I would like to handle this automatically. But much to my chagrin I did not find a way for dealing with the error messages from nls(), e.g. getting them as "control information" or at least just returning from nls() without a

Can anybody of the architects of R enlighten me about the reasoning behind the selection of the IEEE 754 mode for the building of R, please? Special points of interest are: 1. The selection seems to be based on the availability of isnan() and finite(). Is this a common way to test for IEEE 754? 2. The C9x draft introduces isfinite() instead of finite() (and as I see, some systems, e.g. the

I have a problem to understand the way gc_inhibit_torture is defined. In include/Defn.h it is declared as extern, using the __MAIN__ trick this leads to a definition in main/main.c. But in main/memory.c you define (and initialize) gc_inhibit_torture again: int gc_inhibit_torture = 1; /* gets set to zero after initialisations */ Then, in e.g. unix/system.c immediately at the beginning of

I tried to build R-1.1.1 on my HP-UX 10.20 box. For this I installed gcc 2.95.2 and built g77 (version 2.95.2 19991024 (release)). Making R runs to the point, where shared libs are linked. There it reports: make[4]: Entering directory `/disk4/hoffmann/Work/rwork/src/library/ctest/src' gcc -I../../../../include -I/usr/local/include -fPIC -g -O2 -c

After the overwhelming response to my last question (many thanks!), I cannot resist to ask the next maybe trivial question: All sources I have looked into tell me about building data frames by reading from files or by putting vectors of data ("column by column") together. Now I would like to create an (empty) data frame with "headings" for every column (column titles) and then

In Matlab or Octave I can evaluate (evalf) the content of a string. In R I have looked around in the docs, but have not found a way to do something equivalent (which is possible, I think). Can anybody give me a pointer? ============================================================================= Thomas Hoffmann, Institut fuer Halbleiter- und Mikrosystemtechnik, TU Dresden E-mail: hoffmann at

I want to report the progress of my attempts to get R running on HP-UX with gcc: Firstly, I could work around the PIC problem by leaving out all the lib stuff (I rebuilt gcc with --enable-share, but this (alone) did not help), e.g. the line gcc -shared -fPIC -o ctest.sl ansari.o fexact.o kendall.o ks.o prho.o swilk.o -L/usr/local/lib -lg2c -lm -L/lib/pa1.1 -L/usr/lib/pa1.1 -u main

I can''t get my view to display the most recent parameters from an associated object in a view... here is a birds-eye of my app: I have two models, Ferms and Kinetics. Ferm has_many :kinetics, Kinetic belongs_to :ferm. Kinetic has the fields ferm_id, brix and temp. In my ferms/index view I have a table listing the attributes of each ferm instance. I would like to display the most recent

I have dose response data analyzed with the package 'drc'. 'summary(mymodel)' prints my kinetic parameters. I want that text in an ASCII text file. I want to get exactly what I would get if I copied and pasted from the terminal window. I've read the documentation on data export to text files here:

?sink On May 18, 2018 9:47:25 AM PDT, Ed Siefker <ebs15242 at gmail.com> wrote: >I have dose response data analyzed with the package 'drc'. >'summary(mymodel)' prints my kinetic parameters. I want >that text in an ASCII text file. I want to get exactly what I >would get if I copied and pasted from the terminal window. > >I've read the documentation on

I was looking for all the glm-related 'family' functions in stats using the following predicate that returns TRUE for any function whose first argument is called "link". is.family <- function(object) is.function(object) && identical(names(as.list(object))[1], "link") It threw an error when applied to SSfol > is.family(SSfol) Error in

Dear list, I do have a problem with nls. I use the following data: >test time conc dose 0.50 5.40 1 0.75 11.10 1 1.00 8.40 1 1.25 13.80 1 1.50 15.50 1 1.75 18.00 1 2.00 17.00 1 2.50 13.90 1 3.00 11.20 1 3.50 9.90 1 4.00 4.70 1 5.00 5.00 1 6.00 1.90 1 7.00 1.90 1 9.00 1.10 1 12.00 0.95 1 14.00

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

Hi, I was trying to run the example of Indomethacin kinetics from the book: ## From Pinheiro/Bates, Mixed-Effects-Models in S and S-Plus, ## Springer, Second Printing 2001, Section 6.2 library(nlme) plot(Indometh) fm1Indom.nls <- nls(conc~SSbiexp(time,A1,lrc1,A2,lrc2), data=Indometh) summary(fm1Indom.nls) plot(fm1Indom.nls,Subject~resid(.),abline=0) ## .... the last plot command gives me the

Dear all, I am trying to gain some experience with the function gnls from the nlme package. I tried to model the Theophyline data by trying to model the presumed dependency of the clearance on the body weight. This is my function call of gnls: gnls(conc~SSfol(Dose,Time,lKe,lKa,lCl),data=Theoph, params=list(lKe~1,lKa~1,lCl~Wt),start=c(-2.4,0.46,-3.22,0.01)) That's been the result: Error

A couple of weeks ago there was a question regarding apparent convergence in nls when using the SSfol selfStart model for fitting a first-order pharmacokinetic model. I can't manage to find the original message either in my archive or in the list archives but the data were time conc dose 0.50 5.40 1 0.75 11.10 1 1.00 8.40 1 1.25 13.80 1 1.50 15.50 1

Hello, I'm trying to do some differential equation modeling in R using the package 'deSolve.' Briefly, I'm trying to use the law of mass action (the details of which aren't really important) to structure a vector of rate equation which will be passed into an ODE function and solved with associated functions from 'deSolve.' Roughly what this entails is scanning

Hi there LARTC, We are running a set of three systems for semiconductor technology, and would like to finalize our work in getting them to interoperate properly, but have run into some issues that touch on the very fabric of TCP/IP expertise. Iptables has already been used to solve "part of this", can we somehow use Diffserv capabilities to acheive the aims of modfiying packets such

I have never taken a statistics class nor read a statistics text, but I am in dire need of help with a trivial data analysis problem for which I need to write a report in two hours. I have spent 10,000 hours of study in my field of expertise (high frequency noise-making plant biology) but I've always thought that statistics is something that can be mastered on short notice. Briefly, I have

All, I have been able to successfully use the optim( ) function with "L-BFGS-B" to find reasonable parameters for a one-compartment open pharmacokinetic model. My loss function in this case was squared error, and I made no assumptions about the distribution of the plasma values. The model appeared to fit pretty well. Out of curiosity, I decided to try to use nlminb( ) applied to a