flourbeetle.RdThe flourbeetle data frame has 8 rows and 4 columns.
Two columns are explanatory, the other two are responses.
data(flourbeetle)This data frame contains the following columns:
log10 applied to CS2mgL.
a numeric vector, the concentration of gaseous carbon disulphide in mg per litre.
a numeric vector, counts; the number of beetles exposed to the poison.
a numeric vector, counts; the numbers killed.
These data were originally given in Table IV of Bliss (1935) and are the combination of two series of toxicological experiments involving Tribolium confusum, also known as the flour beetle. Groups of such adult beetles were exposed for 5 hours of gaseous carbon disulphide at different concentrations, and their mortality measured.
Bliss, C.I., 1935. The calculation of the dosage-mortality curve. Annals of Applied Biology, 22, 134–167.
fit1 <- vglm(cbind(killed, exposed - killed) ~ logdose,
binomialff(link = probitlink), flourbeetle, trace = TRUE)
#> Iteration 1: deviance = 10.13646
#> Iteration 2: deviance = 10.10674
#> Iteration 3: deviance = 10.10671
#> Iteration 4: deviance = 10.10671
summary(fit1)
#>
#> Call:
#> vglm(formula = cbind(killed, exposed - killed) ~ logdose, family = binomialff(link = probitlink),
#> data = flourbeetle, trace = TRUE)
#>
#> Coefficients:
#> Estimate Std. Error z value Pr(>|z|)
#> (Intercept) -34.939 2.648 -13.20 <2e-16 ***
#> logdose 19.730 1.487 13.27 <2e-16 ***
#> ---
#> Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
#>
#> Name of linear predictor: probitlink(prob)
#>
#> Residual deviance: 10.1067 on 6 degrees of freedom
#>
#> Log-likelihood: -18.1524 on 6 degrees of freedom
#>
#> Number of Fisher scoring iterations: 4
#>
#> Warning: Hauck-Donner effect detected in the following estimate(s):
#> 'logdose'
#>