Displaying 4 results from an estimated 4 matches for "estpar".
Did you mean:
estar
2006 Jul 11
2
Multiple tests on 2 way-ANOVA
...cients of a
# 2x2 factorial design (see Piantadosi 2005, p. 509)
cm1 <- matrix(c(0, 1, 0, 0, 0, 1), nrow = 2, ncol=3, byrow=TRUE,
dimnames = list(c("A", "B"), c("C.1", "C.2", "C.3")))
# cm1 = contrast matrix for main effects
v1 <- csimint(estpar=c(100, 6, 5), df=4*n-3, covm=cov1*sigma^2/n, cmatrix=cm1, conf.level=0.95)
summary(v1)
The adjusted p-values are almost the Bonferroni p-values.
If I understood right: You need not to adjust for multiple testing
on main effects in a 2x2 factorial design
assuming the absence of interaction.
I do...
2003 May 19
1
multcomp and glm
...rmine whether wolves are more likely to cross some
features than others and am using csimtest in the package multcomp to do
so.
x <- coef(model)
var.cov <- vcov(model)
df <- model$df.residual
contrast.matrix <- contrMat(rep(2,length(x), type = "Tukey")
post.hoc <- csimtest(estpar = x, df=df, covm = var.cov, cmatrix =
contrast.matrix)
summary(post.hoc)
My questions are:
(1) Have I entered my parameters for the posthoc analysis correctly and
does it matter which categorical variable I use as my reference category?
(2) How do I omit the continuous variable "null.cross&q...
2006 Jul 25
1
Multiple tests on repeated measurements
...andom= ~ 1|id)
summary(m1)
par4 <- m1$coef$fixed
cov4 <- vcov(m1)
cm4 <- matrix(c(0, 0, 1, 0, 0, 0, 1, 1), nrow = 2, ncol=4, byrow=TRUE,
dimnames = list(c("diff/v=0", "diff/v=1"), c("C.1", "C.2", "C.3", "C.4")))
v4 <- csimint(estpar=par4, df=n-6, # I'm not sure whether I found
# the correct degrees of freedom
covm=cov4,
cmatrix=cm4, conf.level=0.95)
sv4 <- summary(v4)
# 2. Method: I found in Handbook of Statistics Vol 13, p.616,
# same can be found in http://home.clara.net/sisa/bonhlp.htm
# Bonferroni on correl...
2005 Mar 09
1
multiple comparisons for lme using multcomp
...m=~1|box)
> library(multcomp)
Loading required package: mvtnorm
> # first way to do (seem uncorrect)
> summary(csimtest(coef(lm1),vcov(lm1),cmatrix=contrMat(table(treatment),
type="Tukey"),df=59))
Error in csimtest(coef(lm1), vcov(lm1), cmatrix =
contrMat(table(treatment), : estpar not a vector
> #indeed
> coef(lm1)
(Intercept) treatmentAl200 treatmentAl400 treatmentAl600 treatmentAl800
a 1.546679 -0.1648540 -0.4895219 -0.6383375 -0.7066632
b 1.546657 -0.1648540 -0.4895219 -0.6383375 -0.7066632
c 1.546664 -0.1648540...