Reconsideration of r/K Selection Theory Using Stochastic Control Theory and Nonlinear Structured Population Models

Despite the fact that density effects and individual differences in life history are considered to be important for evolution, these factors lead to several difficulties in understanding the evolution of life history, especially when population sizes reach the carrying capacity. r/K selection theory explains what types of life strategies evolve in the presence of density effects and individual differences. However, the relationship between the life schedules of individuals and population size is still unclear, even if the theory can classify life strategies appropriately. To address this issue, we propose a few equations on adaptive life strategies in r/K selection where density effects are absent or present. The equations detail not only the adaptive life history but also the population dynamics. Furthermore, the equations can incorporate temporal individual differences, which are referred to as internal stochasticity. Our framework reveals that maximizing density effects is an evolutionarily stable strategy related to the carrying capacity. A significant consequence of our analysis is that adaptive strategies in both selections maximize an identical function, providing both population growth rate and carrying capacity. We apply our method to an optimal foraging problem in a semelparous species model and demonstrate that the adaptive strategy yields a lower intrinsic growth rate as well as a lower basic reproductive number than those obtained with other strategies. This study proposes that the diversity of life strategies arises due to the effects of density and internal stochasticity.


Introduction
Density effects in life history are considered to be important for evolution. Intraspecific competition has led to the evolution of antlers and peacock feathers in sexual selection, and another type of selection has led to phase variation, depending on the population density, in migratory locusts [1]. Determining what life strategies can evolve in each population size is a difficult question; these questions describe an area of research referred to as r/K-selection theory. This theory has captured the attention of ecologists for several decades. It was proposed in Mac-Arthur and Wilson's book [2], which discussed it in the context of the differences between adaptive reproductive strategies at low and high levels of intraspecific competition. At the lower levels, such as an environment affected by disturbance (r-selection), life histories evolve to maximize the intrinsic rate of natural increase (IRNI), whereas at the higher levels, such as saturated populations (K-selection), life histories evolve to maximize the carrying capacity. These different life histories were classified by Pianka [3]. For instance, species affected by rselection have precocity and prolificacy. Conversely, the other selection involves altricity in addition to producing a smaller number of offspring. According to Reznick et al (2002), the classifications and theory have been discussed in many papers, with both agreement and disagreement, since the theory was first published [4]. The papers that disagree with the theory state that several species do not seem not to match Pianka's classification, e.g., trees with both longevity and prolificacy. To clarify the issue, researchers have developed many precursor theoretical models and approaches, such as optimal life schedule problem (OLSP), adaptive dynamics (AD), structured demographic models (SDM), and so on. The paradigm of r/K selection theory has shifted to a theory of density dependence of mortality in SDMs. Many studies are based on the hypothesis that all intraspecific competition influences fertility and mortality. Furthermore, a recent work extended these models to adapting environmental stochasticity in addition to density effects [5]. Ecologists speculate that adaptive life histories under both selections evolve as the consequence of density and environmental effects.
One reason why r/K selection theory has attracted the attention of ecologists is that there are several difficulties concerning definitions in K-selection; it is not well understood what individuals should maximize (or minimize) in the selection. Because the basic reproductive number R 0 is always one for individuals in saturated populations, it is difficult to find factors augmenting their carrying capacity on the individual scale. AD postulates that adaptive species have an evolutionarily stable strategy (ESS) that cannot be disrupted by any other mutant strategies. The analysis of AD requires appropriate replicator dynamics providing a fitness function (or invasion fitness). The fitness function does not necessarily clarify the relationship with individual life dynamics such as size growth, reproduction, death, and so on. Age-SDMs can be applied to this problem [6,7]. The age-SDMs define carrying capacity as the interior equilibrium generated by the density effects on generations in terms of fertility and mortality. This model uses a characteristic function incorporating fertility and mortality functions with respect to age that includes IRNI and R 0 as the fitness function. Though it is reasonable, the model deals with intraspecific competitions in terms of only fertility and mortality. For instance, competition for food intake may inhibit growth and maturity for each individual. The inhibition will reduce IRNI, even if R 0 is preserved. To handle the wide range of life history evolutions, age-SDMs are not sufficient; they should take into account other ingredients (e.g., size, spatial position, genetic expressions, and individual differences). Given the components of life history, the theory associated with r/K-selection should provide an appropriate definition of adaptive strategy in the context of carrying capacity. Therefore, a theoretical framework is required that describes the population dynamics and life schedule of individuals simultaneously.
These problems mainly lack a theory that connects individuals and the population scale. Recent works have addressed this issue. Oizumi and Takada (2013) derived a characteristic function that incorporates age and other factors to characterize life history using age-size SDMs with diffusion by introducing probability theory [8]. Their analysis showed that diffusion could change the generation time and R 0 even if the life histories were identical. They termed the diffusion "internal stochasticity", which represents individual differences in life schedules, whereas the stochasticity on the population scale, such as environmental or demographic stochasticity, was designated as "external stochasticity". The results of this paper raised the possibility that internal stochasticity leads to life histories based on Pianka's classification. Moreover, Oizumi (2014) extended a previous work to unify a model of linear SDMs and stochastic control theory [9]. This theoretical model consists of two equations that both provide a more general characteristic equation based on Oizumi and Takada's equation. These equations can simultaneously address the analysis of optimal life strategies and population dynamics under internal stochasticity. However, they do not account for density effects. Because the evolutionary diversity of life histories is viewed as the involvement of intraspecific competition at each population scale, Oizumi's method should be extended to address r/Kselection theory.
To this end, this study establishes a more general framework containing both Oizumi's methods and density effects. An appropriate fitness function is proposed that addresses both selections simultaneously. The function uses the Euler-Lotka equation in the absence of density effects and yields IRNI and R 0 . Conversely, the function yields an equation where the density effects are satisfied at equilibrium. These properties of the function originate from Leon's work based on age-SDMs [6]; however, the crucial difference is that the fitness function is generated by life histories, which are formulated using stochastic differential equations. Therefore, this study is able to treat not only age-specific SDMs but also arbitrary state-specific SDMs with diffusion. Applying our framework to specific models concerning semelparous species, we show that with internal stochasticity, we can reduce the carrying capacity. Furthermore, we present a specific example demonstrating that an adaptive risk hedge of internal stochasticity augments the carrying capacity and that it is ESS under K-selection.

