In the past few decades, parameter estimation problem for the stochastic differential equation have been studied by many scholars whose results mostly base on the discrete observation. In order to get more accurate estimators, we should observe through continuous time. Deterministic model, which parameters are deterministic irrespective of environmental fluctuations, are usually used to describe the overall impact of changes between different factors. these models obviously impose limitations in mathematical modeling of whole real systems. However, in the real world, many random factors (i.e. earthquakes, typhoons, car accidents and other unforeseen factors) may make the parameters into random variables. Therefore, it is more reasonable to use the stochastic differential equation with to describe the real systems disturbed by random noises. For example, the stochastic logistic diffusion model has been widely used in the field of social life, application of stochastic logistic diffusion model has been used in the field of applied economics    , biology     , power engineering    and so on. Very recently, considerable research results have been reported on the parameter estimation based on discrete observation. To be special,  used the least squares method to estimate the parameters, also obtained the point estimators and confidence intervals as well as joint confidence regions.  used conditional least squares and weighted conditional least squares method to study the parameter estimation of two-type continuous-state branching processes with immigration based on low frequency observations at equidistant time points.  studied the asymptotic behaviour of parametric estimator for nonstationary reflected Ornstein-Uhlenbeck process by applying maximum likelihood estimation. For the stochastic logistic diffusion model,  worked out the optimization problem with respect to stationary probability density and provide a new equivalent, an ergodic method is used to show the almost surely equivalency between the time averaging yield and sustainable yield.  considered a stochastic logistic growth model involving both birth and death rates in the drift and diffusion coefficients and the associated complete Fokker-Planck equation is also established to describe the law of the process.  focused on stochastic dynamics involve continuous states as well as discrete events, and obtain weak convergence of the underlying system, and utilized the structure of limit system as a bridge to invest stochastic permanence of original system driving by a singular Markov chain with a large number of states.  presented some basic aspects of adequate numerical analysis for the random extensions such as numerical regularity and mean square convergence.  improved two mathematically tractable cases: at the limit of the number of individuals and at the limit of basic reproduction ratio. In the discrete observations, we let the time interval tends to 0 to get more accurate result. Therefore, this means that parametric inference based on continuous time observation is much more accurate in dealing with parameter estimation problem. During the estimation processing, two important theories have been used to estimate parameters based on continuous observation in the existing literature. One is denoting Radon-Nikodym derivative with likelihood ratio and the other one is by using stationary distribution of solution.
The stochastic logistic diffusion model can be described by the following stochastic differential equation:
where represents the population capacity at time t, represents the
natural birth rate, represents mortality, represents load capacity,
usually also represents the largest population that environmental resources can support. represents the dynamic effect of noise on . is a Wiener process modelling the random factor.  studied the existence, uniqueness and global attractively of positive solutions for model (1.1), and established a maximum likelihood estimator for the parameters.  proved that no matter how small , the solution will not explode in a limited time.
In this paper, the continuous observations shall be used to obtain more accurate results than discrete observations, and the likelihood ratio will be employed to get Radon-Nikodym derivative which can be used to solve the parameter estimation problem for logistic diffusion model. As we all know that logistic diffusion model is a diffusion process, for a general diffusion model , the parameter enter into the description of through or or both. However, the nature of diffusions allows us to evaluate exactly under a given continuous record, from the formula a.s. as (a.s. means almost surely) as pointed out in  . This result can be rewritten descriptively as a.s. as . Thus we could consider a parameter involved in as being know a single realization which means we should only consider the estimation of in . Then the likelihood ratio as Radon-Nikodym derivative and the expressions of all estimators would be obtained. As  have studied the stationary distribution of , based on the result and strong law of large numbers of martingales the law of iterated logarithm, the consistency property and the normality of asymptotic shall be proved for the estimation error.
This paper is organized as follows. In Section 2, a new method of estimating parameters is given and estimators are obtained. In Section 3, the strong consistency properties of estimators and estimation of asymptotic normality of error are proved. In Section 4, a numerical example for the estimators and error of estimation between estimators and trues is given to demonstrate the effectiveness of the proposed results. The conclusion is given in Section 5.
In this paper, the parameter estimation problem shall be studied for the logistic diffusion model described by a stochastic differential equation as given in (1.1). In this model , are unknown parameters. We can calculate use the method in  and this means that is also a known parameter. Because of the complication of the transitional density function for this model, it is difficult to obtain commonly used expression for the unknown parameters. Therefore, we will calculate the likelihood ratio (with respect to , and are the true parameters) for a finite set of time points , and then let the number of time points tend to infinity to get the likelihood function. From now on we shall work under the assumptions below.
Assumption 1: and are positive, is positive and independent with .
Assumption 2: , which implies that cannot reach zero.
Assumption 3: is a positive random variable, and there is a such that hold.
Next, the specific steps with respect to derivations of the likelihood function and parameter estimators are given below.
Assume that is known for all . can be observed continuously throughout the time interval . For observations in this detail enable the true diffusion parameter to be determinate through following result
For all , above equation can be rewritten descriptively as follows:
Then, may there be assumed know as:
The parameter and enter into the description of of (1.1), we can get exactly by (2.1). We begin with a class of probability , where the real stochastic process on evolves according to one of probability laws . For each define
with σ-field generated by sets , and B is a Borel set on R and the class N is -null set of .
From the theory of random time transformations, we have
where is a standard Brownian motion(or its measure). (1.1) with the initial distribution . is a Brownian motion with respect to the filtration and is a -measurable random variable with distribution . Then, for the process of option times (2.3), we can create the inverse process that is with filtration , and we can find a Brownian motion with respect to . For more details please see reference  , the existence and uniqueness of will be found. From Kailath and Zakai’s  , we know
by substituting for s,
Similarly, we can get,
Writing for the restriction of to , and we can now define following likelihood function as a Radon-Nikodym derivative:
Let , (the “log()” function has the basement “e”) solving following equation
we can obtain the estimators as follows:
3. Main Results and Proofs
In this section, we shall give the asymptotic distribution of the estimated errors and the corresponding proof. It’s easy to konow the solution of Equation (1.1) has following expression:
Firstly, let us give following four lemmas.
Lemma 1 (The law of iterated logarithm) 
Lemma 2  Assume that there is a positive number C such that
Then, for any given initial value , the solution of equation has following properties
Lemma 3  If the conditions and hold. Then, for any given initial value , the solution of has the property that
Lemma 4 Assume that is a solution to the stochastic differential Equation (1.1) and Assumptions 1 - 3 hold. Then, we have
proof: It is known from lemma 2 that the solution (3.1) obeys
While, from Assumption 2 we have . Then, by the properties in lemma 2, the solution (3.1) satisfies
By the Itô formula, it is easy to know that
Dividing both side by t and then letting , we obtain
The proof is complete.,
Remark: (The one-dimensional Itô formula) Let be an Itô process on with the stochastic differential
where and . Let .Then is again an Itô process with the stochastic differential given by
Theorem 1 Let be a solution to the stochastic differential Equation (1.1) and Assumptions 1-3 hold. Then, we have
proof: It follows from (1.1) that , dividing both sides by t and then letting , one has
is a matingale with zero mean with respect to the σ-algebra . Moreover, according to (1.1), . The equation (among them , ,
) gets maximum when , thus,
Therefore, is bounded. It then follows that
By the strong law of large numbers of martingales, we have
Together with (3.1) and Lemma 1, we obatain
By LHospital rule, we get
According to (3.10), we get
The proof is now complete.,
Remark: (LHospital rule) The general form of LHospital rule covers many cases. Let c and L be extended real numbers.The real valued function f and g are assumed to be differentiable on an open interval with endpoint c, and additionally
on the interval. It is also assumed that . Thus the
rule applies to situations in which the ratio of the derivatives has a finite or infinite limit, and not to situations in which that ratio fluctuates permanenty as x gets closer and closer to c.
Theorem 2 Under Assumptions 1-3, and are strongly consistent.
Proof: Substituting (1.1) into the expression of yields
Letting , and according to Lemma 4 and Theorem 1, we have
It follows from (3.12) that
Substituting (1.1) into the expression of yields
Similarly, we have
It is easy to get from (3.12) that
Thus, and are strongly consistent. The proof is complete.,
Theorem 3 Under Assumptions 1-3, we have
Proof: It follows from (3.15) that
Substituting (3.9) and (3.11) into the above expression and then letting , we have
According to (3.13), one has
Similarly, it follows easily from (3.17) that
Substituting (3.9) and (3.11) into the above expression and then letting yields
The proof is complete.,
In this section, a numerical simulation example shall be presented to demonstrate the effectiveness of the approach results.
The simulation is based on (2.6) (2.7) and (2.12). First according to (2.6) and (2.7), for given values of and such as and , we can get the sample values based on the likelihood ratio estimation and MATLAB. Then, for substituting the sample values into (2.12), the values of can be obtained. Subsequently, we calculate the average values of the estimators. Finally, the average errors between estimators can also be calculated. Simulation results are shown in Table 1. In Table 1, the time is represented as “t” and likelihood ratio estimator is shown as “LR”. Table 1 lists the values of “ “, “ “ and the average errors of “LR”. Table 1 illustrates that the average errors of depended on the size of given value of . But under the hypothesis of normal distribution, obvious difference can not be found between estimators and true values, estimators are good. From the Table 1 we can see clearly that the estimators become more and more close to the true value by increasing the time t. The was of continuous time estimation is better than the discrete observation  . (Those data comes from web of Statistical Data and I use MATLAB to simulate those data to get the result. The confidence intervals is for , for ; for , for ; for , for ; for , for ;)
Table 1. Likelihood ratio estimator simulation results of .
In this paper, parameter estimation problem has been studied for the continuous time stochastic logistic diffusion model by using likelihood ratio. The explicit expressions for the estimation errors have been given and the according asymptotic properties have been proved by applying the he law of iterated logarithm, random time transformations, stationary distribution of solutions of stochastic differential equations and the law of large numbers for martingales. To get more accurate results, we use continuous observation method, and the proposed estimators are closer to the true value that be demonstrated by a simulation example. In the future research, we will consider the state estimation problem for non-linear systems with incomplete observation  and non-linear systems with random disturbance caused by Levy jump or Poisson jump  .
*This work was supported in part by the National Nature Science Foundations of China under Grant No. 61673103 and No.61403248.
 Roman-Roman, P. and Torres-Ruiz, F. (2012) Modelling Logistic Growth by a New Diffusion Process: Application to Biological Systems. Biosystems, 110, 9-21.
 Moller, J., Bergmann, K., Christiansen, L. and Madsen, H. (2012) Development of a Restricted State Space Stochastic Differential Equation Model for Bacterial Growth in Rich Media. Journal of Theoretical Biology, 305, 78-87.
 Chen, J. and Chang, W. (1998) Modeling Differential Diffusion Effects in Turbulent Nonreacting/Reacting Jets with Stochastic Mixing Models. Combustion Science and Technology, 133, 343-375.
 Pan, J., Gray, A., Greenhalgh, D. and Mao, X. (2014) Parameter Estimation for the Stochastic SIS Epidemic Model. Statistical Inference for Stochastic Processes, 17, 75-98.
 Zang, Q. (2016) Asymptotic Behaviour of Parametric Estimation for Nonstationary Reflected Ornstein? CUhlenbeck Processes. Journal of Mathematical Analysis and Applications, 444, 839-851.
 Zou, X. and Wang, K. (2014) Optimal Harvesting for a Stochastic Regime-Switching Logistic Diffusion System with Jumps. Nonlinear Analysis Hybrid Systems, 13, 32-44.
 Li, X. and Yin, G. (2016) Switching Diffusion Logistic Models Involving Singularly Perturbed Markov Chains: Weak Convergence and Stochastic Permanence. Stochastic Analysis and Applications, 35, 364-389.
 Mao, X., Marion, G. and Renshaw, E. (2002) Environmental Brownian Noise Suppresses Explosions in Population Dynamics. Stochastic Processes and Their Applications, 97, 95-110.
 Kailath, T. and Zakai, M. (1971) Absolute Continuity and Radon-Nikodym Derivatives for Certain Measures Relative to Wiener Measure. Annals of Mathematical Statistics, 42, 130-140.
 Liu, M. and Zhang, D. (2015) A Dynamic Logistic Model for Medical Resources Allocation in an Epidemic Control with Demand Forecast Updating. Journal of the Operational Research Society, 67, 841-852.
 Sun, X. and Wang, Y. (2008) Stability Analysis of a Stochastic Logistic Model with Nonlinear Diffusion Term. Applied Mathematical Modelling, 32, 2067-2075.