Minimisation for a dataset
minimise.Rdminimise randomises patients in a data frame using minimisation
Usage
minimise(
data,
groups = 3,
factors,
burnin = 10,
minprob = 0.8,
ratio = rep(1, groups),
group.names = NULL,
seed = NULL
)Arguments
- data
a
data.frameobject with one line per participant and columns for the minimisation factors.- groups
an integer, the number of groups to randomise to, default is 3
- factors
a character vector with the factors for minimisation
- burnin
an integer, the burnin length before minimisation kicks in, default is 10. Must be > 0 and < total sample size.
- minprob
the minimisation probability. The default is to give 0.8 probability to the group which would lead to the least imbalance.
- ratio
a numeric vector of randomisation ratios in ascending order (must be of length equal to the number of groups)
- group.names
optional, a character vector with the group names, must be the same length as
groups.- seed
optional, an integer that is used with
set.seed()for offsetting the random number generator.
Value
(Invisibly) the data.frame with an additional column Group
indicating numerically which group has been allocated.
Examples
# Minimisation to 3 groups, with two factors and a burnin of 15, using the
# patients data from the package
(mini <- minimise(patients, groups = 3, factors = c("sex", "stage"),
burnin = 15))
#> Multi-arm Minimisation
#> --------------------------------------------------------------------------------
#> Groups: A B C
#> Randomisation ratio: 1:1:1
#> Factors: sex, stage
#> Burnin: 15
#> Minimisation probability: 0.8
#> Group sizes: 51, 50, 49
# View data with group info
as.data.frame(mini)
#> sex stage site Group
#> 1 F I 2 A
#> 2 F II 2 B
#> 3 M I 9 A
#> 4 F I 4 B
#> 5 M II 7 C
#> 6 F III 4 B
#> 7 M II 8 C
#> 8 M I 5 B
#> 9 F I 10 C
#> 10 F II 7 A
#> 11 M I 4 A
#> 12 F I 6 A
#> 13 F I 6 A
#> 14 M I 2 C
#> 15 F II 7 A
#> 16 M III 2 C
#> 17 F III 3 A
#> 18 M II 3 B
#> 19 M I 8 B
#> 20 M III 10 A
#> 21 M II 7 B
#> 22 F II 4 C
#> 23 M II 5 A
#> 24 M I 10 C
#> 25 F II 1 C
#> 26 M III 10 B
#> 27 F I 10 C
#> 28 F I 9 B
#> 29 M II 3 A
#> 30 F II 2 B
#> 31 F II 3 C
#> 32 M III 9 A
#> 33 F I 7 C
#> 34 F III 1 A
#> 35 F I 7 B
#> 36 M I 9 B
#> 37 M I 2 C
#> 38 M III 7 C
#> 39 F I 7 A
#> 40 M I 6 A
#> 41 M III 1 C
#> 42 M I 6 B
#> 43 M III 3 B
#> 44 M I 9 C
#> 45 M I 4 B
#> 46 F I 2 B
#> 47 M II 10 C
#> 48 M II 1 A
#> 49 M III 3 A
#> 50 F III 5 C
#> 51 F II 9 B
#> 52 F I 7 C
#> 53 F I 2 A
#> 54 M I 6 A
#> 55 F II 8 B
#> 56 M I 6 C
#> 57 F I 4 C
#> 58 M II 4 A
#> 59 M I 1 A
#> 60 F I 1 C
#> 61 F II 4 A
#> 62 M II 9 B
#> 63 F II 8 B
#> 64 F I 2 B
#> 65 M III 6 B
#> 66 F I 5 C
#> 67 M III 2 C
#> 68 F I 5 A
#> 69 F I 4 B
#> 70 F I 3 A
#> 71 M I 1 B
#> 72 F III 10 B
#> 73 F II 7 A
#> 74 M I 3 B
#> 75 M I 1 C
#> 76 M I 7 A
#> 77 M II 2 C
#> 78 M III 9 C
#> 79 F I 3 C
#> 80 M II 4 C
#> 81 M II 9 C
#> 82 F I 3 C
#> 83 M I 5 B
#> 84 F I 7 A
#> 85 F II 2 A
#> 86 F I 8 B
#> 87 F II 10 B
#> 88 M III 6 A
#> 89 M III 3 B
#> 90 M I 2 A
#> 91 F III 1 C
#> 92 F I 4 A
#> 93 M II 2 B
#> 94 F II 1 C
#> 95 M II 6 C
#> 96 M III 9 A
#> 97 M I 4 B
#> 98 M I 6 C
#> 99 F II 5 B
#> 100 M I 6 A
#> 101 M II 6 A
#> 102 F II 10 A
#> 103 F III 2 B
#> 104 M I 6 B
#> 105 M I 1 C
#> 106 F I 6 C
#> 107 M II 8 A
#> 108 F III 1 A
#> 109 M I 1 B
#> 110 M II 5 B
#> 111 F I 1 A
#> 112 F I 5 C
#> 113 F I 5 B
#> 114 M II 7 C
#> 115 F I 7 A
#> 116 M III 7 B
#> 117 M I 6 A
#> 118 M III 7 A
#> 119 M II 3 A
#> 120 M I 8 C
#> 121 M I 6 B
#> 122 F II 9 C
#> 123 F I 3 B
#> 124 F I 3 C
#> 125 F I 8 A
#> 126 M I 8 B
#> 127 M II 3 B
#> 128 M II 2 A
#> 129 F II 8 B
#> 130 M II 5 C
#> 131 F I 10 C
#> 132 M I 9 A
#> 133 F I 5 B
#> 134 F II 8 A
#> 135 M I 5 C
#> 136 F II 4 B
#> 137 M III 5 C
#> 138 M I 5 B
#> 139 M I 4 B
#> 140 F I 8 C
#> 141 F I 8 C
#> 142 M I 3 A
#> 143 F I 1 A
#> 144 F II 3 A
#> 145 M I 6 A
#> 146 M I 2 C
#> 147 F III 3 A
#> 148 F I 9 B
#> 149 M II 9 C
#> 150 M III 5 B
# Use 1:1:2 ratio
minimise(patients, groups = 3, factors = c("sex", "stage"), burnin = 5,
ratio = c(1,1,2))
#> Multi-arm Minimisation
#> --------------------------------------------------------------------------------
#> Groups: A B C
#> Randomisation ratio: 1:1:2
#> Factors: sex, stage
#> Burnin: 5
#> Minimisation probability: 0.8
#> Group sizes: 37, 38, 75