# \donttest{
require(gkwreg)
require(gkwdist)
data(WeatherTask)
# Example 1: Main effects model
# Probability judgments affected by priming and elicitation format
fit_kw <- gkwreg(
agreement ~ priming + eliciting,
data = WeatherTask,
family = "kw"
)
summary(fit_kw)
# Interpretation:
# - Alpha: Seven-fold priming may shift probability estimates
# Imprecise elicitation may produce different mean estimates
# Example 2: Interaction model with heteroscedasticity
# Priming effects may differ by elicitation format
# Variability may also depend on conditions
fit_kw_interact <- gkwreg(
agreement ~ priming * eliciting |
priming + eliciting,
data = WeatherTask,
family = "kw"
)
summary(fit_kw_interact)
# Interpretation:
# - Alpha: Interaction tests if partition priming works differently
# for precise vs. imprecise probability judgments
# - Beta: Precision varies by experimental condition
# Test interaction
anova(fit_kw, fit_kw_interact)
# Example 3: McDonald distribution for polarized responses
# Probability judgments often show polarization (clustering at extremes)
# particularly under certain priming conditions
fit_mc <- gkwreg(
agreement ~ priming * eliciting | # gamma
priming * eliciting | # delta
priming, # lambda: priming affects polarization
data = WeatherTask,
family = "mc",
control = gkw_control(method = "BFGS", maxit = 1500)
)
summary(fit_mc)
# Interpretation:
# - Lambda varies by priming: Seven-fold priming may produce more
# extreme/polarized probability judgments
# }
Run the code above in your browser using DataLab