[R] confint() on mlm object requires named response matrix
iiori@@@wickeim@ier m@iii@g oii u@i-tuebi@ge@@de
iiori@@@wickeim@ier m@iii@g oii u@i-tuebi@ge@@de
Sat Feb 27 17:36:59 CET 2021
Dear list,
I fit a multivariate model and want to extract univariate confidence
intervals from the mlm object. If the response matrix is unnamed,
confint() just repeats the first interval. If the response matrix is
named, confint() correctly returns CIs for all columns:
set.seed(1143)
Y <- MASS::mvrnorm(10, c(0, 0), diag(2)) # unnamed response matrix
confint(lm(Y ~ 1))
2.5 % 97.5 %
:(Intercept) -0.7069968 0.2783868
:(Intercept) -0.7069968 0.2783868 # incorrect, same as first CI
colnames(Y) <- 1:2 # named response matrix
confint(lm(Y ~ 1))
2.5 % 97.5 %
1:(Intercept) -0.7069968 0.2783868
2:(Intercept) -0.9226622 -0.1208534 # correct
Is this intended behavior of confint()? I think it would be great if the
unnamed case gave a warning. Or am I missing something?
Best, Florian
---
Florian Wickelmaier
Department of Psychology
University of Tuebingen
More information about the R-help
mailing list