The number of cell type 1 ( × ) is plotted against the total cell number in model II. The initial condition of u$_{i}$ is chosen randomly from the interval u$_{i}$ ∈ [0 , 1]. Solid and bro- ken lines indicate N$_{(1)}$ ( N, v ∗ $_{1}$) = − A ( v ∗ $_{1}$) N + B ( v ∗ $_{1}$) and N$_{(1)}$ ( N, v ∗ $_{2}$) = − A ( v ∗ $_{2}$) N + B ( v ∗ $_{2}$), respectively, where A ( v ∗ $_{1}$) = 0 . 011, B ( v ∗ $_{1}$) = 24, A ( v ∗ $_{2}$) = 0 . 0099, and B ( v ∗ $_{2}$) = 93. The parameters are K$_{v}$ = 2 . 0, β = 2 . 0, c$_{2}$ = 0 . 1. 

## Figure 7: 

The slope A ( β ; v ) is plotted as a function of β . Here, A ( β ; v ) for two different values of v , i.e., v ∗ 1 and v ∗ 2 are plotted, which agree within the resolution of the plot in the figure. The parameter values are K$_{v}$ = 2 . 0, c$_{2}$ = 0 . 1. 

## Figure 8: 

The fixed point solutions of model III plotted against the total cell number N . At each N , 100 initial conditions are chosen. The expression levels of u$_{i}$ for a single cluster (+) and two-cluster solutions ( ◦ ) are plotted as a function of N . The value for two-cluster so- lutions is the average over initial conditions. The parameter values are ˜ K$_{v}$ = 0 . 2, β = 2 . 0, c$_{v}$$_{1}$ = c$_{v}$$_{2}$ = 0 . 005. 

## Figure 9: 

The ratio of the number cell type 1 N$_{(1)}$ to the total cell number N is plotted against N for model III. The initial condition of u$_{i}$ is chosen randomly from the interval of u$_{i}$ ∈ [0 , 1]. Solid and broken lines indicate N$_{(1)}$ ( N, v ∗ $_{1}$) /N = ˜ A ( v ∗ $_{1}$) + ˜ B ( v ∗ $_{1}$) /N and N$_{(1)}$ ( N, v ∗ $_{2}$) /N = ˜ A ( v ∗ $_{2}$) + ˜ B ( v ∗ $_{2}$) /N , respectively, where ˜ A ( v ∗ $_{1}$) = 0 . 16, ˜ B ( v ∗ $_{1}$) = 69, ˜ A ( v ∗ $_{2}$) = 0 . 36, and ˜ B ( v ∗ $_{2}$) = 86. The parameter values are ˜ K$_{v}$ = 0 . 2, β = 2 . 0, c$_{v}$$_{1}$ = c$_{v}$$_{2}$ = 0 . 005. 

