Add two cautions to the `simulate_dyad_responses()` help

Open Beginner friendly
#179 0 comments 0 reactions 0 assignees View on GitHub

Maintainers usually reply within 1 day

@asimht is already working on this.

Since Nov 13, 2026.

  • #185 by @asimht — open

Assessment

Difficulty
2/5
Estimated time
1-3 hours
Newbie friendliness
70/100
Issue type
Documentation
Clarity
Clearly specified
Activity status
Active
Tech stack
r
Domain
documentation

Research direction

Start by locating the help entry for simulate_dyad_responses() and reading its existing cautions about dyad variance and lagged responses. Update the wording to cover zero-inflated and strongly overdispersed count models, the partner-check limitation, and the validity of model-centred summaries with lagged outcomes. Done when the help includes both cautions and the package documentation builds.

Written by the indexing model from the issue text.

Description

documentation priority: low

Two cautions from the review. I have not repeated either simulation.

The underestimated dyad variance is not limited to ordinal models. The help says glmmTMB tends to underestimate the dyad variance in ordinal models, because each dyad has only two people. The review found the same for zero-inflated Poisson models (estimated dyad SD 0.85 for a true 1.0, against 0.83 for a five-category ordinal model). A negative binomial model with strong overdispersion gave 0.91. A correct ordinal model was flagged by the partner check with 3,000 dyads, but not once in 60 runs with 500 dyads.

With a lagged outcome as predictor, raw-scale summaries flag a correct model. The help says lagged responses stay fixed. It does not say what follows. Simulations start from the observed lag with a newly drawn person effect, so raw simulated responses vary too little. In the review, the raw response SD was flagged in 30 of 30 runs with the true parameters. Model-centred summaries stayed near 5%.

Suggested fix and checks

  • Widen the ordinal sentence to zero-inflated and strongly overdispersed count models, and say that the partner check is affected only with very many dyads.
  • Say that with lagged outcomes only model-centred summaries are valid.
Dominant language
R
Stars
1
Forks
2
Avg merge
11h 4m
Merged PRs (30d)
42

Getting set up

First steps

  1. Read the whole issue, then the project's contributing guide.
  2. Comment on the issue to say you are picking it up — it saves two people doing the same work.
  3. Fork the repository and make your change on a branch.
  4. Open a pull request that references the issue number.

More from Pascal-Kueng/dyadMLM

All issues in Pascal-Kueng/dyadMLM

Similar issues

More R issues

Get new issues in your inbox

A short digest of beginner-friendly GitHub issues.