Framework for the mathematical methods
To address life histories in the presence of internal stochasticity and density effects, we outline our methodology and examine a general difference in adaptive life histories between r-and Kselection. The life histories of individuals are assumed to consist of four factors: growth, fertility, mortality, and interaction among individuals in our analysis. First, we define the models for these four factors. Second, we describe the population dynamics of the life history. Finally, the difference in adaptive strategies is discussed for both selections.

General life history model
To link population-scale density effects to individual life histories, it is naturally thought that life histories should contain the density effects in the statistics related to the components, e.g., a mean size growth rate, the variance, etc. Moreover, individual life histories containing internal stochasticities are assumed to be diffusion processes of each state transition. Specifically, they include the migration of the population, size growth, and so on [8,[10][11][12]. Accordingly, we construct a general model of the life history influenced by diffusion processes with density effects. We define a set A R d of all states. Suppose that X a (t) 2 A represents a d-dimensional state vector at age a and time t (e.g., size, weight, spatial position, etc.). The initial state is given as X 0 (t) = x, which is independent of time t. Additionally, it is assumed that each state transition is regulated by an ldimensional control vector v at the given age as per the individual life strategy [9] v ¼ v a t ð Þ : where V is a compact convex set of R l that denotes a set of all strategies. Given each state transition that is influenced by the N-dimensional internal stochasticities and parameterized by the Mdimensional density effect coefficient vector G t;x :¼ ðG 1 t;x ; G 2 t;x ; . . . ; G m t;x ; . . . ; G M t;x Þ, which represents the magnitude of each density effect on the statistics of X a (t) at time t, each diffusion process of states X j a ðtÞ is assumed to satisfy Ito's stochastic differential equation as follows: where x j (1 j d) represents an initial state of the j-th state, X j a ðtÞ. The SDE is defined so that the difference in each state at an infinitesimally short time follows a Gaussian distribution. Because time t and age a are considered to have the same temporal development, the term "infinitesimally short time" is used here in the sense of direction of their characteristic line, i.e., This corresponds to the fact that the individual size growth and other transitions are different at each cohort. The first term g j (Á) on the right hand side in Eq (1) denotes the expectation of each state transition rate, and the second term denotes the fluctuation term and provides a covariance between another state provided by the Gaussian distribution, as follows: For the boundary value x 2 @A, each state transition rate and the fluctuation term are assumed to be g j (x, v, Γ t,x ) = σ jk (x, v, Γ t,x ) = 0. This assumption means that several states attribute to the limited region. For instance, the domain of body-size is restricted to positive value. B k a represents an element of N-dimensional Brownian motion generating internal stochasticity. Each element of the density effect coefficient vector Γ t,x can be determined using the current density effects: where α > 0, γ m (a, y) ! 0 and P t ða; x ! yÞ : ½0; 1Þ Â ½0; aÞ Â A Â A ! R þ :¼ ½0; 1Þ denote the maximum age, interaction coefficient of each generation and state, and population density of state y 2 A at age a 2 (0, α) in time t 2 (0,1), respectively. The diffusion process Eq (1) is complicated but includes the transition process in both r-and K-selection.
x is unique. Then, X a (t) becomes a diffusion process with respect to age in extreme environments where Γ t,x is constant. Thus, Eq (1) accounts for the individual differences in the transition process that are characterized by each sample path and have different statistics for each population size. In addition, other components of the life histories are assumed to be functions with respect to X a (t) and Γ t,x .
Let F : A Â R M þ ! R þ be a state-specific fertility, which depends on y and has integrability with respect to it for all G 2 R M þ as follows: This definition of fertility F accounts for a wide class of breeding systems, such as a semelparous fertility with which the Dirichlet boundary conditions can be applied to population dynamics, as will be mentioned later on.
Let m : R d Â R l Â R M ! R þ be a mortality function depending on y, v, and Γ t,x as follows: The density effects on the fertility and the mortality are identical to the assumptions for the density-dependent age-structured models involved in the Allee effect and logistic effect [13]. We refer to Eqs (1), (3), and (4) as the life history.

Dynamics at the population scale
In this subsection, we consider properties of the population dynamics derived from the components and statistics of the life history (Eqs (1), (3), and (4)). When t and a are fixed, the state transition of the population density P t (a, x ! y) at an infinitesimally short time ε follows an integral projection model: K ε, t (ξ!y) represents a projection function surviving and ranging from ξ to y at time interval ε and time t. The projection function is determined by Eqs (1) and (4). Let n t (x) be the total neonatal density per unit time at time t, e.g., the total number of eggs and seeds, as well. The initial population density is written using n t (x): In addition, the neonatal dynamics follow Using the components of the life history Eqs (1), (3), and (4), we derive a non-linear partial differential equation from Eq (5) (Cf. Text A in S1 File) as follows: where c jj 0 (y, v, Γ t,x ) denote Pða; x; yÞ : ½0; aÞ Â A Â A ! R þ is the initial population density. It is then assumed that all individuals have the same initial state x. The differential operator H v y ðG t;x Þ is referred to as the Fokker-Planck Hamiltonian operator, which provides the state transition with diffusion. To analyze the adaptive strategies in r/K-selection, it is important to determine the asymptotic behavior of the population dynamics at the two equilibria representing each selection. Hence, we focus on the population sizes satisfying and @ @a where 0 2 R M þ and G Ã x 2 R M þ denote the zero vector and density effect coefficients for a nontrivial equilibrium, respectively. Eq (9) provides the asymptotic behavior of the population size, which is sufficiently small so that the density effect can be neglected: (G Ã x % 0), i.e., it represents r-selection. On the other hand, Eq (10) implies K-selection, where the size reaches a non-trivial and stable equilibrium. Both conditions can be rewritten using a characteristic function where l 2 R and G 2 R M þ are constant and constant vector, respectively. The expectation E x ½Á is given by The projection function K a (x ! y), then, is the solution of which represents the density surviving and growing from x to y at age a under Eq (1), providing the statistical rule. The IRNI of Eq (9) is provided by λ Ã satisfying 8,9], whereas, G Ã x in Eq (10) satisfies The former implies the characteristic equation of Eq (9) (also called the Euler-Lotka equation), and the latter represents that R 0 is equal to one at equilibrium. Thus, Eq (11) can address various steady states of Eq (8) with different values of Γ. In other words, Eq (11) represents the reproductive contributions of all cohorts to the population dynamics for each population size.
Further, the logarithm of the characteristic function Eq (11) is composed of the cumulant generating function of the breeding age a Ã , which contributes to reproduction, such that where ha Ã i ðkÞ x;G represents the k-th cumulant [8]. These statistics are provided by a breeding age distribution parameterized by Γ This distribution is important for considering the effect of r/K-selection on the breeding system. Therefore, Eq (11) can be used as a common criterion in both selections and plays a pivotal role in this study.
The life histories of r/K-selection are based on the assumption that Eq (1) has a logistic type of behavior; (i) Eq (8) has a unique equilibrium population density except at zero (P Ã (a, x ! y) > 0), and (ii) the vector is stable for any v 2 V when the other (P Ã (a, x ! y) = 0) is unstable, as follows: Under conditions (i) and (ii), the carrying capacity can be derived from the density effect at equilibrium: It is known that the equilibrium population can be decomposed into the total neonatal density n Ã (x) and the projection function K Ã a ðx ! yÞ, such that [9]. Then, the projection function represents a density at which individuals survive and reach y from x at age a under the density effect G Ã x . The projection function is composed of and hence, we have where E Ã x denotes the expectation of the projection function under the non-trivial equilibrium, K Ã a ðx ! yÞ. From this, we obtain for all m. If this equilibrium is stable (the condition referred to Text B in S1 File), the total r/K Selection Theory via Stochastic Analysis population size N Ã is equivalent to the carrying capacity and is Conventional K-selection theory argues that adaptive strategies maximize N Ã . Taking into account the influence of the population structure, the argument should be reconsidered at the individual scale. For instance, when two different species with the same carrying capacity invade identical habitats, conventional K-selection theory cannot assess whether either of the two life histories is adaptive. To address this sort of question, we present an answer to the one mentioned above using an extension of Oizumi's method to incorporate nonlinear systems in the next section.

