In order to compare the behavior of the different splitting formulas, let's consider a simple example:
We have $N$ predictors with value $y=1$ in the left, $N$ predictors with value $y=-1$ in the right, and a single point at the right boundary with value $y^\ast$.
There are two reasonable split points in this scenario, first the point (1) in the middle, and the point (2) at the right boundary.
When $N$ is large, say $N=100$, I would wantprefer the split to be made at the point marked by the red (1). Only if the value of $y^\ast$ becomes very low, I would find a split at (2) reasonable.
Let'sSo, let's calculate how the split position changes in dependence of the value of $y^\ast$: For which $y^\ast$ is the split made at (1), or how low $y^\ast$ might get until the split occurs at (2) -- for which.
For the calculations we first use the criterium the OP posted in the question, that iswhich involves the sum of squares (SS) instead of the variance:
Split at (1):
Mean left: $1$, SS Left: $0$
Mean right: $\mu_r = (y^\ast - N)/(N+1)$, SS right: $N(-1 - \mu_r)^2 + (y^\ast - \mu_r)^2$
Split at (2):
Mean left: $0$, SS Left: $2N$
Mean right: $y^\ast$, SS right: $0$
SettingThus, we get a split at point (1) if
$$N(-1 - \mu_r)^2 + (y^\ast - \mu_r)^2 \leq 2N$$
Inserting $\mu_r$ and setting equal the two sum of squares (SS), Wolfram Alpha gives as solution
$$y^\ast = -1 - \sqrt{2N} \ \underbrace{\sqrt{\frac{(N+1)^2}{(N^2+1)}}}_{\approx 1} \ \approx \ -1 - \sqrt{2N}$$$$y_{SS}^\ast = -1 - \sqrt{2N} \ \underbrace{\sqrt{\frac{(N+1)^2}{(N^2+1)}}}_{\approx 1} \ \approx \ -1 - \sqrt{2N}$$
That means, for a value $y$ above this $y^\ast$$y_{SS}^\ast$, the split is made at (1), below $y^\ast$$y_{SS}^\ast$ the split is made at (2). One sees that for large $N$ the split is preferred in the middle.
Now, doing the same for the version whereof the splitting formula is divided by the number of data points insidewhich includes the respective regionvariance, the result is
$$y^\ast = -1 - \sqrt{\frac{(N+1)^2}{(N^2+1)}}$$$$y_{VAR}^\ast = -1 - \sqrt{\frac{(N+1)^2}{(N^2+1)}} \ \approx \ 1$$
That is, the factor $\sqrt{2N}$ is missing. TheWith this, the split is made at (2) if the value at the boundary is slightly lower than $-1$, regardless of the number of datapoints $N$.
Conclusion: The splitting formula without division ofusing the particle numbersum of squares (SS) leads to the imo more intuitive behaviour that the split is donepreferred at point (1) in the middle.
This is reasonableInstead, as without the factorsformula with the absolute quantities matter. Withvariance does not consider the normalization bynumber of particles in the sample sizerespective regions, one considers average quantities by which -- as the SS value is divided by the number of data points. As seen above -- small regions become similar weight, it easily prefers a split where only one data point is separated to very large regionsa balanced split where the region is halfed.
Summarizing, I would prefer the splitting criterium ascontaining the OP gave it in his questionsum-of-squares to the splitting criterium containing the variance.