Page 306 - Applied Probability
P. 306

13. Sequence Analysis
                              294
                                    (b) Show that the likelihood ratio test rejects the null hypothesis
                                        λ = λ 0 whenever the statistic
                                                                       n

                                                           T
                                                               =2λ 0
                                                                      i=1  X i
                                        is too large or too small.
                                    (c) Prove that T has a χ 2  distribution.
                                                          4n
                                 2. Use the Borel-Cantelli lemma [4] to prove that the pattern SFS of a
                                   success, failure, and success occurs infinitely many times in a sequence
                                   of Bernoulli trials. This result obviously generalizes to more complex
                                   patterns.
                                 3. Renewal theory deals with repeated visits to a special state in a sto-
                                   chastic process [3, 4]. Once the state is entered, the process leaves it
                                   and eventually returns for the first time after n> 0 steps with proba-
                                   bility f n . The return times following different visits are independent.
                                   Define u n to be the probability that the process is in the special state
                                   at epoch n given that it starts in the state at epoch 0. Show that
                                                u n  = f 1u n−1 + f 2 u n−2 + ··· + f n u 0
                                   for n ≥ 1. If we define the generating functions U(s)=    ∞  u n s n
                                                                                       n=0
                                                       n
                                   and F(s)=     ∞  f n s with f 0 = 0 and u 0 = 1, then prove that
                                                 n=0
                                   U(s)= [1 − F(s)] −1 .
                                 4. Repeated visits to a pattern such as GCGC in a DNA sequence consti-
                                   tute a renewal process as noted in the text. Given the assumptions of
                                   Section 13.2, one can calculate the generating functions U(s) and F(s)
                                   defined in Problem 3 for renewals of the pattern R =(r 1 ,... ,r m ). If
                                                        and
                                   we let p R = p r 1  ··· p r m
                                                1                             k =0
                                              &
                                          =                  1                1 ≤ k ≤ m − 1 ,
                                       q k                      (m−k)
                                                p r m−k+1  ··· p r m {R
                                                                    =R (m−k) }
                                                0                             k ≥ m
                                   then show that
                                                                n−1

                                                             =                           (13.11)
                                                         p R        u n−k q k
                                                                k=0
                                   for n ≥ m. Use this to prove that
                                                   p R s m
                                                          =[U(s) − 1]Q(s)                (13.12)
                                                   1 − s
                                                                 m
                                                              p R s +(1 − s)Q(s)
                                                    U(s)  =
                                                                  (1 − s)Q(s)
   301   302   303   304   305   306   307   308   309   310   311