wvs94a | R Documentation |
Between 1990 and 1993, 57,561 adults aged 18 and above from 42 nations were interviewed by local academic institutes in Eastern European nations and by professional survey organizations in other nations. The standardized mean difference (SMD) between males and females on life satisfaction and life control in each country were calculated as the effect sizes. Positive values indicate that males have higher scores than females do.
data(wvs94a)
The variables are:
Country
SMD on life satisfaction
SMD on life control
Sampling variance of lifesat
Sampling covariance between lifesat and lifecon
Sampling variance of lifecon
Gross National Product
World Values Study Group. (1994). World Values Survey, 1981-1984 and 1990-1993 [Computer file]. Ann Arbor, MI: Inter-university Consortium for Political and Social Research.
Au, K., & Cheung, M. W.-L. (2004). Intra-cultural variation and job autonomy in 42 countries. Organization Studies, 25, 1339-1362.
Cheung, M. W.-L. (2013). Multivariate meta-analysis as structural equation models. Structural Equation Modeling, 20, 429-454.
data(wvs94a)
## Random-effects model
random.ma1 <- meta(y=cbind(lifesat, lifecon),
v=cbind(lifesat_var, inter_cov, lifecon_var), data=wvs94a,
model.name="Random effects model")
summary(random.ma1)
## Random-effects model with both population effect sizes fixed at 0
random.ma2 <- meta(y=cbind(lifesat, lifecon),
v=cbind(lifesat_var, inter_cov, lifecon_var), data=wvs94a,
intercept.constraints=matrix(0, nrow=1, ncol=2),
model.name="Effect sizes are fixed at 0")
summary(random.ma2)
## Compare the nested models
anova(random.ma1, random.ma2)
## Fixed-effects model by fixing the variance component at 0
fixed.ma <- meta(y=cbind(lifesat, lifecon),
v=cbind(lifesat_var, inter_cov, lifecon_var), data=wvs94a,
RE.constraints=matrix(0, ncol=2, nrow=2),
model.name="Fixed effects model")
summary(fixed.ma)
## Mixed-effects model
## gnp is divided by 10000 and centered by using
## scale(gnp/10000, scale=FALSE)
mixed.ma1 <- meta(y=cbind(lifesat, lifecon),
v=cbind(lifesat_var, inter_cov, lifecon_var),
x=scale(gnp/10000, scale=FALSE), data=wvs94a,
model.name="GNP as a predictor")
summary(mixed.ma1)
## Mixed-effects model with equal regression coefficients
mixed.ma2 <- meta(y=cbind(lifesat, lifecon),
v=cbind(lifesat_var, inter_cov, lifecon_var),
x=scale(gnp/10000, scale=FALSE), data=wvs94a,
coef.constraints=matrix(c("0.0*Eq_slope",
"0.0*Eq_slope"), nrow=2),
model.name="GNP as a predictor with equal slope")
summary(mixed.ma2)
## Compare the nested models
anova(mixed.ma1, mixed.ma2)
## Plot the multivariate effect sizes
plot(random.ma1, main="Estimated effect sizes and their 95% confidence ellipses",
axis.label=c("Gender difference on life satisfaction",
"Gender difference on life control"))
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.