Problem with multiple factors in nlme
"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 _._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._