[R] Re: Problem with multiple factors in nlme

[Peter Dalgaard BSA]

"Uli Flenker; Raum 704" <uli at biochem.dshs-koeln.de> writes:

> > 	library(nlme)
> > 	data(CO2)
> > 	CO2.func <- deriv(~A*(1-exp(-exp(lB)*(conc-C))),
> > 			c("A","lB","C"),
> > 			function(conc,A,lB,C))
> 
> The correct function is
> 
> CO2.func <- deriv(~A*(1-exp(-exp(lB)*(conc-C))),
>                       c("A","lB","C"),
>                        function(conc,A,lB,C){})
> 
> <snip>
> 
> Sorry once more ...

And log(0,01) should read log(0.01) too...

Seems to boil down to 

Browse[1]> fixed[[nm]]
A ~ Treatment * Type
Browse[1]> fixed[[nm]][[3]]
Treatment * Type

..which is a formula, so fixed[[nm]][[3]] != "1" fails

changing to 

fixed[[nm]][[3]][[1]] != "1"

seems to work and it should be safe since any nonsimple formula will
have an operator and not a "1" in that position.

-- 
   O__  ---- Peter Dalgaard             Blegdamsvej 3  
  c/ /'_ --- Dept. of Biostatistics     2200 Cph. N   
 (*) \(*) -- University of Copenhagen   Denmark      Ph: (+45) 35327918
~~~~~~~~~~ - (p.dalgaard at biostat.ku.dk)             FAX: (+45) 35327907
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !)  To: r-help-request at stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._