search for: factor_contrast

Displaying 10 results from an estimated 10 matches for "factor_contrast".

Did you mean: factor_contrasts
2009 Apr 21
3
broken example: lme() + multcomp() Tukey on repeated measures design
...+ V + Error(B/V), data = oats) > Lme.mod <- lme(Y ~ N + V, random = ~1 | B/V, data = oats) > summary(Aov.mod) > anova(Lme.mod) > summary(Lme.mod) > summary(glht(Lme.mod, linfct=mcp(V="Tukey"))) Error in eval(expr, envir, enclos) : object "Y" not found Error in factor_contrasts(model) : no 'model.matrix' method for 'model' found! It all went ok until the last line, where the actual Tukey multiple comparisons was requested. Thanks very much for any help! Bill
2008 Jan 10
1
general linear hypothesis glht() to work with lme()
...model with interaction dat.lme<-lme(info.index~type, random=~1|batch/type, data=dat) glht(dat.lme, linfct = mcp(type =c("b+t-2*m=0"))) The lme model fit is ok, but I got an error message with glht(): Error in eval(expr, envir, enclos) : object "info.index" not found Error in factor_contrasts(model) : no 'model.matrix' method for 'model' found! according to help page of glht(), it should work with linear mixed model, what is the problem here? Thanks ____________________________________________________________________________________ Be a better friend, newsho...
2007 Jul 12
1
error problem with glht
...vcov(fit_area_1) terms(fit_area_1) #perform multiple comparisons to test null hypothesis that #linear model fits are the same among basins fit1_mc <- glht(anova_area_1,linfct=mcp(basin="Tukey")) #the line above gives ... #Error in terms.default(object) : no terms component #Error in factor_contrasts(model) : no model.matrix method for model found! I appreciate any help offered. Cheers, Pete ------ Peter Furey Research Scientist Northwest Research Associates / CORA 3380 Mitchell Lane Boulder, CO 80301 Also affiliated with: Cooperative Institute for Research in Environmental Sciences (...
2008 Jul 25
1
glht after lmer with "$S4class-" and "missing model.matrix-" errors
Hello everybody. In my case, calculating multiple comparisons (Tukey) after lmer produced the following two errors: > sv.mc <- glht(model.sv,linfct=mcp(comp="Tukey")) Error in x$terms : $ operator not defined for this S4 class Error in factor_contrasts(model) : no 'model.matrix' method for 'model' found! What I have done before: > sv.growth <- groupedData(length~meas|box_id,outer=~comp,data=sv.growth) > model.sv <- lmer(length~comp+(meas|box_id),data=sv.growth) Warning message: In .local(x, ..., value) : Estimat...
2012 Feb 06
1
multiple comparisons in nested design
...trat5 -1.0623118 0.3947496 But it is useless for multiple comparisons: #> Pairs <- glht(AnovaModel.3, linfct = mcp(Tratamiento = "Tukey")) Error in model.matrix.default(model, data = structure(list(method = lm), .Names = "method"), : model structure is invalid Error in factor_contrasts(model) : no 'model.matrix' method for 'model' found! #> TukeyHSD(AnovaModel.3, "Tratamiento") Error in model.tables.aov(x, "means") : this fit does not inherit from "lm" Another way of analysis is through lme4 package: #> AnovaModel.3 &l...
2011 Jan 07
4
Problems with glht function for lme object
...on 'treat', using the multcomp package: summary(glht(m4.feed,linfct=mcp(treat="Tukey"))) However, this does not work, and I get the below error message. Error in if (is.null(pkg) | pkg == "nlme") terms(formula(x)) else slot(x, : argument is of length zero Error in factor_contrasts(model) : no ?model.matrix? method for ?model? found! I suspect this might have quite a straightforward solution, but I'm stuck at this point. Any help would be most appreciated. Sample data below. Kind regards, Andreas Nord Sweden ============== feedings sex site treat year date^...
2011 Feb 08
0
glht{multcomp} : use with lme {nlme}
...on=na.omit) ...runs sweet,..then ....try a multicomparisons approach for the categorical rf > summary(glht(CHOL, linfct=mcp(rf="Tukey"))) * Error in model.frame.default(object, data, xlev = xlev) : l'oggetto non è una matrice * <-- Object isn't a matrix!!* Errore in factor_contrasts(model) : no ‘model.matrix’ method for ‘model’ found! Errore in summary(glht(CHOL, linfct = mcp(rf = "Tukey"))) : error in evaluating the argument 'object' in selecting a method for function 'summary'* data are ok...I have followed the examples on ?glht ..,there is s...
2006 Oct 27
0
glht for aov with Error() term
...ames(amod), names(amod.err[["wool:tension"]])) [1] "contrasts" "xlevels" "call" "model" > > glht(amod[["wool:tension"]], linfct = mcp(tension = "Tukey")) Error in terms.default(object) : no terms component Error in factor_contrasts(model) : no 'model.matrix' method for 'model' found! --------------- end snip ---------------- ---------------------------------------- Michael Zehetleitner, M.Sc. Allgemeine und Experimentelle Psychologie I Ludwig-Maximilian Universität München [[alternative HTML versi...
2011 Apr 11
3
multiple comparisons with generalised least squares
...e factor levels by using the varIdent function. Then I wanted to perform multiple comparisons with the glht function of the multcomp package. glht(M1,linfct = mcp(Season = "Tukey")) However, here I got an error message "Error in terms.default(object) : no terms component. Error in factor_contrasts(model) : no ?model.matrix? method for ?model? found!". Does the glht function work with a gls model? And if not, is there an other way to perform multiple comparisons for a gls model? I've searched this forum for an answer to this question, but I could only found someone with the same q...
2008 Jul 25
0
glht after lmer with "$S4class-" and "missing model.matrix-" errors with DATA
...com>: >> Hello everybody. > In my case, calculating multiple comparisons (Tukey) after lmer > produced the following two errors: > > > sv.mc <- glht(model.sv,linfct=mcp(comp="Tukey")) > Error in x$terms : $ operator not defined for this S4 class > Error in factor_contrasts(model) : > no 'model.matrix' method for 'model' found! > > What I have done before: > > > sv.growth <- groupedData(length~meas|box_id,outer=~comp,data=sv.growth) > > model.sv <- lmer(length~comp+(meas|box_id),data=sv.growth) > Warning message: >...