Description Usage Arguments Value Author(s) Examples
The empirical logit transform allows for a tolerence such that infinity is not returned when the argument is zero or one.
1 | empLogit(x, eps = 1e-3)
|
x |
numerical vector for which the empirical logit transform is desired |
eps |
numerical scalar; a tolerence to prevent infinite values |
the empirical logit transform of x
Rune Haubo B Christensen
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 | ## The function is currently defined as
## function (x, eps = 1e-3) log((eps + x)/(1 - x + eps))
## Lifted from example(predict.glm):
ldose <- rep(0:5, 2)
numdead <- c(1, 4, 9, 13, 18, 20, 0, 2, 6, 10, 12, 16)
sex <- factor(rep(c("M", "F"), c(6, 6)))
SF <- cbind(numdead, numalive=20-numdead)
## budworm.lg <- glm(SF ~ sex*ldose, family=binomial)
## summary(budworm.lg)
empLogit(numdead/20)
## Possible usage:
## Explorative interaction plot:
interaction.plot(ldose, sex, empLogit(numdead/20))
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.