Ecological Archives A016-005-A2

Charles D. Canham and María Uriarte. 2006. Analysis of neighborhood dynamics of forest ecosystems using likelihood methods and modeling. Ecological Applications 16:62–73.

Appendix B. Results of a reanalysis of patterns of seedling dispersion for five tree species originally reported by LePage et al. (2000).

We reanalyzed patterns of seedling dispersion for five tree species originally reported by LePage et al. (2000). The data are from samples in mapped stands in "partial" canopy sites. See LePage et al. (2000) for field methods and descriptions of the study sites and species. We used likelihood methods to fit two different dispersal functions: a lognormal function:

 
and an exponential function:
 

where dbhi and disti are the size and distance to i = 1..n trees within a fixed radius of a seedling quadrat (20 m radius for hemlock and cedar, 30 m for the other three species).  is a normalizer that equals the arcwise integration of the dispersal kernel. We also estimated a set of eight substrate favorability terms as done in the original analysis (see LePage et al. 2000 for details). Maximum-likelihood parameters were estimated using simulated annealing, an algorithm for global optimization. We compared the models using AIC (Table B1):  the best model (lowest AIC) for each species is indicated in bold type. In order to compare our results to the original analyses reported by LePage et al. (2000), we also did a set of analyses with the exponential function, but in which the values of the and   parameters were fixed at 2 and 3, respectively.

TABLE B1.  Parameter estimates, log-likelihood, and AIC for three alternate models for seedling dispersion of five tree species of interior cedar hemlock forests of British Columbia. Data are from LePage et al. (2000).

           
 

Hemlock

Cedar

Spruce

Amabilis Fir

Birch

 

(800 quadrats)

(620 quadrats)

(402 quadrats)

(331 quadrats)

(362 quadrats)

Lognormal Function

       

AIC

19614.5

3319.5

1296.1

1647.4

1371.9

Log-likelihood

-9795.3

-1647.7

-636.0

-811.7

-674.0

r2

0.41

0.46

0.33

0.28

0.32

A

4434.1

1836.1

46154.4

342.0

1292.1

2.06

3.74

1.47

3.57

0.08

B

10.92

13.52

6.59

15.71

5.92

0.179

0.243

1.735

0.229

0.589

Exponential Function

       

AIC

20538.8

3381.1

1295.5

1666.8

1379.8

Log-likelihood

-10257.4

-1678.5

-635.8

-821.4

-677.9

r2

0.35

0.43

0.34

0.25

0.28

A

48806.5

15296.7

35267.2

416.6

1292.7

2.12

3.62

1.18

3.40

0.17

B (*10000)

0.016919

0.011066

103.098247

0.001942

7.723590

3.04

3.31

1.08

4.99

2.61

Exponential Function with fixed a and b

     

AIC

20544.9

3621.9

1294.4

1682.8

1386.4

Log-likelihood

-10262.4

-1801.0

-637.2

-831.4

-683.2

r2

0.34

0.39

0.34

0.26

0.29

A

49621.6

4943.4

39747.3

600.4

2120.0

2

2

2

2

2

B (*10000)

0.022019

0.877780

0.004031

0.929382

2.449535

3

3

3

3

3

 

LITERATURE CITED

LePage, P. T., C. D. Canham, K. D. Coates, and P. Bartemucci. 2000. Seed abundance versus substrate limitation of seedling recruitment in northern temperate forests of British Columbia. Canadian Journal of Forest Research 30:415–427.



[Back to A016-005]