Why Repeat a Random Split?
A single random train-test partition gives one estimate of predictive performance. A different partition of the same data may give a different estimate: some observations that are difficult to predict may land in the test set in one split but not another. Repeated random-split validation makes this sensitivity visible by fitting and evaluating the same learning procedure on several independently generated partitions.
The result has two parts that should be reported separately: the average test error across splits, and the variation in test error from split to split. The average describes the central result under the chosen splitting procedure. The variation describes how sensitive that result is to the random partition. Neither quantity, by itself, is a confidence interval for performance on a new population.
The loss \(L_r\) must match the performance question. For regression it could be mean squared error (MSE) or mean absolute error; for classification it could be misclassification rate. A smaller value is better for these loss measures. Define the metric, the training fraction, and the random-split procedure before comparing results. When test-set sizes differ, also decide whether the target is the average error for a randomly selected split or the error pooled over all test predictions. Those are not generally the same quantity.
What the Mean and Standard Deviation Describe
The mean \(\bar L\) summarizes the errors observed under the chosen split procedure. The standard deviation \(s_L\) summarizes their dispersion in the same units as the loss. For example, if test error is a proportion, its standard deviation is also measured in proportion units; if the metric is an MSE in squared response units, the standard deviation is in those squared units.
Repeated splits reuse the same observed data. Consequently, \(s_L\) is a measure of split sensitivity conditional on the data set and procedure; it is not a measure of all uncertainty about future predictive performance. In particular, repeating the split does not create new independent observations from the population, remove limitations in the observed data, or correct a poorly chosen model procedure. More repetitions can make the estimated split-distribution mean less affected by Monte Carlo randomness, but they do not make the original data set larger.
If the split assignments are independently generated, then, conditional on the fixed data and a deterministic fitting procedure, the split losses are independent draws from the split-randomization distribution. In that setting, \(s_L/\sqrt{R}\) describes the Monte Carlo standard error of the estimated mean over that distribution. It does not measure uncertainty from drawing a different data set from the population. If the fitting procedure itself uses randomness, that randomness must also be specified: either hold it fixed to isolate split variation, or repeat it as part of the procedure whose variability is being summarized.
Worked Examples
Worked Example: Summarizing Repeated Regression MSEs
A hypothetical insurer uses a fixed regression procedure to predict claim-handling time. It makes four random splits, each with four test observations. The sum of squared test residuals (SSE) is 16, 36, 4, and 24 minutes squared, respectively. Calculate each test MSE, the repeated-split mean MSE, and the sample standard deviation of the four MSEs.
Method. For each split, divide its test SSE by its test-set size, \(n_{\text{test}}=4\). Then calculate the mean and sample standard deviation of the four split MSEs.
Split 1: \(L_1=16/4=4\) minutes squared.
Split 2: \(L_2=36/4=9\) minutes squared.
Split 3: \(L_3=4/4=1\) minute squared.
Split 4: \(L_4=24/4=6\) minutes squared.
The mean is \(\bar L=(4+9+1+6)/4=20/4=5\) minutes squared. The deviations from 5 are \(-1,4,-4,1\); their squared values are \(1,16,16,1\), with sum 34. Thus the sample variance is \(34/(4-1)=34/3=11.3333\) minutes to the fourth power, and the sample standard deviation is \(\sqrt{34/3}=3.3665\) minutes squared.
The four test sets have equal size, so pooling all test residuals gives the same MSE: total SSE is \(16+36+4+24=80\), and the total number of test predictions is \(4(4)=16\), giving \(80/16=5\). Check: this matches the mean of the four split MSEs. The answer is \(5\) minutes squared with a split-to-split sample standard deviation of \(3.3665\) minutes squared (rounded to four decimals). The common mistake to avoid is dividing an individual split's SSE by the number of repetitions rather than by that split's test-set size.
Worked Example: Unequal Test Sizes and Two Different Averages
A hypothetical classifier predicts whether an equipment inspection requires urgent follow-up. Across three random splits, test-set sizes are 8, 10, and 12, with 1, 3, and 2 misclassifications, respectively. Find the mean of the split error rates and the pooled error rate. Explain why they differ.
Method. Calculate each split's error rate by dividing its number of errors by its own test-set size. Average those three rates with equal weight to obtain the mean split error. Separately, divide the total number of errors by the total number of test predictions to obtain the pooled rate.
The split rates are \(1/8=0.1250\), \(3/10=0.3000\), and \(2/12=0.1667\) (the last displayed to four decimals). Their equal-weight mean, using the unrounded fractions, is \((0.125+0.300+0.166666\ldots)/3=0.197222\ldots\), or \(0.1972\).
The pooled rate is \((1+3+2)/(8+10+12)=6/30=0.2000\). These answers differ because the split mean gives each split equal weight, whereas the pooled rate gives each test prediction equal weight. Check: the largest split rate is 0.3000, and the pooled rate lies between the smallest and largest split rates, as expected. The common exam trap is to call the pooled rate the “average split error” without checking whether test-set sizes are equal. Here the requested split-average is 0.1972; the pooled error rate is 0.2000.
Worked Example: Average RMSE Is Not the Square Root of Average MSE
A hypothetical regression procedure is evaluated on four test sets, each of size four. Its test SSEs are 16, 64, 100, and 4. Calculate the RMSE for each split, the mean split RMSE, the sample standard deviation of split RMSE, and the square root of the mean split MSE.
Method. For each split, calculate \(\mathrm{MSE}=\mathrm{SSE}/4\) and \(\mathrm{RMSE}=\sqrt{\mathrm{MSE}}\). Then summarize the RMSE values separately from the MSE values.
The MSEs are \(16/4=4\), \(64/4=16\), \(100/4=25\), and \(4/4=1\). The RMSEs are therefore \(2,4,5,1\). Their mean is \((2+4+5+1)/4=12/4=3\). Deviations from 3 are \(-1,1,2,-2\); the squared deviations sum to \(1+1+4+4=10\). The sample standard deviation of RMSE is \(\sqrt{10/(4-1)}=\sqrt{10/3}=1.8257\).
The mean MSE is \((4+16+25+1)/4=46/4=11.5\), so the square root of the mean MSE is \(\sqrt{11.5}=3.3912\). This is not the mean RMSE of 3. Check: the two calculations apply the square root at different stages, and a nonlinear transformation does not generally commute with averaging. The requested mean split RMSE is 3.0000, its sample standard deviation is 1.8257, and the square root of the mean MSE is 3.3912 (rounded to four decimals). The common mistake is to report \(\sqrt{\text{mean MSE}}\) when asked for the average of split RMSEs.
Why the Sample Variance Uses \(R-1\)
The sample variance formula uses \(R-1\), rather than \(R\), because the same data are used to calculate both the sample mean and the deviations from it. The following justification shows the adjustment when the split losses can be treated as independent draws from a distribution with mean \(\mu_L\) and variance \(\sigma_L^2\). This is a statement about repeated split randomization under that setup; it does not assert that the split losses describe all uncertainty in predictive performance.
To verify the identity, write \(L_r-\bar L=(L_r-\mu_L)-(\bar L-\mu_L)\), square, and sum over \(r\). The cross-product is zero because \(\sum_{r=1}^R(L_r-\mu_L)=R(\bar L-\mu_L)\). Taking expectations, the first sum on the right has expectation \(R\sigma_L^2\). The second term has expectation \(R\operatorname{Var}(\bar L)=R(\sigma_L^2/R)=\sigma_L^2\), using independence and equal variances. Therefore:
Dividing by \(R-1\) thus makes the sample variance unbiased for the split-randomization variance under these assumptions. Dividing by \(R\) instead would systematically underestimate that variance on average. This justification does not turn \(s_L\) into a population uncertainty measure; it explains the denominator used to summarize variation across the repeated splits.
Common Traps and Practical Interpretation
- Calling split variation population uncertainty. A small \(s_L\) means the observed data set produced similar results under the chosen partitions. It does not establish that the same error will hold in a new population.
- Ignoring test-set size. The unweighted mean of split rates and the pooled rate answer different questions when test sizes differ. State which one is being reported.
- Confusing averaging scales. Mean RMSE is the average of the split-specific square roots of MSE. It generally differs from the square root of average MSE.
- Assuming more repetitions fix a weak evaluation design. Increasing \(R\) gives a more stable summary of the selected split procedure. It does not undo an unsuitable split fraction or make the available observations more representative.
- Leaving the procedure unspecified. The reported result depends on the metric, training fraction, split mechanism, and fitting procedure. If the goal is to assess a full modeling procedure, each split must evaluate that procedure as it would actually be applied using training data.
Repeated random-split validation is most informative when its summary is read as a pair: the average error indicates typical performance under the specified partitioning rule, while the standard deviation indicates how much that estimate changes across partitions. A low mean with a high standard deviation signals that the central estimate alone hides substantial split sensitivity. A low standard deviation with a high mean indicates stable but poor performance under the chosen metric.
Check Your Understanding
For each problem, identify whether the requested quantity is a split average, a pooled error, or a measure of split variation.
- A regression procedure has test MSEs of 3, 5, 8, and 4 across four equally sized test sets. Calculate the mean MSE and the sample standard deviation of the MSEs. (Answer: mean MSE \(=5\); sample standard deviation \(=\sqrt{14/3}\approx2.1602\).)
- A classifier makes 2 errors among 20 test cases in one split and 3 errors among 10 test cases in another. Find the equal-weight mean of the two split error rates and the pooled error rate. (Answer: split mean \(=(0.10+0.30)/2=0.2000\); pooled rate \(=5/30\approx0.1667\).)
- Three repeated splits give loss values 0.12, 0.18, and 0.16. Calculate the mean loss and sample standard deviation, rounding both to four decimals. (Answer: mean \(=0.1533\); sample standard deviation \(\approx0.0306\).)
- Under independently generated split assignments, a practitioner increases the number of repetitions from \(R=25\) to \(R=100\). What happens approximately to the Monte Carlo standard error of the mean split error, and what does not change? (Answer: it is approximately halved because it scales as \(1/\sqrt{R}\); this does not eliminate uncertainty from the original data set or change the chosen split procedure.)