arXiv:1706.07034v1 [math.ST] 21 Jun 2017

LOCAL BANDWIDTH SELECTION FOR KERNEL DENSITY ESTIMATION
IN BIFURCATING MARKOV CHAIN MODEL
S. VALÈRE BITSEKI PENDA AND ANGELINA ROCHE

Abstract. We propose an adaptive estimator for the stationary distribution of a bifurcating
Markov Chain on Rd . Bifurcating Markov chains (BMC for short) are a class of stochastic
processes indexed by regular binary trees. A kernel estimator is proposed whose bandwidth is
selected by a method inspired by the works of Goldenshluger and Lepski [18]. Drawing inspiration
from dimension jump methods for model selection, we also provide an algorithm to select the
best constant in the penalty.

Keywords: Bandwidth selection ; Bifurcating autoregressive process ; Bifurcating Markov
chains ; Binary trees ; Nonparametric kernel estimation.
Mathematics Subject Classification (2010): 62G05, 62G10, 62G20, 60J80, 60F05, 60F10,
60J20, 92D25.
1. Introduction
First introduced by Basawa and Zhou (2004) [2], Bifurcating Markov chains models (BMCM
for short) has recently received particular attention for its application to cell lineage study. Guyon
(2007) [21], have proposed such a model to detect cellular aging in Escherichia Coli and proved
laws of large numbers and central limit theorem for this class of stochastic process. BitsekiPenda et al. (2014) [7] have then completed these asymptotic results and proved concentration
inequalities.
To the best of our knowledge, kernel density estimation for the BMCM were considered first
by Doumic & al. [17], where they estimate the division rate of population of cells reproducing by
symmetric division, i.e. cell reproduction where each fission produces two equal daughter cells.
After this work, Bitseki & al. [8], have used the wavelets methodology to study the nonparametric
estimation of the density of BMCM. They propose an adaptive estimator in dimension 1. Recently,
Bitseki and Olivier [9] have studied the Nadaraya-Watson type estimators of a BMCM that they
called nonlinear bifurcating autoregressive process. The latter model can be seen as an adaptation
of nonlinear autoregressive process on binary regular tree. We mention that, except in [8], all the
estimations done in the previous works are non adaptive. In particular, the question of data-driven
bandwidth selection was not addressed in [17] and [9]. The main objective of this work is then to
propose a data-driven method for choosing the bandwidth for the kernel estimator of the invariant
measure in the multi-dimensional BMCM, following ideas from the works of Goldenshluger and
Lepski (2011) [18].
The idea of the method is to select the bandwidth minimizing an empirical criterion imitating
the bias-variance decomposition of the risk of the kernel estimator. More precisely, let H be a
collection of bandwidths and let (b
νh )h∈H be a family of kernel estimators of an unknown density
ν. Then, we select the bandwidth b
h as
b
b
h = arg min{A(h)
+ bV (h)}
h∈H

1

2

S. VALÈRE BITSEKI PENDA AND ANGELINA ROCHE

with A(h) an empirical version of the bias of the estimator νbh and V (h) a penalty term with the
same order than the variance. The now so-called Goldenshuger-Lepski methodology ([18], but also
[19, 26, 27, 24, 22]), initially developed for density estimation, has been applied in many contexts
such as deconvolution problems [14], conditional cumulative distribution function estimation [11,
10], regression problems [10, 12], conditional density estimation [10, 5], hazard rate estimation [10],
white noise model [23], kernel empirical risk minimization (including robust regression) [13], Lévy
processes [3], Cox model [20], stochastic differential model [16]. Under suitable assumptions on the
kernel, it is shown in [18] that this selection rule leads to a minimax adaptive estimator on a general
class of regular functions, for a general class of Ls -risks. Pointwise versions of the GoldenshlugerLepski selection rule have been less considered (Comte et Lacour (2013) [14], Rebelles (2015) [27]
and Chagny and Roche (2016) [12]). The interest of such an approach is that the bandwidth is
selected with a local criterion which realizes the best bias-variance compromise at the point where
the estimator is calculated. On the contrary, integrated versions of the Goldenshluger-Lepski
selection rule select the same bandwidth at all points. Let us mention that, in all the articles
cited above, the theoretical results rely on concentration inequalities for sums of i.i.d. random
variables such as Bernstein Inequality or Talagrand Inequality. In our context, such results are
not applicable. Hence, we prove a Bernstein-type Inequality for functionals of BMC, where the
functions are kernels and convolution of kernels. Compare to those obtained in [8], our inequalities
are more complete in the sense that the deviation parameter can take all the positive values. More
precisely, its values do not depend on the size of the samples, with is essential for our theoretical
results.
Ideally, the penalty term V , called “the majorant” by Goldenshluger and Lepski, depends
entirely on the kernel and the observations (it does not depend on the density ν). The selected
estimator is then νbbh . However, for the BMCM, the variance term V in the previous selection rule
contains a term which may depend on the unknown density ν. Moreover, this term is generally not
estimable from observations. To resolve this problem, we propose a modification of GoldenshlugerLepski rule’s selection, inspired by the works of Lacour and Massart (2016) [22]. As suggested in
[22], the constant term in bV (x, h) is then selected automatically from the data with an algorithm
inspired by the works of Arlot and Massart (2009) [1].
The paper is organized as follows. The model is defined in Section 2. Section 3 is devoted to
the definition of the estimator. In Section 4, we provide a numerical study of our estimator. The
proofs are given in Section 5.
2. Definitions
We are now going to give a precise definition of a BMC. First we note that this class of stochastic
processes has been introduced by Guyon [21] in order to understand the mechanisms of cell division.
Indeed, these stochastic processes are well adapted to study a population (or more generally, any
dynamic system) where each individual (or more generally, each particle) in one generation gives
birth to two individuals in the next one. In the sequel, we will then use the language of the
population dynamic to define the sets of interest.
Let (Ω, F, (Fm , m ∈ N), P) be a filtered probability space. Let (Xu , u ∈ T) be a sequence of
random variablesSdefined on (Ω, P), taking values in Rd , where d ≥ 1, and indexed by the infinite
∞
binary tree T = m=0 {0, 1}m , with the convention that {0, 1}0 = ∅. We equip Rd with its usual
Borel σ-field. Now, we will see T as a given population. Then each individual u of this population
is represented by a sequence of 0’s and 1’s, and has two descendants, u0 and u1. The initial
individual of the population is ∅. For all m ∈ N, let Gm be the set of individuals belonging to the
m-th generation, and Tm the set of individuals belonging to the m first generations. We have:

BANDWIDTH SELECTION FOR KERNEL DENSITY ESTIMATION IN A BMCM

Gm = {0, 1}m , Tm =

