##載入套件
library(haven)
library(dplyr)
library(lme4)
library(lmerTest)
library(knitr)##讀取 .sav 檔案
data <- read_sav("C:/Users/isisg/OneDrive/Desktop/社會變遷基本調查/2010Hofstade個人開車.sav")##定義summerize函數
# 定义 summarize_model 函数
summarize_model <- function(model, anova_result) {
full_model_summary <- summary(model)
# 提取并格式化固定效应结果
fixed_effects <- as.data.frame(full_model_summary$coefficients)
fixed_effects_rounded <- round(fixed_effects, 2)
# 提取并格式化随机效应结果
random_effects <- as.data.frame(VarCorr(model))
numeric_random_effects <- random_effects[sapply(random_effects, is.numeric)]
random_effects_rounded <- cbind(random_effects[, !sapply(random_effects, is.numeric)], round(numeric_random_effects, 2))
# 提取并格式化 ANOVA 结果
anova_result_rounded <- round(as.data.frame(anova_result), 2)
# 禁用科学记数法
options(scipen = 999)
# 自定义列名，避免乱码
colnames(fixed_effects_rounded) <- c("Estimate", "Std. Error", "df", "t value", "p-value")
# 打印固定效应结果并避免科学记数法
cat("Fixed Effects:\n")
print(kable(fixed_effects_rounded, format = "markdown", digits = 2))
# 打印随机效应结果并避免科学记数法
cat("\nRandom Effects:\n")
print(kable(random_effects_rounded, format = "markdown", digits = 2))
# 打印 ANOVA 结果并避免科学记数法
cat("\nANOVA Result:\n")
print(kable(anova_result_rounded, format = "markdown", digits = 2))
}##直接低成本行動 ##直接效果
full_model <- lmer(Zdirect_low ~ Z環境態度2 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 +
(1 + Z人均GDP | country), data = data)
reduced_model <- lmer(Zdirect_low ~ Z環境態度2 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 +
(1 | country), data = data)
anova_result <- anova(reduced_model, full_model)
summarize_model(full_model, anova_result)## Fixed Effects:
##
##
## | | Estimate| Std. Error| df| t value| p-value|
## |:-----------|--------:|----------:|--------:|-------:|-------:|
## |(Intercept) | 0.08| 0.05| 8.29| 1.81| 0.11|
## |Z環境態度2 | 0.24| 0.01| 27114.87| 44.88| 0.00|
## |Z年齡 | 0.16| 0.01| 27121.11| 27.86| 0.00|
## |Z性別 | 0.06| 0.01| 27102.99| 11.78| 0.00|
## |Z教育程度 | 0.03| 0.01| 27077.66| 5.14| 0.00|
## |Z社經地位 | 0.00| 0.01| 27127.25| -0.52| 0.60|
##
## Random Effects:
##
##
## |grp |var1 |var2 | vcov| sdcor|
## |:--------|:-----------|:--------|-----:|-----:|
## |country |(Intercept) |NA | 0.22| 0.46|
## |country |Z人均GDP |NA | 0.06| 0.24|
## |country |(Intercept) |Z人均GDP | -0.11| -1.00|
## |Residual |NA |NA | 0.69| 0.83|
##
## ANOVA Result:
##
##
## | | npar| AIC| BIC| logLik| deviance| Chisq| Df| Pr(>Chisq)|
## |:-------------|----:|--------:|--------:|---------:|--------:|-----:|--:|----------:|
## |reduced_model | 8| 67167.02| 67232.69| -33575.51| 67151.02| NA| NA| NA|
## |full_model | 10| 67162.37| 67244.46| -33571.18| 67142.37| 8.65| 2| 0.01|##政府信任
full_model <- lmer(Zdirect_low ~ Z環境態度2 + Z政府信任 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度2 * Z政府信任 + (1 + Z人均GDP | country), data = data)
reduced_model <- lmer(Zdirect_low ~ Z環境態度2 + Z政府信任 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度2 * Z政府信任 + (1 | country), data = data)
anova_result <- anova(reduced_model, full_model)
summarize_model(full_model, anova_result)## Fixed Effects:
##
##
## | | Estimate| Std. Error| df| t value| p-value|
## |:--------------------|--------:|----------:|--------:|-------:|-------:|
## |(Intercept) | 0.09| 0.05| 8.16| 1.93| 0.09|
## |Z環境態度2 | 0.24| 0.01| 27111.76| 44.80| 0.00|
## |Z政府信任 | -0.02| 0.01| 26960.72| -2.80| 0.01|
## |Z年齡 | 0.16| 0.01| 27118.24| 27.92| 0.00|
## |Z性別 | 0.06| 0.01| 27101.26| 11.83| 0.00|
## |Z教育程度 | 0.03| 0.01| 27079.76| 5.20| 0.00|
## |Z社經地位 | 0.00| 0.01| 27124.17| -0.25| 0.80|
## |Z環境態度2:Z政府信任 | 0.00| 0.00| 27104.23| -0.23| 0.81|
##
## Random Effects:
##
##
## |grp |var1 |var2 | vcov| sdcor|
## |:--------|:-----------|:--------|-----:|-----:|
## |country |(Intercept) |NA | 0.22| 0.47|
## |country |Z人均GDP |NA | 0.06| 0.24|
## |country |(Intercept) |Z人均GDP | -0.11| -1.00|
## |Residual |NA |NA | 0.69| 0.83|
##
## ANOVA Result:
##
##
## | | npar| AIC| BIC| logLik| deviance| Chisq| Df| Pr(>Chisq)|
## |:-------------|----:|--------:|--------:|---------:|--------:|-----:|--:|----------:|
## |reduced_model | 10| 67162.95| 67245.04| -33571.47| 67142.95| NA| NA| NA|
## |full_model | 12| 67158.38| 67256.88| -33567.19| 67134.38| 8.57| 2| 0.01|##權力距離
full_model <- lmer(Zdirect_low ~ Z環境態度2 + Z權力距離 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度2 * Z權力距離 + (1 + Z人均GDP + Z權力距離 | country), data = data)
reduced_model <- lmer(Zdirect_low ~ Z環境態度2 + Z權力距離 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度2 * Z權力距離 + (1 | country), data = data)
anova_result <- anova(reduced_model, full_model)
summarize_model(full_model, anova_result)## Fixed Effects:
##
##
## | | Estimate| Std. Error| df| t value| p-value|
## |:--------------------|--------:|----------:|--------:|-------:|-------:|
## |(Intercept) | 0.05| 0.06| 16.94| 0.94| 0.36|
## |Z環境態度2 | 0.24| 0.01| 26651.69| 44.21| 0.00|
## |Z權力距離 | -0.11| 0.06| 12.21| -1.80| 0.10|
## |Z年齡 | 0.16| 0.01| 26653.29| 27.35| 0.00|
## |Z性別 | 0.06| 0.01| 26638.07| 11.15| 0.00|
## |Z教育程度 | 0.03| 0.01| 26653.20| 4.92| 0.00|
## |Z社經地位 | 0.00| 0.01| 26653.08| -0.73| 0.46|
## |Z環境態度2:Z權力距離 | -0.03| 0.01| 26649.45| -5.09| 0.00|
##
## Random Effects:
##
##
## |grp |var1 |var2 | vcov| sdcor|
## |:--------|:-----------|:---------|-----:|-----:|
## |country |(Intercept) |NA | 0.07| 0.27|
## |country |Z人均GDP |NA | 0.21| 0.46|
## |country |Z權力距離 |NA | 0.14| 0.38|
## |country |(Intercept) |Z人均GDP | -0.07| -0.53|
## |country |(Intercept) |Z權力距離 | -0.01| -0.10|
## |country |Z人均GDP |Z權力距離 | 0.16| 0.90|
## |Residual |NA |NA | 0.69| 0.83|
##
## ANOVA Result:
##
##
## | | npar| AIC| BIC| logLik| deviance| Chisq| Df| Pr(>Chisq)|
## |:-------------|----:|--------:|--------:|---------:|--------:|-----:|--:|----------:|
## |reduced_model | 10| 66059.57| 66141.49| -33019.79| 66039.57| NA| NA| NA|
## |full_model | 15| 66054.79| 66177.66| -33012.39| 66024.79| 14.79| 5| 0.01|##個人主義
full_model <- lmer(Zdirect_low ~ Z環境態度2 + Z個人主義 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度2 * Z個人主義 + (1 + Z人均GDP + Z個人主義 | country), data = data)
reduced_model <- lmer(Zdirect_low ~ Z環境態度2 + Z個人主義 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度2 * Z個人主義 + (1 | country), data = data)
anova_result <- anova(reduced_model, full_model)
summarize_model(full_model, anova_result)## Fixed Effects:
##
##
## | | Estimate| Std. Error| df| t value| p-value|
## |:--------------------|--------:|----------:|--------:|-------:|-------:|
## |(Intercept) | -0.06| 0.06| 19.67| -0.96| 0.35|
## |Z環境態度2 | 0.24| 0.01| 27092.84| 44.22| 0.00|
## |Z個人主義 | 0.46| 0.09| 21.31| 5.04| 0.00|
## |Z年齡 | 0.16| 0.01| 26863.45| 27.88| 0.00|
## |Z性別 | 0.06| 0.01| 27103.35| 11.64| 0.00|
## |Z教育程度 | 0.03| 0.01| 26839.03| 5.17| 0.00|
## |Z社經地位 | 0.00| 0.01| 26820.05| -0.59| 0.56|
## |Z環境態度2:Z個人主義 | 0.01| 0.01| 27111.56| 2.10| 0.04|
##
## Random Effects:
##
##
## |grp |var1 |var2 | vcov| sdcor|
## |:--------|:-----------|:---------|-----:|-----:|
## |country |(Intercept) |NA | 0.06| 0.24|
## |country |Z人均GDP |NA | 0.00| 0.06|
## |country |Z個人主義 |NA | 0.22| 0.47|
## |country |(Intercept) |Z人均GDP | 0.00| -0.22|
## |country |(Intercept) |Z個人主義 | -0.09| -0.77|
## |country |Z人均GDP |Z個人主義 | -0.01| -0.46|
## |Residual |NA |NA | 0.69| 0.83|
##
## ANOVA Result:
##
##
## | | npar| AIC| BIC| logLik| deviance| Chisq| Df| Pr(>Chisq)|
## |:-------------|----:|--------:|--------:|---------:|--------:|-----:|--:|----------:|
## |reduced_model | 10| 67151.19| 67233.28| -33565.60| 67131.19| NA| NA| NA|
## |full_model | 15| 67152.34| 67275.47| -33561.17| 67122.34| 8.85| 5| 0.12|##成就動機
full_model <- lmer(Zdirect_low ~ Z環境態度2 + Z成就動機 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度2 * Z成就動機 + (1 + Z人均GDP + Z成就動機 | country), data = data)
reduced_model <- lmer(Zdirect_low ~ Z環境態度2 + Z成就動機 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度2 * Z成就動機 + (1 | country), data = data)
anova_result <- anova(reduced_model, full_model)
summarize_model(full_model, anova_result)## Fixed Effects:
##
##
## | | Estimate| Std. Error| df| t value| p-value|
## |:--------------------|--------:|----------:|--------:|-------:|-------:|
## |(Intercept) | 0.14| 0.06| 10.98| 2.23| 0.05|
## |Z環境態度2 | 0.24| 0.01| 26640.41| 44.50| 0.00|
## |Z成就動機 | 0.01| 0.04| 7.53| 0.31| 0.76|
## |Z年齡 | 0.16| 0.01| 26653.37| 27.52| 0.00|
## |Z性別 | 0.06| 0.01| 26638.91| 11.25| 0.00|
## |Z教育程度 | 0.03| 0.01| 26252.91| 4.96| 0.00|
## |Z社經地位 | 0.00| 0.01| 26659.73| -0.69| 0.49|
## |Z環境態度2:Z成就動機 | -0.02| 0.01| 26652.12| -4.21| 0.00|
##
## Random Effects:
##
##
## |grp |var1 |var2 | vcov| sdcor|
## |:--------|:-----------|:---------|-----:|-----:|
## |country |(Intercept) |NA | 0.21| 0.45|
## |country |Z人均GDP |NA | 0.08| 0.28|
## |country |Z成就動機 |NA | 0.01| 0.11|
## |country |(Intercept) |Z人均GDP | -0.12| -0.96|
## |country |(Intercept) |Z成就動機 | 0.05| 0.91|
## |country |Z人均GDP |Z成就動機 | -0.02| -0.76|
## |Residual |NA |NA | 0.69| 0.83|
##
## ANOVA Result:
##
##
## | | npar| AIC| BIC| logLik| deviance| Chisq| Df| Pr(>Chisq)|
## |:-------------|----:|--------:|--------:|---------:|--------:|-----:|--:|----------:|
## |reduced_model | 10| 66068.29| 66150.21| -33024.15| 66048.29| NA| NA| NA|
## |full_model | 15| 66066.11| 66188.98| -33018.05| 66036.11| 12.18| 5| 0.03|##不確定規避
full_model <- lmer(Zdirect_low ~ Z環境態度2 + Z不確定規避 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度2 * Z不確定規避 + (1 + Z人均GDP + Z不確定規避 | country), data = data)
reduced_model <- lmer(Zdirect_low ~ Z環境態度2 + Z不確定規避 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度2 * Z不確定規避 + (1 | country), data = data)
anova_result <- anova(reduced_model, full_model)
summarize_model(full_model, anova_result)## Fixed Effects:
##
##
## | | Estimate| Std. Error| df| t value| p-value|
## |:----------------------|--------:|----------:|--------:|-------:|-------:|
## |(Intercept) | -0.01| 0.08| 27.76| -0.09| 0.93|
## |Z環境態度2 | 0.24| 0.01| 26646.74| 44.50| 0.00|
## |Z不確定規避 | -0.06| 0.04| 18.11| -1.70| 0.11|
## |Z年齡 | 0.16| 0.01| 26651.73| 27.37| 0.00|
## |Z性別 | 0.06| 0.01| 26638.83| 11.28| 0.00|
## |Z教育程度 | 0.03| 0.01| 26510.99| 5.02| 0.00|
## |Z社經地位 | 0.00| 0.01| 26656.30| -0.71| 0.48|
## |Z環境態度2:Z不確定規避 | -0.02| 0.01| 26643.47| -3.14| 0.00|
##
## Random Effects:
##
##
## |grp |var1 |var2 | vcov| sdcor|
## |:--------|:-----------|:-----------|-----:|-----:|
## |country |(Intercept) |NA | 0.17| 0.42|
## |country |Z人均GDP |NA | 0.04| 0.19|
## |country |Z不確定規避 |NA | 0.02| 0.16|
## |country |(Intercept) |Z人均GDP | -0.05| -0.58|
## |country |(Intercept) |Z不確定規避 | 0.06| 0.88|
## |country |Z人均GDP |Z不確定規避 | 0.00| -0.11|
## |Residual |NA |NA | 0.69| 0.83|
##
## ANOVA Result:
##
##
## | | npar| AIC| BIC| logLik| deviance| Chisq| Df| Pr(>Chisq)|
## |:-------------|----:|--------:|--------:|---------:|--------:|-----:|--:|----------:|
## |reduced_model | 10| 66076.30| 66158.21| -33028.15| 66056.30| NA| NA| NA|
## |full_model | 15| 66073.22| 66196.10| -33021.61| 66043.22| 13.07| 5| 0.02|##長期導向
full_model <- lmer(Zdirect_low ~ Z環境態度2 + Z長期導向 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度2 * Z長期導向 + (1 + Z人均GDP + Z長期導向 | country), data = data)
reduced_model <- lmer(Zdirect_low ~ Z環境態度2 + Z長期導向 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度2 * Z長期導向 + (1 | country), data = data)
anova_result <- anova(reduced_model, full_model)
summarize_model(full_model, anova_result)## Fixed Effects:
##
##
## | | Estimate| Std. Error| df| t value| p-value|
## |:--------------------|--------:|----------:|--------:|-------:|-------:|
## |(Intercept) | 0.01| 0.06| 12.07| 0.10| 0.92|
## |Z環境態度2 | 0.24| 0.01| 27117.14| 44.88| 0.00|
## |Z長期導向 | 0.35| 0.07| 2.53| 4.77| 0.03|
## |Z年齡 | 0.16| 0.01| 27115.36| 27.84| 0.00|
## |Z性別 | 0.06| 0.01| 27101.22| 11.78| 0.00|
## |Z教育程度 | 0.03| 0.01| 27104.06| 5.11| 0.00|
## |Z社經地位 | 0.00| 0.01| 27124.40| -0.52| 0.61|
## |Z環境態度2:Z長期導向 | 0.00| 0.01| 27094.58| -0.16| 0.87|
##
## Random Effects:
##
##
## |grp |var1 |var2 | vcov| sdcor|
## |:--------|:-----------|:---------|-----:|-----:|
## |country |(Intercept) |NA | 0.15| 0.38|
## |country |Z人均GDP |NA | 0.04| 0.21|
## |country |Z長期導向 |NA | 0.04| 0.21|
## |country |(Intercept) |Z人均GDP | -0.07| -0.87|
## |country |(Intercept) |Z長期導向 | 0.07| 0.81|
## |country |Z人均GDP |Z長期導向 | -0.04| -0.99|
## |Residual |NA |NA | 0.69| 0.83|
##
## ANOVA Result:
##
##
## | | npar| AIC| BIC| logLik| deviance| Chisq| Df| Pr(>Chisq)|
## |:-------------|----:|--------:|--------:|---------:|--------:|-----:|--:|----------:|
## |reduced_model | 10| 67164.34| 67246.43| -33572.17| 67144.34| NA| NA| NA|
## |full_model | 15| 67162.27| 67285.40| -33566.13| 67132.27| 12.07| 5| 0.03|##放任約束
full_model <- lmer(Zdirect_low ~ Z環境態度2 + Z放任約束 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度2 * Z放任約束 + (1 + Z人均GDP + Z放任約束 | country), data = data)
reduced_model <- lmer(Zdirect_low ~ Z環境態度2 + Z放任約束 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度2 * Z放任約束 + (1 | country), data = data)
anova_result <- anova(reduced_model, full_model)
summarize_model(full_model, anova_result)## Fixed Effects:
##
##
## | | Estimate| Std. Error| df| t value| p-value|
## |:--------------------|--------:|----------:|--------:|-------:|-------:|
## |(Intercept) | 0.13| 0.06| 25.47| 2.34| 0.03|
## |Z環境態度2 | 0.24| 0.01| 26482.44| 44.31| 0.00|
## |Z放任約束 | 0.21| 0.07| 11.38| 3.19| 0.01|
## |Z年齡 | 0.16| 0.01| 26487.49| 27.52| 0.00|
## |Z性別 | 0.06| 0.01| 26471.53| 12.11| 0.00|
## |Z教育程度 | 0.03| 0.01| 26496.32| 4.60| 0.00|
## |Z社經地位 | 0.00| 0.01| 26493.63| 0.37| 0.71|
## |Z環境態度2:Z放任約束 | -0.02| 0.01| 26479.73| -3.39| 0.00|
##
## Random Effects:
##
##
## |grp |var1 |var2 | vcov| sdcor|
## |:--------|:-----------|:---------|-----:|-----:|
## |country |(Intercept) |NA | 0.00| 0.00|
## |country |Z人均GDP |NA | 0.30| 0.55|
## |country |Z放任約束 |NA | 0.15| 0.38|
## |country |(Intercept) |Z人均GDP | 0.00| NaN|
## |country |(Intercept) |Z放任約束 | 0.00| NaN|
## |country |Z人均GDP |Z放任約束 | -0.17| -0.80|
## |Residual |NA |NA | 0.68| 0.83|
##
## ANOVA Result:
##
##
## | | npar| AIC| BIC| logLik| deviance| Chisq| Df| Pr(>Chisq)|
## |:-------------|----:|--------:|--------:|---------:|--------:|-----:|--:|----------:|
## |reduced_model | 10| 65365.16| 65447.01| -32672.58| 65345.16| NA| NA| NA|
## |full_model | 15| 65367.61| 65490.38| -32668.80| 65337.61| 7.56| 5| 0.18|##直接高成本行動 ##直接效果
full_model <- lmer(Zdirect_high ~ Z環境態度2 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 +
(1 + Z人均GDP | country), data = data)
reduced_model <- lmer(Zdirect_high ~ Z環境態度2 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 +
(1 | country), data = data)
anova_result <- anova(reduced_model, full_model)
summarize_model(full_model, anova_result)## Fixed Effects:
##
##
## | | Estimate| Std. Error| df| t value| p-value|
## |:-----------|--------:|----------:|--------:|-------:|-------:|
## |(Intercept) | -0.01| 0.04| 32.14| -0.20| 0.84|
## |Z環境態度2 | 0.27| 0.01| 27127.80| 45.23| 0.00|
## |Z年齡 | 0.11| 0.01| 27131.39| 18.25| 0.00|
## |Z性別 | 0.07| 0.01| 27107.67| 11.74| 0.00|
## |Z教育程度 | 0.04| 0.01| 27091.60| 5.98| 0.00|
## |Z社經地位 | -0.02| 0.01| 27044.12| -2.86| 0.00|
##
## Random Effects:
##
##
## |grp |var1 |var2 | vcov| sdcor|
## |:--------|:-----------|:--------|-----:|-----:|
## |country |(Intercept) |NA | 0.05| 0.22|
## |country |Z人均GDP |NA | 0.01| 0.09|
## |country |(Intercept) |Z人均GDP | -0.01| -0.28|
## |Residual |NA |NA | 0.85| 0.92|
##
## ANOVA Result:
##
##
## | | npar| AIC| BIC| logLik| deviance| Chisq| Df| Pr(>Chisq)|
## |:-------------|----:|--------:|--------:|---------:|--------:|-----:|--:|----------:|
## |reduced_model | 8| 72651.42| 72717.09| -36317.71| 72635.42| NA| NA| NA|
## |full_model | 10| 72654.13| 72736.22| -36317.07| 72634.13| 1.29| 2| 0.53|##政府信任
full_model <- lmer(Zdirect_high ~ Z環境態度2 + Z政府信任 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度2 * Z政府信任 + (1 + Z人均GDP | country), data = data)
reduced_model <- lmer(Zdirect_high ~ Z環境態度2 + Z政府信任 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度2 * Z政府信任 + (1 | country), data = data)
anova_result <- anova(reduced_model, full_model)
summarize_model(full_model, anova_result)## Fixed Effects:
##
##
## | | Estimate| Std. Error| df| t value| p-value|
## |:--------------------|--------:|----------:|--------:|-------:|-------:|
## |(Intercept) | -0.01| 0.04| 32.13| -0.19| 0.85|
## |Z環境態度2 | 0.27| 0.01| 27125.36| 45.01| 0.00|
## |Z政府信任 | 0.00| 0.01| 27056.21| -0.36| 0.72|
## |Z年齡 | 0.11| 0.01| 27129.84| 18.25| 0.00|
## |Z性別 | 0.07| 0.01| 27105.77| 11.75| 0.00|
## |Z教育程度 | 0.04| 0.01| 27090.31| 6.00| 0.00|
## |Z社經地位 | -0.02| 0.01| 27056.07| -2.81| 0.00|
## |Z環境態度2:Z政府信任 | 0.00| 0.01| 27115.86| -0.45| 0.66|
##
## Random Effects:
##
##
## |grp |var1 |var2 | vcov| sdcor|
## |:--------|:-----------|:--------|-----:|-----:|
## |country |(Intercept) |NA | 0.05| 0.22|
## |country |Z人均GDP |NA | 0.01| 0.09|
## |country |(Intercept) |Z人均GDP | -0.01| -0.28|
## |Residual |NA |NA | 0.85| 0.92|
##
## ANOVA Result:
##
##
## | | npar| AIC| BIC| logLik| deviance| Chisq| Df| Pr(>Chisq)|
## |:-------------|----:|--------:|--------:|---------:|--------:|-----:|--:|----------:|
## |reduced_model | 10| 72655.08| 72737.17| -36317.54| 72635.08| NA| NA| NA|
## |full_model | 12| 72657.79| 72756.29| -36316.89| 72633.79| 1.3| 2| 0.52|##權力距離
full_model <- lmer(Zdirect_high ~ Z環境態度2 + Z權力距離 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度2 * Z權力距離 + (1 + Z人均GDP + Z權力距離 | country), data = data)
reduced_model <- lmer(Zdirect_high ~ Z環境態度2 + Z權力距離 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度2 * Z權力距離 + (1 | country), data = data)
anova_result <- anova(reduced_model, full_model)
summarize_model(full_model, anova_result)## Fixed Effects:
##
##
## | | Estimate| Std. Error| df| t value| p-value|
## |:--------------------|--------:|----------:|--------:|-------:|-------:|
## |(Intercept) | -0.01| 0.04| 19.42| -0.38| 0.71|
## |Z環境態度2 | 0.27| 0.01| 26659.91| 44.28| 0.00|
## |Z權力距離 | 0.05| 0.04| 12.12| 1.14| 0.28|
## |Z年齡 | 0.11| 0.01| 26636.48| 18.01| 0.00|
## |Z性別 | 0.06| 0.01| 26644.63| 11.01| 0.00|
## |Z教育程度 | 0.04| 0.01| 26487.00| 5.73| 0.00|
## |Z社經地位 | -0.02| 0.01| 26514.77| -2.92| 0.00|
## |Z環境態度2:Z權力距離 | -0.06| 0.01| 26509.31| -9.64| 0.00|
##
## Random Effects:
##
##
## |grp |var1 |var2 | vcov| sdcor|
## |:--------|:-----------|:---------|-----:|-----:|
## |country |(Intercept) |NA | 0.02| 0.16|
## |country |Z人均GDP |NA | 0.03| 0.18|
## |country |Z權力距離 |NA | 0.05| 0.22|
## |country |(Intercept) |Z人均GDP | -0.01| -0.35|
## |country |(Intercept) |Z權力距離 | 0.00| -0.04|
## |country |Z人均GDP |Z權力距離 | 0.04| 0.92|
## |Residual |NA |NA | 0.85| 0.92|
##
## ANOVA Result:
##
##
## | | npar| AIC| BIC| logLik| deviance| Chisq| Df| Pr(>Chisq)|
## |:-------------|----:|--------:|--------:|---------:|--------:|-----:|--:|----------:|
## |reduced_model | 10| 71417.93| 71499.84| -35698.96| 71397.93| NA| NA| NA|
## |full_model | 15| 71421.65| 71544.52| -35695.82| 71391.65| 6.28| 5| 0.28|##個人主義
full_model <- lmer(Zdirect_high ~ Z環境態度2 + Z個人主義 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度2 * Z個人主義 + (1 + Z人均GDP + Z個人主義 | country), data = data)
reduced_model <- lmer(Zdirect_high ~ Z環境態度2 + Z個人主義 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度2 * Z個人主義 + (1 | country), data = data)
anova_result <- anova(reduced_model, full_model)
summarize_model(full_model, anova_result)## Fixed Effects:
##
##
## | | Estimate| Std. Error| df| t value| p-value|
## |:--------------------|--------:|----------:|--------:|-------:|-------:|
## |(Intercept) | -0.08| 0.03| 22.05| -2.56| 0.02|
## |Z環境態度2 | 0.26| 0.01| 27090.75| 43.87| 0.00|
## |Z個人主義 | -0.03| 0.03| 4.32| -1.23| 0.28|
## |Z年齡 | 0.11| 0.01| 27104.98| 18.30| 0.00|
## |Z性別 | 0.06| 0.01| 27106.04| 11.22| 0.00|
## |Z教育程度 | 0.04| 0.01| 26692.39| 5.97| 0.00|
## |Z社經地位 | -0.02| 0.01| 27009.90| -2.85| 0.00|
## |Z環境態度2:Z個人主義 | 0.05| 0.01| 26921.86| 8.47| 0.00|
##
## Random Effects:
##
##
## |grp |var1 |var2 | vcov| sdcor|
## |:--------|:-----------|:---------|-----:|-----:|
## |country |(Intercept) |NA | 0.07| 0.26|
## |country |Z人均GDP |NA | 0.03| 0.18|
## |country |Z個人主義 |NA | 0.04| 0.19|
## |country |(Intercept) |Z人均GDP | 0.05| 1.00|
## |country |(Intercept) |Z個人主義 | -0.05| -1.00|
## |country |Z人均GDP |Z個人主義 | -0.03| -1.00|
## |Residual |NA |NA | 0.84| 0.92|
##
## ANOVA Result:
##
##
## | | npar| AIC| BIC| logLik| deviance| Chisq| Df| Pr(>Chisq)|
## |:-------------|----:|--------:|--------:|---------:|--------:|-----:|--:|----------:|
## |reduced_model | 10| 72582.52| 72664.60| -36281.26| 72562.52| NA| NA| NA|
## |full_model | 15| 72584.25| 72707.38| -36277.12| 72554.25| 8.27| 5| 0.14|##成就動機
full_model <- lmer(Zdirect_high ~ Z環境態度2 + Z成就動機 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度2 * Z成就動機 + (1 + Z人均GDP + Z成就動機 | country), data = data)
reduced_model <- lmer(Zdirect_high ~ Z環境態度2 + Z成就動機 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度2 * Z成就動機 + (1 | country), data = data)
anova_result <- anova(reduced_model, full_model)
summarize_model(full_model, anova_result)## Fixed Effects:
##
##
## | | Estimate| Std. Error| df| t value| p-value|
## |:--------------------|--------:|----------:|--------:|-------:|-------:|
## |(Intercept) | 0.02| 0.04| 26.02| 0.48| 0.63|
## |Z環境態度2 | 0.27| 0.01| 25480.54| 44.77| 0.00|
## |Z成就動機 | 0.10| 0.03| 22.02| 3.74| 0.00|
## |Z年齡 | 0.12| 0.01| 25654.68| 18.26| 0.00|
## |Z性別 | 0.06| 0.01| 26649.29| 11.29| 0.00|
## |Z教育程度 | 0.04| 0.01| 26258.57| 5.80| 0.00|
## |Z社經地位 | -0.02| 0.01| 25248.45| -2.83| 0.00|
## |Z環境態度2:Z成就動機 | -0.02| 0.01| 19974.21| -4.25| 0.00|
##
## Random Effects:
##
##
## |grp |var1 |var2 | vcov| sdcor|
## |:--------|:-----------|:---------|-----:|-----:|
## |country |(Intercept) |NA | 0.04| 0.21|
## |country |Z人均GDP |NA | 0.00| 0.03|
## |country |Z成就動機 |NA | 0.01| 0.11|
## |country |(Intercept) |Z人均GDP | -0.01| -1.00|
## |country |(Intercept) |Z成就動機 | 0.02| 1.00|
## |country |Z人均GDP |Z成就動機 | 0.00| -1.00|
## |Residual |NA |NA | 0.85| 0.92|
##
## ANOVA Result:
##
##
## | | npar| AIC| BIC| logLik| deviance| Chisq| Df| Pr(>Chisq)|
## |:-------------|----:|--------:|--------:|---------:|--------:|-----:|--:|----------:|
## |reduced_model | 10| 71488.36| 71570.27| -35734.18| 71468.36| NA| NA| NA|
## |full_model | 15| 71480.06| 71602.94| -35725.03| 71450.06| 18.29| 5| 0|##不確定規避
full_model <- lmer(Zdirect_high ~ Z環境態度2 + Z不確定規避 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度2 * Z不確定規避 + (1 + Z人均GDP + Z不確定規避 | country), data = data)
reduced_model <- lmer(Zdirect_high ~ Z環境態度2 + Z不確定規避 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度2 * Z不確定規避 + (1 | country), data = data)
anova_result <- anova(reduced_model, full_model)
summarize_model(full_model, anova_result)## Fixed Effects:
##
##
## | | Estimate| Std. Error| df| t value| p-value|
## |:----------------------|--------:|----------:|--------:|-------:|-------:|
## |(Intercept) | -0.02| 0.04| 26.00| -0.55| 0.59|
## |Z環境態度2 | 0.27| 0.01| 26646.45| 44.77| 0.00|
## |Z不確定規避 | 0.05| 0.03| 2.77| 1.58| 0.22|
## |Z年齡 | 0.11| 0.01| 26661.43| 18.03| 0.00|
## |Z性別 | 0.06| 0.01| 26644.42| 11.24| 0.00|
## |Z教育程度 | 0.04| 0.01| 26563.32| 5.91| 0.00|
## |Z社經地位 | -0.02| 0.01| 26591.51| -2.88| 0.00|
## |Z環境態度2:Z不確定規避 | -0.04| 0.01| 26602.88| -6.77| 0.00|
##
## Random Effects:
##
##
## |grp |var1 |var2 | vcov| sdcor|
## |:--------|:-----------|:-----------|----:|-----:|
## |country |(Intercept) |NA | 0.03| 0.18|
## |country |Z人均GDP |NA | 0.03| 0.16|
## |country |Z不確定規避 |NA | 0.02| 0.13|
## |country |(Intercept) |Z人均GDP | 0.00| 0.11|
## |country |(Intercept) |Z不確定規避 | 0.01| 0.52|
## |country |Z人均GDP |Z不確定規避 | 0.02| 0.91|
## |Residual |NA |NA | 0.85| 0.92|
##
## ANOVA Result:
##
##
## | | npar| AIC| BIC| logLik| deviance| Chisq| Df| Pr(>Chisq)|
## |:-------------|----:|--------:|--------:|---------:|--------:|-----:|--:|----------:|
## |reduced_model | 10| 71467.03| 71548.95| -35723.52| 71447.03| NA| NA| NA|
## |full_model | 15| 71469.28| 71592.15| -35719.64| 71439.28| 7.75| 5| 0.17|##長期導向
full_model <- lmer(Zdirect_high ~ Z環境態度2 + Z長期導向 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度2 * Z長期導向 + (1 + Z人均GDP + Z長期導向 | country), data = data)
reduced_model <- lmer(Zdirect_high ~ Z環境態度2 + Z長期導向 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度2 * Z長期導向 + (1 | country), data = data)
anova_result <- anova(reduced_model, full_model)
summarize_model(full_model, anova_result)## Fixed Effects:
##
##
## | | Estimate| Std. Error| df| t value| p-value|
## |:--------------------|--------:|----------:|--------:|-------:|-------:|
## |(Intercept) | -0.04| 0.04| 23.47| -1.04| 0.31|
## |Z環境態度2 | 0.27| 0.01| 26984.35| 45.17| 0.00|
## |Z長期導向 | 0.14| 0.11| 11.06| 1.30| 0.22|
## |Z年齡 | 0.11| 0.01| 27113.27| 18.28| 0.00|
## |Z性別 | 0.07| 0.01| 27101.69| 11.64| 0.00|
## |Z教育程度 | 0.04| 0.01| 27066.76| 5.97| 0.00|
## |Z社經地位 | -0.02| 0.01| 27038.36| -2.84| 0.00|
## |Z環境態度2:Z長期導向 | 0.02| 0.01| 27104.91| 3.38| 0.00|
##
## Random Effects:
##
##
## |grp |var1 |var2 | vcov| sdcor|
## |:--------|:-----------|:---------|-----:|-----:|
## |country |(Intercept) |NA | 0.00| 0.00|
## |country |Z人均GDP |NA | 0.04| 0.21|
## |country |Z長期導向 |NA | 0.27| 0.52|
## |country |(Intercept) |Z人均GDP | 0.00| NaN|
## |country |(Intercept) |Z長期導向 | 0.00| NaN|
## |country |Z人均GDP |Z長期導向 | -0.04| -0.37|
## |Residual |NA |NA | 0.85| 0.92|
##
## ANOVA Result:
##
##
## | | npar| AIC| BIC| logLik| deviance| Chisq| Df| Pr(>Chisq)|
## |:-------------|----:|--------:|--------:|---------:|--------:|-----:|--:|----------:|
## |reduced_model | 10| 72644.19| 72726.27| -36312.09| 72624.19| NA| NA| NA|
## |full_model | 15| 72670.86| 72793.99| -36320.43| 72640.86| 0| 5| 1|##放任約束
full_model <- lmer(Zdirect_high ~ Z環境態度2 + Z放任約束 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度2 * Z放任約束 + (1 + Z人均GDP + Z放任約束 | country), data = data)
reduced_model <- lmer(Zdirect_high ~ Z環境態度2 + Z放任約束 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度2 * Z放任約束 + (1 | country), data = data)
anova_result <- anova(reduced_model, full_model)
summarize_model(full_model, anova_result)## Fixed Effects:
##
##
## | | Estimate| Std. Error| df| t value| p-value|
## |:--------------------|--------:|----------:|--------:|-------:|-------:|
## |(Intercept) | -0.03| 0.04| 15.94| -0.75| 0.47|
## |Z環境態度2 | 0.27| 0.01| 26434.90| 44.97| 0.00|
## |Z放任約束 | 0.08| 0.03| 21.66| 3.03| 0.01|
## |Z年齡 | 0.12| 0.01| 26471.69| 18.54| 0.00|
## |Z性別 | 0.07| 0.01| 26478.47| 11.76| 0.00|
## |Z教育程度 | 0.04| 0.01| 25832.87| 6.05| 0.00|
## |Z社經地位 | -0.02| 0.01| 26187.94| -2.51| 0.01|
## |Z環境態度2:Z放任約束 | 0.01| 0.01| 26022.69| 1.07| 0.29|
##
## Random Effects:
##
##
## |grp |var1 |var2 | vcov| sdcor|
## |:--------|:-----------|:---------|-----:|-----:|
## |country |(Intercept) |NA | 0.01| 0.11|
## |country |Z人均GDP |NA | 0.11| 0.33|
## |country |Z放任約束 |NA | 0.04| 0.19|
## |country |(Intercept) |Z人均GDP | 0.01| 0.16|
## |country |(Intercept) |Z放任約束 | 0.00| 0.19|
## |country |Z人均GDP |Z放任約束 | -0.06| -0.94|
## |Residual |NA |NA | 0.85| 0.92|
##
## ANOVA Result:
##
##
## | | npar| AIC| BIC| logLik| deviance| Chisq| Df| Pr(>Chisq)|
## |:-------------|----:|--------:|--------:|---------:|--------:|-----:|--:|----------:|
## |reduced_model | 10| 70924.45| 71006.30| -35452.23| 70904.45| NA| NA| NA|
## |full_model | 15| 70921.51| 71044.28| -35445.75| 70891.51| 12.94| 5| 0.02|##間接低成本行動 ##直接效果
full_model <- lmer(Zindirect_low ~ Z環境態度2 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 +
(1 + Z人均GDP | country), data = data)
reduced_model <- lmer(Zindirect_low ~ Z環境態度2 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 +
(1 | country), data = data)
anova_result <- anova(reduced_model, full_model)
summarize_model(full_model, anova_result)## Fixed Effects:
##
##
## | | Estimate| Std. Error| df| t value| p-value|
## |:-----------|--------:|----------:|--------:|-------:|-------:|
## |(Intercept) | -0.15| 0.02| 11.84| -7.01| 0.00|
## |Z環境態度2 | 0.22| 0.01| 21093.32| 36.72| 0.00|
## |Z年齡 | 0.01| 0.01| 25482.69| 1.83| 0.07|
## |Z性別 | 0.00| 0.01| 27063.28| -0.71| 0.48|
## |Z教育程度 | 0.13| 0.01| 16675.93| 19.51| 0.00|
## |Z社經地位 | 0.02| 0.01| 25238.71| 2.78| 0.01|
##
## Random Effects:
##
##
## |grp |var1 |var2 | vcov| sdcor|
## |:--------|:-----------|:--------|----:|-----:|
## |country |(Intercept) |NA | 0.03| 0.16|
## |country |Z人均GDP |NA | 0.01| 0.12|
## |country |(Intercept) |Z人均GDP | 0.02| 0.88|
## |Residual |NA |NA | 0.84| 0.92|
##
## ANOVA Result:
##
##
## | | npar| AIC| BIC| logLik| deviance| Chisq| Df| Pr(>Chisq)|
## |:-------------|----:|--------:|--------:|---------:|--------:|-----:|--:|----------:|
## |reduced_model | 8| 72389.75| 72455.42| -36186.88| 72373.75| NA| NA| NA|
## |full_model | 10| 72382.78| 72464.87| -36181.39| 72362.78| 10.97| 2| 0|##政府信任
full_model <- lmer(Zindirect_low ~ Z環境態度2 + Z政府信任 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度2 * Z政府信任 + (1 + Z人均GDP | country), data = data)
reduced_model <- lmer(Zindirect_low ~ Z環境態度2 + Z政府信任 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度2 * Z政府信任 + (1 | country), data = data)
anova_result <- anova(reduced_model, full_model)
summarize_model(full_model, anova_result)## Fixed Effects:
##
##
## | | Estimate| Std. Error| df| t value| p-value|
## |:--------------------|--------:|----------:|--------:|-------:|-------:|
## |(Intercept) | -0.15| 0.02| 12.02| -6.83| 0.00|
## |Z環境態度2 | 0.22| 0.01| 19326.89| 36.74| 0.00|
## |Z政府信任 | 0.01| 0.01| 22485.37| 1.09| 0.28|
## |Z年齡 | 0.01| 0.01| 24660.25| 1.79| 0.07|
## |Z性別 | 0.00| 0.01| 27042.29| -0.76| 0.45|
## |Z教育程度 | 0.13| 0.01| 14791.22| 19.43| 0.00|
## |Z社經地位 | 0.02| 0.01| 25386.61| 2.67| 0.01|
## |Z環境態度2:Z政府信任 | 0.02| 0.01| 26442.47| 2.84| 0.00|
##
## Random Effects:
##
##
## |grp |var1 |var2 | vcov| sdcor|
## |:--------|:-----------|:--------|----:|-----:|
## |country |(Intercept) |NA | 0.02| 0.16|
## |country |Z人均GDP |NA | 0.01| 0.12|
## |country |(Intercept) |Z人均GDP | 0.02| 0.86|
## |Residual |NA |NA | 0.84| 0.92|
##
## ANOVA Result:
##
##
## | | npar| AIC| BIC| logLik| deviance| Chisq| Df| Pr(>Chisq)|
## |:-------------|----:|--------:|--------:|---------:|--------:|-----:|--:|----------:|
## |reduced_model | 10| 72383.84| 72465.92| -36181.92| 72363.84| NA| NA| NA|
## |full_model | 12| 72377.22| 72475.73| -36176.61| 72353.22| 10.61| 2| 0|##權力距離
full_model <- lmer(Zindirect_low ~ Z環境態度2 + Z權力距離 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度2 * Z權力距離 + (1 + Z人均GDP + Z權力距離 | country), data = data)
reduced_model <- lmer(Zindirect_low ~ Z環境態度2 + Z權力距離 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度2 * Z權力距離 + (1 | country), data = data)
anova_result <- anova(reduced_model, full_model)
summarize_model(full_model, anova_result)## Fixed Effects:
##
##
## | | Estimate| Std. Error| df| t value| p-value|
## |:--------------------|--------:|----------:|--------:|-------:|-------:|
## |(Intercept) | -0.06| 0.02| 18.09| -3.45| 0.00|
## |Z環境態度2 | 0.22| 0.01| 25790.59| 36.13| 0.00|
## |Z權力距離 | -0.12| 0.02| 10.54| -5.20| 0.00|
## |Z年齡 | 0.01| 0.01| 26568.12| 1.52| 0.13|
## |Z性別 | -0.01| 0.01| 26641.80| -1.21| 0.23|
## |Z教育程度 | 0.13| 0.01| 18955.72| 19.22| 0.00|
## |Z社經地位 | 0.02| 0.01| 21977.33| 2.58| 0.01|
## |Z環境態度2:Z權力距離 | -0.07| 0.01| 24222.98| -12.40| 0.00|
##
## Random Effects:
##
##
## |grp |var1 |var2 | vcov| sdcor|
## |:--------|:-----------|:---------|----:|-----:|
## |country |(Intercept) |NA | 0.00| 0.00|
## |country |Z人均GDP |NA | 0.08| 0.27|
## |country |Z權力距離 |NA | 0.03| 0.16|
## |country |(Intercept) |Z人均GDP | 0.00| NaN|
## |country |(Intercept) |Z權力距離 | 0.00| NaN|
## |country |Z人均GDP |Z權力距離 | 0.04| 0.92|
## |Residual |NA |NA | 0.83| 0.91|
##
## ANOVA Result:
##
##
## | | npar| AIC| BIC| logLik| deviance| Chisq| Df| Pr(>Chisq)|
## |:-------------|----:|--------:|--------:|---------:|--------:|-----:|--:|----------:|
## |reduced_model | 10| 70994.16| 71076.08| -35487.08| 70974.16| NA| NA| NA|
## |full_model | 15| 71002.00| 71124.87| -35486.00| 70972.00| 2.16| 5| 0.83|##個人主義
full_model <- lmer(Zindirect_low ~ Z環境態度2 + Z個人主義 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度2 * Z個人主義 + (1 + Z人均GDP + Z個人主義 | country), data = data)
reduced_model <- lmer(Zindirect_low ~ Z環境態度2 + Z個人主義 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度2 * Z個人主義 + (1 | country), data = data)
anova_result <- anova(reduced_model, full_model)
summarize_model(full_model, anova_result)## Fixed Effects:
##
##
## | | Estimate| Std. Error| df| t value| p-value|
## |:--------------------|--------:|----------:|--------:|-------:|-------:|
## |(Intercept) | -0.05| 0.01| 38.55| -3.47| 0.00|
## |Z環境態度2 | 0.21| 0.01| 26085.73| 34.54| 0.00|
## |Z個人主義 | 0.12| 0.01| 25.80| 10.37| 0.00|
## |Z年齡 | 0.01| 0.01| 23795.34| 1.85| 0.06|
## |Z性別 | -0.01| 0.01| 27101.77| -1.61| 0.11|
## |Z教育程度 | 0.13| 0.01| 19631.55| 19.53| 0.00|
## |Z社經地位 | 0.02| 0.01| 19692.32| 2.69| 0.01|
## |Z環境態度2:Z個人主義 | 0.09| 0.01| 8810.54| 16.01| 0.00|
##
## Random Effects:
##
##
## |grp |var1 |var2 | vcov| sdcor|
## |:--------|:-----------|:---------|-----:|-----:|
## |country |(Intercept) |NA | 0.00| 0.05|
## |country |Z人均GDP |NA | 0.08| 0.28|
## |country |Z個人主義 |NA | 0.03| 0.17|
## |country |(Intercept) |Z人均GDP | 0.01| 1.00|
## |country |(Intercept) |Z個人主義 | -0.01| -1.00|
## |country |Z人均GDP |Z個人主義 | -0.05| -1.00|
## |Residual |NA |NA | 0.83| 0.91|
##
## ANOVA Result:
##
##
## | | npar| AIC| BIC| logLik| deviance| Chisq| Df| Pr(>Chisq)|
## |:-------------|----:|--------:|--------:|---------:|--------:|-----:|--:|----------:|
## |reduced_model | 10| 72135.24| 72217.32| -36057.62| 72115.24| NA| NA| NA|
## |full_model | 15| 72130.81| 72253.94| -36050.41| 72100.81| 14.43| 5| 0.01|##成就動機
full_model <- lmer(Zindirect_low ~ Z環境態度2 + Z成就動機 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度2 * Z成就動機 + (1 + Z人均GDP + Z成就動機 | country), data = data)
reduced_model <- lmer(Zindirect_low ~ Z環境態度2 + Z成就動機 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度2 * Z成就動機 + (1 | country), data = data)
anova_result <- anova(reduced_model, full_model)
summarize_model(full_model, anova_result)## Fixed Effects:
##
##
## | | Estimate| Std. Error| df| t value| p-value|
## |:--------------------|--------:|----------:|--------:|-------:|-------:|
## |(Intercept) | -0.04| 0.02| 8.63| -2.67| 0.03|
## |Z環境態度2 | 0.22| 0.01| 24523.53| 36.68| 0.00|
## |Z成就動機 | 0.15| 0.03| 1.88| 5.98| 0.03|
## |Z年齡 | 0.01| 0.01| 26402.94| 1.89| 0.06|
## |Z性別 | 0.00| 0.01| 26582.89| -0.83| 0.41|
## |Z教育程度 | 0.13| 0.01| 13164.59| 19.30| 0.00|
## |Z社經地位 | 0.02| 0.01| 11822.75| 2.61| 0.01|
## |Z環境態度2:Z成就動機 | -0.02| 0.01| 23080.29| -3.98| 0.00|
##
## Random Effects:
##
##
## |grp |var1 |var2 | vcov| sdcor|
## |:--------|:-----------|:---------|-----:|-----:|
## |country |(Intercept) |NA | 0.00| 0.00|
## |country |Z人均GDP |NA | 0.05| 0.22|
## |country |Z成就動機 |NA | 0.02| 0.14|
## |country |(Intercept) |Z人均GDP | 0.00| NaN|
## |country |(Intercept) |Z成就動機 | 0.00| NaN|
## |country |Z人均GDP |Z成就動機 | -0.03| -0.92|
## |Residual |NA |NA | 0.84| 0.92|
##
## ANOVA Result:
##
##
## | | npar| AIC| BIC| logLik| deviance| Chisq| Df| Pr(>Chisq)|
## |:-------------|----:|--------:|--------:|---------:|--------:|-----:|--:|----------:|
## |reduced_model | 10| 71138.31| 71220.23| -35559.16| 71118.31| NA| NA| NA|
## |full_model | 15| 71139.92| 71262.79| -35554.96| 71109.92| 8.39| 5| 0.14|##不確定規避
full_model <- lmer(Zindirect_low ~ Z環境態度2 + Z不確定規避 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度2 * Z不確定規避 + (1 + Z人均GDP + Z不確定規避 | country), data = data)
reduced_model <- lmer(Zindirect_low ~ Z環境態度2 + Z不確定規避 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度2 * Z不確定規避 + (1 | country), data = data)
anova_result <- anova(reduced_model, full_model)
summarize_model(full_model, anova_result)## Fixed Effects:
##
##
## | | Estimate| Std. Error| df| t value| p-value|
## |:----------------------|--------:|----------:|--------:|-------:|-------:|
## |(Intercept) | -0.03| 0.02| 24.09| -1.08| 0.29|
## |Z環境態度2 | 0.22| 0.01| 26208.87| 36.68| 0.00|
## |Z不確定規避 | -0.10| 0.03| 13.08| -3.99| 0.00|
## |Z年齡 | 0.01| 0.01| 26457.73| 1.69| 0.09|
## |Z性別 | -0.01| 0.01| 26648.36| -0.89| 0.37|
## |Z教育程度 | 0.13| 0.01| 23999.19| 19.41| 0.00|
## |Z社經地位 | 0.02| 0.01| 24396.95| 2.54| 0.01|
## |Z環境態度2:Z不確定規避 | -0.04| 0.01| 26540.69| -7.03| 0.00|
##
## Random Effects:
##
##
## |grp |var1 |var2 | vcov| sdcor|
## |:--------|:-----------|:-----------|----:|-----:|
## |country |(Intercept) |NA | 0.00| 0.00|
## |country |Z人均GDP |NA | 0.04| 0.20|
## |country |Z不確定規避 |NA | 0.01| 0.08|
## |country |(Intercept) |Z人均GDP | 0.00| NaN|
## |country |(Intercept) |Z不確定規避 | 0.00| NaN|
## |country |Z人均GDP |Z不確定規避 | 0.01| 0.68|
## |Residual |NA |NA | 0.84| 0.92|
##
## ANOVA Result:
##
##
## | | npar| AIC| BIC| logLik| deviance| Chisq| Df| Pr(>Chisq)|
## |:-------------|----:|--------:|--------:|---------:|--------:|-----:|--:|----------:|
## |reduced_model | 10| 71096.92| 71178.84| -35538.46| 71076.92| NA| NA| NA|
## |full_model | 15| 71103.41| 71226.28| -35536.70| 71073.41| 3.51| 5| 0.62|##長期導向
full_model <- lmer(Zindirect_low ~ Z環境態度2 + Z長期導向 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度2 * Z長期導向 + (1 + Z人均GDP + Z長期導向 | country), data = data)
reduced_model <- lmer(Zindirect_low ~ Z環境態度2 + Z長期導向 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度2 * Z長期導向 + (1 | country), data = data)
anova_result <- anova(reduced_model, full_model)
summarize_model(full_model, anova_result)## Fixed Effects:
##
##
## | | Estimate| Std. Error| df| t value| p-value|
## |:--------------------|--------:|----------:|--------:|-------:|-------:|
## |(Intercept) | -0.02| 0.02| 12.17| -0.95| 0.36|
## |Z環境態度2 | 0.22| 0.01| 27012.61| 36.56| 0.00|
## |Z長期導向 | 0.14| 0.04| 7.05| 3.64| 0.01|
## |Z年齡 | 0.01| 0.01| 26126.10| 1.88| 0.06|
## |Z性別 | 0.00| 0.01| 27097.14| -0.78| 0.43|
## |Z教育程度 | 0.13| 0.01| 26027.95| 19.55| 0.00|
## |Z社經地位 | 0.02| 0.01| 24612.24| 2.62| 0.01|
## |Z環境態度2:Z長期導向 | 0.02| 0.01| 26995.98| 3.27| 0.00|
##
## Random Effects:
##
##
## |grp |var1 |var2 | vcov| sdcor|
## |:--------|:-----------|:---------|-----:|-----:|
## |country |(Intercept) |NA | 0.00| 0.00|
## |country |Z人均GDP |NA | 0.04| 0.21|
## |country |Z長期導向 |NA | 0.03| 0.17|
## |country |(Intercept) |Z人均GDP | 0.00| NaN|
## |country |(Intercept) |Z長期導向 | 0.00| NaN|
## |country |Z人均GDP |Z長期導向 | -0.03| -0.72|
## |Residual |NA |NA | 0.84| 0.92|
##
## ANOVA Result:
##
##
## | | npar| AIC| BIC| logLik| deviance| Chisq| Df| Pr(>Chisq)|
## |:-------------|----:|--------:|--------:|---------:|--------:|-----:|--:|----------:|
## |reduced_model | 10| 72383.37| 72465.45| -36181.68| 72363.37| NA| NA| NA|
## |full_model | 15| 72388.92| 72512.05| -36179.46| 72358.92| 4.45| 5| 0.49|##放任約束
full_model <- lmer(Zindirect_low ~ Z環境態度2 + Z放任約束 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度2 * Z放任約束 + (1 + Z人均GDP + Z放任約束 | country), data = data)
reduced_model <- lmer(Zindirect_low ~ Z環境態度2 + Z放任約束 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度2 * Z放任約束 + (1 | country), data = data)
anova_result <- anova(reduced_model, full_model)
summarize_model(full_model, anova_result)## Fixed Effects:
##
##
## | | Estimate| Std. Error| df| t value| p-value|
## |:--------------------|--------:|----------:|--------:|-------:|-------:|
## |(Intercept) | -0.06| 0.03| 20.90| -2.10| 0.05|
## |Z環境態度2 | 0.22| 0.01| 25467.09| 36.13| 0.00|
## |Z放任約束 | 0.10| 0.02| 15.24| 4.57| 0.00|
## |Z年齡 | 0.01| 0.01| 26204.06| 1.93| 0.05|
## |Z性別 | 0.00| 0.01| 26459.29| -0.86| 0.39|
## |Z教育程度 | 0.13| 0.01| 25049.30| 18.97| 0.00|
## |Z社經地位 | 0.02| 0.01| 25484.79| 2.95| 0.00|
## |Z環境態度2:Z放任約束 | 0.05| 0.01| 26195.04| 8.06| 0.00|
##
## Random Effects:
##
##
## |grp |var1 |var2 | vcov| sdcor|
## |:--------|:-----------|:---------|----:|-----:|
## |country |(Intercept) |NA | 0.02| 0.15|
## |country |Z人均GDP |NA | 0.00| 0.02|
## |country |Z放任約束 |NA | 0.00| 0.05|
## |country |(Intercept) |Z人均GDP | 0.00| 1.00|
## |country |(Intercept) |Z放任約束 | 0.01| 0.80|
## |country |Z人均GDP |Z放任約束 | 0.00| 0.80|
## |Residual |NA |NA | 0.84| 0.92|
##
## ANOVA Result:
##
##
## | | npar| AIC| BIC| logLik| deviance| Chisq| Df| Pr(>Chisq)|
## |:-------------|----:|--------:|--------:|---------:|--------:|-----:|--:|----------:|
## |reduced_model | 10| 70668.44| 70750.29| -35324.22| 70648.44| NA| NA| NA|
## |full_model | 15| 70671.16| 70793.93| -35320.58| 70641.16| 7.28| 5| 0.2|##間接高成本行動 ##直接效果
full_model <- lmer(Zindirect_high ~ Z環境態度2 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 +
(1 + Z人均GDP | country), data = data)
reduced_model <- lmer(Zindirect_high ~ Z環境態度2 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 +
(1 | country), data = data)
anova_result <- anova(reduced_model, full_model)
summarize_model(full_model, anova_result)## Fixed Effects:
##
##
## | | Estimate| Std. Error| df| t value| p-value|
## |:-----------|--------:|----------:|--------:|-------:|-------:|
## |(Intercept) | -0.01| 0.02| 34.07| -0.58| 0.57|
## |Z環境態度2 | 0.10| 0.01| 26506.69| 16.25| 0.00|
## |Z年齡 | -0.01| 0.01| 26751.82| -1.11| 0.27|
## |Z性別 | -0.02| 0.01| 27114.18| -3.25| 0.00|
## |Z教育程度 | 0.05| 0.01| 25622.26| 7.00| 0.00|
## |Z社經地位 | 0.00| 0.01| 25116.89| 0.63| 0.53|
##
## Random Effects:
##
##
## |grp |var1 |var2 | vcov| sdcor|
## |:--------|:-----------|:--------|----:|-----:|
## |country |(Intercept) |NA | 0.02| 0.13|
## |country |Z人均GDP |NA | 0.00| 0.00|
## |country |(Intercept) |Z人均GDP | 0.00| 1.00|
## |Residual |NA |NA | 0.97| 0.99|
##
## ANOVA Result:
##
##
## | | npar| AIC| BIC| logLik| deviance| Chisq| Df| Pr(>Chisq)|
## |:-------------|----:|--------:|--------:|---------:|--------:|-----:|--:|----------:|
## |reduced_model | 8| 76368.85| 76434.52| -38176.43| 76352.85| NA| NA| NA|
## |full_model | 10| 76372.85| 76454.94| -38176.43| 76352.85| 0| 2| 1|##政府信任
full_model <- lmer(Zindirect_high ~ Z環境態度2 + Z政府信任 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度2 * Z政府信任 + (1 + Z人均GDP | country), data = data)
reduced_model <- lmer(Zindirect_high ~ Z環境態度2 + Z政府信任 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度2 * Z政府信任 + (1 | country), data = data)
anova_result <- anova(reduced_model, full_model)
summarize_model(full_model, anova_result)## Fixed Effects:
##
##
## | | Estimate| Std. Error| df| t value| p-value|
## |:--------------------|--------:|----------:|--------:|-------:|-------:|
## |(Intercept) | -0.01| 0.02| 33.95| -0.49| 0.63|
## |Z環境態度2 | 0.10| 0.01| 26385.99| 16.01| 0.00|
## |Z政府信任 | 0.00| 0.01| 24114.48| 0.77| 0.44|
## |Z年齡 | -0.01| 0.01| 26756.71| -1.12| 0.26|
## |Z性別 | -0.02| 0.01| 27112.11| -3.24| 0.00|
## |Z教育程度 | 0.05| 0.01| 25561.21| 7.01| 0.00|
## |Z社經地位 | 0.00| 0.01| 25290.46| 0.55| 0.58|
## |Z環境態度2:Z政府信任 | -0.01| 0.01| 27103.62| -1.64| 0.10|
##
## Random Effects:
##
##
## |grp |var1 |var2 | vcov| sdcor|
## |:--------|:-----------|:--------|----:|-----:|
## |country |(Intercept) |NA | 0.02| 0.13|
## |country |Z人均GDP |NA | 0.00| 0.00|
## |country |(Intercept) |Z人均GDP | 0.00| 1.00|
## |Residual |NA |NA | 0.97| 0.99|
##
## ANOVA Result:
##
##
## | | npar| AIC| BIC| logLik| deviance| Chisq| Df| Pr(>Chisq)|
## |:-------------|----:|--------:|--------:|---------:|--------:|-----:|--:|----------:|
## |reduced_model | 10| 76369.68| 76451.77| -38174.84| 76349.68| NA| NA| NA|
## |full_model | 12| 76373.68| 76472.19| -38174.84| 76349.68| 0| 2| 1|##權力距離
full_model <- lmer(Zindirect_high ~ Z環境態度2 + Z權力距離 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度2 * Z權力距離 + (1 + Z人均GDP + Z權力距離 | country), data = data)
reduced_model <- lmer(Zindirect_high ~ Z環境態度2 + Z權力距離 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度2 * Z權力距離 + (1 | country), data = data)
anova_result <- anova(reduced_model, full_model)
summarize_model(full_model, anova_result)## Fixed Effects:
##
##
## | | Estimate| Std. Error| df| t value| p-value|
## |:--------------------|--------:|----------:|--------:|-------:|-------:|
## |(Intercept) | -0.03| 0.02| 24.36| -1.13| 0.27|
## |Z環境態度2 | 0.10| 0.01| 24070.27| 15.81| 0.00|
## |Z權力距離 | 0.05| 0.02| 13.17| 2.08| 0.06|
## |Z年齡 | -0.01| 0.01| 26237.23| -1.05| 0.29|
## |Z性別 | -0.02| 0.01| 26645.87| -3.28| 0.00|
## |Z教育程度 | 0.05| 0.01| 14117.72| 6.73| 0.00|
## |Z社經地位 | 0.00| 0.01| 25839.23| 0.70| 0.48|
## |Z環境態度2:Z權力距離 | -0.01| 0.01| 25580.19| -2.16| 0.03|
##
## Random Effects:
##
##
## |grp |var1 |var2 | vcov| sdcor|
## |:--------|:-----------|:---------|----:|-----:|
## |country |(Intercept) |NA | 0.02| 0.12|
## |country |Z人均GDP |NA | 0.00| 0.03|
## |country |Z權力距離 |NA | 0.00| 0.07|
## |country |(Intercept) |Z人均GDP | 0.00| -0.70|
## |country |(Intercept) |Z權力距離 | 0.00| 0.46|
## |country |Z人均GDP |Z權力距離 | 0.00| 0.31|
## |Residual |NA |NA | 0.97| 0.98|
##
## ANOVA Result:
##
##
## | | npar| AIC| BIC| logLik| deviance| Chisq| Df| Pr(>Chisq)|
## |:-------------|----:|--------:|--------:|---------:|--------:|-----:|--:|----------:|
## |reduced_model | 10| 74966.51| 75048.42| -37473.25| 74946.51| NA| NA| NA|
## |full_model | 15| 74974.35| 75097.22| -37472.17| 74944.35| 2.16| 5| 0.83|##個人主義
full_model <- lmer(Zindirect_high ~ Z環境態度2 + Z個人主義 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度2 * Z個人主義 + (1 + Z人均GDP + Z個人主義 | country), data = data)
reduced_model <- lmer(Zindirect_high ~ Z環境態度2 + Z個人主義 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度2 * Z個人主義 + (1 | country), data = data)
anova_result <- anova(reduced_model, full_model)
summarize_model(full_model, anova_result)## Fixed Effects:
##
##
## | | Estimate| Std. Error| df| t value| p-value|
## |:--------------------|--------:|----------:|--------:|-------:|-------:|
## |(Intercept) | -0.03| 0.02| 15.17| -1.52| 0.15|
## |Z環境態度2 | 0.10| 0.01| 26612.45| 15.58| 0.00|
## |Z個人主義 | -0.04| 0.03| 11.06| -1.13| 0.28|
## |Z年齡 | -0.01| 0.01| 26179.04| -1.04| 0.30|
## |Z性別 | -0.02| 0.01| 27109.88| -3.53| 0.00|
## |Z教育程度 | 0.05| 0.01| 19232.67| 6.96| 0.00|
## |Z社經地位 | 0.01| 0.01| 19851.13| 0.75| 0.45|
## |Z環境態度2:Z個人主義 | 0.03| 0.01| 27019.39| 5.10| 0.00|
##
## Random Effects:
##
##
## |grp |var1 |var2 | vcov| sdcor|
## |:--------|:-----------|:---------|-----:|-----:|
## |country |(Intercept) |NA | 0.00| 0.00|
## |country |Z人均GDP |NA | 0.06| 0.24|
## |country |Z個人主義 |NA | 0.12| 0.34|
## |country |(Intercept) |Z人均GDP | 0.00| NaN|
## |country |(Intercept) |Z個人主義 | 0.00| NaN|
## |country |Z人均GDP |Z個人主義 | -0.08| -0.97|
## |Residual |NA |NA | 0.97| 0.99|
##
## ANOVA Result:
##
##
## | | npar| AIC| BIC| logLik| deviance| Chisq| Df| Pr(>Chisq)|
## |:-------------|----:|--------:|--------:|---------:|--------:|-----:|--:|----------:|
## |reduced_model | 10| 76344.59| 76426.67| -38162.29| 76324.59| NA| NA| NA|
## |full_model | 15| 76372.99| 76496.12| -38171.49| 76342.99| 0| 5| 1|##成就動機
full_model <- lmer(Zindirect_high ~ Z環境態度2 + Z成就動機 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度2 * Z成就動機 + (1 + Z人均GDP + Z成就動機 | country), data = data)
reduced_model <- lmer(Zindirect_high ~ Z環境態度2 + Z成就動機 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度2 * Z成就動機 + (1 | country), data = data)
anova_result <- anova(reduced_model, full_model)
summarize_model(full_model, anova_result)## Fixed Effects:
##
##
## | | Estimate| Std. Error| df| t value| p-value|
## |:--------------------|--------:|----------:|--------:|-------:|-------:|
## |(Intercept) | -0.01| 0.02| 32.07| -0.23| 0.82|
## |Z環境態度2 | 0.10| 0.01| 24387.85| 15.94| 0.00|
## |Z成就動機 | -0.04| 0.02| 5.87| -2.20| 0.07|
## |Z年齡 | -0.01| 0.01| 25706.16| -1.05| 0.30|
## |Z性別 | -0.02| 0.01| 26653.99| -3.24| 0.00|
## |Z教育程度 | 0.05| 0.01| 17275.41| 6.74| 0.00|
## |Z社經地位 | 0.00| 0.01| 19874.73| 0.68| 0.49|
## |Z環境態度2:Z成就動機 | -0.01| 0.01| 22724.73| -1.88| 0.06|
##
## Random Effects:
##
##
## |grp |var1 |var2 | vcov| sdcor|
## |:--------|:-----------|:---------|-----:|-----:|
## |country |(Intercept) |NA | 0.02| 0.13|
## |country |Z人均GDP |NA | 0.00| 0.02|
## |country |Z成就動機 |NA | 0.01| 0.07|
## |country |(Intercept) |Z人均GDP | 0.00| -1.00|
## |country |(Intercept) |Z成就動機 | -0.01| -1.00|
## |country |Z人均GDP |Z成就動機 | 0.00| 1.00|
## |Residual |NA |NA | 0.97| 0.98|
##
## ANOVA Result:
##
##
## | | npar| AIC| BIC| logLik| deviance| Chisq| Df| Pr(>Chisq)|
## |:-------------|----:|--------:|--------:|---------:|--------:|-----:|--:|----------:|
## |reduced_model | 10| 74969.10| 75051.02| -37474.55| 74949.10| NA| NA| NA|
## |full_model | 15| 74974.73| 75097.61| -37472.37| 74944.73| 4.37| 5| 0.5|##不確定規避
full_model <- lmer(Zindirect_high ~ Z環境態度2 + Z不確定規避 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度2 * Z不確定規避 + (1 + Z人均GDP + Z不確定規避 | country), data = data)
reduced_model <- lmer(Zindirect_high ~ Z環境態度2 + Z不確定規避 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度2 * Z不確定規避 + (1 | country), data = data)
anova_result <- anova(reduced_model, full_model)
summarize_model(full_model, anova_result)## Fixed Effects:
##
##
## | | Estimate| Std. Error| df| t value| p-value|
## |:----------------------|--------:|----------:|--------:|-------:|-------:|
## |(Intercept) | -0.01| 0.02| 20.46| -0.75| 0.46|
## |Z環境態度2 | 0.10| 0.01| 24567.73| 15.81| 0.00|
## |Z不確定規避 | 0.05| 0.01| 16.04| 3.86| 0.00|
## |Z年齡 | -0.01| 0.01| 25786.05| -0.93| 0.35|
## |Z性別 | -0.02| 0.01| 26571.74| -3.06| 0.00|
## |Z教育程度 | 0.05| 0.01| 22617.01| 6.74| 0.00|
## |Z社經地位 | 0.01| 0.01| 21963.58| 0.74| 0.46|
## |Z環境態度2:Z不確定規避 | 0.02| 0.01| 21379.89| 3.64| 0.00|
##
## Random Effects:
##
##
## |grp |var1 |var2 | vcov| sdcor|
## |:--------|:-----------|:-----------|----:|-----:|
## |country |(Intercept) |NA | 0.01| 0.09|
## |country |Z人均GDP |NA | 0.01| 0.07|
## |country |Z不確定規避 |NA | 0.01| 0.08|
## |country |(Intercept) |Z人均GDP | 0.00| 0.72|
## |country |(Intercept) |Z不確定規避 | 0.01| 0.93|
## |country |Z人均GDP |Z不確定規避 | 0.01| 0.92|
## |Residual |NA |NA | 0.97| 0.98|
##
## ANOVA Result:
##
##
## | | npar| AIC| BIC| logLik| deviance| Chisq| Df| Pr(>Chisq)|
## |:-------------|----:|--------:|--------:|---------:|--------:|-----:|--:|----------:|
## |reduced_model | 10| 74956.66| 75038.57| -37468.33| 74936.66| NA| NA| NA|
## |full_model | 15| 74957.53| 75080.40| -37463.76| 74927.53| 9.13| 5| 0.1|##長期導向
full_model <- lmer(Zindirect_high ~ Z環境態度2 + Z長期導向 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度2 * Z長期導向 + (1 + Z人均GDP + Z長期導向 | country), data = data)
reduced_model <- lmer(Zindirect_high ~ Z環境態度2 + Z長期導向 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度2 * Z長期導向 + (1 | country), data = data)
anova_result <- anova(reduced_model, full_model)
summarize_model(full_model, anova_result)## Fixed Effects:
##
##
## | | Estimate| Std. Error| df| t value| p-value|
## |:--------------------|--------:|----------:|--------:|-------:|-------:|
## |(Intercept) | -0.04| 0.02| 50.06| -1.85| 0.07|
## |Z環境態度2 | 0.10| 0.01| 26444.57| 16.25| 0.00|
## |Z長期導向 | -0.03| 0.06| 21.62| -0.51| 0.61|
## |Z年齡 | -0.01| 0.01| 26844.23| -1.06| 0.29|
## |Z性別 | -0.02| 0.01| 27115.53| -3.25| 0.00|
## |Z教育程度 | 0.05| 0.01| 25182.38| 7.11| 0.00|
## |Z社經地位 | 0.00| 0.01| 25426.21| 0.61| 0.54|
## |Z環境態度2:Z長期導向 | -0.01| 0.01| 27068.16| -0.82| 0.41|
##
## Random Effects:
##
##
## |grp |var1 |var2 | vcov| sdcor|
## |:--------|:-----------|:---------|----:|-----:|
## |country |(Intercept) |NA | 0.00| 0.00|
## |country |Z人均GDP |NA | 0.01| 0.09|
## |country |Z長期導向 |NA | 0.11| 0.34|
## |country |(Intercept) |Z人均GDP | 0.00| NaN|
## |country |(Intercept) |Z長期導向 | 0.00| NaN|
## |country |Z人均GDP |Z長期導向 | 0.03| 1.00|
## |Residual |NA |NA | 0.97| 0.99|
##
## ANOVA Result:
##
##
## | | npar| AIC| BIC| logLik| deviance| Chisq| Df| Pr(>Chisq)|
## |:-------------|----:|--------:|--------:|---------:|--------:|-----:|--:|----------:|
## |reduced_model | 10| 76368.85| 76450.93| -38174.42| 76348.85| NA| NA| NA|
## |full_model | 15| 76408.68| 76531.81| -38189.34| 76378.68| 0| 5| 1|##放任約束
full_model <- lmer(Zindirect_high ~ Z環境態度2 + Z放任約束 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度2 * Z放任約束 + (1 + Z人均GDP + Z放任約束 | country), data = data)
reduced_model <- lmer(Zindirect_high ~ Z環境態度2 + Z放任約束 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度2 * Z放任約束 + (1 | country), data = data)
anova_result <- anova(reduced_model, full_model)
summarize_model(full_model, anova_result)## Fixed Effects:
##
##
## | | Estimate| Std. Error| df| t value| p-value|
## |:--------------------|--------:|----------:|--------:|-------:|-------:|
## |(Intercept) | -0.03| 0.02| 21.02| -1.45| 0.16|
## |Z環境態度2 | 0.10| 0.01| 25869.61| 15.79| 0.00|
## |Z放任約束 | -0.03| 0.02| 26.41| -1.15| 0.26|
## |Z年齡 | -0.01| 0.01| 25543.36| -1.14| 0.26|
## |Z性別 | -0.02| 0.01| 26486.88| -3.44| 0.00|
## |Z教育程度 | 0.05| 0.01| 19514.11| 6.85| 0.00|
## |Z社經地位 | 0.00| 0.01| 24138.02| 0.45| 0.65|
## |Z環境態度2:Z放任約束 | 0.00| 0.01| 26293.33| -0.59| 0.56|
##
## Random Effects:
##
##
## |grp |var1 |var2 | vcov| sdcor|
## |:--------|:-----------|:---------|----:|-----:|
## |country |(Intercept) |NA | 0.02| 0.13|
## |country |Z人均GDP |NA | 0.00| 0.03|
## |country |Z放任約束 |NA | 0.00| 0.01|
## |country |(Intercept) |Z人均GDP | 0.00| -1.00|
## |country |(Intercept) |Z放任約束 | 0.00| -1.00|
## |country |Z人均GDP |Z放任約束 | 0.00| 1.00|
## |Residual |NA |NA | 0.96| 0.98|
##
## ANOVA Result:
##
##
## | | npar| AIC| BIC| logLik| deviance| Chisq| Df| Pr(>Chisq)|
## |:-------------|----:|--------:|--------:|---------:|--------:|-----:|--:|----------:|
## |reduced_model | 10| 74301.86| 74383.71| -37140.93| 74281.86| NA| NA| NA|
## |full_model | 15| 74310.85| 74433.63| -37140.43| 74280.85| 1.01| 5| 0.96|