Adaptive life histories
At the beginning of the analysis, we introduce an important theorem related to Eq (11): and c v l ðx; GÞ be given by Eq (11) and c v l Ã;v ðx; GÞ ¼ 1, such that Definel byl :¼ l Ã;ṽ . Then, we have This theorem describes that a strategyṽ for maximizing Eq (11) is equivalent to maximizing λ Ã . This theorem is easily verified because of the monotonicity of Eq (11) with respect to λ. Theorem 1 provides thatṽ is defined as the adaptive strategy at Γ. Accordingly, we adopt Eq (11) as the fitness function in this paper. To findṽ, the maximized characteristic function is to be determined. Hence, we introduce Oizumi's equations parameterized by Γ such that Eq (25) is one of the Hamilton-Jacobi-Bellman (HJB) equations in Control Theory [14], the solution of which,w l;a ðx; GÞ, is the following function described using the age in descending order (a ! α − a); w l;a x; G ð Þ The reason behind the variable change in age originated from Bellman's principle with regard to derivation of the HJB equation (Cf. Text C in S1 File and [9]).
H v x ðGÞ represents the formal adjoint operator of H v y ðGÞ in terms of the variable x in Eq (8). Then, Eq (27) represents the maximized reproductive contribution of cohort at age a, which includes a maximized characteristic function, such thatc Eq (26) where the strategy is expressed in terms of the original age a. We then add a new assumption as follows: (iii) Eq (29) exists and is unique. Substituting Γ = 0 into Eq (29),lðx; 0Þ represents the IRNI of the adaptive species andṽðG ¼ 0Þ denotes the adaptive strategy in the absence of the density effect, which is referred to as the r-strategyṽ r . In contrast, we show that in the presence of the density effect,G satisfyingl ¼ 0 provides the adaptive strategy in K-selection (K-strategy:ṽ K ¼ṽðG ¼GÞ).

