# Ordinal Model================
ord<- clmm(Satisfaction ~ Self.Reported.Health+Education+Income+Wealth+Age+Private.Health.Insurance+Some.Diff.ADL+Not.Happy.Depress+Restless.Sleep+(1|Wave), data = HRSF)
summary(ord)## Cumulative Link Mixed Model fitted with the Laplace approximation
##
## formula:
## Satisfaction ~ Self.Reported.Health + Education + Income + Wealth +
## Age + Private.Health.Insurance + Some.Diff.ADL + Not.Happy.Depress +
## Restless.Sleep + (1 | Wave)
## data: HRSF
##
## link threshold nobs logLik AIC niter max.grad cond.H
## logit flexible 67328 -51066.02 102194.04 8755(26268) 3.94e-01 4.3e+03
##
## Random effects:
## Groups Name Variance Std.Dev.
## Wave (Intercept) 0.05555 0.2357
## Number of groups: Wave 10
##
## Coefficients:
## Estimate Std. Error z value Pr(>|z|)
## Self.Reported.Health2.Very good 0.32058 0.03904 8.211 < 2e-16 ***
## Self.Reported.Health3.Good 0.78129 0.03801 20.554 < 2e-16 ***
## Self.Reported.Health4.Fair 1.15726 0.03978 29.092 < 2e-16 ***
## Self.Reported.Health5.Poor 1.60299 0.04680 34.253 < 2e-16 ***
## Education1.GED 0.14718 0.03946 3.729 0.000192 ***
## Education2.HS -0.05350 0.02323 -2.303 0.021260 *
## Education3.HS/GED 0.05261 0.02708 1.943 0.052013 .
## Education4.AA/ Lt BA 0.08027 0.04553 1.763 0.077918 .
## Education5.BA -0.11743 0.03368 -3.486 0.000490 ***
## Education6.MA/MBA -0.22668 0.04240 -5.346 8.99e-08 ***
## Education7.Law/MD/PhD 0.14491 0.07217 2.008 0.044663 *
## Education8.Other -0.26521 0.68278 -0.388 0.697703
## Income1e+05+ -0.45278 0.03941 -11.490 < 2e-16 ***
## Income25000-49999 -0.27407 0.02085 -13.145 < 2e-16 ***
## Income50000-74999 -0.47296 0.02960 -15.977 < 2e-16 ***
## Income75000-99999 -0.51928 0.04158 -12.490 < 2e-16 ***
## Wealth100K-250K -0.40569 0.02747 -14.769 < 2e-16 ***
## Wealth1M+ -0.84777 0.04017 -21.102 < 2e-16 ***
## Wealth250K-500K -0.57932 0.03070 -18.870 < 2e-16 ***
## Wealth25K-100K -0.13727 0.02701 -5.083 3.72e-07 ***
## Wealth500K-1M -0.81745 0.03458 -23.638 < 2e-16 ***
## Age60-69 -0.43834 0.03341 -13.120 < 2e-16 ***
## Age70-79 -0.68434 0.03379 -20.255 < 2e-16 ***
## Age80+ -1.00485 0.03607 -27.858 < 2e-16 ***
## Private.Health.Insurance1 -0.29333 0.01780 -16.475 < 2e-16 ***
## Some.Diff.ADL1 0.38289 0.02158 17.745 < 2e-16 ***
## Not.Happy.Depress1 0.72661 0.01938 37.488 < 2e-16 ***
## Restless.Sleep1.Yes 0.30325 0.01856 16.336 < 2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Threshold coefficients:
## Estimate Std. Error z value
## 1.Very|2.Moderately -0.08927 0.09085 -0.983
## 2.Moderately|3.Not at all 2.70513 0.09149 29.569
## (353202 observations deleted due to missingness)frd<- plm(Satisfaction~ Some.Diff.ADL+Restless.Sleep+Not.Happy.Depress+Private.Health.Insurance+Excellent+Very.Good+Good+Fair+Poor-1, data= HRSFP, model='fd')
frd1<- plm(Satisfaction~ Some.Diff.ADL+Restless.Sleep+Not.Happy.Depress+Private.Health.Insurance+Self.Reported.Health.Dummy-1, data= HRSFP, model='fd')
summary(frd)## Oneway (individual) effect First-Difference Model
##
## Call:
## plm(formula = Satisfaction ~ Some.Diff.ADL + Restless.Sleep +
## Not.Happy.Depress + Private.Health.Insurance + Excellent +
## Very.Good + Good + Fair + Poor - 1, data = HRSFP, model = "fd")
##
## Unbalanced Panel: n = 20762, T = 1-10, N = 77004
## Observations used in estimation: 56242
##
## Residuals:
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## -2.2243 -0.0619 0.0000 0.0013 0.0569 2.2346
##
## Coefficients:
## Estimate Std. Error t-value Pr(>|t|)
## Some.Diff.ADL 0.0506472 0.0060928 8.3126 < 2.2e-16 ***
## Restless.Sleep 0.0234154 0.0049194 4.7598 1.943e-06 ***
## Not.Happy.Depress 0.0683778 0.0051018 13.4028 < 2.2e-16 ***
## Private.Health.Insurance -0.0114628 0.0053329 -2.1495 0.0316 *
## Excellent -0.1037045 0.0678974 -1.5274 0.1267
## Very.Good -0.0993915 0.0675161 -1.4721 0.1410
## Good -0.0618292 0.0674321 -0.9169 0.3592
## Fair -0.0184339 0.0674427 -0.2733 0.7846
## Poor 0.0434334 0.0677161 0.6414 0.5213
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Total Sum of Squares: 19862
## Residual Sum of Squares: 19637
## R-Squared: 0.0113
## Adj. R-Squared: 0.011159
## F-statistic: 71.3773 on 9 and 56233 DF, p-value: < 2.22e-16summary(frd1)## Oneway (individual) effect First-Difference Model
##
## Call:
## plm(formula = Satisfaction ~ Some.Diff.ADL + Restless.Sleep +
## Not.Happy.Depress + Private.Health.Insurance + Self.Reported.Health.Dummy -
## 1, data = HRSFP, model = "fd")
##
## Unbalanced Panel: n = 20757, T = 1-10, N = 76946
## Observations used in estimation: 56189
##
## Residuals:
## Min. 1st Qu. Median Mean 3rd Qu. Max.
## -2.22348 -0.06225 0.00000 0.00142 0.05659 2.20871
##
## Coefficients:
## Estimate Std. Error t-value Pr(>|t|)
## Some.Diff.ADL 0.0523271 0.0060830 8.6022 < 2.2e-16 ***
## Restless.Sleep 0.0237403 0.0049208 4.8245 1.407e-06 ***
## Not.Happy.Depress 0.0694456 0.0050980 13.6221 < 2.2e-16 ***
## Private.Health.Insurance -0.0108834 0.0053342 -2.0403 0.04133 *
## Self.Reported.Health.Dummy 0.0385075 0.0027231 14.1412 < 2.2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Total Sum of Squares: 19832
## Residual Sum of Squares: 19616
## R-Squared: 0.010866
## Adj. R-Squared: 0.010796
## F-statistic: 123.377 on 5 and 56184 DF, p-value: < 2.22e-16```