我正在使用
在下面的示例中,估计的方差为:
我希望你能把这些(例如)修正为:
lme4
R软件包,使用lmer()
函数创建线性混合模型。在此模型中,我有四个随机效应和一个固定效应(截距)。我的问题是关于随机效应的估计方差。是否可以像在SAS
中使用PARMS
参数一样指定协方差参数的初始值。在下面的示例中,估计的方差为:
c(0.00000, 0.03716, 0.00000, 0.02306)
我希望你能把这些(例如)修正为:
c(0.09902947, 0.02460464, 0.05848691, 0.06093686)
因此,它们没有被估计。
> summary(mod1)
Linear mixed model fit by maximum likelihood ['lmerMod']
Formula: log_cumcover_mod ~ (1 | kildestationsnavn) + (1 | year) + (1 |
kildestationsnavn:year) + (1 | proevetager)
Data: res
AIC BIC logLik deviance df.resid
109.9 122.9 -48.9 97.9 59
Scaled residuals:
Min 1Q Median 3Q Max
-2.1056 -0.6831 0.2094 0.8204 1.7574
Random effects:
Groups Name Variance Std.Dev.
kildestationsnavn:year (Intercept) 0.00000 0.0000
kildestationsnavn (Intercept) 0.03716 0.1928
proevetager (Intercept) 0.00000 0.0000
year (Intercept) 0.02306 0.1518
Residual 0.23975 0.4896
Number of obs: 65, groups:
kildestationsnavn:year, 6; kildestationsnavn, 3; proevetager, 2; year, 2
Fixed effects:
Estimate Std. Error t value
(Intercept) 4.9379 0.1672 29.54