K-strategy
Substituting G ¼G into Eq (26) givesc Theorem 1 indicates that the adaptive species under the density effectsG has the maximum IRNIl equal to zero: the strategyṽðG ¼GÞ is adaptive in terms of the maximum density effect. Suppose that individuals with a strategy v reach an equilibrium given by Γ Ã . If for Γ Ã , lðx; G Ã Þ > 0 in Eq (29), the population can be invaded because there exists a strategy that allows individuals to bear more offspring with the given population size. No strategy, however, can invadeṽðG ¼GÞ at equilibrium because for all v. Accordingly,ṽðG ¼GÞ is consistent with the definition of ESS; hence, we call this strategy K-strategyṽ K . Eq (31), then, has a comparable meaning to the evolutionary invasion analysis in AD. Focusing on what K-strategy modifies on the population scale, the possibility is raised thatṽ K generates complicated state structures of the population depending on each interaction coefficient γ m (a, y). It is significantly different from traditional K-selection theory becauseṽ K does not maximize the carrying capacity. This difference is accounted for by considering the structure of the intraspecific interaction as the density effect. If and only if the m 0 -th interaction coefficient in Eq (2) is constant, g m 0 ða; yÞ g m 0 0 ,G is equivalent to the maximum carrying capacity as follows: This may not be necessarily true otherwise. Taken together,ṽ K contains classical K-selection theory as a particular case. Eqs (25) and (26) can possibly be applied to OLSP influenced by various effects of intraspecific competition and internal stochasticity. In order to show the utility of both equations, in the next section, we illustrate that there exists a specific model matching our framework and examine how both effects influence the life history using the model.

