Hi there,
I can't see the SEs in $output... I guess that means they are computed by summary()?
But I also don't see how to extract them programatically from the output of summary(fit)
none of these work:
x = summary(fit)
x$SE; x$standard_error; x$error_estimate
any clues?
x below is the captured result of the summary function.
As an aside, in other R parameter summary output tables, the "parameter estimate" column is called "Estimate" and "error estimate" is called "Std.Error".
I think it would be useful to change to these conventions before too many scripts start using the current non-standard names for the parameters.
Column names changed.
answering own question:
summary(fit)->x
x$parameters->y
names(y)
y$"error estimate"
y$"parameter estimate"
CI95 = y$"error estimate" *1.96
out = data.frame(cbind(y$"parameter estimate" - CI95,y$"parameter estimate" + CI95))
names(out)<-c("lower95", "upper95")
cbind(y,out)
<
pre>
name matrix row col parameter estimate error estimate lower95 upper95
1 z.a 1 1 0.70073451 0.02057023 0.66041686 0.7410522
2 z.a 2 1 0.24912664 0.02591953 0.19832436 0.2999289
3 z.a 3 1 0.20698145 0.02510899 0.15776784 0.2561951
4 z.a 2 2 0.51334363 0.02391689 0.46646651 0.5602207
5 z.a 3 2 0.14779040 0.03128508 0.08647163 0.2091092
6 z.a 3 3 0.46045070 0.02237860 0.41658864 0.5043128
7 z.e 1 1 0.56046335 0.01486462 0.53132870 0.5895980
8 z.e 2 1 0.09410643 0.02322942 0.04857677 0.1396361
9 z.e 3 1 0.15132024 0.02320856 0.10583146 0.1968090
10 z.e 2 2 0.57513631 0.01714868 0.54152490 0.6087477
11 z.e 3 2 0.23929379 0.02285542 0.19449718 0.2840904
12 z.e 3 3 0.57438104 0.01486940 0.54523701 0.6035251
13 Trait1mean z.expMean 1 1 2.49599883 0.02137079 2.45411208 2.5378856
14 Trait2mean z.expMean 1 2 2.97599306 0.01945248 2.93786620 3.0141199
15 Trait3mean z.expMean 1 3 2.08150033 0.01837959 2.04547632 2.1175243