代码:
test.lme1 <- lme(logit.y1X1+X2+X3+X4+X5+X6+X7+X8+X10+X11+X12+X13+X14+X15+X16+X17,
random = 1|X1,data = mid1, method = 'ML') #We set this to maximum-likelihood as the default is REML (restricted maximum likelihood)
summary(test.lme1)
结果:
Linear mixed-effects model fit by maximum likelihood
Data: mid1
AIC BIC logLik
-592 -471 325
Random effects:
Formula: 1 | X1
(Intercept) Residual
StdDev: 0.115 0.0432
Fixed effects: logit.y1 ~ X1 + X2 + X3 + X4 + X5 + X6 + X7 + X8 + X10 + X11 + X12 + X13 + X14 + X15 + X16 + X17
Value Std.Error DF t-value p-value
(Intercept) -32.7 31.9 453 -1.02 0.3068
X1 0.0 0.0 453 -1.51 0.1322
X2 0.0 0.0 453 1.56 0.1188
X3aug 0.2 0.0 453 6.46 0.0000
X3dec -0.1 0.0 453 -4.04 0.0001
X3feb -0.3 0.0 453 -8.99 0.0000
X3jan -0.2 0.0 453 -8.57 0.0000
X3jul 0.2 0.0 453 6.81 0.0000
X3jun 0.1 0.0 453 4.84 0.0000
X3mat -0.1 0.0 453 -4.82 0.0000
X3may 0.2 0.0 453 5.45 0.0000
X3nov 0.0 0.0 453 1.12 0.2619
X3oct 0.2 0.0 453 5.30 0.0000
X3sep 0.1 0.0 453 3.25 0.0012
X4 0.2 0.1 453 2.33 0.0203
X5 0.1 0.1 453 1.36 0.1742
X6 0.0 0.3 453 -0.10 0.9197
X7 1.1 0.4 453 3.03 0.0026
X8 0.0 0.1 453 -0.41 0.6806
X10 0.0 0.0 453 1.54 0.1236
X11 1.7 1.8 453 0.97 0.3319
X12 29.9 31.8 453 0.94 0.3468
X13 -0.1 0.1 453 -1.70 0.0893
X14 149.3 164.7 453 0.91 0.3650
X15 -0.1 0.1 453 -1.90 0.0580
X16 0.0 0.0 453 0.59 0.5551
X17 -0.2 0.3 453 -0.82 0.4124
Correlation:
(Intr) X1 X2 X3aug X3dec X3feb X3jan X3jul X3jun X3mat X3may X3nov X3oct
X1 0.634
X2 -0.477 -0.751
X3aug -0.011 -0.017 0.012
X3dec 0.010 0.017 -0.012 0.500
X3feb 0.005 0.008 -0.006 0.500 0.500
X3jan 0.007 0.012 -0.009 0.500 0.500 0.500
X3jul -0.008 -0.012 0.009 0.500 0.500 0.500 0.500
X3jun -0.006 -0.008 0.006 0.500 0.500 0.500 0.500 0.500
X3mat 0.002 0.004 -0.003 0.500 0.500 0.500 0.500 0.500 0.500
X3may -0.003 -0.004 0.003 0.500 0.500 0.500 0.500 0.500 0.500 0.500
X3nov 0.013 0.021 -0.016 0.499 0.500 0.500 0.500 0.500 0.500 0.500 0.500
X3oct 0.015 0.025 -0.019 0.499 0.500 0.500 0.500 0.500 0.500 0.500 0.500 0.500
X3sep 0.018 0.029 -0.022 0.499 0.500 0.500 0.500 0.499 0.500 0.500 0.500 0.500 0.500
X4 0.413 0.477 -0.149 -0.008 0.008 0.004 0.006 -0.006 -0.004 0.002 -0.002 0.010 0.012
X5 0.439 0.372 -0.300 -0.006 0.006 0.003 0.005 -0.005 -0.003 0.002 -0.002 0.008 0.009
X6 0.194 0.388 -0.422 -0.006 0.006 0.003 0.005 -0.005 -0.003 0.002 -0.002 0.008 0.010
X7 0.399 0.430 -0.327 -0.007 0.007 0.004 0.005 -0.005 -0.004 0.002 -0.002 0.009 0.011
X8 0.214 -0.105 0.052 0.002 -0.002 -0.001 -0.001 0.001 0.001 0.000 0.000 -0.002 -0.003
X10 0.649 0.570 -0.553 -0.009 0.009 0.005 0.007 -0.007 -0.005 0.002 -0.002 0.012 0.014
X11 0.670 0.483 -0.517 -0.008 0.008 0.004 0.006 -0.006 -0.004 0.002 -0.002 0.010 0.012
X12 -1.000 -0.634 0.475 0.010 -0.010 -0.005 -0.008 0.008 0.005 -0.003 0.003 -0.013 -0.016
X13 -0.120 0.069 -0.115 -0.001 0.001 0.001 0.001 -0.001 -0.001 0.000 0.000 0.001 0.002
X14 -1.000 -0.634 0.477 0.010 -0.010 -0.005 -0.008 0.008 0.005 -0.003 0.003 -0.013 -0.016
X15 -0.661 -0.420 0.460 0.007 -0.007 -0.003 -0.005 0.005 0.003 -0.002 0.002 -0.009 -0.010
X16 -0.070 0.151 -0.100 -0.003 0.003 0.001 0.002 -0.002 -0.001 0.001 -0.001 0.003 0.004
X17 0.094 -0.191 0.253 0.003 -0.003 -0.002 -0.002 0.002 0.002 -0.001 0.001 -0.004 -0.005
X3sep X4 X5 X6 X7 X8 X10 X11 X12 X13 X14 X15 X16
X1
X2
X3aug
X3dec
X3feb
X3jan
X3jul
X3jun
X3mat
X3may
X3nov
X3oct
X3sep
X4 0.014
X5 0.011 0.701
X6 0.011 0.081 0.095
X7 0.012 0.412 0.393 0.212
X8 -0.003 0.184 0.711 -0.256 0.215
X10 0.016 0.295 0.465 0.115 0.477 0.337
X11 0.014 0.072 0.198 0.288 0.516 0.126 0.868
X12 -0.018 -0.416 -0.442 -0.192 -0.396 -0.217 -0.648 -0.666
X13 0.002 -0.391 -0.179 -0.101 -0.292 0.045 0.172 0.042 0.116
X14 -0.018 -0.413 -0.441 -0.194 -0.401 -0.217 -0.654 -0.675 1.000 0.117
X15 -0.012 -0.126 -0.214 -0.234 -0.561 -0.105 -0.801 -0.943 0.656 0.130 0.665
X16 0.004 -0.269 -0.355 -0.020 -0.188 -0.464 -0.173 -0.122 0.072 0.186 0.073 0.088
X17 -0.006 0.309 0.306 -0.052 0.174 0.441 0.162 0.104 -0.098 -0.161 -0.098 -0.083 -0.929
Standardized Within-Group Residuals:
Min Q1 Med Q3 Max
-1.023 -0.216 0.014 0.243 0.916
Number of Observations: 480
Number of Groups: 480
请问:(1)怎么描述其固定效应解释了多少变量,随机效应解释了多少变量
(2)这个模型总共解释了多少变量
(3)响应变量和解释变量的关系怎么描述
(4)要求以采样点为随机截距(高斯误差分布,同一性链接) 我这个代码高出这些没有
谢谢大佬了