Specific example
This section addresses a specific model following the assumptions (i)-(iii) to examine how the density effect and internal stochasticity influence the adaptive strategy. As a mathematical example, the model which we would handle is suggestive on the examination because it has the exact solution and several solvable indices. By use of the example, we then focus on these effects for alternations of generations. According to Pianka's classification mentioned in the Introduction [3], r-strategy has a shorter generation time than other strategies. To examine this, we demonstrate in the next subsection that internal stochasticity affects the breeding age density (see Eq (17)) and the carrying capacity. Second, we analyze the evolution of adaptive life history by using two-resource utilization models following the former results. Then, the adaptive strategy of the model is analyzed in terms of the statistics of the breeding age compared with Pianka's table.
Resource acquisition competition model. The role of internal stochasticity at the population scale should be clarified before analyzing adaptive control. Accordingly, we start by analyzing the effect of stochasticity on the life history without control parameters. Suppose that the size of the individuals X 1 a 2 A & R þ grows in accordance with a geometric Brownian motion and that the average size growth rate is reduced by intraspecific competition for foraging Γ t as follows: where b 1 ; G; s 1 2 R þ . If Γ t does not exist, the size of the individuals grows exponentially with fluctuation. We call this model the resource acquisition competition model. The breeding system is assumed to be semelparity. Let x Ã > x be the mature size, and let the breeding age a Ã be defined as the first age when a size reaches the mature size as follows: The breeding age is, therefore, not unique but is represented using a statistic because of the noise. We refer to the age as the mature age in this model. Therefore, the mature age is a statistic because of the noise. The size-specific fertility F S (y), then, is defined by where the fertility rate function ðyÞ : R þ ! R þ denotes the fecundity per individual, which satisfies ϕ(x Ã ) > 1 at the mature size.
The mortality in Eq (4) is assumed to be where μ 0 > 0, which does not have control parameters. Because semelparous species die upon reproduction, we can write the survivorship as The interaction coefficient γ (a, y) is assumed to be constant as follows; Eq (38) seems to be oversimplified, but it is not important for the adaptive strategy because the interaction coefficients of our framework are only related to the stability and population size at equilibrium. These assumptions of semelparous life history are reasonable for comparing differences between the absence and the presence of the density effect because each individual bears an identical number of offspring. Accordingly, the density effect is only involved in the statistics of the mature age; in other words, it is sufficient to pay attention to the effect on generation time and the size of Γ Ã .
respectively. ρ Ã denotes If the IRNI is positive, Eq (42) becomes the carrying capacity, i.e., it is unique and stable (cf. Text D in S1 File). To examine the relationship between the mature age structure and internal stochasticity, we derive Eq (17) from Eq (40) (cf. Text E in S1 File). Fig 1 illustrates the effect of internal stochasticity on the mature age density, implying that the stochasticity promotes maturity in both r-and K-selection. According to the previous paper [8], the bias of the density is due to low survivorship of late-mature individuals.
Focusing on the effect of internal stochasticity on the IRNI and the carrying capacity, it was shown in the previous paper that the sensitivity of the IRNI with respect to the magnitude of the stochasticity σ 1 is determined by ρ Ã being greater or less than one (cf. [8]). If ρ Ã is greater than one, the sensitivity is positive; otherwise, it is zero or negative. Similarly, the carrying capacity has an identical condition: Furthermore, because the sensitivities of the IRNI and N Ã with respect to ρ Ã are always positive in the persistent domain (λ Ã > 0 i.e., ρ Ã > ρ 0 (0)), the adaptive mature sizex is obtained using the chain rule in calculus as follows: That is,r x is, thus, ESS in both selections, which means that nature selects a mature size maximizing ρ Ã . Because ρ Ã plays a key role in the risk hedging of the life history, we term it "Efficiency Exponent of Semelparity" (EES). Maximizing EES evokes two different interpretations; one is that the size maximizes the number of the individual's own offspring in the numerator, and the other reduces the proportion of the initial size to the mature size in the denominator. The latter accelerates alternations of generations. If the EES is equal to one, the influence of internal stochasticity disappears, which is interpreted as the remaining valance between the length of the life span and basic reproductive number. This suggests that species with an EES of one seem to follow the deterministic rule σ 1 = 0.
Internal stochasticity becomes a negative risk for species with EES less than one ρ Ã < 1 because of increasing the risk of extinction by reduction of the carrying capacity and delaying recovery of the population from large disturbances. An optimal risk hedging is required for the life history. We consider the optimal risk hedging of resource utilization as an application of Eqs (25) and (26) in the next section.

