rm(list=ls())
library(lavaan)## Warning: package 'lavaan' was built under R version 3.1.3## This is lavaan 0.5-20
## lavaan is BETA software! Please report any bugs.setwd("~/Google Drive/work/cost benefit UG/new design with Max/analysis/data")
sub_temp = read.csv("SEM_moderated_by_EG_deciders.csv")Set model vars
Run model
for (i in 1 : length(Emotion_name)){
model1 = out_model[i]
model2 = med_model[i]
model3 = '# effets
indirect_dg_recipient := a*b
direct_dg_recipient := c
total_dg_recipient := direct_dg_recipient + indirect_dg_recipient
#### moderation by EG
mod_indirect_pi := a1pi*b
mod_direct_pi := c1pi
mod_total_pi := mod_direct_pi + mod_indirect_pi
mod_indirect_ig := a1ig*b
mod_direct_ig := c1ig
mod_total_ig := mod_direct_ig + mod_indirect_ig
mod_indirect_ug := a1ug*b
mod_direct_ug := c1ug
mod_total_ug := mod_direct_ug + mod_indirect_ug
### moderation by deciders
mod_indirect_pi_decider := a2pi*b
mod_direct_pi_decider := c2pi
mod_total_pi_decider := mod_direct_pi_decider + mod_indirect_pi_decider
mod_indirect_ig_decider := a2ig*b
mod_direct_ig_decider := c2ig
mod_total_ig_decider := mod_direct_ig_decider + mod_indirect_ig_decider
mod_indirect_ug_decider := a2ug*b
mod_direct_ug_decider := c2ug
mod_total_ug_decider := mod_direct_ug_decider + mod_indirect_ug_decider
'
model = paste (model1, model2, model3, sep = " ")
fit <- sem(model, data = sub_temp)
print(Emotion_name[i])
summary (fit)
}## [1] "happy"
## lavaan (0.5-20) converged normally after 117 iterations
##
## Number of observations 1008
##
## Estimator ML
## Minimum Function Test Statistic 0.000
## Degrees of freedom 0
##
## Parameter Estimates:
##
## Information Expected
## Standard Errors Standard
##
## Regressions:
## Estimate Std.Err Z-value P(>|z|)
## eWTR ~
## gift (c) 0.002 0.000 5.569 0.000
## gift_PI (c1pi) 0.000 0.000 0.788 0.431
## gift_IG (c1ig) 0.000 0.000 0.895 0.371
## gift_UG (c1ug) -0.000 0.000 -0.802 0.423
## gf_PI_D (c2pi) 0.000 0.000 0.154 0.878
## gf_IG_D (c2ig) -0.001 0.001 -1.357 0.175
## gf_UG_D (c2ug) 0.000 0.000 0.949 0.343
## happy (b) 0.054 0.006 8.644 0.000
## happy ~
## gift (a) 0.022 0.002 13.719 0.000
## gift_PI (a1pi) 0.000 0.002 0.096 0.924
## gift_IG (a1ig) 0.001 0.002 0.232 0.817
## gift_UG (a1ug) -0.005 0.002 -2.342 0.019
## gf_PI_D (a2pi) -0.004 0.003 -1.421 0.155
## gf_IG_D (a2ig) -0.004 0.003 -1.671 0.095
## gf_UG_D (a2ug) -0.001 0.002 -0.601 0.548
##
## Variances:
## Estimate Std.Err Z-value P(>|z|)
## eWTR 0.029 0.001 22.450 0.000
## happy 0.747 0.033 22.450 0.000
##
## Defined Parameters:
## Estimate Std.Err Z-value P(>|z|)
## indrct_dg_rcpn 0.001 0.000 7.313 0.000
## dirct_dg_rcpnt 0.002 0.000 5.569 0.000
## total_dg_rcpnt 0.003 0.000 9.457 0.000
## mod_indirect_p 0.000 0.000 0.096 0.924
## mod_direct_pi 0.000 0.000 0.788 0.431
## mod_total_pi 0.000 0.000 0.786 0.432
## mod_indirect_g 0.000 0.000 0.232 0.817
## mod_direct_ig 0.000 0.000 0.895 0.371
## mod_total_ig 0.000 0.000 0.924 0.355
## mod_indirect_g -0.000 0.000 -2.260 0.024
## mod_direct_ug -0.000 0.000 -0.802 0.423
## mod_total_ug -0.001 0.000 -1.391 0.164
## md_ndrct_p_dcd -0.000 0.000 -1.402 0.161
## md_drct_p_dcdr 0.000 0.000 0.154 0.878
## mod_ttl_p_dcdr -0.000 0.001 -0.225 0.822
## md_ndrct_g_dcd -0.000 0.000 -1.641 0.101
## md_drct_g_dcdr -0.001 0.001 -1.357 0.175
## mod_ttl_g_dcdr -0.001 0.001 -1.750 0.080
## md_ndrct_g_dcd -0.000 0.000 -0.599 0.549
## md_drct_g_dcdr 0.000 0.000 0.949 0.343
## mod_ttl_g_dcdr 0.000 0.000 0.758 0.449
##
## [1] "grateful"
## lavaan (0.5-20) converged normally after 116 iterations
##
## Number of observations 1008
##
## Estimator ML
## Minimum Function Test Statistic 0.000
## Degrees of freedom 0
##
## Parameter Estimates:
##
## Information Expected
## Standard Errors Standard
##
## Regressions:
## Estimate Std.Err Z-value P(>|z|)
## eWTR ~
## gift (c) 0.002 0.000 5.513 0.000
## gift_PI (c1pi) 0.000 0.000 0.557 0.577
## gift_IG (c1ig) 0.000 0.000 0.542 0.588
## gift_UG (c1ug) -0.000 0.000 -0.916 0.359
## gf_PI_D (c2pi) 0.000 0.000 0.613 0.540
## gf_IG_D (c2ig) -0.000 0.000 -0.697 0.486
## gf_UG_D (c2ug) 0.001 0.000 1.343 0.179
## gratefl (b) 0.059 0.006 9.769 0.000
## grateful ~
## gift (a) 0.021 0.002 12.834 0.000
## gift_PI (a1pi) 0.002 0.002 0.860 0.390
## gift_IG (a1ig) 0.003 0.002 1.379 0.168
## gift_UG (a1ug) -0.004 0.002 -1.746 0.081
## gf_PI_D (a2pi) -0.007 0.003 -2.765 0.006
## gf_IG_D (a2ig) -0.009 0.003 -3.670 0.000
## gf_UG_D (a2ug) -0.004 0.002 -1.795 0.073
##
## Variances:
## Estimate Std.Err Z-value P(>|z|)
## eWTR 0.028 0.001 22.450 0.000
## grateful 0.761 0.034 22.450 0.000
##
## Defined Parameters:
## Estimate Std.Err Z-value P(>|z|)
## indrct_dg_rcpn 0.001 0.000 7.773 0.000
## dirct_dg_rcpnt 0.002 0.000 5.513 0.000
## total_dg_rcpnt 0.003 0.000 9.457 0.000
## mod_indirect_p 0.000 0.000 0.857 0.392
## mod_direct_pi 0.000 0.000 0.557 0.577
## mod_total_pi 0.000 0.000 0.786 0.432
## mod_indirect_g 0.000 0.000 1.366 0.172
## mod_direct_ig 0.000 0.000 0.542 0.588
## mod_total_ig 0.000 0.000 0.924 0.355
## mod_indirect_g -0.000 0.000 -1.719 0.086
## mod_direct_ug -0.000 0.000 -0.916 0.359
## mod_total_ug -0.001 0.000 -1.391 0.164
## md_ndrct_p_dcd -0.000 0.000 -2.660 0.008
## md_drct_p_dcdr 0.000 0.000 0.613 0.540
## mod_ttl_p_dcdr -0.000 0.001 -0.225 0.822
## md_ndrct_g_dcd -0.001 0.000 -3.436 0.001
## md_drct_g_dcdr -0.000 0.000 -0.697 0.486
## mod_ttl_g_dcdr -0.001 0.001 -1.750 0.080
## md_ndrct_g_dcd -0.000 0.000 -1.765 0.078
## md_drct_g_dcdr 0.001 0.000 1.343 0.179
## mod_ttl_g_dcdr 0.000 0.000 0.758 0.449
##
## [1] "surprised"
## lavaan (0.5-20) converged normally after 111 iterations
##
## Number of observations 1008
##
## Estimator ML
## Minimum Function Test Statistic 0.000
## Degrees of freedom 0
##
## Parameter Estimates:
##
## Information Expected
## Standard Errors Standard
##
## Regressions:
## Estimate Std.Err Z-value P(>|z|)
## eWTR ~
## gift (c) 0.003 0.000 8.728 0.000
## gift_PI (c1pi) 0.000 0.000 0.724 0.469
## gift_IG (c1ig) 0.000 0.000 0.840 0.401
## gift_UG (c1ug) -0.001 0.000 -1.378 0.168
## gf_PI_D (c2pi) -0.000 0.001 -0.109 0.913
## gf_IG_D (c2ig) -0.001 0.001 -1.441 0.150
## gf_UG_D (c2ug) 0.000 0.000 0.960 0.337
## surprsd (b) 0.018 0.006 3.092 0.002
## surprised ~
## gift (a) 0.011 0.002 6.237 0.000
## gift_PI (a1pi) 0.002 0.002 0.674 0.500
## gift_IG (a1ig) 0.002 0.002 0.904 0.366
## gift_UG (a1ug) -0.000 0.002 -0.202 0.840
## gf_PI_D (a2pi) -0.003 0.003 -1.198 0.231
## gf_IG_D (a2ig) -0.009 0.003 -3.184 0.001
## gf_UG_D (a2ug) -0.005 0.003 -2.062 0.039
##
## Variances:
## Estimate Std.Err Z-value P(>|z|)
## eWTR 0.031 0.001 22.450 0.000
## surprised 0.920 0.041 22.450 0.000
##
## Defined Parameters:
## Estimate Std.Err Z-value P(>|z|)
## indrct_dg_rcpn 0.000 0.000 2.770 0.006
## dirct_dg_rcpnt 0.003 0.000 8.728 0.000
## total_dg_rcpnt 0.003 0.000 9.457 0.000
## mod_indirect_p 0.000 0.000 0.658 0.510
## mod_direct_pi 0.000 0.000 0.724 0.469
## mod_total_pi 0.000 0.000 0.786 0.432
## mod_indirect_g 0.000 0.000 0.868 0.385
## mod_direct_ig 0.000 0.000 0.840 0.401
## mod_total_ig 0.000 0.000 0.924 0.355
## mod_indirect_g -0.000 0.000 -0.202 0.840
## mod_direct_ug -0.001 0.000 -1.378 0.168
## mod_total_ug -0.001 0.000 -1.391 0.164
## md_ndrct_p_dcd -0.000 0.000 -1.117 0.264
## md_drct_p_dcdr -0.000 0.001 -0.109 0.913
## mod_ttl_p_dcdr -0.000 0.001 -0.225 0.822
## md_ndrct_g_dcd -0.000 0.000 -2.218 0.027
## md_drct_g_dcdr -0.001 0.001 -1.441 0.150
## mod_ttl_g_dcdr -0.001 0.001 -1.750 0.080
## md_ndrct_g_dcd -0.000 0.000 -1.716 0.086
## md_drct_g_dcdr 0.000 0.000 0.960 0.337
## mod_ttl_g_dcdr 0.000 0.000 0.758 0.449
##
## [1] "angry"
## lavaan (0.5-20) converged normally after 113 iterations
##
## Number of observations 1008
##
## Estimator ML
## Minimum Function Test Statistic 0.000
## Degrees of freedom 0
##
## Parameter Estimates:
##
## Information Expected
## Standard Errors Standard
##
## Regressions:
## Estimate Std.Err Z-value P(>|z|)
## eWTR ~
## gift (c) 0.003 0.000 7.417 0.000
## gift_PI (c1pi) 0.000 0.000 0.711 0.477
## gift_IG (c1ig) 0.000 0.000 0.923 0.356
## gift_UG (c1ug) -0.000 0.000 -1.161 0.246
## gf_PI_D (c2pi) -0.000 0.001 -0.024 0.981
## gf_IG_D (c2ig) -0.001 0.001 -1.658 0.097
## gf_UG_D (c2ug) 0.000 0.000 0.905 0.365
## angry (b) -0.038 0.006 -6.485 0.000
## angry ~
## gift (a) -0.016 0.002 -9.389 0.000
## gift_PI (a1pi) -0.001 0.002 -0.443 0.658
## gift_IG (a1ig) -0.000 0.002 -0.099 0.921
## gift_UG (a1ug) 0.003 0.002 1.262 0.207
## gf_PI_D (a2pi) 0.003 0.003 1.008 0.314
## gf_IG_D (a2ig) 0.002 0.003 0.624 0.533
## gf_UG_D (a2ug) 0.002 0.003 0.646 0.519
##
## Variances:
## Estimate Std.Err Z-value P(>|z|)
## eWTR 0.030 0.001 22.450 0.000
## angry 0.855 0.038 22.450 0.000
##
## Defined Parameters:
## Estimate Std.Err Z-value P(>|z|)
## indrct_dg_rcpn 0.001 0.000 5.336 0.000
## dirct_dg_rcpnt 0.003 0.000 7.417 0.000
## total_dg_rcpnt 0.003 0.000 9.457 0.000
## mod_indirect_p 0.000 0.000 0.442 0.659
## mod_direct_pi 0.000 0.000 0.711 0.477
## mod_total_pi 0.000 0.000 0.786 0.432
## mod_indirect_g 0.000 0.000 0.099 0.921
## mod_direct_ig 0.000 0.000 0.923 0.356
## mod_total_ig 0.000 0.000 0.924 0.355
## mod_indirect_g -0.000 0.000 -1.239 0.215
## mod_direct_ug -0.000 0.000 -1.161 0.246
## mod_total_ug -0.001 0.000 -1.391 0.164
## md_ndrct_p_dcd -0.000 0.000 -0.996 0.319
## md_drct_p_dcdr -0.000 0.001 -0.024 0.981
## mod_ttl_p_dcdr -0.000 0.001 -0.225 0.822
## md_ndrct_g_dcd -0.000 0.000 -0.621 0.535
## md_drct_g_dcdr -0.001 0.001 -1.658 0.097
## mod_ttl_g_dcdr -0.001 0.001 -1.750 0.080
## md_ndrct_g_dcd -0.000 0.000 -0.642 0.521
## md_drct_g_dcdr 0.000 0.000 0.905 0.365
## mod_ttl_g_dcdr 0.000 0.000 0.758 0.449
##
## [1] "annoyed"
## lavaan (0.5-20) converged normally after 124 iterations
##
## Number of observations 1008
##
## Estimator ML
## Minimum Function Test Statistic 0.000
## Degrees of freedom 0
##
## Parameter Estimates:
##
## Information Expected
## Standard Errors Standard
##
## Regressions:
## Estimate Std.Err Z-value P(>|z|)
## eWTR ~
## gift (c) 0.002 0.000 6.327 0.000
## gift_PI (c1pi) 0.000 0.000 0.795 0.427
## gift_IG (c1ig) 0.000 0.000 0.906 0.365
## gift_UG (c1ug) -0.000 0.000 -0.909 0.363
## gf_PI_D (c2pi) -0.000 0.000 -0.086 0.931
## gf_IG_D (c2ig) -0.001 0.001 -1.471 0.141
## gf_UG_D (c2ug) 0.000 0.000 0.727 0.467
## annoyed (b) -0.049 0.006 -8.222 0.000
## annoyed ~
## gift (a) -0.020 0.002 -11.689 0.000
## gift_PI (a1pi) -0.000 0.002 -0.065 0.948
## gift_IG (a1ig) -0.000 0.002 -0.189 0.850
## gift_UG (a1ug) 0.004 0.002 2.029 0.042
## gf_PI_D (a2pi) 0.001 0.003 0.564 0.573
## gf_IG_D (a2ig) 0.003 0.003 1.293 0.196
## gf_UG_D (a2ug) -0.001 0.002 -0.214 0.830
##
## Variances:
## Estimate Std.Err Z-value P(>|z|)
## eWTR 0.029 0.001 22.450 0.000
## annoyed 0.796 0.035 22.450 0.000
##
## Defined Parameters:
## Estimate Std.Err Z-value P(>|z|)
## indrct_dg_rcpn 0.001 0.000 6.725 0.000
## dirct_dg_rcpnt 0.002 0.000 6.327 0.000
## total_dg_rcpnt 0.003 0.000 9.457 0.000
## mod_indirect_p 0.000 0.000 0.065 0.948
## mod_direct_pi 0.000 0.000 0.795 0.427
## mod_total_pi 0.000 0.000 0.786 0.432
## mod_indirect_g 0.000 0.000 0.189 0.850
## mod_direct_ig 0.000 0.000 0.906 0.365
## mod_total_ig 0.000 0.000 0.924 0.355
## mod_indirect_g -0.000 0.000 -1.970 0.049
## mod_direct_ug -0.000 0.000 -0.909 0.363
## mod_total_ug -0.001 0.000 -1.391 0.164
## md_ndrct_p_dcd -0.000 0.000 -0.563 0.574
## md_drct_p_dcdr -0.000 0.000 -0.086 0.931
## mod_ttl_p_dcdr -0.000 0.001 -0.225 0.822
## md_ndrct_g_dcd -0.000 0.000 -1.278 0.201
## md_drct_g_dcdr -0.001 0.001 -1.471 0.141
## mod_ttl_g_dcdr -0.001 0.001 -1.750 0.080
## md_ndrct_g_dcd 0.000 0.000 0.214 0.831
## md_drct_g_dcdr 0.000 0.000 0.727 0.467
## mod_ttl_g_dcdr 0.000 0.000 0.758 0.449
##
## [1] "disgust"
## lavaan (0.5-20) converged normally after 113 iterations
##
## Number of observations 1008
##
## Estimator ML
## Minimum Function Test Statistic 0.000
## Degrees of freedom 0
##
## Parameter Estimates:
##
## Information Expected
## Standard Errors Standard
##
## Regressions:
## Estimate Std.Err Z-value P(>|z|)
## eWTR ~
## gift (c) 0.002 0.000 7.170 0.000
## gift_PI (c1pi) 0.000 0.000 0.819 0.413
## gift_IG (c1ig) 0.000 0.000 0.962 0.336
## gift_UG (c1ug) -0.000 0.000 -1.045 0.296
## gf_PI_D (c2pi) -0.000 0.001 -0.137 0.891
## gf_IG_D (c2ig) -0.001 0.001 -1.633 0.103
## gf_UG_D (c2ug) 0.000 0.000 0.781 0.435
## disgust (b) -0.039 0.006 -6.493 0.000
## disgust ~
## gift (a) -0.018 0.002 -10.332 0.000
## gift_PI (a1pi) 0.000 0.002 0.082 0.935
## gift_IG (a1ig) 0.000 0.002 0.092 0.926
## gift_UG (a1ug) 0.004 0.002 1.824 0.068
## gf_PI_D (a2pi) 0.001 0.003 0.454 0.650
## gf_IG_D (a2ig) 0.002 0.003 0.747 0.455
## gf_UG_D (a2ug) 0.000 0.003 0.036 0.971
##
## Variances:
## Estimate Std.Err Z-value P(>|z|)
## eWTR 0.030 0.001 22.450 0.000
## disgust 0.836 0.037 22.450 0.000
##
## Defined Parameters:
## Estimate Std.Err Z-value P(>|z|)
## indrct_dg_rcpn 0.001 0.000 5.498 0.000
## dirct_dg_rcpnt 0.002 0.000 7.170 0.000
## total_dg_rcpnt 0.003 0.000 9.457 0.000
## mod_indirect_p -0.000 0.000 -0.082 0.935
## mod_direct_pi 0.000 0.000 0.819 0.413
## mod_total_pi 0.000 0.000 0.786 0.432
## mod_indirect_g -0.000 0.000 -0.092 0.926
## mod_direct_ig 0.000 0.000 0.962 0.336
## mod_total_ig 0.000 0.000 0.924 0.355
## mod_indirect_g -0.000 0.000 -1.756 0.079
## mod_direct_ug -0.000 0.000 -1.045 0.296
## mod_total_ug -0.001 0.000 -1.391 0.164
## md_ndrct_p_dcd -0.000 0.000 -0.453 0.650
## md_drct_p_dcdr -0.000 0.001 -0.137 0.891
## mod_ttl_p_dcdr -0.000 0.001 -0.225 0.822
## md_ndrct_g_dcd -0.000 0.000 -0.742 0.458
## md_drct_g_dcdr -0.001 0.001 -1.633 0.103
## mod_ttl_g_dcdr -0.001 0.001 -1.750 0.080
## md_ndrct_g_dcd -0.000 0.000 -0.036 0.971
## md_drct_g_dcdr 0.000 0.000 0.781 0.435
## mod_ttl_g_dcdr 0.000 0.000 0.758 0.449