# \donttest{
require(gkwreg)
require(gkwdist)
data(LossAversion)
# Control bounds
LossAversion$invest <- with(
LossAversion,
ifelse(invest <= 0, 0.000001,
ifelse(invest >= 1, 0.999999, invest)
)
)
# Example 1: Test for myopic loss aversion
# Do short-term players invest less? (They shouldn't, per Sutter et al.)
fit_kw <- gkwreg(
invest ~ arrangement + age + male + grade |
arrangement + male,
data = LossAversion,
family = "kw"
)
summary(fit_kw)
# Interpretation:
# - Alpha: Effect of investment horizon (arrangement) on mean investment
# Age and gender effects on risk-taking
# - Beta: Precision varies by horizon and gender
# (some groups more consistent than others)
# Example 2: Interaction effects
# Does the horizon effect differ by age/grade?
fit_kw_interact <- gkwreg(
invest ~ grade * (arrangement + age) + male |
arrangement + male + grade,
data = LossAversion,
family = "kw"
)
summary(fit_kw_interact)
# Interpretation:
# - Grade × arrangement interaction tests if myopic loss aversion
# emerges differently at different developmental stages
# Example 3: Extended-support for boundary observations
# Some students invest 0% or 100% of tokens
# Original 'invest' variable may include exact 0 and 1 values
fit_xbx <- gkwreg(
invest ~ grade * (arrangement + age) + male |
arrangement + male + grade,
data = LossAversion,
family = "kw" # Note: for true [0,1] support, use extended-support models
)
summary(fit_xbx)
# Interpretation:
# - Model accommodates extreme risk-taking (all-in or all-out strategies)
# Compare models
anova(fit_kw, fit_kw_interact)
# Visualization: Investment by horizon
boxplot(invest ~ arrangement,
data = LossAversion,
xlab = "Investment Horizon", ylab = "Proportion Invested",
main = "No Myopic Loss Aversion in Adolescents",
col = c("lightblue", "lightgreen", "lightyellow")
)
# }
Run the code above in your browser using DataLab