Density-dependent two-resource utilization model
Based on Eqs (33), (35), (36), and (38), we consider a species utilizing different resources (we call them R 1 and R 2 ). The size of a specialist utilizing R 1 is the following growth rate and that of utilizing the other resource, R 2 , is B 1 a and B 2 a are Brownian motions that are independent of each other. We then assume that with b 1 2 R þ > b 2 2 R (b 2 could be negative), σ 1 > σ 2 ! 0, i.e., if Γ t % 0, choosing R 1 implies a higher risk and expectation of growth rate than choosing R 2 . Conversely, choosing R 2 under the same conditions confers another risk that individuals have a lower survivorship until they reach a mature age than choosing R 1 because of the slower growth rate, on average. A constant κ 2 (0,1) denotes the ratio of strength in intraspecific competition, which alters the statistically significant difference between these resources. Therefore, individuals should find their optimal risk hedging,ṽ a 2 ½0; 1, in accordance with each population size under the following growth rate: where y 1 :¼ b 1 À kG; y 2 :¼ b 2 À 1 À k ð ÞG: ð50Þ r/K Selection Theory via Stochastic Analysis Eq (49) is referred as to the density-dependent two-resource utilization model, satisfying the conditions (i-iii) when EES is less than one (the reason for the model satisfying the condition (iii) is mentioned in a later section). This model generates the following adjoint Fokker-Planck Hamiltonian, Adaptive utilization and density effect Because our model is α = 1, the characteristic function of the adaptive species satisfies (Cf. Text C in S1 File). Substituting Eqs (35), (37) and (51) into Eq (52), we obtain the optimal risk hedgingṽ S under (see Text F in S1 File). Eq (53) implies that the K-strategyṽ K ðr Ã Þ is biased towards either one of the two strategies of utilization at the border κ = 0.5. Yet, in κ = 0.5, due to a counterbalance between competition rates in both resources, the r and K strategies have an identical utilization ratioṽ. In another situation (κ 6 ¼ 0.5), the K-strategy can evolve all patterns of utilization ratios in accordance with the value of κ, which is different from the r-strategy (see Fig 2); for instance, the evolution of R 2 -specialists never occurs in r-selection [9]. Accordingly, it is also possible that the K-strategy potentially has smaller values of IRNI and R 0 than the otherṽ r ðrÞ (see Figs 3 and 4). This result suggests that the use of IRNI and R 0 as the fitness must be restricted. As shown in Eq (32), the K-strategy always has a larger carrying capacity than any other strategy even if the IRNI and R 0 are potentially smaller than those obtained with other strategies (cf. Fig 5). Note that those simulations are collected to illustrate shapes of functions and properties (cf. ; these are exactly given: Eqs (32), (53), and Eq (F.15) at Text F in S1 File. For this reason, comparison with the exact values of those parameters in the nature is not so important to the purpose of analysis. Because augmenting EES reduces the effect of internal stochasticity on the carrying capacity (see Eq (44)), Eq (53) has the following limit lim r Ã "1ṽ If κ = κ 0 , an adaptive strategy does not exist because all strategies have identical carrying capacity. It is also shown that at the individual scale, Eq (31) becomes equal for all v 2 [0, 1]. On the other hand, when σ 2 = 0, we obtain an identical result from the deterministic case, such as lim s 1 #0ṽ ðr Ã ; GÞ. Generalists do not evolve in higher EES and in the absence of internal stochasticity. Therefore, risk hedging is important for the persistence of the species with smaller EES than those under internal stochasticity.

Discussion
The central point of r/K-selection theory is K-selection. This study proposes an interpretation of the benefits of K-strategies by adapting concepts of ESS to structured population models. Natural selection operates essentially through an identical parameter Eq (27) in both r-and Kselection. This parameter represents the reproductive contribution of all cohorts to the population dynamics. Our specific example raises the possibility that the K-strategy leads to opposite traits than that obtained using the r-strategy as Pianka suggested [3]. However, it also indicates that internal stochasticity is able to induce strategies that are different from Pianka's classification. The example is a K-strategy with a reproductive number identical to that of the r-strategy and with an IRNI lower than that of the other strategies, which is not included in Pianka's classification (Cf . Figs 3 and 4). The benefit of the K-strategy for individuals is that the population is not invaded by mutants with any other strategies. On the other hand, our study suggests that several traits of the K-strategies in Pianka's classification are not essentially adaptive strategies but are by-products of adaptation. Although recent studies do not focus on examining Pianka's classifications, it is still significant in any ecological research that considers adaptive traits of life history in the presence of density effects [4]. Various life histories of organisms are thought to result from adaptation to intraspecific competition. Both r-and K-strategies arise from identical consequences, of which individuals maximize their reproductive contribution to the population dynamics (Eq (11)) under different magnitudes of density effects. Density effects are, Fig 2. The K-strategy of the two-resource utilization model. This figure shows the shape of adaptive utilization under K-selection with respect to the ratio of the strength in the intraspecific competition κ. As the ratio increases, the adaptive strategy shifts to R 1 -specialists, generalists, and R 2 -specialists in turn. It is known that R 2 -specialists do not evolve in the absence of density effects [9]. The parameters are b 1 = 0.15, b 2 = 0.05, σ 1 = 0.6, σ 2 = 0.02, x = 0.01, x* = 1.5, ϕ(x*) = 10, and μ 0 = 0.01. doi:10.1371/journal.pone.0157715.g002 r/K Selection Theory via Stochastic Analysis however, unable to impair the r-strategy as a result of intraspecific competition. Consequently, species evolve the strategy that seems to be intuitively disadvantageous (e.g., altricity or low basic reproductive number). Though it seems that it is difficult for the strategy to evolve under r-selection, these disadvantages are not essential for the adaptive strategy in K-selection. Our results present a different viewpoint in r/K selection theory. As mentioned above, the reproductive traits of K-strategies are not always the purpose of adaptation, but they are by-products caused by at least the adaptive strategies. A significant criterion for the K-strategy is whether a strategy becomes ESS in terms of the carrying capacity or not; it should have the highest level of intraspecific competitionG. Then, the individual lifespan, reproductive number, maturity, and other vital rates are determined by what the population density affects in their life history. Density effects, thus, are able to influence not only mortality but also the other components of life history, and they alter the statistical characteristics of the total contributions of all cohorts Eq (28). Various traits of life histories, besides precocity and prolificacy, might arise due to the complicated system of density effects. In addition, this study proposes how we analyze the evolution of life histories under multiple density effects with internal stochasticity.
By contrast, conventional r-selection theory has thought that the evolution of the r-strategy occurs with high natural mortality as long as the intraspecific interaction is negligible. The rstrategyṽ r maximizes the IRNI at a sufficiently small population density; nevertheless, Eq (31) indicates that the population is eventually dominated by the K-strategyṽ K . Accordingly, the  15) at Text F in S1 File. The IRNI of the r-strategy is always greater than that of the others based on its definition. The domains of both the specialists and generalists represent the K-strategy because the r-strategy is independent of intraspecific competition. The parameters used are the same as those used for Fig 2. doi:10.1371/journal.pone.0157715.g003 r/K Selection Theory via Stochastic Analysis evolution of the r-strategy is unlikely to occur under assumptions (i-iii). For example, if an adaptive strategy is independent of the density effects, the r-strategy might remain as ESS at all population sizes because it combines both strategies (ṽ r ¼ṽ K ), such as κ = 0.5 in Eq (49). In contrast, the r-strategy appears to confer an advantage over a significant disturbance because of its higher population growth rate. As discussed in previous studies on this subject, the rstrategy possibly serves as plasticity with respect to the reduction of the population size by external stochasticity. Eqs (25) and (26) can only provide local adaptive control at each density effect level Γ Ã , which does not mean prohibiting a species from evolving to have plastic control at each population size. Therefore, the r-strategy is thought to be incorporated into plastic control, which switches from the r-strategy to the K-strategy associated with an increase in population size (such as the phase variation in migratory locusts [1]).

