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

From: David Airey <david.airey_at_vanderbilt.edu>
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)".



R-help_at_r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html 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 https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.

list of date sections of archive