Forest modelling: the gamma shape mixture model and simulation of tree diameter distributions

Table 2 The p values for the Anderson-Darling k-sample test comparing DBH distributions within drawn DBH samples

Plot	Samples of 100 DBHs		Samples of 250 DBHs		Samples of 500 DBHs
Plot	MDGR procedure^a	MH method^b	MDGR procedure	MH method	MDGR procedure	MH method
RS01	0.9253	0.0297	0.8078	0.1976	0.3301	0.0009
RS02	0.0582	0.0038	0.9798	0.0037	0.6010	0.0291
RS03	0.6648	0.0001	0.7463	0.0001	0.0924	0.0001
RS04	0.8603	0.0305	0.1516	0.0577	0.8856	0.0047
RS05	0.4721	0.0031	0.8919	0.0001	0.0263	0.0006
RS06	0.4708	0.5851	0.8513	0.4444	0.0935	0.2412
RS07	0.4452	0.1675	0.6376	0.0951	0.5768	0.8390
RS08	0.3950	0.0810	0.2444	0.2489	0.7798	0.8061
RS09	0.1200	0.0070	0.0719	0.0118	0.0330	0.2469
RS10	0.4200	0.2767	0.6434	0.0179	0.1769	0.8099
BMS01	0.6946	0.8743	0.9177	0.1956	0.1090	0.4082
BMS02	0.5359	0.0172	0.7526	0.1888	0.7488	0.5131
BMS03	0.5717	0.6526	0.4623	0.2471	0.5513	0.5977
BMS04	0.5778	0.8791	0.3161	0.5320	0.7599	0.2769
BMS05	0.2923	0.6314	0.9241	0.7821	0.8939	0.0861
UID01	0.2080	0.7176	0.1280	0.5187	0.2531	0.0930
UID02	0.1098	0.7715	0.5172	0.0355	0.7183	0.1302
UID03	0.2167	0.0469	0.8524	0.4696	0.2260	0.8897
UID04	0.3961	0.8530	0.2067	0.3511	0.3385	0.4795
UID05	0.0109	0.1558	0.9657	0.1492	0.0172	0.5158

^aGamma random numbers were generated in multinomial distribution cells using the acceptance-rejection principle with proper choice of the majorisation function (when the shape parameter was less than 1) or as the sum of two independent gamma variates (when the shape parameter was greater than or equal to 1) (Ahrens and Dieter methods; for detailed information, see Ahrens and Dieter 1974, 1982)
^bThe standard Metropolis–Hastings algorithm with jumping normal distribution was used (Robert and Casella 2004)

ISSN: 1297-966X