[R-sig-ME] Multiple binary responses per ID and time

Thierry Onkelinx thierry.onkelinx at inbo.be
Mon Mar 21 16:29:51 CET 2016

Dear Lize,

So the response for individual i at time i is (x_it, n_it). Does the
predictor has the same amount of trials as the responses (y_it, n_it)? If
so, do you have information on the n_it Bernouilli trials of both the
response and the predictor? If that is the case then you can model the
individual Bernoulli trials.

If you don't have the information at that detail, then you have to turn the
binomial predictor into a proportion. With a Bayesian hierarchical model
you can first model the predictor and then uses this modelled proportion as
a predictor for the response. There is an example in the INLA FAQ:

Best regards,


ir. Thierry Onkelinx
Instituut voor natuur- en bosonderzoek / Research Institute for Nature and
team Biometrie & Kwaliteitszorg / team Biometrics & Quality Assurance
Kliniekstraat 25
1070 Anderlecht

To call in the statistician after the experiment is done may be no more
than asking him to perform a post-mortem examination: he may be able to say
what the experiment died of. ~ Sir Ronald Aylmer Fisher
The plural of anecdote is not data. ~ Roger Brinner
The combination of some data and an aching desire for an answer does not
ensure that a reasonable answer can be extracted from a given body of data.
~ John Tukey

2016-03-21 15:33 GMT+01:00 Lize van der merwe <lizestats op gmail.com>:

> Dear List,
> Please advise.    I cannot get my head around modelling this data.
> Study involves 200 individuals with several (not always the same number)
> dichotomous outcomes, at 10 different times.  The predictor also has
> several
> (not the same as each other, nor the same as what the individal has at that
> time-point) dichotomous outcomes for the same individuals at the the same
> timepoints.  There are time-level covariates and also individual level
> covariates to include.
> How do I model these?   Not even sure how to lay out the data.
> Binomial pair (x,n) outcome,  for each individual and each time and another
> binomial pair for the predictor?
> Regards
> Lize van der Merwe
> _______________________________________________
> R-sig-mixed-models op r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models

	[[alternative HTML version deleted]]

More information about the R-sig-mixed-models mailing list