##載入套件
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環境態度1 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 +
(1 + Z人均GDP | country), data = data)
reduced_model <- lmer(Zdirect_low ~ Z環境態度1 + 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.06| 0.05| 0.23| 1.18| 0.70|
## |Z環境態度1 | 0.26| 0.01| 27104.23| 48.00| 0.00|
## |Z年齡 | 0.15| 0.01| 27091.59| 26.50| 0.00|
## |Z性別 | 0.06| 0.01| 27098.59| 11.58| 0.00|
## |Z教育程度 | 0.04| 0.01| 25632.50| 6.84| 0.00|
## |Z社經地位 | -0.01| 0.01| 27115.52| -1.16| 0.25|
##
## Random Effects:
##
##
## |grp |var1 |var2 | vcov| sdcor|
## |:--------|:-----------|:--------|-----:|-----:|
## |country |(Intercept) |NA | 0.20| 0.44|
## |country |Z人均GDP |NA | 0.05| 0.23|
## |country |(Intercept) |Z人均GDP | -0.10| -0.99|
## |Residual |NA |NA | 0.68| 0.83|
##
## ANOVA Result:
##
##
## | | npar| AIC| BIC| logLik| deviance| Chisq| Df| Pr(>Chisq)|
## |:-------------|----:|--------:|--------:|---------:|--------:|-----:|--:|----------:|
## |reduced_model | 8| 66898.48| 66964.15| -33441.24| 66882.48| NA| NA| NA|
## |full_model | 10| 66894.25| 66976.33| -33437.12| 66874.25| 8.23| 2| 0.02|##政府信任
full_model <- lmer(Zdirect_low ~ Z環境態度1 + Z政府信任 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度1 * Z政府信任 + (1 + Z人均GDP | country), data = data)
reduced_model <- lmer(Zdirect_low ~ Z環境態度1 + Z政府信任 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度1 * 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.05| 9.24| 1.34| 0.21|
## |Z環境態度1 | 0.26| 0.01| 27115.98| 47.95| 0.00|
## |Z政府信任 | -0.01| 0.01| 26975.98| -2.51| 0.01|
## |Z年齡 | 0.15| 0.01| 27118.54| 26.53| 0.00|
## |Z性別 | 0.06| 0.01| 27101.65| 11.62| 0.00|
## |Z教育程度 | 0.04| 0.01| 27088.65| 6.88| 0.00|
## |Z社經地位 | -0.01| 0.01| 27124.97| -0.92| 0.36|
## |Z環境態度1:Z政府信任 | 0.01| 0.00| 27103.86| 1.30| 0.19|
##
## Random Effects:
##
##
## |grp |var1 |var2 | vcov| sdcor|
## |:--------|:-----------|:--------|-----:|-----:|
## |country |(Intercept) |NA | 0.20| 0.45|
## |country |Z人均GDP |NA | 0.05| 0.23|
## |country |(Intercept) |Z人均GDP | -0.10| -1.00|
## |Residual |NA |NA | 0.68| 0.83|
##
## ANOVA Result:
##
##
## | | npar| AIC| BIC| logLik| deviance| Chisq| Df| Pr(>Chisq)|
## |:-------------|----:|--------:|--------:|---------:|--------:|-----:|--:|----------:|
## |reduced_model | 10| 66895.07| 66977.16| -33437.54| 66875.07| NA| NA| NA|
## |full_model | 12| 66890.87| 66989.38| -33433.44| 66866.87| 8.2| 2| 0.02|##權力距離
full_model <- lmer(Zdirect_low ~ Z環境態度1 + Z權力距離 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度1 * Z權力距離 + (1 + Z人均GDP + Z權力距離 | country), data = data)
reduced_model <- lmer(Zdirect_low ~ Z環境態度1 + Z權力距離 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度1 * 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.06| 14.35| 0.69| 0.50|
## |Z環境態度1 | 0.26| 0.01| 26652.92| 47.46| 0.00|
## |Z權力距離 | -0.13| 0.06| 14.58| -2.14| 0.05|
## |Z年齡 | 0.15| 0.01| 26652.71| 25.99| 0.00|
## |Z性別 | 0.06| 0.01| 26638.38| 11.05| 0.00|
## |Z教育程度 | 0.04| 0.01| 26660.47| 6.62| 0.00|
## |Z社經地位 | -0.01| 0.01| 26643.80| -1.37| 0.17|
## |Z環境態度1:Z權力距離 | -0.02| 0.01| 26650.29| -3.25| 0.00|
##
## Random Effects:
##
##
## |grp |var1 |var2 | vcov| sdcor|
## |:--------|:-----------|:---------|-----:|-----:|
## |country |(Intercept) |NA | 0.06| 0.25|
## |country |Z人均GDP |NA | 0.21| 0.46|
## |country |Z權力距離 |NA | 0.16| 0.40|
## |country |(Intercept) |Z人均GDP | -0.06| -0.53|
## |country |(Intercept) |Z權力距離 | -0.02| -0.15|
## |country |Z人均GDP |Z權力距離 | 0.17| 0.92|
## |Residual |NA |NA | 0.69| 0.83|
##
## ANOVA Result:
##
##
## | | npar| AIC| BIC| logLik| deviance| Chisq| Df| Pr(>Chisq)|
## |:-------------|----:|--------:|--------:|---------:|--------:|-----:|--:|----------:|
## |reduced_model | 10| 65802.58| 65884.49| -32891.29| 65782.58| NA| NA| NA|
## |full_model | 15| 65798.37| 65921.25| -32884.19| 65768.37| 14.21| 5| 0.01|##個人主義
full_model <- lmer(Zdirect_low ~ Z環境態度1 + Z個人主義 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度1 * Z個人主義 + (1 + Z人均GDP + Z個人主義 | country), data = data)
reduced_model <- lmer(Zdirect_low ~ Z環境態度1 + Z個人主義 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度1 * 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.07| 0.06| 20.41| -1.20| 0.24|
## |Z環境態度1 | 0.26| 0.01| 27078.55| 47.79| 0.00|
## |Z個人主義 | 0.45| 0.09| 23.00| 5.23| 0.00|
## |Z年齡 | 0.15| 0.01| 26810.62| 26.51| 0.00|
## |Z性別 | 0.06| 0.01| 27104.31| 11.58| 0.00|
## |Z教育程度 | 0.04| 0.01| 26844.47| 6.87| 0.00|
## |Z社經地位 | -0.01| 0.01| 26814.79| -1.22| 0.22|
## |Z環境態度1:Z個人主義 | 0.00| 0.01| 27110.52| -0.33| 0.74|
##
## Random Effects:
##
##
## |grp |var1 |var2 | vcov| sdcor|
## |:--------|:-----------|:---------|-----:|-----:|
## |country |(Intercept) |NA | 0.05| 0.22|
## |country |Z人均GDP |NA | 0.00| 0.07|
## |country |Z個人主義 |NA | 0.20| 0.45|
## |country |(Intercept) |Z人均GDP | 0.00| -0.35|
## |country |(Intercept) |Z個人主義 | -0.07| -0.75|
## |country |Z人均GDP |Z個人主義 | -0.01| -0.36|
## |Residual |NA |NA | 0.68| 0.83|
##
## ANOVA Result:
##
##
## | | npar| AIC| BIC| logLik| deviance| Chisq| Df| Pr(>Chisq)|
## |:-------------|----:|--------:|--------:|---------:|--------:|-----:|--:|----------:|
## |reduced_model | 10| 66885.17| 66967.25| -33432.58| 66865.17| NA| NA| NA|
## |full_model | 15| 66886.90| 67010.03| -33428.45| 66856.90| 8.27| 5| 0.14|##成就動機
full_model <- lmer(Zdirect_low ~ Z環境態度1 + Z成就動機 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度1 * Z成就動機 + (1 + Z人均GDP + Z成就動機 | country), data = data)
reduced_model <- lmer(Zdirect_low ~ Z環境態度1 + Z成就動機 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度1 * 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.11| 0.06| 12.11| 1.82| 0.09|
## |Z環境態度1 | 0.26| 0.01| 26652.43| 47.66| 0.00|
## |Z成就動機 | 0.01| 0.04| 7.96| 0.14| 0.90|
## |Z年齡 | 0.15| 0.01| 26653.19| 26.15| 0.00|
## |Z性別 | 0.06| 0.01| 26638.98| 11.04| 0.00|
## |Z教育程度 | 0.04| 0.01| 26344.73| 6.65| 0.00|
## |Z社經地位 | -0.01| 0.01| 26661.76| -1.35| 0.18|
## |Z環境態度1:Z成就動機 | -0.02| 0.01| 26653.82| -4.01| 0.00|
##
## Random Effects:
##
##
## |grp |var1 |var2 | vcov| sdcor|
## |:--------|:-----------|:---------|-----:|-----:|
## |country |(Intercept) |NA | 0.20| 0.44|
## |country |Z人均GDP |NA | 0.07| 0.27|
## |country |Z成就動機 |NA | 0.01| 0.11|
## |country |(Intercept) |Z人均GDP | -0.11| -0.97|
## |country |(Intercept) |Z成就動機 | 0.04| 0.91|
## |country |Z人均GDP |Z成就動機 | -0.02| -0.77|
## |Residual |NA |NA | 0.68| 0.83|
##
## ANOVA Result:
##
##
## | | npar| AIC| BIC| logLik| deviance| Chisq| Df| Pr(>Chisq)|
## |:-------------|----:|--------:|--------:|---------:|--------:|-----:|--:|----------:|
## |reduced_model | 10| 65797.88| 65879.80| -32888.94| 65777.88| NA| NA| NA|
## |full_model | 15| 65795.82| 65918.69| -32882.91| 65765.82| 12.06| 5| 0.03|##不確定規避
full_model <- lmer(Zdirect_low ~ Z環境態度1 + Z不確定規避 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度1 * Z不確定規避 + (1 + Z人均GDP + Z不確定規避 | country), data = data)
reduced_model <- lmer(Zdirect_low ~ Z環境態度1 + Z不確定規避 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度1 * 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.07| 26.23| -0.19| 0.85|
## |Z環境態度1 | 0.26| 0.01| 26653.35| 47.72| 0.00|
## |Z不確定規避 | -0.05| 0.04| 20.03| -1.48| 0.15|
## |Z年齡 | 0.15| 0.01| 26651.10| 26.03| 0.00|
## |Z性別 | 0.06| 0.01| 26638.87| 11.08| 0.00|
## |Z教育程度 | 0.04| 0.01| 26465.12| 6.66| 0.00|
## |Z社經地位 | -0.01| 0.01| 26653.52| -1.38| 0.17|
## |Z環境態度1:Z不確定規避 | -0.02| 0.01| 26652.62| -3.17| 0.00|
##
## Random Effects:
##
##
## |grp |var1 |var2 | vcov| sdcor|
## |:--------|:-----------|:-----------|-----:|-----:|
## |country |(Intercept) |NA | 0.16| 0.40|
## |country |Z人均GDP |NA | 0.04| 0.19|
## |country |Z不確定規避 |NA | 0.02| 0.15|
## |country |(Intercept) |Z人均GDP | -0.04| -0.56|
## |country |(Intercept) |Z不確定規避 | 0.05| 0.86|
## |country |Z人均GDP |Z不確定規避 | 0.00| -0.05|
## |Residual |NA |NA | 0.69| 0.83|
##
## ANOVA Result:
##
##
## | | npar| AIC| BIC| logLik| deviance| Chisq| Df| Pr(>Chisq)|
## |:-------------|----:|--------:|--------:|---------:|--------:|-----:|--:|----------:|
## |reduced_model | 10| 65804.07| 65885.99| -32892.04| 65784.07| NA| NA| NA|
## |full_model | 15| 65801.58| 65924.45| -32885.79| 65771.58| 12.49| 5| 0.03|##長期導向
full_model <- lmer(Zdirect_low ~ Z環境態度1 + Z長期導向 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度1 * Z長期導向 + (1 + Z人均GDP + Z長期導向 | country), data = data)
reduced_model <- lmer(Zdirect_low ~ Z環境態度1 + Z長期導向 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度1 * 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| 13.58| -0.16| 0.87|
## |Z環境態度1 | 0.26| 0.01| 27119.39| 48.03| 0.00|
## |Z長期導向 | 0.34| 0.08| 3.34| 4.48| 0.02|
## |Z年齡 | 0.15| 0.01| 27115.33| 26.48| 0.00|
## |Z性別 | 0.06| 0.01| 27101.39| 11.64| 0.00|
## |Z教育程度 | 0.04| 0.01| 27111.23| 6.80| 0.00|
## |Z社經地位 | -0.01| 0.01| 27125.13| -1.14| 0.25|
## |Z環境態度1:Z長期導向 | -0.01| 0.01| 27107.21| -1.66| 0.10|
##
## Random Effects:
##
##
## |grp |var1 |var2 | vcov| sdcor|
## |:--------|:-----------|:---------|-----:|-----:|
## |country |(Intercept) |NA | 0.13| 0.37|
## |country |Z人均GDP |NA | 0.04| 0.20|
## |country |Z長期導向 |NA | 0.04| 0.21|
## |country |(Intercept) |Z人均GDP | -0.06| -0.85|
## |country |(Intercept) |Z長期導向 | 0.06| 0.74|
## |country |Z人均GDP |Z長期導向 | -0.04| -0.98|
## |Residual |NA |NA | 0.68| 0.83|
##
## ANOVA Result:
##
##
## | | npar| AIC| BIC| logLik| deviance| Chisq| Df| Pr(>Chisq)|
## |:-------------|----:|--------:|--------:|---------:|--------:|-----:|--:|----------:|
## |reduced_model | 10| 66893.35| 66975.43| -33436.67| 66873.35| NA| NA| NA|
## |full_model | 15| 66892.35| 67015.48| -33431.17| 66862.35| 11| 5| 0.05|##放任約束
full_model <- lmer(Zdirect_low ~ Z環境態度1 + Z放任約束 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度1 * Z放任約束 + (1 + Z人均GDP + Z放任約束 | country), data = data)
reduced_model <- lmer(Zdirect_low ~ Z環境態度1 + Z放任約束 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度1 * 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.12| 0.05| 26.15| 2.15| 0.04|
## |Z環境態度1 | 0.26| 0.01| 26479.93| 47.30| 0.00|
## |Z放任約束 | 0.20| 0.06| 11.58| 3.07| 0.01|
## |Z年齡 | 0.15| 0.01| 26488.10| 26.22| 0.00|
## |Z性別 | 0.06| 0.01| 26472.22| 11.89| 0.00|
## |Z教育程度 | 0.04| 0.01| 26495.64| 6.28| 0.00|
## |Z社經地位 | 0.00| 0.01| 26492.61| -0.22| 0.82|
## |Z環境態度1:Z放任約束 | -0.01| 0.01| 26480.04| -2.53| 0.01|
##
## Random Effects:
##
##
## |grp |var1 |var2 | vcov| sdcor|
## |:--------|:-----------|:---------|-----:|-----:|
## |country |(Intercept) |NA | 0.00| 0.00|
## |country |Z人均GDP |NA | 0.28| 0.53|
## |country |Z放任約束 |NA | 0.14| 0.37|
## |country |(Intercept) |Z人均GDP | 0.00| NaN|
## |country |(Intercept) |Z放任約束 | 0.00| NaN|
## |country |Z人均GDP |Z放任約束 | -0.16| -0.79|
## |Residual |NA |NA | 0.68| 0.82|
##
## ANOVA Result:
##
##
## | | npar| AIC| BIC| logLik| deviance| Chisq| Df| Pr(>Chisq)|
## |:-------------|----:|--------:|--------:|---------:|--------:|-----:|--:|----------:|
## |reduced_model | 10| 65110.48| 65192.33| -32545.24| 65090.48| NA| NA| NA|
## |full_model | 15| 65112.93| 65235.71| -32541.47| 65082.93| 7.54| 5| 0.18|##直接高成本行動 ##直接效果
full_model <- lmer(Zdirect_high ~ Z環境態度1 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 +
(1 + Z人均GDP | country), data = data)
reduced_model <- lmer(Zdirect_high ~ Z環境態度1 + 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.02| 0.04| 32.45| -0.55| 0.58|
## |Z環境態度1 | 0.29| 0.01| 27115.50| 47.31| 0.00|
## |Z年齡 | 0.11| 0.01| 27128.95| 16.86| 0.00|
## |Z性別 | 0.07| 0.01| 27110.35| 11.61| 0.00|
## |Z教育程度 | 0.05| 0.01| 27060.13| 7.84| 0.00|
## |Z社經地位 | -0.02| 0.01| 26963.89| -3.45| 0.00|
##
## Random Effects:
##
##
## |grp |var1 |var2 | vcov| sdcor|
## |:--------|:-----------|:--------|----:|-----:|
## |country |(Intercept) |NA | 0.04| 0.20|
## |country |Z人均GDP |NA | 0.01| 0.10|
## |country |(Intercept) |Z人均GDP | 0.00| -0.10|
## |Residual |NA |NA | 0.84| 0.92|
##
## ANOVA Result:
##
##
## | | npar| AIC| BIC| logLik| deviance| Chisq| Df| Pr(>Chisq)|
## |:-------------|----:|--------:|--------:|---------:|--------:|-----:|--:|----------:|
## |reduced_model | 8| 72474.79| 72540.46| -36229.39| 72458.79| NA| NA| NA|
## |full_model | 10| 72477.35| 72559.44| -36228.68| 72457.35| 1.44| 2| 0.49|##政府信任
full_model <- lmer(Zdirect_high ~ Z環境態度1 + Z政府信任 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度1 * Z政府信任 + (1 + Z人均GDP | country), data = data)
reduced_model <- lmer(Zdirect_high ~ Z環境態度1 + Z政府信任 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度1 * 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| 32.47| -0.55| 0.59|
## |Z環境態度1 | 0.29| 0.01| 27116.49| 46.99| 0.00|
## |Z政府信任 | 0.00| 0.01| 27035.45| 0.15| 0.88|
## |Z年齡 | 0.11| 0.01| 27127.79| 16.85| 0.00|
## |Z性別 | 0.07| 0.01| 27108.47| 11.61| 0.00|
## |Z教育程度 | 0.05| 0.01| 27058.45| 7.84| 0.00|
## |Z社經地位 | -0.02| 0.01| 26981.00| -3.44| 0.00|
## |Z環境態度1:Z政府信任 | 0.00| 0.01| 27115.62| -0.43| 0.67|
##
## Random Effects:
##
##
## |grp |var1 |var2 | vcov| sdcor|
## |:--------|:-----------|:--------|----:|-----:|
## |country |(Intercept) |NA | 0.04| 0.20|
## |country |Z人均GDP |NA | 0.01| 0.10|
## |country |(Intercept) |Z人均GDP | 0.00| -0.10|
## |Residual |NA |NA | 0.84| 0.92|
##
## ANOVA Result:
##
##
## | | npar| AIC| BIC| logLik| deviance| Chisq| Df| Pr(>Chisq)|
## |:-------------|----:|--------:|--------:|---------:|--------:|-----:|--:|----------:|
## |reduced_model | 10| 72478.59| 72560.68| -36229.30| 72458.59| NA| NA| NA|
## |full_model | 12| 72481.16| 72579.66| -36228.58| 72457.16| 1.43| 2| 0.49|##權力距離
full_model <- lmer(Zdirect_high ~ Z環境態度1 + Z權力距離 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度1 * Z權力距離 + (1 + Z人均GDP + Z權力距離 | country), data = data)
reduced_model <- lmer(Zdirect_high ~ Z環境態度1 + Z權力距離 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度1 * 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.03| 20.79| -0.69| 0.50|
## |Z環境態度1 | 0.28| 0.01| 26619.92| 46.51| 0.00|
## |Z權力距離 | 0.04| 0.04| 13.77| 0.97| 0.35|
## |Z年齡 | 0.11| 0.01| 26637.89| 16.62| 0.00|
## |Z性別 | 0.06| 0.01| 26647.15| 11.07| 0.00|
## |Z教育程度 | 0.05| 0.01| 26352.99| 7.60| 0.00|
## |Z社經地位 | -0.02| 0.01| 26241.89| -3.46| 0.00|
## |Z環境態度1:Z權力距離 | -0.04| 0.01| 26605.07| -6.83| 0.00|
##
## Random Effects:
##
##
## |grp |var1 |var2 | vcov| sdcor|
## |:--------|:-----------|:---------|-----:|-----:|
## |country |(Intercept) |NA | 0.01| 0.12|
## |country |Z人均GDP |NA | 0.03| 0.16|
## |country |Z權力距離 |NA | 0.04| 0.20|
## |country |(Intercept) |Z人均GDP | -0.01| -0.39|
## |country |(Intercept) |Z權力距離 | 0.00| -0.16|
## |country |Z人均GDP |Z權力距離 | 0.02| 0.75|
## |Residual |NA |NA | 0.84| 0.92|
##
## ANOVA Result:
##
##
## | | npar| AIC| BIC| logLik| deviance| Chisq| Df| Pr(>Chisq)|
## |:-------------|----:|--------:|--------:|---------:|--------:|-----:|--:|----------:|
## |reduced_model | 10| 71282.96| 71364.88| -35631.48| 71262.96| NA| NA| NA|
## |full_model | 15| 71285.80| 71408.67| -35627.90| 71255.80| 7.16| 5| 0.21|##個人主義
full_model <- lmer(Zdirect_high ~ Z環境態度1 + Z個人主義 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度1 * Z個人主義 + (1 + Z人均GDP + Z個人主義 | country), data = data)
reduced_model <- lmer(Zdirect_high ~ Z環境態度1 + Z個人主義 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度1 * 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| 12.41| -1.70| 0.11|
## |Z環境態度1 | 0.28| 0.01| 26949.41| 46.56| 0.00|
## |Z個人主義 | -0.04| 0.03| 8.61| -1.11| 0.30|
## |Z年齡 | 0.11| 0.01| 27013.65| 16.81| 0.00|
## |Z性別 | 0.06| 0.01| 27109.99| 11.32| 0.00|
## |Z教育程度 | 0.05| 0.01| 25987.09| 7.88| 0.00|
## |Z社經地位 | -0.02| 0.01| 27008.28| -3.34| 0.00|
## |Z環境態度1:Z個人主義 | 0.03| 0.01| 27091.48| 5.59| 0.00|
##
## Random Effects:
##
##
## |grp |var1 |var2 | vcov| sdcor|
## |:--------|:-----------|:---------|-----:|-----:|
## |country |(Intercept) |NA | 0.05| 0.23|
## |country |Z人均GDP |NA | 0.01| 0.12|
## |country |Z個人主義 |NA | 0.02| 0.13|
## |country |(Intercept) |Z人均GDP | 0.03| 1.00|
## |country |(Intercept) |Z個人主義 | -0.03| -1.00|
## |country |Z人均GDP |Z個人主義 | -0.01| -1.00|
## |Residual |NA |NA | 0.84| 0.92|
##
## ANOVA Result:
##
##
## | | npar| AIC| BIC| logLik| deviance| Chisq| Df| Pr(>Chisq)|
## |:-------------|----:|--------:|--------:|---------:|--------:|-----:|--:|----------:|
## |reduced_model | 10| 72446.94| 72529.03| -36213.47| 72426.94| NA| NA| NA|
## |full_model | 15| 72448.32| 72571.46| -36209.16| 72418.32| 8.62| 5| 0.13|##成就動機
full_model <- lmer(Zdirect_high ~ Z環境態度1 + Z成就動機 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度1 * Z成就動機 + (1 + Z人均GDP + Z成就動機 | country), data = data)
reduced_model <- lmer(Zdirect_high ~ Z環境態度1 + Z成就動機 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度1 * 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.00| 0.04| 27.02| -0.05| 0.96|
## |Z環境態度1 | 0.29| 0.01| 18162.59| 46.96| 0.00|
## |Z成就動機 | 0.10| 0.03| 17.01| 3.43| 0.00|
## |Z年齡 | 0.11| 0.01| 26140.34| 16.87| 0.00|
## |Z性別 | 0.06| 0.01| 26065.22| 11.19| 0.00|
## |Z教育程度 | 0.05| 0.01| 14958.93| 7.65| 0.00|
## |Z社經地位 | -0.02| 0.01| 14485.68| -3.41| 0.00|
## |Z環境態度1:Z成就動機 | -0.02| 0.01| 12200.66| -2.72| 0.01|
##
## Random Effects:
##
##
## |grp |var1 |var2 | vcov| sdcor|
## |:--------|:-----------|:---------|----:|-----:|
## |country |(Intercept) |NA | 0.04| 0.19|
## |country |Z人均GDP |NA | 0.00| 0.03|
## |country |Z成就動機 |NA | 0.02| 0.14|
## |country |(Intercept) |Z人均GDP | 0.00| -0.09|
## |country |(Intercept) |Z成就動機 | 0.03| 1.00|
## |country |Z人均GDP |Z成就動機 | 0.00| -0.01|
## |Residual |NA |NA | 0.84| 0.92|
##
## ANOVA Result:
##
##
## | | npar| AIC| BIC| logLik| deviance| Chisq| Df| Pr(>Chisq)|
## |:-------------|----:|--------:|--------:|---------:|--------:|-----:|--:|----------:|
## |reduced_model | 10| 71318.25| 71400.16| -35649.12| 71298.25| NA| NA| NA|
## |full_model | 15| 71310.14| 71433.02| -35640.07| 71280.14| 18.11| 5| 0|##不確定規避
full_model <- lmer(Zdirect_high ~ Z環境態度1 + Z不確定規避 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度1 * Z不確定規避 + (1 + Z人均GDP + Z不確定規避 | country), data = data)
reduced_model <- lmer(Zdirect_high ~ Z環境態度1 + Z不確定規避 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度1 * 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.03| 23.90| -0.95| 0.35|
## |Z環境態度1 | 0.29| 0.01| 26613.10| 46.95| 0.00|
## |Z不確定規避 | 0.06| 0.02| 2.96| 2.51| 0.09|
## |Z年齡 | 0.11| 0.01| 26661.55| 16.67| 0.00|
## |Z性別 | 0.06| 0.01| 26646.59| 11.17| 0.00|
## |Z教育程度 | 0.05| 0.01| 26359.69| 7.64| 0.00|
## |Z社經地位 | -0.02| 0.01| 26521.76| -3.45| 0.00|
## |Z環境態度1:Z不確定規避 | -0.03| 0.01| 26517.33| -4.98| 0.00|
##
## Random Effects:
##
##
## |grp |var1 |var2 | vcov| sdcor|
## |:--------|:-----------|:-----------|----:|-----:|
## |country |(Intercept) |NA | 0.03| 0.16|
## |country |Z人均GDP |NA | 0.02| 0.16|
## |country |Z不確定規避 |NA | 0.01| 0.12|
## |country |(Intercept) |Z人均GDP | 0.01| 0.32|
## |country |(Intercept) |Z不確定規避 | 0.01| 0.69|
## |country |Z人均GDP |Z不確定規避 | 0.02| 0.91|
## |Residual |NA |NA | 0.84| 0.92|
##
## ANOVA Result:
##
##
## | | npar| AIC| BIC| logLik| deviance| Chisq| Df| Pr(>Chisq)|
## |:-------------|----:|-------:|--------:|---------:|--------:|-----:|--:|----------:|
## |reduced_model | 10| 71306.3| 71388.21| -35643.15| 71286.3| NA| NA| NA|
## |full_model | 15| 71307.9| 71430.77| -35638.95| 71277.9| 8.4| 5| 0.14|##長期導向
full_model <- lmer(Zdirect_high ~ Z環境態度1 + Z長期導向 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度1 * Z長期導向 + (1 + Z人均GDP + Z長期導向 | country), data = data)
reduced_model <- lmer(Zdirect_high ~ Z環境態度1 + Z長期導向 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度1 * 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| 28.91| -0.60| 0.55|
## |Z環境態度1 | 0.29| 0.01| 26843.40| 47.29| 0.00|
## |Z長期導向 | 0.03| 0.03| 1.38| 0.76| 0.56|
## |Z年齡 | 0.11| 0.01| 27047.94| 16.86| 0.00|
## |Z性別 | 0.06| 0.01| 27104.29| 11.53| 0.00|
## |Z教育程度 | 0.05| 0.01| 27001.97| 7.84| 0.00|
## |Z社經地位 | -0.02| 0.01| 26722.67| -3.46| 0.00|
## |Z環境態度1:Z長期導向 | 0.01| 0.01| 26390.37| 1.86| 0.06|
##
## Random Effects:
##
##
## |grp |var1 |var2 | vcov| sdcor|
## |:--------|:-----------|:---------|-----:|-----:|
## |country |(Intercept) |NA | 0.04| 0.19|
## |country |Z人均GDP |NA | 0.01| 0.09|
## |country |Z長期導向 |NA | 0.01| 0.11|
## |country |(Intercept) |Z人均GDP | 0.00| 0.16|
## |country |(Intercept) |Z長期導向 | -0.02| -0.95|
## |country |Z人均GDP |Z長期導向 | 0.00| -0.47|
## |Residual |NA |NA | 0.84| 0.92|
##
## ANOVA Result:
##
##
## | | npar| AIC| BIC| logLik| deviance| Chisq| Df| Pr(>Chisq)|
## |:-------------|----:|--------:|--------:|---------:|--------:|-----:|--:|----------:|
## |reduced_model | 10| 72475.58| 72557.67| -36227.79| 72455.58| NA| NA| NA|
## |full_model | 15| 72478.98| 72602.11| -36224.49| 72448.98| 6.6| 5| 0.25|##放任約束
full_model <- lmer(Zdirect_high ~ Z環境態度1 + Z放任約束 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度1 * Z放任約束 + (1 + Z人均GDP + Z放任約束 | country), data = data)
reduced_model <- lmer(Zdirect_high ~ Z環境態度1 + Z放任約束 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度1 * 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.03| 19.13| -1.17| 0.26|
## |Z環境態度1 | 0.29| 0.01| 26482.48| 47.09| 0.00|
## |Z放任約束 | 0.06| 0.02| 23.78| 2.60| 0.02|
## |Z年齡 | 0.11| 0.01| 26423.79| 17.16| 0.00|
## |Z性別 | 0.07| 0.01| 26467.84| 11.64| 0.00|
## |Z教育程度 | 0.05| 0.01| 25745.71| 7.90| 0.00|
## |Z社經地位 | -0.02| 0.01| 25687.69| -3.02| 0.00|
## |Z環境態度1:Z放任約束 | 0.01| 0.01| 25479.11| 1.39| 0.16|
##
## Random Effects:
##
##
## |grp |var1 |var2 | vcov| sdcor|
## |:--------|:-----------|:---------|-----:|-----:|
## |country |(Intercept) |NA | 0.01| 0.09|
## |country |Z人均GDP |NA | 0.08| 0.29|
## |country |Z放任約束 |NA | 0.02| 0.15|
## |country |(Intercept) |Z人均GDP | 0.01| 0.40|
## |country |(Intercept) |Z放任約束 | 0.00| -0.04|
## |country |Z人均GDP |Z放任約束 | -0.04| -0.93|
## |Residual |NA |NA | 0.84| 0.92|
##
## ANOVA Result:
##
##
## | | npar| AIC| BIC| logLik| deviance| Chisq| Df| Pr(>Chisq)|
## |:-------------|----:|--------:|--------:|---------:|--------:|-----:|--:|----------:|
## |reduced_model | 10| 70747.54| 70829.39| -35363.77| 70727.54| NA| NA| NA|
## |full_model | 15| 70740.51| 70863.29| -35355.26| 70710.51| 17.03| 5| 0|##間接低成本行動 ##直接效果
full_model <- lmer(Zindirect_low ~ Z環境態度1 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 +
(1 + Z人均GDP | country), data = data)
reduced_model <- lmer(Zindirect_low ~ Z環境態度1 + 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.18| 0.02| 12.52| -10.31| 0.00|
## |Z環境態度1 | 0.20| 0.01| 26311.21| 32.28| 0.00|
## |Z年齡 | 0.01| 0.01| 26060.28| 0.99| 0.32|
## |Z性別 | 0.00| 0.01| 27006.25| -0.26| 0.79|
## |Z教育程度 | 0.14| 0.01| 15930.95| 21.94| 0.00|
## |Z社經地位 | 0.02| 0.01| 23576.39| 2.68| 0.01|
##
## Random Effects:
##
##
## |grp |var1 |var2 | vcov| sdcor|
## |:--------|:-----------|:--------|----:|-----:|
## |country |(Intercept) |NA | 0.03| 0.19|
## |country |Z人均GDP |NA | 0.01| 0.12|
## |country |(Intercept) |Z人均GDP | 0.02| 1.00|
## |Residual |NA |NA | 0.85| 0.92|
##
## ANOVA Result:
##
##
## | | npar| AIC| BIC| logLik| deviance| Chisq| Df| Pr(>Chisq)|
## |:-------------|----:|--------:|--------:|---------:|--------:|-----:|--:|----------:|
## |reduced_model | 8| 72684.91| 72750.58| -36334.46| 72668.91| NA| NA| NA|
## |full_model | 10| 72675.34| 72757.43| -36327.67| 72655.34| 13.58| 2| 0|##政府信任
full_model <- lmer(Zindirect_low ~ Z環境態度1 + Z政府信任 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度1 * Z政府信任 + (1 + Z人均GDP | country), data = data)
reduced_model <- lmer(Zindirect_low ~ Z環境態度1 + Z政府信任 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度1 * 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.18| 0.02| 12.19| -10.11| 0.00|
## |Z環境態度1 | 0.20| 0.01| 26393.87| 32.27| 0.00|
## |Z政府信任 | 0.01| 0.01| 23165.78| 1.47| 0.14|
## |Z年齡 | 0.01| 0.01| 26037.03| 0.93| 0.35|
## |Z性別 | 0.00| 0.01| 27002.79| -0.31| 0.76|
## |Z教育程度 | 0.14| 0.01| 16058.01| 21.85| 0.00|
## |Z社經地位 | 0.02| 0.01| 23902.34| 2.53| 0.01|
## |Z環境態度1:Z政府信任 | 0.01| 0.01| 26710.32| 2.19| 0.03|
##
## Random Effects:
##
##
## |grp |var1 |var2 | vcov| sdcor|
## |:--------|:-----------|:--------|----:|-----:|
## |country |(Intercept) |NA | 0.03| 0.18|
## |country |Z人均GDP |NA | 0.01| 0.12|
## |country |(Intercept) |Z人均GDP | 0.02| 1.00|
## |Residual |NA |NA | 0.85| 0.92|
##
## ANOVA Result:
##
##
## | | npar| AIC| BIC| logLik| deviance| Chisq| Df| Pr(>Chisq)|
## |:-------------|----:|--------:|--------:|---------:|--------:|-----:|--:|----------:|
## |reduced_model | 10| 72680.65| 72762.73| -36330.32| 72660.65| NA| NA| NA|
## |full_model | 12| 72671.67| 72770.18| -36323.84| 72647.67| 12.97| 2| 0|##權力距離
full_model <- lmer(Zindirect_low ~ Z環境態度1 + Z權力距離 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度1 * Z權力距離 + (1 + Z人均GDP + Z權力距離 | country), data = data)
reduced_model <- lmer(Zindirect_low ~ Z環境態度1 + Z權力距離 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度1 * 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.12| 0.03| 19.21| -4.58| 0.00|
## |Z環境態度1 | 0.19| 0.01| 25576.55| 31.66| 0.00|
## |Z權力距離 | -0.07| 0.02| 15.87| -3.52| 0.00|
## |Z年齡 | 0.00| 0.01| 26351.08| 0.69| 0.49|
## |Z性別 | 0.00| 0.01| 26636.05| -0.52| 0.61|
## |Z教育程度 | 0.14| 0.01| 23430.97| 21.63| 0.00|
## |Z社經地位 | 0.02| 0.01| 20748.02| 2.58| 0.01|
## |Z環境態度1:Z權力距離 | -0.06| 0.01| 23259.55| -9.74| 0.00|
##
## Random Effects:
##
##
## |grp |var1 |var2 | vcov| sdcor|
## |:--------|:-----------|:---------|----:|-----:|
## |country |(Intercept) |NA | 0.02| 0.13|
## |country |Z人均GDP |NA | 0.02| 0.15|
## |country |Z權力距離 |NA | 0.01| 0.09|
## |country |(Intercept) |Z人均GDP | 0.01| 0.62|
## |country |(Intercept) |Z權力距離 | 0.00| 0.12|
## |country |Z人均GDP |Z權力距離 | 0.01| 0.85|
## |Residual |NA |NA | 0.85| 0.92|
##
## ANOVA Result:
##
##
## | | npar| AIC| BIC| logLik| deviance| Chisq| Df| Pr(>Chisq)|
## |:-------------|----:|--------:|--------:|---------:|--------:|-----:|--:|----------:|
## |reduced_model | 10| 71345.72| 71427.64| -35662.86| 71325.72| NA| NA| NA|
## |full_model | 15| 71344.10| 71466.97| -35657.05| 71314.10| 11.62| 5| 0.04|##個人主義
full_model <- lmer(Zindirect_low ~ Z環境態度1 + Z個人主義 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度1 * Z個人主義 + (1 + Z人均GDP + Z個人主義 | country), data = data)
reduced_model <- lmer(Zindirect_low ~ Z環境態度1 + Z個人主義 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度1 * 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.07| 0.02| 5.11| -4.30| 0.01|
## |Z環境態度1 | 0.19| 0.01| 25246.66| 30.85| 0.00|
## |Z個人主義 | 0.10| 0.01| 5.90| 7.05| 0.00|
## |Z年齡 | 0.00| 0.01| 23554.38| 0.73| 0.46|
## |Z性別 | 0.00| 0.01| 27105.25| -0.88| 0.38|
## |Z教育程度 | 0.14| 0.01| 21480.62| 21.91| 0.00|
## |Z社經地位 | 0.02| 0.01| 11111.51| 2.87| 0.00|
## |Z環境態度1:Z個人主義 | 0.08| 0.01| 18329.18| 12.97| 0.00|
##
## Random Effects:
##
##
## |grp |var1 |var2 | vcov| sdcor|
## |:--------|:-----------|:---------|-----:|-----:|
## |country |(Intercept) |NA | 0.00| 0.06|
## |country |Z人均GDP |NA | 0.06| 0.24|
## |country |Z個人主義 |NA | 0.02| 0.16|
## |country |(Intercept) |Z人均GDP | 0.01| 0.72|
## |country |(Intercept) |Z個人主義 | -0.01| -0.63|
## |country |Z人均GDP |Z個人主義 | -0.04| -0.99|
## |Residual |NA |NA | 0.84| 0.92|
##
## ANOVA Result:
##
##
## | | npar| AIC| BIC| logLik| deviance| Chisq| Df| Pr(>Chisq)|
## |:-------------|----:|-------:|--------:|--------:|--------:|-----:|--:|----------:|
## |reduced_model | 10| 72506.2| 72588.28| -36243.1| 72486.2| NA| NA| NA|
## |full_model | 15| 72504.4| 72627.53| -36237.2| 72474.4| 11.8| 5| 0.04|##成就動機
full_model <- lmer(Zindirect_low ~ Z環境態度1 + Z成就動機 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度1 * Z成就動機 + (1 + Z人均GDP + Z成就動機 | country), data = data)
reduced_model <- lmer(Zindirect_low ~ Z環境態度1 + Z成就動機 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度1 * 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.01| 6.36| -2.84| 0.03|
## |Z環境態度1 | 0.20| 0.01| 23509.66| 32.24| 0.00|
## |Z成就動機 | 0.15| 0.02| 1.37| 8.54| 0.04|
## |Z年齡 | 0.01| 0.01| 26242.84| 1.00| 0.32|
## |Z性別 | 0.00| 0.01| 26563.03| -0.37| 0.71|
## |Z教育程度 | 0.14| 0.01| 15654.05| 21.70| 0.00|
## |Z社經地位 | 0.02| 0.01| 14590.71| 2.57| 0.01|
## |Z環境態度1:Z成就動機 | -0.02| 0.01| 13288.56| -3.82| 0.00|
##
## Random Effects:
##
##
## |grp |var1 |var2 | vcov| sdcor|
## |:--------|:-----------|:---------|-----:|-----:|
## |country |(Intercept) |NA | 0.00| 0.01|
## |country |Z人均GDP |NA | 0.06| 0.24|
## |country |Z成就動機 |NA | 0.02| 0.14|
## |country |(Intercept) |Z人均GDP | 0.00| -1.00|
## |country |(Intercept) |Z成就動機 | 0.00| 1.00|
## |country |Z人均GDP |Z成就動機 | -0.03| -1.00|
## |Residual |NA |NA | 0.85| 0.92|
##
## ANOVA Result:
##
##
## | | npar| AIC| BIC| logLik| deviance| Chisq| Df| Pr(>Chisq)|
## |:-------------|----:|--------:|--------:|---------:|--------:|-----:|--:|----------:|
## |reduced_model | 10| 71433.09| 71515.00| -35706.54| 71413.09| NA| NA| NA|
## |full_model | 15| 71432.62| 71555.49| -35701.31| 71402.62| 10.47| 5| 0.06|##不確定規避
full_model <- lmer(Zindirect_low ~ Z環境態度1 + Z不確定規避 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度1 * Z不確定規避 + (1 + Z人均GDP + Z不確定規避 | country), data = data)
reduced_model <- lmer(Zindirect_low ~ Z環境態度1 + Z不確定規避 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度1 * 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.03| 18.10| -3.20| 0.00|
## |Z環境態度1 | 0.20| 0.01| 25879.00| 32.27| 0.00|
## |Z不確定規避 | -0.08| 0.03| 21.71| -2.73| 0.01|
## |Z年齡 | 0.01| 0.01| 26211.14| 0.86| 0.39|
## |Z性別 | 0.00| 0.01| 26618.76| -0.38| 0.70|
## |Z教育程度 | 0.14| 0.01| 22569.44| 21.68| 0.00|
## |Z社經地位 | 0.02| 0.01| 25359.87| 2.48| 0.01|
## |Z環境態度1:Z不確定規避 | -0.04| 0.01| 26017.01| -5.85| 0.00|
##
## Random Effects:
##
##
## |grp |var1 |var2 | vcov| sdcor|
## |:--------|:-----------|:-----------|----:|-----:|
## |country |(Intercept) |NA | 0.02| 0.14|
## |country |Z人均GDP |NA | 0.01| 0.12|
## |country |Z不確定規避 |NA | 0.00| 0.06|
## |country |(Intercept) |Z人均GDP | 0.01| 0.44|
## |country |(Intercept) |Z不確定規避 | 0.00| -0.46|
## |country |Z人均GDP |Z不確定規避 | 0.00| 0.59|
## |Residual |NA |NA | 0.85| 0.92|
##
## ANOVA Result:
##
##
## | | npar| AIC| BIC| logLik| deviance| Chisq| Df| Pr(>Chisq)|
## |:-------------|----:|--------:|--------:|---------:|--------:|-----:|--:|----------:|
## |reduced_model | 10| 71406.89| 71488.81| -35693.45| 71386.89| NA| NA| NA|
## |full_model | 15| 71407.80| 71530.68| -35688.90| 71377.80| 9.09| 5| 0.11|##長期導向
full_model <- lmer(Zindirect_low ~ Z環境態度1 + Z長期導向 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度1 * Z長期導向 + (1 + Z人均GDP + Z長期導向 | country), data = data)
reduced_model <- lmer(Zindirect_low ~ Z環境態度1 + Z長期導向 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度1 * 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| 11.86| -1.10| 0.29|
## |Z環境態度1 | 0.20| 0.01| 26896.46| 32.15| 0.00|
## |Z長期導向 | 0.12| 0.03| 2.34| 3.45| 0.06|
## |Z年齡 | 0.01| 0.01| 23895.51| 1.01| 0.31|
## |Z性別 | 0.00| 0.01| 27094.58| -0.30| 0.77|
## |Z教育程度 | 0.14| 0.01| 25824.30| 21.90| 0.00|
## |Z社經地位 | 0.02| 0.01| 23448.38| 2.57| 0.01|
## |Z環境態度1:Z長期導向 | 0.01| 0.01| 26865.81| 2.21| 0.03|
##
## 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.15|
## |country |(Intercept) |Z人均GDP | 0.00| NaN|
## |country |(Intercept) |Z長期導向 | 0.00| NaN|
## |country |Z人均GDP |Z長期導向 | -0.03| -0.80|
## |Residual |NA |NA | 0.85| 0.92|
##
## ANOVA Result:
##
##
## | | npar| AIC| BIC| logLik| deviance| Chisq| Df| Pr(>Chisq)|
## |:-------------|----:|--------:|--------:|---------:|--------:|-----:|--:|----------:|
## |reduced_model | 10| 72684.14| 72766.23| -36332.07| 72664.14| NA| NA| NA|
## |full_model | 15| 72684.86| 72807.99| -36327.43| 72654.86| 9.28| 5| 0.1|##放任約束
full_model <- lmer(Zindirect_low ~ Z環境態度1 + Z放任約束 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度1 * Z放任約束 + (1 + Z人均GDP + Z放任約束 | country), data = data)
reduced_model <- lmer(Zindirect_low ~ Z環境態度1 + Z放任約束 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度1 * 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.03| 16.20| -1.25| 0.23|
## |Z環境態度1 | 0.20| 0.01| 26479.79| 31.77| 0.00|
## |Z放任約束 | 0.12| 0.02| 22.38| 5.62| 0.00|
## |Z年齡 | 0.01| 0.01| 26395.92| 1.02| 0.31|
## |Z性別 | 0.00| 0.01| 26469.30| -0.34| 0.74|
## |Z教育程度 | 0.14| 0.01| 24826.07| 21.30| 0.00|
## |Z社經地位 | 0.02| 0.01| 25714.45| 2.89| 0.00|
## |Z環境態度1:Z放任約束 | 0.04| 0.01| 25426.03| 7.10| 0.00|
##
## Random Effects:
##
##
## |grp |var1 |var2 | vcov| sdcor|
## |:--------|:-----------|:---------|-----:|-----:|
## |country |(Intercept) |NA | 0.01| 0.12|
## |country |Z人均GDP |NA | 0.02| 0.13|
## |country |Z放任約束 |NA | 0.01| 0.10|
## |country |(Intercept) |Z人均GDP | 0.00| -0.24|
## |country |(Intercept) |Z放任約束 | 0.01| 0.47|
## |country |Z人均GDP |Z放任約束 | -0.01| -0.97|
## |Residual |NA |NA | 0.85| 0.92|
##
## ANOVA Result:
##
##
## | | npar| AIC| BIC| logLik| deviance| Chisq| Df| Pr(>Chisq)|
## |:-------------|----:|--------:|--------:|---------:|--------:|-----:|--:|----------:|
## |reduced_model | 10| 70974.09| 71055.94| -35477.05| 70954.09| NA| NA| NA|
## |full_model | 15| 70974.70| 71097.48| -35472.35| 70944.70| 9.39| 5| 0.09|##間接高成本行動 ##直接效果
full_model <- lmer(Zindirect_high ~ Z環境態度1 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 +
(1 + Z人均GDP | country), data = data)
reduced_model <- lmer(Zindirect_high ~ Z環境態度1 + 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.02| 0.02| 33.52| -0.71| 0.48|
## |Z環境態度1 | 0.10| 0.01| 26501.68| 14.81| 0.00|
## |Z年齡 | -0.01| 0.01| 26758.31| -1.50| 0.13|
## |Z性別 | -0.02| 0.01| 27106.23| -3.09| 0.00|
## |Z教育程度 | 0.06| 0.01| 25384.96| 8.00| 0.00|
## |Z社經地位 | 0.00| 0.01| 25072.88| 0.58| 0.56|
##
## 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| 76413.03| 76478.7| -38198.51| 76397.03| NA| NA| NA|
## |full_model | 10| 76417.02| 76499.1| -38198.51| 76397.02| 0.01| 2| 1|##政府信任
full_model <- lmer(Zindirect_high ~ Z環境態度1 + Z政府信任 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度1 * Z政府信任 + (1 + Z人均GDP | country), data = data)
reduced_model <- lmer(Zindirect_high ~ Z環境態度1 + Z政府信任 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度1 * 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.40| -0.59| 0.56|
## |Z環境態度1 | 0.10| 0.01| 26531.60| 14.57| 0.00|
## |Z政府信任 | 0.01| 0.01| 24864.51| 1.04| 0.30|
## |Z年齡 | -0.01| 0.01| 26757.32| -1.51| 0.13|
## |Z性別 | -0.02| 0.01| 27103.62| -3.10| 0.00|
## |Z教育程度 | 0.06| 0.01| 25173.07| 8.00| 0.00|
## |Z社經地位 | 0.00| 0.01| 25231.40| 0.48| 0.63|
## |Z環境態度1:Z政府信任 | -0.01| 0.01| 27042.59| -1.19| 0.23|
##
## Random Effects:
##
##
## |grp |var1 |var2 | vcov| sdcor|
## |:--------|:-----------|:--------|----:|-----:|
## |country |(Intercept) |NA | 0.02| 0.13|
## |country |Z人均GDP |NA | 0.00| 0.01|
## |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| 76414.79| 76496.87| -38197.39| 76394.79| NA| NA| NA|
## |full_model | 12| 76418.76| 76517.26| -38197.38| 76394.76| 0.03| 2| 0.99|##權力距離
full_model <- lmer(Zindirect_high ~ Z環境態度1 + Z權力距離 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度1 * Z權力距離 + (1 + Z人均GDP + Z權力距離 | country), data = data)
reduced_model <- lmer(Zindirect_high ~ Z環境態度1 + Z權力距離 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度1 * 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.87| -1.38| 0.18|
## |Z環境態度1 | 0.09| 0.01| 24024.42| 14.46| 0.00|
## |Z權力距離 | 0.05| 0.02| 10.99| 2.01| 0.07|
## |Z年齡 | -0.01| 0.01| 26281.12| -1.45| 0.15|
## |Z性別 | -0.02| 0.01| 26643.89| -3.10| 0.00|
## |Z教育程度 | 0.05| 0.01| 15401.29| 7.70| 0.00|
## |Z社經地位 | 0.00| 0.01| 25854.89| 0.64| 0.52|
## |Z環境態度1:Z權力距離 | -0.01| 0.01| 24762.61| -1.65| 0.10|
##
## 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.67|
## |country |(Intercept) |Z權力距離 | 0.00| 0.43|
## |country |Z人均GDP |Z權力距離 | 0.00| 0.38|
## |Residual |NA |NA | 0.97| 0.99|
##
## ANOVA Result:
##
##
## | | npar| AIC| BIC| logLik| deviance| Chisq| Df| Pr(>Chisq)|
## |:-------------|----:|--------:|--------:|---------:|--------:|-----:|--:|----------:|
## |reduced_model | 10| 75009.60| 75091.52| -37494.80| 74989.60| NA| NA| NA|
## |full_model | 15| 75017.89| 75140.76| -37493.94| 74987.89| 1.72| 5| 0.89|##個人主義
full_model <- lmer(Zindirect_high ~ Z環境態度1 + Z個人主義 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度1 * Z個人主義 + (1 + Z人均GDP + Z個人主義 | country), data = data)
reduced_model <- lmer(Zindirect_high ~ Z環境態度1 + Z個人主義 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度1 * 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| 16.09| -1.81| 0.09|
## |Z環境態度1 | 0.09| 0.01| 26614.22| 14.43| 0.00|
## |Z個人主義 | -0.04| 0.03| 12.46| -1.20| 0.25|
## |Z年齡 | -0.01| 0.01| 26400.82| -1.51| 0.13|
## |Z性別 | -0.02| 0.01| 27104.65| -3.28| 0.00|
## |Z教育程度 | 0.06| 0.01| 21114.80| 7.96| 0.00|
## |Z社經地位 | 0.01| 0.01| 19696.41| 0.78| 0.44|
## |Z環境態度1:Z個人主義 | 0.03| 0.01| 27062.88| 4.13| 0.00|
##
## Random Effects:
##
##
## |grp |var1 |var2 | vcov| sdcor|
## |:--------|:-----------|:---------|-----:|-----:|
## |country |(Intercept) |NA | 0.00| 0.00|
## |country |Z人均GDP |NA | 0.07| 0.26|
## |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.09| -0.97|
## |Residual |NA |NA | 0.97| 0.99|
##
## ANOVA Result:
##
##
## | | npar| AIC| BIC| logLik| deviance| Chisq| Df| Pr(>Chisq)|
## |:-------------|----:|--------:|--------:|---------:|--------:|-----:|--:|----------:|
## |reduced_model | 10| 76398.18| 76480.27| -38189.09| 76378.18| NA| NA| NA|
## |full_model | 15| 76427.35| 76550.48| -38198.68| 76397.35| 0| 5| 1|##成就動機
full_model <- lmer(Zindirect_high ~ Z環境態度1 + Z成就動機 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度1 * Z成就動機 + (1 + Z人均GDP + Z成就動機 | country), data = data)
reduced_model <- lmer(Zindirect_high ~ Z環境態度1 + Z成就動機 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度1 * 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.00| 0.02| 31.87| -0.13| 0.89|
## |Z環境態度1 | 0.09| 0.01| 24027.60| 14.50| 0.00|
## |Z成就動機 | -0.04| 0.02| 12.30| -1.99| 0.07|
## |Z年齡 | -0.01| 0.01| 25511.89| -1.47| 0.14|
## |Z性別 | -0.02| 0.01| 26659.19| -3.08| 0.00|
## |Z教育程度 | 0.05| 0.01| 17685.68| 7.71| 0.00|
## |Z社經地位 | 0.00| 0.01| 20756.51| 0.58| 0.56|
## |Z環境態度1:Z成就動機 | -0.01| 0.01| 21848.02| -1.38| 0.17|
##
## Random Effects:
##
##
## |grp |var1 |var2 | vcov| sdcor|
## |:--------|:-----------|:---------|-----:|-----:|
## |country |(Intercept) |NA | 0.02| 0.13|
## |country |Z人均GDP |NA | 0.00| 0.01|
## |country |Z成就動機 |NA | 0.01| 0.08|
## |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.99|
##
## ANOVA Result:
##
##
## | | npar| AIC| BIC| logLik| deviance| Chisq| Df| Pr(>Chisq)|
## |:-------------|----:|--------:|--------:|--------:|--------:|-----:|--:|----------:|
## |reduced_model | 10| 75011.81| 75093.72| -37495.9| 74991.81| NA| NA| NA|
## |full_model | 15| 75018.00| 75140.88| -37494.0| 74988.00| 3.8| 5| 0.58|##不確定規避
full_model <- lmer(Zindirect_high ~ Z環境態度1 + Z不確定規避 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度1 * Z不確定規避 + (1 + Z人均GDP + Z不確定規避 | country), data = data)
reduced_model <- lmer(Zindirect_high ~ Z環境態度1 + Z不確定規避 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度1 * 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| 25.98| -1.31| 0.20|
## |Z環境態度1 | 0.09| 0.01| 24902.35| 14.47| 0.00|
## |Z不確定規避 | 0.05| 0.01| 63.69| 3.61| 0.00|
## |Z年齡 | -0.01| 0.01| 25709.21| -1.34| 0.18|
## |Z性別 | -0.02| 0.01| 26657.97| -2.99| 0.00|
## |Z教育程度 | 0.05| 0.01| 19366.16| 7.66| 0.00|
## |Z社經地位 | 0.00| 0.01| 24877.00| 0.65| 0.52|
## |Z環境態度1:Z不確定規避 | 0.02| 0.01| 22610.77| 2.74| 0.01|
##
## Random Effects:
##
##
## |grp |var1 |var2 | vcov| sdcor|
## |:--------|:-----------|:-----------|----:|-----:|
## |country |(Intercept) |NA | 0.01| 0.11|
## |country |Z人均GDP |NA | 0.00| 0.02|
## |country |Z不確定規避 |NA | 0.00| 0.04|
## |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.97| 0.99|
##
## ANOVA Result:
##
##
## | | npar| AIC| BIC| logLik| deviance| Chisq| Df| Pr(>Chisq)|
## |:-------------|----:|--------:|--------:|---------:|--------:|-----:|--:|----------:|
## |reduced_model | 10| 75002.63| 75084.54| -37491.31| 74982.63| NA| NA| NA|
## |full_model | 15| 75003.46| 75126.34| -37486.73| 74973.46| 9.17| 5| 0.1|##長期導向
full_model <- lmer(Zindirect_high ~ Z環境態度1 + Z長期導向 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度1 * Z長期導向 + (1 + Z人均GDP + Z長期導向 | country), data = data)
reduced_model <- lmer(Zindirect_high ~ Z環境態度1 + Z長期導向 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度1 * 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| 61.51| -1.80| 0.08|
## |Z環境態度1 | 0.10| 0.01| 26058.73| 14.82| 0.00|
## |Z長期導向 | -0.04| 0.06| 28.37| -0.60| 0.55|
## |Z年齡 | -0.01| 0.01| 26908.16| -1.45| 0.15|
## |Z性別 | -0.02| 0.01| 27114.54| -3.08| 0.00|
## |Z教育程度 | 0.06| 0.01| 25608.76| 8.11| 0.00|
## |Z社經地位 | 0.00| 0.01| 25788.90| 0.54| 0.59|
## |Z環境態度1:Z長期導向 | -0.01| 0.01| 27051.94| -1.44| 0.15|
##
## 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.33|
## |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| 76410.62| 76492.71| -38195.31| 76390.62| NA| NA| NA|
## |full_model | 15| 76451.43| 76574.56| -38210.72| 76421.43| 0| 5| 1|##放任約束
full_model <- lmer(Zindirect_high ~ Z環境態度1 + Z放任約束 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度1 * Z放任約束 + (1 + Z人均GDP + Z放任約束 | country), data = data)
reduced_model <- lmer(Zindirect_high ~ Z環境態度1 + Z放任約束 + Z年齡 + Z性別 + Z教育程度 + Z社經地位 + Z環境態度1 * 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| 23.06| -1.71| 0.10|
## |Z環境態度1 | 0.09| 0.01| 25782.83| 14.25| 0.00|
## |Z放任約束 | -0.03| 0.02| 26.14| -1.16| 0.25|
## |Z年齡 | -0.01| 0.01| 25542.75| -1.51| 0.13|
## |Z性別 | -0.02| 0.01| 26489.31| -3.28| 0.00|
## |Z教育程度 | 0.06| 0.01| 20563.35| 7.83| 0.00|
## |Z社經地位 | 0.00| 0.01| 23957.84| 0.41| 0.68|
## |Z環境態度1:Z放任約束 | 0.00| 0.01| 25850.85| -0.22| 0.82|
##
## 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.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| 74347.78| 74429.64| -37163.89| 74327.78| NA| NA| NA|
## |full_model | 15| 74356.76| 74479.54| -37163.38| 74326.76| 1.02| 5| 0.96|