# \donttest{
require(gkwreg)
require(gkwdist)
data(ImpreciseTask)
# Example 1: Basic model with task effects
# Probability location varies by task type and uncertainty level
fit_kw <- gkwreg(location ~ task * difference,
data = ImpreciseTask,
family = "kw"
)
summary(fit_kw)
# Interpretation:
# - Alpha: Task type and uncertainty (difference) interact to affect
# probability estimates
# - Different tasks may have different baseline probability assessments
# Example 2: Heteroscedastic model
# Precision of estimates may vary by task and uncertainty
fit_kw_hetero <- gkwreg(
location ~ task * difference |
task + difference,
data = ImpreciseTask,
family = "kw"
)
summary(fit_kw_hetero)
# Interpretation:
# - Beta: Variability in estimates differs between tasks
# Higher uncertainty (difference) may lead to less precise estimates
# Example 3: McDonald distribution for extreme uncertainty
# Some participants may show very extreme probability assessments
fit_mc <- gkwreg(
location ~ task * difference | # gamma: full interaction
task * difference | # delta: full interaction
task, # lambda: task affects extremity
data = ImpreciseTask,
family = "mc",
control = gkw_control(
method = "BFGS",
maxit = 1500
)
)
summary(fit_mc)
# Interpretation:
# - Lambda varies by task: Weather vs. stock may produce
# different patterns of extreme probability assessments
# }
Run the code above in your browser using DataLab