Displaying 2 results from an estimated 2 matches for "minrty".
2006 Nov 28
1
Slight discrepancy between predict.lm() and all.effects()
In the course of exploring response prediction, I stumbled upon a
small discrepancy between the CIs produced by predict.lm() and
all.effects()
require(mlmRev)
require(effects)
hsb.lm <- lm(mAch ~ minrty * sector, Hsb82)
hsb.new <- data.frame(
minrty = rep(c('No', 'Yes'), 2),
sector = rep(c('Public', 'Catholic'), each = 2))
hsb.eff <- all.effects(hsb.lm)
cbind(
hsb.new,
predict(hsb.lm, hsb.new, interval = 'confidence', type =
'...
2006 Nov 27
1
Help with response CIs for lme
Hi,
Can someone please offer a procedure for going from CIs produced by
intervals.lme() to fixed-effects response CIs.
Here's a simple example:
library(mlmRev)
library(nlme)
hsb.lme <- lme(mAch ~ minrty * sector, random = ~ 1 | cses, Hsb82)
(intervals(hsb.lme))
(hsb.new <- data.frame
minrty = rep(c('No', 'Yes'), 2),
sector = rep(c('Public', 'Catholic'), each = 2)))
cbind(hsb.new, predict(hsb.lme, hsb.new, level = 0))
Is the following correct (I know fr...