search for: lmod

Hello, is it possible to obtain type III sums of squares for a nested model as in the following: lmod <- lm(resp ~ A * B + (C %in% A), mydata)) I have tried library(car) Anova(lmod, type="III") but this gives me an error (and I also understand from the documentation of Anova as well as from a previous request (http://finzi.psych.upenn.edu/R/Rhelp02a/archive/64477.html) that it is no...

...hat attaches the 'coef' attribute to symbols. Since there is only one symbol object in a session with a given name, this means that this attaching has a global effect. Should this be quietly allowed or should there be a warning or an error? E.g., str(quote(Education)) # symbol Education lmod <- stats::lm(Fertility ~ ., data = datasets::swiss) glmod <- multcomp::glht(lmod, c("Agriculture=0", "Education=0")) str(quote(Education)) # symbol Education # - attr(*, "coef")= num 1 Bill Dunlap TIBCO Software wdunlap tibco.com [[alternative HTML version del...

melt.updn <- structure(list(date = structure(c(11808, 11869, 11961, 11992, 12084, 12173, 12265, 12418, 12600, 12631, 12753, 12996, 13057, 13149, 11808, 11869, 11961, 11992, 12084, 12173, 12265, 12418, 12600, 12631, 12753, 12996, 13057, 13149), class = "Date"), variable = structure(c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L,

...mbol object in a session with a >> given name, this means that >> this attaching has a global effect. Should this be quietly allowed >> or should there be a warning or >> an error? >> E.g., >> >> str(quote(Education)) >> # symbol Education >> lmod <- stats::lm(Fertility ~ ., data = datasets::swiss) >> glmod <- multcomp::glht(lmod, c("Agriculture=0", "Education=0")) >> str(quote(Education)) >> # symbol Education >> # - attr(*, "coef")= num 1 >> >> Bill Dunlap >> TIB...

...te to > symbols.? Since there is only one symbol object in a session with a given name, this means that > this attaching has a global effect.? Should this be quietly allowed or should there be a warning or > an error? > E.g., > > str(quote(Education)) > # symbol Education > lmod <- stats::lm(Fertility ~ ., data = datasets::swiss) > glmod <- multcomp::glht(lmod, c("Agriculture=0", "Education=0")) > str(quote(Education)) > # symbol Education > # - attr(*, "coef")= num 1 > > Bill Dunlap > TIBCO Software > wdunlap ti...

...cant or not. >From the help archives and the P/B book on mixed models i think the correct procedure is to build a model without the random effect and compare it to the model with the random effect using anova() Thus i created a model without the random effect using the code and compared lmm1 to lmod lmod <- lm(asinh(density_recruit) ~ region, data=datag) 3. my printout from anova() was Model df AIC BIC logLik Test L.Ratio p-value lmm1 1 5 1458.319 1477.708 -724.1596 lm1 2 4 1677.033 1692.544 - 834.5165 1 vs 2 220.7138 <.0001 Does this mean that lm...

-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 The lmer() function in the lme4 package has some code of the form mc <- match.call() mc[[1]] <- as.name("lFormula") lmod <- eval(mc, parent.frame(1L)) this is a fairly common idiom in R, found e.g. in lm(), used when one wants to pass all of the arguments of a function to a different function (in the case of lm() it's model.frame()). This idiom fails ("error in eval(...) : could not find function ...

...atches to previous arguments that I can see. Anybody know the reason or fix for this? calcMin(pvec=data.frame(val=par,min=rep(.000001,length(par)),max=rep(100 0000,length(par)),active=rep(TRUE,length(par))),func=optimwrap2, resid=resid,caa=caa,na=na,vcode=vcode,maa=maa,ny=ny,nind=nind,qage=qage, selmod=selmod,oldagei=oldagei,vpaflag=vpaflag) Error in f(x, ...) : unused argument(s) (resid = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0...