m
[

Gq and T =

q=0

[

3

Gm .

m≥0

For an individual u ∈ Gm , we set |u| := m its length (i.e. the generation to which it belongs).
2.1. Bifurcating Markov chain [21].
Definition 1 (T-transition probability). Let P : Rd × B((Rd )2 ) → [0, 1] (with B(Rd )2 = B(Rd ) ⊗
B(Rd ) the usual product σ-field on (Rd )2 ). Then P is a T-transition probability if
• x → P(x, A) is measurable for all A ∈ B(Rd )2 .
• A → P(x, A) is a probability measure on ((Rd )2 , B(Rd )2 ) for all x ∈ Rd .
For a B(Rd )3 -measurable function f : (Rd )3 → R, we denote (when it is defined) by Pf the
B(Rd )-measurable function
Z
x ∈ Rd 7→
f (x, y, z)P(x, dy, dz).
(Rd )2

Definition 2 (Bifurcating Markov chain). Let µ be a probability measure on (Rd , B(Rd )) and P a
T-transition probability. We say that (Xu , u ∈ T) is a (Fm )-bifurcating Markov chain with initial
distribution µ and T-transition probability P (denoted T-BMC in the sequel) if
• Xu is Fm measurable for all u ∈ Gm .
• X∅ has distribution µ.
• For all m ∈ N, and for all family (fu , u ∈ Gm ) of B(Rd )-measurable functions from (Rd )3
to R,
"
#
Y
Y
E
fu (Xu , Xu0 , Xu1 )|Fm =
Pfu (Xu ),
u∈Gm

u∈Gm

where u0 := (u, 0) ∈ Gm+1 and u1 := (u, 1) ∈ Gm+1 .
2.2. Tagged-branched chain [21, 9]. Let (Xu , u ∈ T) be a T-BMC with initial distribution µ and
T-transition probability P. We denote by P0 and P1 respectively the first and second marginals
of P. More precisely
P0 (x, B) = P(x, B × Rd ) and P1 (x, B) = P(x, Rd × B),
for all x ∈ Rd and all B ∈ B(Rd ). Let Q be the mixture of P0 and P1 with equal weights
1
1
P0 + P1 .
2
2
on Rd with initial value Y0 := X∅ and transition probability
Q=

The Markov chain Y := (Ym )m∈N
Q is called the tagged-branch chain.
In all the paper, we will denote by Qm the mth iterated of Q recursively defined by the formulas
Z
Q0 (x, ·) = δx and Qm+1 (x, B) =
Q(x, dy)Qm (y, B) ∀B ∈ B(Rd ).
S
m

It is well known that Q is a transition probability in (Rd , B(Rd )). In particular, we have
E[f (Ym )] = µQm f for all measurable function f : Rd → R.
In the sequel, we will assume that the Markov chain Y is ergodic- that is say, that there exists
a unique distribution ν on (Rd , B(Rd ) such that, for all measurable function f : Rd → R,

4

S. VALÈRE BITSEKI PENDA AND ANGELINA ROCHE

Z
lim E[f (Ym )] =

m→∞

f dν.
Rd

We will also assume that the distribution ν has a density, that we also denote by ν, with respect
to the Lebesgue measure.
As previous works have shown (see for example [21, 7]), the analysis of a BMC (Xu , u ∈ T)
is strongly related to the asymptotic behavior of the tagged-branched chain (Ym , m ∈ N), and
therefore to the knowledge of the invariant distribution ν. We stress that this distribution is
unknown and it is not directly observable, in such a way that its estimation from the data is of
great interest. The aim is to estimate ν from the observation of a subpopulation (Xu , u ∈ Tn ).
3. Estimation of the stationary distribution ν
3.1. Definition of the estimator. We suppose that we observe the process (Xu , u ∈ T) up to
the n-th generation. We denote by |Tn | = 2n+1 − 1 the cardinality of Tn . Based on the observation
of (Xu , u ∈ Tn ), we propose the following estimator of ν
(1)

νbh (x) =

1 X
Kh (x − Xu )
|Tn |

∀x ∈ Rd ,

u∈Tn

t

d

where for all x = (x1 , ..., xd ) ∈ R and for all h = (h1 , ..., hd )t ∈]0, +∞[d
Kh (x) = K(x1 /h1 , ..., xd /hd )/

d
Y

hj

j=1

R
and K is a kernel, that is to say a function K : Rd → R which verifies Rd K(u)du = 1.
The vector h is usually called a bandwidth. In the sequel, |h| will denote the product h1 ×· · ·×hd .
It is well known in the kernel estimation theory that the choice of h is of great interest. Indeed
the amount of smoothing is controlled by a judicious choice of the bandwidth. Now, in order to
tackle the issue of the choice of this bandwidth, we will need the following assumptions.
Assumption 3 (Assumptions on the kernel).
Z
Z
kKk1 =
|K(t)|dt < +∞, kKk22 =
|K(t)|2 dt < +∞
Rd

Rd

and

kKk∞ = sup |K(t)| < +∞.
t∈Rd

Assumption 4 (Uniform geometric ergodicity condition). There exists two constants ρ ∈ (0, 1/2)
and M > 0 such that for all bounded ν-integrable function g, for all x ∈ Rd and m ≥ 0
Z
|Qm g(x) −
gdν| ≤ M kgk∞ ρm .
Rd

This assumption is verified for severals models of BMC. We refer for example to [9] where
Bitseki and Olivier have shown it for NBAR processes (see [9] Lemma 20). For a precise definition
of NBAR processes, we refer to section 4.
Assumption 5 (Assumption on P, P0 , P1 , Q and ν). We assume that the transitions P, P0 , P1
and Q admit densities with respect to the Lebesgue measure that we denote with the same notations.
Moreover, we assume that
kPk∞ < +∞, kP0 k∞ < +∞, kP1 k∞ < +∞, kνk∞ < +∞, and kQk∞ < +∞.

BANDWIDTH SELECTION FOR KERNEL DENSITY ESTIMATION IN A BMCM

5

3.2. Bias-variance decomposition. We consider a pointwise quadratic risk
E[(b
νh (x) − ν(x))2 ],
where x ∈ Rd is a given point.
The results of [9] (see the proof of Proposition 21) allows to obtain the following upper-bound
on the risk, where we have make explicit the constants that appear.
Proposition 6. Under Assumption 3 to Assumption 5, we have
E[(b
νh (x) − ν(x))2 ] ≤ 2(Kh ∗ ν(x) − ν(x))2 + 2

(2)

where ∗ denotes the convolution product f ∗g(x) =
over Rd and

R
Rd

C(P, ν)
,
|Tn ||h|

f (x−t)g(t)dt for all functions f, g integrable

CI
C(P, ν) = √
( 2 − 1)2
with
CI = (1 +

1
)(kQk∞ + kνk∞ )2 + M 2 + CP
1 − 2ρ2

and
CP = 2kKk22 (kQk∞ + kνk∞ ) + kPk∞ + kνk∞ (kKk22 + kP0 k∞ + kP1 k∞ )
Inequality (2) can be seen as a bias-variance decomposition in the sense that
Z
1 X
Eν [b
νh (x)] =
Eν [Kh (x − Xu )] =
Kh (x − t)dν(t) = Kh ∗ ν(x),
|Tn |
Rd
u∈Tn

where Eν is the expectation
with respect to the

 measure ν and C(P, µ)/|Tn ||h| is an upper-bound
on the variance term E (b
νh (x) − Kh ∗ ν(x))2 .
3.3. Selection rule. Given a family of bandwidths Hn ⊂ [0, +∞[d , we have a family of estimators
(b
νh )h∈Hn and the aim is to select an estimator in this family with risk close to the unknown oracle
risk




E (b
νh∗ (x) − ν(x))2 with h∗ = arg minh∈Hn E (b
νh (x) − ν(x))2 .
Imitating the decomposition given in Equation (2) we could consider the following selection rule
inspired by the work of [18]
b
h = arg minh∈Hn {A(x, h) + bV (x, h)}

(3)
where

• Hn ⊂ [0, +∞[d is a finite collection of bandwidths;

• A(x, h) = maxh0 ∈Hn (b
νh0 (x) − Kh ∗ νbh0 (x))2 − aV (x, h0 ) + with b ≥ a ≥ 1;
• V (x, h) = C(P, µ) log(|Tn |)/|Tn ||h|.
The constants a and b can be different, as suggested in Section 5 of [22]. In Section 3.4, we
provide a method to choose both a and b, as well as the quantity C(P, µ) appearing in the variance
term V (x, h). The lower bound a = 1 is a critical minimal value, in the sense that the procedure
fails for lower values of a (see [22]).
We prove the following oracle-inequality on the selected estimator νbbh .

6

S. VALÈRE BITSEKI PENDA AND ANGELINA ROCHE

Theorem 7. Under Assumption 4, if minh∈Hn |h| ≥ log(|Tn |)/|Tn |, then there exists a minimal
value amin > 0 independent of n, such that, for all a > amin ,




log(|Tn |)
C2
(4)
E (b
νbh (x) − ν(x))2 ≤ C1 min Bh (x) +
+
,
h∈Hn
|Tn ||h|
|Tn |
where C1 , C2 > 0 do not depend on n nor x and
Bh (x) = max
(Kh0 ∗ ν(x) − Kh ∗ Kh0 ∗ ν(x))2 .
0
h ∈Hn

Remark 8. The form of the bias term Bh (x) in Inequality (4) is very similar to the one obtained
for pointwise adaptive kernel density estimation in [27, Theorem 1]. It can be replaced by an
upper-bound, e.g. kKk1 kν − Kh ∗ νk∞ , coming from the Young’s inequality, as in [14, Theorem 2].
Once the previous theorem is proved, an immediate corollary follows
Corollary 9. Suppose that the assumptions of Theorem 7 are verified and that ν ∈ Σ(β, L, D)
where β = (β1 , ..., βd ) ∈]0, +∞[d and Σ(β, L, D) is the set of all β-Hölder densities on the open set
D ⊆ Rd i.e. the set of all functions f : D → R which admits, for all j = 1, ..., d, partial derivatives
with respect to xj up to the order bβj c and verifies, for all x = (x1 , ..., xd ) ∈ D, x0j ∈ R such that
(x1 , ..., xj−1 , x0j , xj+1 , ..., xd ) ∈ D,
∂ bβj c f
bβ c
∂xj j

(x1 , , ..., xj−1 , x0j , xj+1 , ..., xd ) −

∂ bβj c f
bβ c
∂xj j

(x) ≤ L|xj − x0j |βj −bβj c .

RMoreover, suppose
R jthat K is a kernel of order ` ∈ N (with ` ≥ maxj=1,...,d {βj }) that is to say,
K(t)dt
=
1,
x K(x) = 0 for all j = 1, ..., ` and that
R
R
Z
|x|` |K(u)| < +∞.
R

Suppose also that, for all n, there exists h∗ = (h∗1 , ..., h∗d ) ∈ Hn such that h∗j 



log(|Tn |)
|Tn |

β̄/(βj (2β̄+d))

where β̄ = d/(1/β1 + ... + 1/βd ) is the harmonic mean of β. Then, there exists a constant C > 0
such that

−2β̄/(2β̄+1)


|Tn |
sup E (b
νbh (x) − ν(x))2 ≤ C
.
log(|Tn |)
x∈D
This corollary is a direct consequence of Theorem 7 and [28] Proposition 2.1.
Remark that the assumption on the bandwidth collection is verified e.g. by
n
od
Hn = hmax k −α , k = 1, ...., (|Tn |hmax / log(|Tn |))1/α
for all constant hmax > 0 and all α > 1.
3.4. Estimation of the constant appearing in the variance. Due to the term C(P, µ) which is
hardly estimable, the variance term V (x, h) is not calculable in practice. We propose an algorithm
to estimate it, based on slope estimation, as developed in [1, 22].
As suggested in [22], we take b = 2a, hence in order to calculate the estimator, it is sufficient to
have a good value of κ = aC(P, µ). The following algorithm is inspired by the procedure described
in [1], for model selection purposes. However, note that, in our case, both terms A(x, h) and
V (x, h) depends on the constant κ, which is not the case in model selection contexts where only
the penalty term depends on the constant. The selection of the grid of κ is then different here.

BANDWIDTH SELECTION FOR KERNEL DENSITY ESTIMATION IN A BMCM

7

Algorithm 10.
1. Initialization : set
(5)

κm =

|Tn |
max |h0 |(b
νh0 (x) − Kh ∗ νbh0 (x))2 and κ1 = 0.
log(|Tn |) h,h0 ∈Hn

2. While s ≤ smax
j−1
(κm − κ1 ), for all j =
(i) Generate a sequence (κj )1≤j≤m such that and κj = κ1 + m−1
1, ..., m.
(ii) Calculate b
hj := b
h(κj ) as the minimizer of the criterion (3) with a = κ/C(P, µ) and
b = 2a.
(iii) Set
jjump = arg maxj=1,...,m−1

1
1
−
,
b
b
|hj | |hj+1 |

κ1 = κjjump , κ1 = κjjump +1 and s = s + 1.
3. Return b
hjjump +1 .
The algorithm search the value of κ for which the variance of the estimator increases significantly
and select a slightly larger value. This allows to select an estimator with a reasonable variance. The
same reasoning has given rise to the so-called dimension jump method for model selection purposes
(see [1]). The chosen value for the initialization of κm comes from the following reasoning. Setting
κ ≥ κm as suggested: by definition, for all h, h0 ∈ Hn ,
|Tn |h0
(b
νh0 (x) − Kh ∗ νbh0 (x))2 ≤ κ
log(|Tn |)
which implies that A(x, h) = 0 for all h ∈ Hn and that the criterion (3) will select the smaller
bandwidth in Hn . On the contrary, if κ > κm , let h, h0 ∈ Hn for which the minimum is attained
in (5), we have
|Tn |
(b
νh0 (x) − Kh ∗ νbh0 (x))2 > κ
log(|Tn |h0 )
which implies that A(x, h) > 0 and the criterion may select a bandwidth which is not the smaller
one. Hence, the values of κ for which the criterion (3) may select suitable values of the bandwidth
can not be greater than κm . That is the reason why we consider an initial grid in the interval
[0, κm ]. However, the initial value κm may be very large compared to the optimal value of κ. The
loop 2. allows to search among small values of κ while avoiding the choice of a too large m which
could be very expensive in terms of computation time. In practice smax = 2 and m = 20 seems to
be a reasonable choice.
Remark 11. The choice b = 2a is arbitrary since, in theory, any value b such that b ≥ a might
work as soon as a sufficiently large. Other choices may be made such as b = a which is the usual
choice of Lepski’s method, b = a + ε, with ε > 0 or b = (1 + ε)a.
4. Simulations
We shall now illustrate the previous results on simulated data coming from various models of
bifurcating Markov chains.

8

S. VALÈRE BITSEKI PENDA AND ANGELINA ROCHE

4.1. Bifurcating autoregressive processes.
Bifurcating autoregressive processes (BAR, for short) were first introduced by Cowan and
Staudte [15] in order to study the data from cell division, where each individual in one generation gives birth to two children in the next generation. This model has been widely studied over
the last thirty years (see for e.g. [9] and references therein). Recently, Bitseki and Olivier in [9]
have proposed an extension of BAR process initially introduced by Cowan and Staudte. Their
model is defined as follows.
Let Xu ∈ R be a quantitative data associated to the cell u ∈ T, for example the growth rate of
E. Coli. Then the quantities Xu0 and Xu1 associated to u0 and u1 the two children of u are linked
to Xu through the following autoregressive equations

 Xu0 = f0 (Xu ) + εu0 ,
(6)
L(X∅ ) = µ, and for u ∈ T,

Xu1 = f1 (Xu ) + εu1 ,

where µ is a distribution probability on R and f0 , f1 : R ; R. The noise (εu0 , εu1 ), u ∈ T forms
a sequence of independent and identically distributed bivariate centered random variables with
common density gε on R × R. The process (Xu , u ∈ T) defined by (6) is a bifurcating Markov
chain with T-transition probability
P(x, dy, dz) = gε (y − f0 (x), z − f1 (x))dydz.
Under some assumptions on µ, f0 , f1 and gε , it has been shown in [9] that the process X satisfies
all the good properties needed for our theoretical results (we refer to [9] for more details). We note
that the previous model can be seen as an adaptation of nonlinear autoregressive model when the
data have a binary tree structure. Furthermore, the original BAR process in [15] is defined for
linear link functions f0 and f1 with f0 = f1 .
Now, for our numerical illustrations, we build a BAR process living in S := [0, 1] as follows. First,
we choose X∅ such that L(X∅ ) = Beta(2, 2), where Beta(2, 2) is the standard Beta distribution
with shape parameters (2, 2). Then for u ∈ T and conditionally on Xu = x, we construct Xu0
and Xu1 independently in such a way that P(Xu0 ∈ dy, Xu1 ∈ dz) = P(x, y, z)dydz, where
P(x, ·, ·) := P(x, ·) ⊗ P(x, ·) and
P(x, y) := (1 − x)

y 2 (1 − y)
y(1 − y)2
+x
,
B(2, 3)
B(3, 2)

x, y ∈ [0, 1]

with B(α, β) the normalizing constant of a standard Beta distribution with shape parameters α
and β. Now, one can prove that this process is stationary, it has an explicit invariant density, which
is crucial to evaluate the quality of estimation of our method: this is a standard Beta distribution
with shape parameters (2, 2). One can also prove that
E [Xu0 |Xu ] = E [Xu1 |Xu ] = 1/5Xu + 2/5,
in such a way that the equations (6) are satisfied with f0 (x) = f1 (x) = 1/5x + 2/5 (for more
details, we refer for e.g. to [25]). Now, it is no hard to verify that this process satisfies our required
assumptions.
We simulate the n first generations of the process X, with n = 10 (hence the size of Tn is
2
|Tn | = 210 = 1024). We consider the Gaussian kernel K(t) = √12π e−t /2 . The results are given in
Figure 1.

BANDWIDTH SELECTION FOR KERNEL DENSITY ESTIMATION IN A BMCM

9

Figure 1. Plot of 10 estimators (red dashed lines) obtained from independent
copies of X. The black solid line represents the function ν to estimate.
4.2. Estimation of splitting rate in a growth-fragmentation model.
We are now interested in growth-fragmentation models. Theses models describes (for e.g.) the
evolution of cells which grow and divide randomly over time (see for e.g. [17] and references
therein). The model we are going to study is a simplification of the one studied in [17]; it is
defined as follows. Let S be a subset of [0, ∞) and let B : S 7→ [0, ∞) be a continuous function
(the splitting rate). Each cell u ∈ T grows exponentially with a common rate τ > 0 and when it
reaches a certain size x, it splits at rate B(x), and gives birth to two offspring (u0 and u1) of size
x/2. Next, this two offspring, u0 and u1, start a new life independently of each other. Clearly,
the process (Xu , u ∈ T), where Xu is the size of the cell u at birth, is a BMC. It is proved in [17]
that the T-transition probability of this BMC is P(x, ·, ·) := P(x, ·) ⊗ P(x, ·), where the density of
P(x, ·) is given by
!
Z y
B(2y)
B(2z)
P(x, y) :=
exp −
dz 1{y≥x/2} ,
τy
τz
x/2
for x ∈ S and y ∈ S/2. It can also be seen that the probability transition of the tagged-branch
chain is Q = P.
Doumic et al. [17] have proved that P admits an invariant probability measure ν having a
density, that we still denote by ν(·), with respect to the Lebesgue measure. It is also known from
[17] that the rate function B(·) and the invariant density ν(·) verify
B(x) =

τ x νB (x/2)
Rx
2 x/2
νB (z)dz

in such a way that a natural estimator for B(x), based on the observation of (Xu , u ∈ Tn ) is
bn (x) = τ x
B
2

1
|Tn |

νbh (x/2)

,
u∈Tn 1{x/2≤Xu <x} ∨ $n

P

where νbh is the kernel estimator defined in (1) and $n is a threshold which ideally go to 0.
Moreover, Bitseki et al. [8] have proved that under suitable assumptions on the splitting rate
B(·), the process X satisfies all the good properties needed for our theoretical results. For all the
previous assertions, we refer to [17, 8] for more details. The strategy is then to use our results for
the estimation of νB (x/2) for all x ∈ S.

10

S. VALÈRE BITSEKI PENDA AND ANGELINA ROCHE

30

25

20

15

10

5

0
1.5

2

2.5

3

3.5

4

4.5

5

Figure 2. Plot of 1 estimator (red solid line) obtained from a copy of X. The
yellow solid line represents the function B to estimate.
Now, for our numerical illustrations, we will work with the splitting rate using in [8]. We choose
τ = 2, S = (0, 5) and for all x ∈ S, B has the form

x
B(x) =
+ 3T 2(x − 27 )
x−5
where T (x) = (1 + x)1{−1≤x<0} + (1 − x)1{0≤x≤1} is a tent shaped function. With this choice of
B, the required assumptions for our theoretical results are satisfied (see [8]).
We simulate the n first generations of the process X, with n = 15. The results are given in
Figure 2.
5. Proofs
The proof rely on the lemma below, which is a Bernstein-type inequality.
Lemma 12. Let (Xu , u ∈ Tn ) be a bifurcating Markov chain on Rd with initial distribution µ and
T-transition probability P . Under the assumption of uniform geometric ergodicity, we have for all
δ>0

(7)

P

!
X
1
(Kh ∗ Kh0 (x − Xu ) − Eν [Kh ∗ Kh0 (x − Xu )]) > δ
|Tn |
u∈Tn


!
2
0
δcK,Q,ν,M c0ρ
δ |Tn ||h |
 ,
≤ 2 exp 4c kKk kKk δ
exp − 
1
∞
ρ,M
4c
kKk kKk∞ δ
0
+ cK,Q,ν,M cρ
2 cK,Q,ν,M c0ρ + ρ,M 3 1
3

where
M (1 + ρ) 0
2
, cρ = 3 +
,
1 − 2ρ
1 − 2ρ
cK,Q,ν,M = 8 max{2kKk21 kKk22 (kQk∞,∞ + kνk∞ ); max{kQk∞,∞ + kνk∞ ; M kKk1 kKk∞ }2 }

cρ,M =

We also have for all δ > 0,

BANDWIDTH SELECTION FOR KERNEL DENSITY ESTIMATION IN A BMCM

(8)

P

11

!
X
1
(Kh (x − Xu ) − Eν [Kh (x − Xu )]) > δ
|Tn |
u∈Tn


!
2
δc0K,Q,ν,M c0ρ
δ
|T
||h|
n

≤ 2 exp 4c kKk δ
exp − 
∞
ρ,M
4c
kKk δ
0
0
0
0
+
c
c
2 cK,Q,ν,M cρ + ρ,M 3 ∞
K,Q,ν,M ρ
3

where c0K,Q,ν,M = 8 max{M kKk∞ ; (kQk∞,∞ + kνk∞ )kKk1 ; (kQk∞,∞ + kνk∞ )kKk22 }.
Remark 13. As mentioned above, these inequalities are more complete than those obtained in [8],
since the deviation parameter δ does not depend on the size of the data. We stress that this fact is
essential for our theoretical results.
Proof. We will do the proof of (7). The proof of (8) follows the same lines.
Let λ > 0 and δ > 0. By Chernoff inequality, we have
!
1 X
(Kh ∗ Kh0 (x − Xu ) − Eν [Kh ∗ Kh0 (x − Xu )]) > δ
P
|Tn |
u∈Tn
"
≤ exp (−δλ|Tn |) E exp λ

!#
X

g(Xu )

,

u∈Tn

where the function g is defined by
g(y) = Kh ∗ Kh0 (x − y) − Eν [Kh ∗ Kh0 (x − X∅ )] .
For all u ∈ Gn−1 , we have, on the one hand
|g(Xu0 ) + g(Xu1 ) − 2Qg(Xu )| ≤ 2M (1 + ρ)kKh ∗ Kh0 k∞ .
Using the Young’s inequality, we have
kKh ∗ Kh0 k∞ ≤ kKk1 kKk∞ /|h0 |
and therefore
|g(Xu0 ) + g(Xu1 ) − 2Qg(Xu ))| ≤ 2cρ,M kKk1 kKk∞ /|h0 |.
On the other hand, we have
h
i
2
E (g(Xu0 ) + g(Xu1 ) − 2Qg(Xu )) |Fn−1 ≤ 4kQk∞,∞ kKh ∗ Kh0 k22 ≤ cK,Q,ν,M /|h0 |).
Now for all λ ∈ (0, 3/(2cρ,M kKk1 kKk∞ (|h0 |)−1 )), we have from Bennett inequality


cK,Q,ν,M (|h0 |)−1 λ2

E [exp (λ (g (Xu0 ) + g (Xu1 ) − 2Qg (Xu ))) |Fn−1 ] ≤ exp  
2 1−
and using the Markov property, this leads us to

2cρ,M kKk1 kKk∞ (|h0 |)−1 λ
3




12

S. VALÈRE BITSEKI PENDA AND ANGELINA ROCHE




0 −1 2
c
(|h
|)
λ
|G
|
K,Q,ν,M
n−1

g(Xu )
≤ exp  
E exp λ
2cρ,M kKk1 kKk∞ (|h0 |)−1 λ
2 1−
u∈Tn
3





X
X
I (2) (Xu )
g(Xu ) × exp λ
× E exp λ
!#

"

X

u∈Gn−2

u∈Tn−3

Y

×


i

E exp J (2) (Xu , Xu0 , Xu1 ) Fn−2 
h



u∈Gn−2

where
I (2) (Xu ) = (g + 2Qg + 22 Q2 g)(Xu )
J (2) (Xu , Xu0 , Xu1 ) = (g + 2Qg)(Xu0 ) + (g + 2Qg)(Xu1 ) − 2(Qg + 22 Q2 g)(Xu ).
Now for the second step, we will do the same thing with J (2) (Xu , Xu0 , Xu1 ) instead of g(Xu0 ) +
g(Xu1 ) − 2Qg(Xu ). For all u ∈ Gn−2 we have
|J (2) (Xu , Xu0 , Xu1 )| ≤ 2cρ,M kKk1 kKk∞ (|h0 |)−1 .
We also have
E[(J (2) (Xu , Xu0 , Xu1 ))2 |Fn−2 ] ≤ 4Q((g + 2Qg)2 )(Xu ).
Now we will control the right hand side of the previous inequality. On the one hand, we have
for all y ∈ Rd ,
Z Z
Qg(y) =
Kh (t) Kh0 (x − z − t) Q(y, z)dtdz − Eν [Kh ∗ Kh0 (x − X∅ ].
Using the change of variables
t = hu and z = x − hu − h0 v,
where for x = (x1 , ..., xd )t , y = (y1 , ..., yd )t ∈ Rd , xy denotes the vector (x1 y1 , ..., xd yd )t , we
obtain
Z Z
Kh (t) Kh0 (x − z − t) Q(y, z)dtdz| ≤ kQk∞,∞ .
In the same way we prove that
|Eν [Kh ∗ Kh0 (x − X∅ ]| ≤ kνk∞ .
This leads us to
|Qg(y)| ≤ kQk∞,∞ + kνk∞
and therefore, for all m ≥ 1
|Qm g(y)| ≤ kQk∞,∞ + kνk∞ .

BANDWIDTH SELECTION FOR KERNEL DENSITY ESTIMATION IN A BMCM

13

On the other hand, uniform geometric ergodicity assumption implies that for all m ≥ 1 and for
all y ∈ Rd ,
Qm g(y) ≤ M ρm kKh ∗ Kh0 k ≤ M kKk1 kKk∞ ρm (|h0 |)−1 ,
and therefore, for all m ≥ 1 and for all y ∈ Rd , we have
Qm g(y) ≤ max{kQk∞,∞ + kνk∞ , M kKk1 kKk∞ } × inf{1, ρm (|h0 |)−1 }.
For all y ∈ Rd , we also have
Qg 2 (y) ≤ 2Q(Kh ∗ Kh0 ) + 2(E[Kh ∗ Kh0 (x − X∅ )])2
≤ 2kKk21 kKk22 (kQk∞,∞ + kνk∞ ),
in such a way that for all u ∈ Gn−2 ,
4Q((g + 2Qg)2 )(Xu ) ≤ 8Qg 2 (Xu ) + 8Q((2Qg)2 )(Xu )
≤ cK,Q,ν,M ((|h0 |)−1 + 2 inf{1, ρ(|h0 |)−1 }2 )
and thus, we obtain
E[J (2) (Xu , Xu0 , Xu1 )|Fn−2 ] ≤ cK,Q,ν,M ((|h0 |)−1 + 2 inf{1, ρ(|h0 |)−1 }2 ).
Once again, for all λ ∈ (0, 3/(2cρ,M kKk1 kKk∞ (|h0 |)−1 )) and for all u ∈ Gn−2 , we have from
Bennett’s inequality


0 −1
0 −1 2 2
c
((|h
|)
+
2
inf{1,
ρ(|h
|)
}
)λ
K,Q,ν,M
.


E[J (2) (Xu , Xu0 , Xu1 )|Fn−2 ] ≤ exp 
2c
kKk1 kKk∞ (|h0 |)−1 λ
2 1 − ρ,M
3
It then follows that

0 −1 2
c
(|h
|)
λ
|G
|
K,Q,ν,M
n−1

E exp λ
g(Xu )
≤ exp  
2cρ,M kKk1 kKk∞ (|h0 |)−1 λ
2
1
−
u∈Tn
3


0 −1
0 −1 2 2
cK,Q,ν,M ((|h |) + 2 inf{1, ρ(|h |) } )λ |Gn−2 | 


× exp 
2c
kKk1 kKk∞ (|h0 |)−1 λ
2 1 − ρ,M
3





X
X
(3)
× E exp λ
g(Xu ) × exp λ
I (Xu )
"

!#



X

u∈Tn−4

u∈Gn−3

×

Y



i
E exp J (3) (Xu , Xu0 , Xu1 ) Fn−2 
h



u∈Gn−3

where
I (3) (Xu ) = (g + 2Qg + 22 Q2 g + 23 Q3 g)(Xu );
J (3) (Xu , Xu0 , Xu1 ) = (g + 2Qg + 22 Q2 g)(Xu0 ) + (g + 2Qg + 22 Q2 g)(Xu1 ) − 2(Qg + 2Q2 g + 22 Q3 g)(Xu ).

14

S. VALÈRE BITSEKI PENDA AND ANGELINA ROCHE

Now, iterating this method, we are led to
"

!#
X

E exp λ

g(Xu )

u∈Tn




2 !
2
|Gn−m| λ 
2 ((|h |) ρ ∧ 1)
(|h |) +

l=1




2cρ,M kKk1 kKk∞ (|h0 |)−1 λ

2 1−
3

n
P

0

cK,Q,ν,M

 m=1

≤ exp 


m−1
P

−1

l

× exp cK,Q,ν,M

0

−1 l

0

(|h |)

−1

+

n
X

!!
2

m

0

m

−1

1 ∧ ρ (|h |)



.

m=1

Set m∗ = blog |h0 |/ log ρc. Then we have
n
X


(|h0 |)−1 +

m=1
n
X

=

m−1
X

!2 
2l ((|h0 |)−1 ρl ∧ 1)

 |Gn−m|

l=1


(|h0 |)−1 +

m=m∗ +1

∗

m
X
l=1

2l +

m−1
X

!2 
2l h−1 ρl

l=m∗

∗

 |Gn−m| +

m
X


(|h0 |)−1 +

m=1


≤

m−1
X

!2 
2l

 |Gn−m|

l=1


1
2
6 + (1 +
) (|h0 |)−1 |Tn |.
1 − 2ρ

We also have
(|h0 |)−1 +

n
X


2m 1 ∧ ρm (|h0 |)−1 ≤ c0ρ (|h0 |)−1 .

m=1

In view of the above, for all λ ∈ (0, 3/(2cρ,M kKk1 kKk∞ (|h0 |)−1 )) we have

P

!
1 X
(Kh ∗ Kh0 (x − Xu ) − Eν [Kh ∗ Kh0 (x − Xu )]) > δ
|Tn |
u∈Tn



cK,Q,ν,M c0ρ (|h0 |)−1 |Tn |λ2
  × exp λcK,Q,ν,M c0ρ (|h0 |)−1 .
≤ exp −λδ|Tn | + 
2cρ,M kKk1 kKk∞ (|h0 |)−1 λ
2 1−
3

Taking λ = (δ(|h0 |))/(cK,Q,ν,M c0ρ + (4cρ,M kKk1 kKk∞ δ)/3), we obtain

P

!
1 X
(Kh ∗ Kh0 (x − Xu ) − Eν [Kh ∗ Kh0 (x − Xu )]) > δ
|Tn |
u∈Tn


!
2
0
δcK,Q,ν,M c0ρ
δ
|T
|(|h
|)
n
 .
≤ exp 4c kKk kKk δ
exp − 
1
∞
ρ,M
0
0 + 4cρ,M kKk1 kKk∞ δ
+
c
c
2
c
c
K,Q,ν,M
ρ
K,Q,ν,M ρ
3
3

BANDWIDTH SELECTION FOR KERNEL DENSITY ESTIMATION IN A BMCM

15

The result follows since we can do the same thing for −g instead of g. Now, the proof of (8)
follows the same lines and this ends the proof.


Proof of Theorem 7. We start from the following decomposition, true for all h ∈ Hn
(b
νbh (x) − ν(x))2 ≤ 3(b
νbh (x) − Kh ∗ νbbh (x))2 + 3(Kh ∗ νbbh (x) − νbh (x))2 + 3(b
νh (x) − ν(x))2 .
Hence, since Kh ∗ νbbh (x) = Kbh ∗ νbh (x), by definition of A(x, h) and then by definition of b
h and the
fact that a ≤ b
(b
νbh (x) − ν(x))2

≤ 3((b
νbh (x) − Kh ∗ νbbh (x))2 − aV (x, b
h) + aV (x, b
h))
+3((Kbh ∗ νbh (x) − νbh (x))2 − aV (x, h) + aV (x, h)) + 3(b
νh (x) − ν(x))2
≤ 3(A(x, h) + bV (x, b
h))
b
+3(A(x, h) + bV (x, h)) + 3(b
νh (x) − ν(x))2
≤ 6(A(x, h) + bV (x, h)) + 3(b
νh (x) − ν(x))2 .

(9)

Now it remains to upper-bound E[A(x, h)]. We have
(b
νh0 (x) − Kh ∗ νbh0 (x))2

≤ 3(b
νh0 (x) − Kh0 ∗ ν(x))2 + 3(Kh ∗ νbh0 (x) − Kh ∗ Kh0 ∗ ν(x))2
+3(Kh0 ∗ ν(x) − Kh ∗ Kh0 ∗ ν(x))2 .

With a rough upper-bound of the maxh∈Hn by the
"
E[A(x, h)]

≤ 3E


max

h0 ∈Hn

(b
νh0 (x) − Kh0

P

h∈Hn

we get

V (x, h0 )
∗ ν(x)) − a
6
2

 #
+

"


 #
0
V
(x,
h
)
+ 3 E max
Kh ∗ νbh0 (x) − Kh ∗ Kh0 ∗ ν(x))2 − a
h0 ∈Hn
6
+
+ 3 max
(Kh0 ∗ ν(x) − Kh ∗ Kh0 ∗ ν(x))2
h0 ∈Hn

 
!2
0
X
X
1
V
(x,
h
)
 
≤ 3
E 
Kh0 (x − Xu ) − Eν [Kh0 (x − Xu )] − a
|Tn |
6
u∈Tn
h0 ∈Hn
+

 
!2
X
1 X
V (x, h0 )  
+3
E 
Kh ∗ Kh0 (x − Xu ) − Eν [Kh ∗ Kh0 (x − Xu )] − a
|Tn |
6
0
h ∈Hn

u∈Tn

+
2

+3 max
(Kh0 ∗ ν(x) − Kh ∗ Kh0 ∗ ν(x))
0
h ∈Hn

≤

T1 + T2 + Bh (x).

16

S. VALÈRE BITSEKI PENDA AND ANGELINA ROCHE

We first give an upper-bound for T1 . Let h0 ∈ Hn fixed, now remark that, by Lemma 12,



!2
Z +∞
0
X
1
V
(x,
h
)
 ≥ t dt
P 
T1 =
Kh0 (x − Xu ) − Eν [Kh0 (x − Xu )] − a
|Tn |
6
0
u∈Tn
+
!
r
Z +∞
1 X
V (x, h0 )
≤
P
Kh0 (x − Xu ) − Eν [Kh0 (x − Xu )] ≥ t + a
dt
|Tn |
6
0
u∈Tn
!
Z +∞
√
1 X
Kh0 (x − Xu ) − Eν [Kh0 (x − Xu )] ≥ u du
≤
P
|Tn |
aV (x,h0 )/6
u∈Tn
!
√


Z +∞
ucK c0ρ
u|Tn ||h|
√
≤
exp 4c kKk √u
exp −
du
ρ
∞
2(cK c0ρ + (4/3)cρ kKk∞ u)
aV (x,h0 )/6
+ cK c0ρ
3
≤ I1 + I2 ,
where
I1

I2

√

aC(P,µ)/6

!

ucK c0ρ


exp −


u|Tn ||h|
√
du
=
exp 4c kKk √u
ρ
∞
2(cK c0ρ + (4/3)cρ kKk∞ u)
aV (x,h0 )/6
+ cK c0ρ
3
!
√


Z +∞
ucK c0ρ
u|Tn ||h|
√
=
du,
exp 4c kKk √u
exp −
ρ
∞
2(cK c0ρ + (4/3)cρ kKk∞ u)
aC(P,µ)/6
+ cK c0
Z

ρ

3

where we recall that, since for all h0 ∈ Hn , |h0 | ≥
We first upper-bound I1 ,

log(|Tn |
|Tn | ,

we have V (x, h0 ) ≤ C(P, µ).

!
u|Tn ||h0 |
p
√
√
I1 ≤ C 1
exp −
du
2(cK c0ρ + (2 2/3 3)cρ kKk∞ aC(P, µ))
aV (x,h0 )/6
√


√ ∗
√ ∗
C10
0
0
0
− ac∗
1,
exp
−
ac
|T
||h
|
≤
C
log(|T
|)
exp
−
ac
log(|T
|)
=
C
|T
|
(10)
≤
n
n
n
n
1
1
1
1
|Tn ||h0 |
p
p
√
√
with C1 = exp( aC(P, µ)/6), C10 = C1 2(cK c0ρ +(2 2/3 3)cρ kKk∞ aC(P, µ)), c∗1 = C(P, µ)/(12(cK c0ρ +
p
√
√
(2 2/3 3)cρ kKk∞ C(P, µ))), using the fact that a ≥ 1 (remark that
c∗ does not depend on a).
√ 1
√
4c kKk
u
We turn now to I2 , remark that the function u 7→ ucK c0ρ /( ρ 3 ∞ +cK c0ρ ) is non decreasing
0
and converges to 3cK cρ /(4cρ kKk∞ ) when u → ∞, hence it is bounded by this quantity. We have
then, using again a ≥ 1,


√
Z +∞


u|Tn ||h0 |
 du


I2 ≤ C2
exp 
−
√


6cK c0ρ
aC(P,µ)/6
2 √
+ (4/3)cρ kKk∞
C(P,µ)


√
C20
C200
≤
+
exp(− ac∗2 |Tn ||h0 |)
0
0
2
|Tn ||h | (|Tn ||h |)
Z

(11)
with

aC(P,µ)/6

≤ C20 log(|Tn |)|Tn |−

√

ac∗
2

+ C200 log2 (|Tn |)|Tn |−

√

ac∗
2

,

BANDWIDTH SELECTION FOR KERNEL DENSITY ESTIMATION IN A BMCM

17

√
16 p
C(P, µ)cρ kKk∞ ),
C20 = C2 a(4cK c0ρ + √
3 6
p
C(P, µ)
∗
.
and c2 = √ √
p
2 6
6cK c0ρ / C(P, µ) + 43 cρ kKk∞

C2 = exp(3cK c0ρ /(4cρ kKk∞ ),
C200


= 8C2

√

p

6cK c0ρ /

4
C(P, µ) + cρ kKk∞
3

2

Hence, gathering (10) and (11), there exists C 0 > 0 depending only on C(P, µ), K, cK and ρ,
such that
T1

≤ C card(Hn ) log2 (|Tn |)|Tn |−

√

ac∗

≤ C 0 |Tn |−1

with C = 3 max{C10 , C20 , C200 }, c∗ = min{c∗1 , c∗2 } as soon as a > 4/(c∗ )2 .



6. Acknowledgement
The authors want to thank Claire Lacour for her helpful advices on constant calibration.
References
[1] S. Arlot, and P. Massart. Data-driven calibration of penalties for least-squares regression. Journal of Machine
Learning Research, 10 (2009), 245–279.
[2] I. V. Basawa and J. Zhou. Non-Gaussian bifurcating models and quasi-likelihood estimation. Journal of Applied
Probability, 41 (2004), 55-64.
[3] M. Bec and C. Lacour. Adaptive pointwise estimation for pure jump Lévy processes. Statistical Inference and
Stochastic Processes, 18 (2015), 229–256.
[4] I. Benjamini and Y. Peres. Markov chains indexed by trees. The Annals of Probability, 22 (1994), 219–243.
[5] K. Bertin, C. Lacour and V. Rivoirard. Adaptive pointwise estimation of conditional density function. Annales
de l’Institut Henri Poincaré Probabilités et Statistiques, 52 (2016), 939–980.
[6] L. Birgé and P. Massart. Minimum contrast estimators on sieves: exponential bounds and rates of convergence.
Bernoulli, 4 (1998), 329–375.
[7] S. V. Bitseki Penda, H. Djellout and A. Guillin. Deviation inequalities, moderate deviations and some limit
theorems for bifurcating Markov chains with application. The Annals of Applied Probability, 24 (2014), 235–291.
[8] S. V. Bitseki Penda, M. Hoffmann and A. Olivier. Adaptive estimation for bifurcating Markov chains. To appear
in Bernoulli (2016).
[9] S. V. Bitseki Penda., and A. Olivier. Autoregressive Functions Estimation in Nonlinear Bifurcating Autoregressive Models. Statistical Inference for Stochastic Processes, 20 (2017), 179–210.
[10] G. Chagny. Adaptive warped kernel estimators. Scandinavian Journal of Statistics, 42 (2015), 336–360.
[11] G. Chagny and A. Roche. Adaptive and minimax estimation of the cumulative distribution function given a
functional covariate. Electronic Journal of Statistics, 8 (2014), 2352–2404.
[12] G. Chagny and A. Roche. Adaptive estimation in the functional nonparametric regression model. Journal of
Multivariate Analysis, 146 (2016), 105–118.
[13] M. Chichignoud and S. Loustau. Bandwidth selection in kernel empirical risk minimization via the gradient.
The Annals of Statistics, 43 (2015), 1617–1646.
[14] F. Comte, and C. Lacour. Anisotropic adaptive kernel deconvolution. Annales de l’Institut Henri Poincaré
Probabilités et Statistiques, 49 (2013), 569–609.
[15] R. Cowan and R. G. Staudte. The bifurcating autoregressive model in cell lineage studies. Biometrics, 42 (1986),
769–783.
[16] C. Dion and V. Genon-Catalot. Bidimensional random effect estimation in mixed stochastic differential model.
Statistical Inference for Stochastic Processes, 19 (2016), 131–158.
[17] M. Doumic, M. Hoffmann, N. Krell and L. Robert. Statistical estimation of a growth-fragmentation model
observed on a genealogical tree. Bernoulli, 21 (2015), 1760–1799.
[18] A. Goldenshluger., and O. Lepski. Bandwidth selection in kernel density estimation: oracle inequalities and
adaptive minimax optimality. The Annals of Statistics, 3 (2011), 1608–1632.
[19] A. Goldenshluger and O. Lepski. On adaptive minimax density estimation on Rd . Probability Theory and
Related Fields, 159 (2014), 479–543.

18

S. VALÈRE BITSEKI PENDA AND ANGELINA ROCHE

[20] A. Guilloux, S. Lemler and M.-L. Taupin. Adaptive kernel estimation of the baseline function in the Cox model
with high-dimensional covariates, Journal of Multivariate Analysis, 148 (2016), 141–159.
[21] J. Guyon. Limit theorems for bifurcating Markov chains. Application to the detection of cellular aging. The
Annals of Applied Probability, 17 (2007), 1538–1569.
[22] C. Lacour, and P. Massart. Minimal penalty for the Goldenshluger-Lepski method. Stochastic Processes and
their Applications, 126 (2016), 3774–3789.
[23] O. Lepski. Adaptive estimation over anisotropic functional classes via oracle approach. The Annals of Statistics,
43 (2015), 1178–1242.
[24] T. Patschkowski and A. Rohde. Adaptation to lowest density regions with application to support recovery. The
Annals of Statistics, 44 (2016), 255–287.
[25] M. K. Pitt, C. Chatfield and S. G. Walker Constructing First Order Stationary Autoregressive Models via
Latent Processes. Scandinavian Journal of Statistics, 29, 4 (2002), 657–663
[26] G. Rebelles. Lp adaptive estimation of an anisotropic density under independence hypothesis. Electronic Journal
of Statistics, 9 (2015), 106–134.
[27] G. Rebelles Pointwise adaptive estimation of a multivariate density under independence hypothesis. Bernoulli,
20 (2015), 1984–2023.
[28] A. Tsybakov. Introduction to Nonparametric Estimation. Springer series in statistics, Springer-Verlag, NewYork, 2009.
S. Valère Bitseki Penda, IMB, CNRS-UMR 5584, Université Bourgogne Franche-Comté, 9 avenue
Alain Savary, 21078 Dijon Cedex, France.
E-mail address: simeon-valere.bitseki-penda@u-bourgogne.fr
Angelina Roche, CEREMADE, CNRS-UMR 7534, Université Paris-Dauphine, Place du maréchal de
Lattre de Tassigny, 75775 Paris Cedex 16, France.
E-mail address: roche@ceremade.dauphine.Fr

