biv.binom {repeated} R Documentation

## Marginal Bivariate Binomial Regression Models

### Description

`biv.binom` fits (logit) linear regression models to a marginal bivariate binomial distribution. The covariates must be of length K, that is the number of 2x2 tables.

### Usage

```biv.binom(freq, marg1=~1, marg2=~1, interaction=~1, pmarg1=1,
pmarg2=1,pinteraction=1, print.level=0, typsiz=abs(p),
ndigit=10, gradtol=0.00001, stepmax=10*sqrt(p%*%p),
steptol=0.00001, iterlim=100, fscale=1)
```

### Arguments

 `freq` A four-column matrix containing K 2x2 frequency tables. `marg1` The model formula for the first margin. `marg2` The model formula for the second margin. `interaction` The model formula for the interaction. `pmarg1` Initial parameter estimates for the first margin regression. `pmarg2` Initial parameter estimates for the second margin regression. `pinteraction` Initial parameter estimates for the interaction regression. `other` Arguments for nlm.

### Value

A list of class `bivbinom` is returned.

J.K. Lindsey

### Examples

```# 5 2x2 tables
Freq <- matrix(rpois(20,10),ncol=4)
x <- c(6,8,10,12,14)
print(z <- biv.binom(Freq,marg1=~x,marg2=~x,inter=~x,pmarg1=c(-2,0.08),
pmarg2=c(-2,0.1),pinter=c(3,0)))
```

[Package repeated version 1.0 Index]