search for: getresponseformula

...he following work fine: > res.same <- lm(SFQ.LOCF.ENDPOINT ~ SFQ.BASELINE + GROUP,land) > res.diff <- lm(SFQ.LOCF.ENDPOINT ~ SFQ.BASELINE + GROUP + SFQ.BASELINE*GROUP,land) > anova(res.same,res.diff) I try: > lmList(SFQ.LOCF.ENDPOINT ~ SFQ.BASELINE | GROUP, land) Call: Error in getResponseFormula(form) : "Form" must be a two sided formula I'm puzzled. That looks like a two sided formula very like the one in the help for lme4 (which had been loaded) and the data look OK: > table(land$SFQ.LOCF.ENDPOINT) 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 23 1 1 2...

...with the library lme4. When I load both packages I get the following message and the getCovariateFormula function no longer works: library(nlme) library(lme4) Attaching package 'lme4': The following object(s) are masked from package:nlme : contr.SAS getCovariateFormula getResponseFormula groupedData lmeControl I have tried removing the package after using it, with: detach(package:lme4) library(nlme) but I still get an error message when I try to use getCovariateFormula. The line I use it in, and the error message is below: rownames(table)<-c((getCovariateFormula(model1)),(ge...

...rix package with apt-get. and both can be found in /usr/lib/R/site-library. When I tried an example for lmer, R could not find the function lmer(), > library(lme4) Attaching package: 'lme4' The following object(s) are masked from package:nlme : getCovariateFormula getResponseFormula groupedData Error in autoloader(name = "confint", package = "MASS") : autoloader did not find 'confint' in 'MASS' > (fm1 <- lmer(decrease ~ treatment + (1|rowpos) + (1|colpos), + OrchardSprays)) Error: couldn't find function &quot...

...{ if (any(naV <- is.na(match(allV, names(data))))) { stop(paste(allV[naV], "not found in data")) } } } else data <- NULL if (inherits(data, "groupedData")) { ff <- formula(data) rF <- deparse(getResponseFormula(ff)[[2]]) cF <- deparse(getCovariateFormula(ff)[[2]]) lbs <- attr(data, "labels") unts <- attr(data, "units") if (!is.null(lbs$x)) cL <- paste(lbs$x, unts$x) else cF <- NULL if (!is.null(lbs$y))...

I'm looking at Pinheiro & Bates "Mixed-Effects Models in S and S-PLUS" at the moment. Several datasets are used, one of which is called "PBIB" (a partially balanced incomplete block design). All the other datasets can be found somewhere or other in R. However, I cannot locate PBIB, and it does not seem to be mentioned in the latest edition of the R Full Reference

...me)) model.data <- groupedData(y~X|group, data=data.frame(X, y)) fit <- lme(y~-1+X, data=model.data, random=pdBlocked(list( pdBlocked(Z.block,pdClass="pdIdent"), pdIdent(~-1+ Replication) )) ,control=list(maxIter=1000, msMaxIter=1000, niterEM=1000)) It gives errror: "Error: getResponseFormula(el) : "Form" must be a two sided formula" Does anybody help how can I write random part? Thanks.. -- View this message in context: http://r.789695.n4.nabble.com/lme-random-effects-in-additive-models-with-interaction-tp4634067.html Sent from the R help mailing list archive at Nabb...