CorrBin (version 1.4)

mc.est.CMData: Distribution of the number of responses assuming marginal compatibility.

Description

The mc.est function estimates the distribution of the number of responses in a cluster under the assumption of marginal compatibility: information from all cluster sizes is pooled. The estimation is performed independently for each treatment group.

Usage

## S3 method for class 'CMData':
mc.est(object, eps = 1e-06, ...)

mc.est(object, ...)

## S3 method for class 'CBData':
mc.est(object, ...)

Arguments

eps

numeric; EM iterations proceed until the sum of squared changes fall below eps

object

a CBData or CMData object

...

other potential arguments; not currently used

Value

• A data frame giving the estimated pdf for each treatment and clustersize. The probabilities add up to 1 for each Trt/ClusterSize combination. It has the following columns:
• Prob
  numeric, the probability of NResp responses in a cluster of size ClusterSize in group Trt
• Trt
  factor, the treatment group
• ClusterSize
  numeric, the cluster size
• NResp
  numeric, the number of responses

Details

The EM algorithm given by Stefanescu and Turnbull (2003) is used.

References

Stefanescu, C. & Turnbull, B. W. (2003) Likelihood inference for exchangeable binary data with varying cluster sizes. Biometrics, 59, 18-24

Examples

data(shelltox)
sh.mc <- mc.est(shelltox)

library(lattice)
xyplot(Prob~NResp|factor(ClusterSize), groups=Trt, data=sh.mc, subset=ClusterSize>0,
   type="l", as.table=TRUE, auto.key=list(columns=4, lines=TRUE, points=FALSE),
   xlab="Number of responses", ylab="Probability P(R=r|N=n)")