splitbin.RdThe function splits grouped data and optional covariates into individual data.
Two types of grouped data are managed by splitbin:
Grouped data with weights;
Grouped data of form cbind(success, failure).
When weights, successes or failures involve non-integer numbers, these numbers are rounded before splitting.
splitbin(formula, data, id = "idbin")A formula. The left-hand side represents grouped data. The right-hand side defines the covariates. See examples for syntax.
A data frame where all the variables described in formula are found.
An optional character string naming the identifier (= grouping factor). Default to “idbin”.
A data frame built according to the formula and function used in the call.
# Grouped data with weights
mydata <- data.frame(
success = c(0, 1, 0, 1),
f1 = c("A", "A", "B", "B"),
f2 = c("C", "D", "C", "D"),
n = c(4, 2, 1, 3)
)
mydata
#> success f1 f2 n
#> 1 0 A C 4
#> 2 1 A D 2
#> 3 0 B C 1
#> 4 1 B D 3
splitbin(formula = n ~ f1, data = mydata)$tab
#> idbin f1
#> 1 1 A
#> 2 1 A
#> 3 1 A
#> 4 1 A
#> 5 2 A
#> 6 2 A
#> 7 3 B
#> 8 4 B
#> 9 4 B
#> 10 4 B
splitbin(formula = n ~ f1 + f2 + success , data = mydata)$tab
#> idbin f1 f2 success
#> 1 1 A C 0
#> 2 1 A C 0
#> 3 1 A C 0
#> 4 1 A C 0
#> 5 2 A D 1
#> 6 2 A D 1
#> 7 3 B C 0
#> 8 4 B D 1
#> 9 4 B D 1
#> 10 4 B D 1
# Grouped data of form "cbind(success, failure)"
mydata <- data.frame(
success = c(4, 1),
failure = c(1, 2),
f1 = c("A", "B"),
f2 = c("C", "D")
)
mydata$n <- mydata$success + mydata$failure
mydata
#> success failure f1 f2 n
#> 1 4 1 A C 5
#> 2 1 2 B D 3
splitbin(formula = cbind(success, failure) ~ 1, data = mydata)$tab
#> idbin success
#> 1 1 0
#> 2 1 1
#> 3 1 1
#> 4 1 1
#> 5 1 1
#> 6 2 0
#> 7 2 0
#> 8 2 1
splitbin(formula = cbind(success, failure) ~ f1 + f2, data = mydata)$tab
#> idbin success f1 f2
#> 1 1 0 A C
#> 2 1 1 A C
#> 3 1 1 A C
#> 4 1 1 A C
#> 5 1 1 A C
#> 6 2 0 B D
#> 7 2 0 B D
#> 8 2 1 B D
splitbin(formula = cbind(success, n - success) ~ f1 + f2, data = mydata)$tab
#> idbin success f1 f2
#> 1 1 0 A C
#> 2 1 1 A C
#> 3 1 1 A C
#> 4 1 1 A C
#> 5 1 1 A C
#> 6 2 0 B D
#> 7 2 0 B D
#> 8 2 1 B D
splitbin(formula = cbind(success, n - 0.5 * failure - success) ~ f1 + f2,
data = mydata)$tab
#> idbin success f1 f2
#> 1 1 1 A C
#> 2 1 1 A C
#> 3 1 1 A C
#> 4 1 1 A C
#> 5 2 0 B D
#> 6 2 1 B D