Page 416 - Contributed Paper Session (CPS) - Volume 6
P. 416

CPS2007 Jai-Hua Yen et al.
                  which is formed by singleton species at the same time. Accordingly, we have
                  the expected sampled doubleton richness:
                        ( ) ≈ { [(1 − ) +  × (1 − )] × (− × )}
                            2
                                     2
                                                             1      1
                                      +  { ×  ×  × (1 − ) ×    × (− × )}.
                                             1
                                                             
                                                                   ,
                  where  denotes the number of sampling unit. By solving the two equations
                  above, we have the singleton and doubleton richness adjustment:
                                                      1
                                         1  =  (1−̅ ̂ ×̂)(−̅ ̂ ×̂)′              (4)
                  and
                                                     1
                                       2 − 1 ×̅ ̂ ×̂×(1− )×   1  ×(−̅ ̂ ×̂)
                                                     
                                                         ,
                                 2  =     (1−̅ ̂ ×̂)×(−̅ ̂ ×̂)   .          (5)
                      However, the estimation of traditional Chao2 estimator will be inaccurate
                  even though   and   are asymptoticly unbiased. It causes the value of
                                         2
                                 1
                    2
                   1  overestimated. Hence, we choose first-order Jackknife and Chao2 richness
                  2 2
                  estimator  as  the  theoretical  foundation  of  deriving  the  adjusted  richness
                  estimator.  We  propose  an  adjusted  richness  estimator  by  Taylor  series
                                    2
                  expansion  of  (   1  ) by  the  mean   and  .  Then  we  get  the  difference
                                                       1
                                                              2
                                   2 2
                  between  [( 2 )] 2  and  (   1 2 ) to have the adjust term:
                           (2 2 )    2 2
                                                                             2
                       1 2  [( )] 2  ̂( )  ( )̂( ,  )  [( )] ̂( )
                                              1
                                                                 2
                                                                           1
                                                              1
                                                      1
                                                                                    2
                                 1
                   (   ) ≈          +         −                   +                 ,
                     2 2    (2 )  2( )       [( )] 2       2[( )] 3
                                             2
                                                           2
                                                                               2
                                   2
                                                                 
                  where ̂( ,  ) = −   1  2  , ̂( ) =  (1 − ). Therefore, we have the
                                                    
                               1
                                                          
                                  2
                                                                 ̂
                                          ̂
                  adjusted richness estimator:
                                                           2
                                                           1
                                                                        1
                                                                 1
                                  ̂   =  ,  +  −1  {( 2 2  −  2 2  −  2 2 2  ) , 0}.       (6)
                                                
                                                                        2
                  When 0 ≤ 2 ≤ 1, by simulation studies, the adjusted richness estimator will
                  be:
                                         ̂   =  ,  +  −1  .             (7)
                                                            1
                                                        

                  3.  Result
                  a.  Simulation Results
                      To test the performance of the adjusted richness estimator, we presented
                  the  simulation  results  by  several  species  detection  models  and  different
                  settings of number of sampling units. We fixed    = 40 and  = 100. 500
                  simulation  data  sets  were  generated  and  200  bootstrapping  trials  were
                  conducted  by  each  simulation  data.  The  bootstrapping  method  is
                  regenerating  , ,  1′  and   by  binomial  distribution  independently  in
                                                2
                  order  to  increase  the  estimated  standard  error  while  the  traditional
                  bootstrapping method usually underestimates the standard error in this case.
                  In true method, the estimation of species richness used the traditional Chao2
                                                                     405 | I S I   W S C   2 0 1 9
   411   412   413   414   415   416   417   418   419   420   421