Running examples in 'pairwiseCI-Ex.R' failed
The error most likely occurred in:

> ### Name: overdispersed.binomial.ratio
> ### Title: Confidence intervals for risk ratios of overdispersed binomial
> ###   data
> ### Aliases: Quasibin.ratio Betabin.ratio ODbin.ratio
> ### Keywords: htest
> 
> ### ** Examples
> 
> 
> # Toxicologoical data: Number of Pups alive four days after birth (16 litters of rats)
> # Original source: Weil, 1970: Selection of valid number[...].
> #  Food and Cosmetics, Toxicology, 8, 177-182.
> # Cited from Zaihra and Paul(2010): Interval Estimation  of Some Epidemiological
> # Measures of Association. Int. J Biostatistics, 6(1), Article 35.
> 
> 
> mchem = c(12, 11, 10, 9, 11, 10, 10, 9, 9, 5, 9, 7, 10, 6, 10, 7)
> xchem = c(12, 11, 10, 9, 10,  9,  9, 8, 8, 4, 7, 4,  5, 3,  3, 0)
> dchem <- cbind("alive"=xchem, "dead"=mchem-xchem)
> 
> mcon = c(13, 12, 9, 9, 8, 8, 13, 12, 10, 10, 9, 13, 5, 7, 10, 10)
> xcon = c(13, 12, 9, 9, 8, 8, 12, 11,  9,  9, 8, 11, 4, 5,  7,  7)
> dcon <- cbind("alive"=xcon, "dead"=mcon-xcon)
> 
> 
> # Zaihra and Paul report: MR2: [0.714; 1.034]
> ODbin.ratio(x=dchem, y=dcon, CImethod="LOD", resbin=FALSE) 
$conf.int
    lower     upper 
0.7141806 1.0342595 
attr(,"methodname")
[1] "Delta-method, log-scale, overdispersed binomial, residual variance."

$estimate
[1] 0.8594463

$estOD
[1] 0.0047984677 0.0007118983

> 
> # Zaihra and Paul report: MR4: [0.710; 1.029]
> ODbin.ratio(x=dchem, y=dcon, CImethod="FOD", resbin=FALSE) 
$conf.int
    lower     upper 
0.7107622 1.0297294 
attr(,"methodname")
[1] "Fieller-Bailey-method, overdispersed binomial, residual variance."

$estimate
[1] 0.8594463

$estOD
[1] 0.0047984677 0.0007118983

> 
> 
> 
> Betabin.ratio(x=dchem, y=dcon, CImethod="FBB", iccpool=TRUE, resbin=TRUE) 
$conf.int
    lower     upper 
0.7248458 1.0146844 
attr(,"methodname")
[1] "Fieller-Bailey-method, betabinomial, icc estimated pooled over samples, variance limited to binomial variance."

$estimate
[1] 0.8594463

$esticc
[1] 0.1932541 0.1932541

> 
> Quasibin.ratio(x=dchem, y=dcon)
Error in `[<-.data.frame`(`*tmp*`, is.na(oci[, 1]), 1, value = 0) : 
  replacement has 1 row, data has 0
Calls: Quasibin.ratio -> QBmover -> [<- -> [<-.data.frame
Execution halted