Re: [R] extension to nlme self start SSmicmen?

From: David Airey <>
Date: Tue, 22 Jan 2008 14:29:16 -0600

> Dear list,
> Has anyone created a version of SSmicmen that allows testing for
> group differences? The basic Michaelis-Menten equation is:
> (Bmax * X) / (Kd + X).
> The nlme package allows modeling of random effects for Bmax and Kd
> as needed, but I curious how I can build in group differences? I
> have receptor binding data for strains of mice, and following
> Pinheiro and Bates' lead in their book, a random effort for animal
> differences in Bmax seems reasonable, within strain. However, I'd
> like a self-start function, or even just the equation above,
> modified to allow testing for strain differences in Bmax and or Kd.
> Thank you.
> -Dave

Reading more closely, I now think I don't have to change the self start function SSmicmen to test for treatment group fixed effects. I just need to get the list of fixed effects right, to include the strain factor. Something like "fixed = list(bmax + kd ~ strain), random = list(bmax + kd ~ 1)". mailing list PLEASE do read the posting guide and provide commented, minimal, self-contained, reproducible code. Received on Tue 22 Jan 2008 - 20:31:04 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Tue 22 Jan 2008 - 21:30:08 GMT.

Mailing list information is available at Please read the posting guide before posting to the list.

list of date sections of archive