Conclusion
This study proposes unifying the theory on optimal life schedule and nonlinear structured population models using a function representing the reproductive contribution of each cohort. The function provides a solution to two equations (Eqs (25) and (26)). Applying these equations to r/K-selection theory, we can clarify the relationship between the individual scale and population scale in evolution. r/K-selection theory consists of a specific environment, where the population growth converges to a unique carrying capacity (Cf. assumptions (i)-(iii)). In nature, all The potential basic reproductive numbers in the rand K-strategies. This figure illustrates R 0 for both strategies in the absence of density effects. The horizontal axis shows the ratio of the strength in the intraspecific competition (κ), and the vertical axis represents R 0 under Γ t = 0. The solid and the dashed lines represent the same entities as those in Fig 3, and the parameters used are also the same. It is shown that the K-strategy of R 0 is not always higher than the other. Therefore, the basic reproductive number of the Kstrategy potentially depends on the ratio of the strength in the intraspecific competition.
doi:10.1371/journal.pone.0157715.g004 r/K Selection Theory via Stochastic Analysis habitats are not necessarily identical to such kinds of environment. Theoretical studies based on dynamical systems predict that density effects can generate multiple equilibria and periodic, or chaotic behaviors of population growth [15]. Furthermore, population dynamics are also thought to be affected by external stochasticity [16]. Because Eqs (25) and (26) can address the adaptive strategy surrounding each stable equilibrium, we can apply them to models that contain multiple equilibria. In this case, each stable equilibrium is supposed to give rise to the adaptive strategy that can be regarded as speciation such as Adaptive Dynamics theory. The difference between our method and AD is that it treats not the trait value but the dynamical controls of life history as the object in natural selection. This means that our method releases the specific adaptive phenotype from the concept of trait value in terms of lumping various life strategies in quantitative genetics. Moreover, we can derive a novel fitness function from three ingredients of life history consisting of state transition rate dX a , fertility F(Á), and mortality μ(Á) via Eq (25). Depending on density effect of each state (e.g., γ m (Á)), the population growth do not always have stable equilibria generally as which theoreticians predicted. In this case, we have few methods to handle it. The similar matter is also argued in the presence of external stochasticity [16][17][18]. Many species are thought to have various population dynamics besides the logistic type and to also be affected by external stochasticity. A further extension of this framework will be required that can account for adaptive strategies under unstable or fluctuating population dynamics in future work. Consequently, we hope that this study will become a basis for systematization in evolutionary theory addressing the effects described above.