Non parametric mixture priors based on an exponential random scheme

1 downloads 0 Views 291KB Size Report
(1999) for a recent ... 1999), defined for exchangeable data in [0, 1]. ...... The problem consists in estimating the batting averages ( 1,...,26) of all 26 major league baseball teams in USA, in the 1986 season, starting from the first 300 at bats.
                

         Ý

   



                                

  

                                                        

                                      

                      

            

                          

       

                   !    " "      #    

         

     

               

                                                                   

                         !"""# $ %  &     '( ))"*+,,'*"#  - .//  0   Ý    $        $ 1     2/ $ &   !$ !**""$ &$    2 3  4         &$    - ./0 / $  $      5    $    6 !7$ !"*87  $   9 20  4  0//  £



                         !            "         

                                           #      $  % 

  &'''(    ) *  )   !                +               +

   ,              &'-.(  / &'-0(  )                                     

    1   1      "                          1   1            1

1 2   +                              "        #3            3                    +                      

 +                                *    

                 &  '''( !         45 6  

     

!              7            45 6     +                        

 +                        8   &955(          )  +        7       !        ! 1

1                          

 &  .( 7    

          :  &955(       7         

    &  9( 7    ;         ,)       )                    

      )     *  )                                       ) &!  !(   ) +     0                 

                 

     

 &'?   :( !  )      )         &    ( )    7    

 /  @           

         )        A  9        &( A   &(  5      7  

   

       &( 

     &(      

            &$   ( A  & &(( 7 

       ;    

B                      A  9             &(    &(  5        

  &$   ( A  &(



                               &(    &(  5  

7 "    )        A  9     

  

    &*    %  ''0   9-.( C     

                 & (       *           "    

  &(   )           )   & ( A   ! )           A     !  "      &( A 

 &(  5         

                                 )                      

                 7          :  &955(        )                       /  &( A &  &(( &(       !       &( A  &( &(

   D A  @  &(  E 7                 &(  

   ,    &,(   "

               /       .

 

+             )     & ( A   &(       ! ) )                 ;             &(

 & ( A   &&((        ) )     &(  

           ) 

    & F    )  ( 

         &  (         +                  

             "   )               &(        &( A  7    

      A    

            )  &( A   &(  5      /             ;          )     &(  !  5 

  & ( &$   ( A    &(&(  A    &9(     7 

      ;      )     :  &955(

                      & $   ( !  "             #  & $   (                 !     & (               



&$   ( A   &$ ( &(   &  (   & ( &( &(           & $  (         $                   & (         $                        %      &(        

& $   (           &  (   &(       & ( A 5 & $   (             2  &(  &(         &$   (      & $ (           &)(              &$   (       &&G((  G     & (   G  A  9               & $   (   

         7     )        

   &$  ( A



     )  



0

    &*      ( &         '       (          & $   (                    '                  & $   (     )     & $ (        (      

 & $ (   )     & (   J & $  (  # 5 &8 

    " & $  (   J &   (( 7 

        )  



&$   ( A J &$  ( &( *         )    = ) . &    (         -                  $ & $ $   (   

+  7 +   & $ (  $ & $   (     

          



&$   ( A $ &$   ( &(           $               J   J   2  &'''(

          

              

       

           F         +       *     )    

       ) +  /            +           "    

     7        

             



   *         .   

          

    

   8                  

             /        "       )    )       7    B          

         



% &          ( A   &(&& ( 

&.(

 &       

  

         8          

                  *               A  9       &     

  

     

 A  & $   ( @  A  9    $    =

    

 

        )                  A  9    C       &       &  (    " 

      

  7    

 &      & $   (   

  

            

                #        ' & (     &' (  '       

                        &  '/       & ' (   

  

    & & ( A 5( A             & $   (           & $   (       

               A  9     &.(  ")    

     &(    (            

    



&$   ( A  &$ ( &(   &$ (   +       )     &( &  (  ) <    

  %  

   

                   +   & $ (            

    7               +  7    ;          &'''(      45 6 %            

     ) %%             "                    

 )         )          )     1

 1    

)  &(                 <     &$ ($ &(              $ &(     ) <    

    

             

    

    

#            

      7     

                 

         A A A  A A     +          ! ) )                  7        

      +       

 &     ( 7  

 

      = 7    ?



          #      *  + '?0$   J  2   ''' I         &      (  

#         %% I  )      

    

  #         %%         J             )              +           

    )             !  C 

      +        J         

              *  

       

!      I B &''0(      #   

  !     

       0101      

 #          = * 3       

           =9

     8 )      %%      

           

   7       

+    ;     #      !  &I  ''0(   &' (               &( A  Æ  &(  &0(



 Æ& (           * &     (         A   A &  (&  ( & (  ! 9                  & '(       )       

               &( A  Æ  &(







&>(

   &0(  ! $             $          7 

 ''-      K   9555 8    

           K      

        

   * J           &       (       )   )           & ( &           ( !    !          &+    + (  

!   )                A   & ( +    +                     &(              !    !  -     

7 

        & (      

, &      !    ! (

 & ( & $  (    







&=(

 ,                 )  !    !           #           !     !  )  ; #         &            ( 

A  E )   A ' E )  E E )  )      !     !   "  ! ! A     $  I          &=(           (        #                  )       K   &9555(   

    3                        )          

 I      &5 (  !      

    J         .    .           

     )    

   

#     ' A 9   A " &95 ( 7             55              .  8       .    )     &  & (   . A 555 8     0101       )     )  ! $ A 55  9   

#        

    8         >  =              .                )   )              I)                 

  )          F        * F    )  

)      "      

                 

     

               I!

            3           3      &     (    "   I+ +          +    +           

         )          +     )

      3                             3 +     &J  '-=(                +   5

  3              &      (            )   +      )      %    :  &''>(  )

 )   "    Æ           " ) *  "                 /      &      (                  

#               

#    7 )              

 7               = ,   #  8  &''>(            3 &   &   >??>-? .



 8 &'?( &              :   N  8    H +   7 I &'-.(               2      &  2;)  N 2  ( *     H + 9-?.59 J  , &'-=( %    I+  ,   3 &   &   0.?00> J   I J  NP  2  2: &'''(       #        &    0.>- J  2   &'''(  

        #      M)    @QQ   +Q Q      K    &9555(  + )     %      #              4 (  .?.'5 /  *H &'-0( C 

           #   &    .>.>?     7 

 / &''-( *            #    1 3    9-.9'?   I &'''( 2         3    .?..'.   I  :   &955( 

         

            O  ) M)R /    

    I  8   / &955( %       & )    3 2     I(  +  2  /I &'-.( 5    )   *    C

  2  P &''5( #      !    !     )    

 3 &   &   =?=90 I  N &''0( *  ) !   #        =.'=>5 7  * &''>( 5       4     M)S )

 2  C S+ J  

0

(b)

200

100

400

600

sampled k

300 200

sampled k

400

800

500

1000

(a)

0

200

400

600

800

0

iterations

200

400

600

800

iterations

 @ %%            . A >55   . A 555

>

(b)

0.0

0.0

0.05

0.05

0.10

0.15

density

0.10

density

0.20

0.15

0.25

0.30

0.20

(a)

10

15

20

25

30

35

10

Rel. vel.(km/sec x 10^-3)

15

20

25

30

35

Rel. vel.(km/sec x 10^-3)

 9@ #     

) )       + Q    -9 

 &2 ''5( & ( %%       )   & ( 7                

=

0.20 0.15 0.10 0.0

0.05

probability

6

7

8

9

10

11

12

13

14

82

number of clusters

 .@ %%                    )     

    7 )          

?

15 0

5

10

density

20

25

30

(a)

0.0

0.2

0.4

0.6

0.8

1.0

theta

0.26 0.24

o

oo

oo

o oo

o

o

o

oo

o ooo ooo

o

0.22

o

0.18

0.20

our estimates

0.28

0.30

(b)

0.18

0.20

0.22

0.24

0.26

0.28

0.30

MLE

 0@ & (,         

   & ( C                /, 7  ;     )

 

-

Suggest Documents