BANOVA.Bernoulli {BANOVA} | R Documentation |

## Estimation of BANOVA with a Bernoulli dependent variable

### Description

`BANOVA.Bernoulli`

implements a Bayesian ANOVA for binary dependent variable, using a logit link and a normal heterogeneity distribution.

### Usage

```
BANOVA.Bernoulli(l1_formula = "NA", l2_formula = "NA", data,
id, l2_hyper = c(1, 1, 0.0001), burnin = 5000, sample = 2000, thin = 10,
adapt = 0, conv_speedup = F, jags = runjags.getOption('jagspath'))
## S3 method for class 'BANOVA.Bernoulli'
summary(object, ...)
## S3 method for class 'BANOVA.Bernoulli'
predict(object, newdata = NULL,...)
## S3 method for class 'BANOVA.Bernoulli'
print(x, ...)
```

### Arguments

`l1_formula` |
formula for level 1 e.g. 'Y~X1+X2' |

`l2_formula` |
formula for level 2 e.g. '~Z1+Z2', response variable must not be included |

`data` |
a data.frame in long format including all features in level 1 and level 2(covariates and categorical factors) and responses |

`id` |
subject ID of each response unit |

`l2_hyper` |
level 2 hyperparameters, c(a, b, |

`burnin` |
the number of burn in draws in the MCMC algorithm, default 5000 |

`sample` |
target samples in the MCMC algorithm after thinning, default 2000 |

`thin` |
the number of samples in the MCMC algorithm that needs to be thinned, default 10 |

`adapt` |
the number of adaptive iterations, default 0 (see run.jags) |

`conv_speedup` |
whether to speedup convergence, default F |

`jags` |
the system call or path for activating 'JAGS'. Default calls findjags() to attempt to locate 'JAGS' on your system |

`object` |
object of class |

`newdata` |
test data, either a matrix, vector or a data.frame. It must have the same format with the original data (the same number of features and the same data classes) |

`x` |
object of class |

`...` |
additional arguments,currently ignored |

### Details

Level 1 model:

`y_i`

~ `Binomial(1,p_i)`

, `p_i = logit^{-1}(\eta_i)`

where `\eta_i = \sum_{p = 0}^{P}\sum_{j=1}^{J_p}X_{i,j}^p\beta_{j,s_i}^p`

, `s_i`

is the subject id of data record `i`

. see `BANOVA-package`

### Value

`BANOVA.Bernoulli`

returns an object of class `"BANOVA.Bernoulli"`

. The returned object is a list containing:

`anova.table` |
table of effect sizes |

`coef.tables` |
table of estimated coefficients |

`pvalue.table` |
table of p-values |

`dMatrice` |
design matrices at level 1 and level 2 |

`samples_l2_param` |
posterior samples of level 2 parameters |

`data` |
original data.frame |

`mf1` |
model.frame of level 1 |

`mf2` |
model.frame of level 2 |

`JAGSmodel` |
'JAGS' model |

### Examples

```
data(bernlogtime)
# model with the dependent variable : response
res <- BANOVA.Bernoulli(response~typical, ~blur + color, bernlogtime,
bernlogtime$subject, burnin = 5000, sample = 2000, thin = 10)
summary(res)
```

*BANOVA*version 1.2.1 Index]