C O N T I N U O U S U N E Q U A L F O R E S T I N V E N T O R Y P R O B A B I L I T Y USING M U L T I S T A G E S A M P L I N G W I T H PARTIAL R E P L A C E M E N T By Shongming Huang B. Sc. F. Nanjing Forestry University, 1985 A THESIS SUBMITTED T H E IN PARTIAL REQUIREMENTS M A S T E R FOR OF FULFILLMENT T H E D E G R E E OF FORESTRY in T H E FACULTY OF GRADUATE DEPARTMENT OF STUDIES FORESTRY We accept this thesis as conforming to the required standard T H E UNIVERSITY OF BRITISH COLUMBIA November 1988 © Shongming Huang, 1988 OF In presenting degree freely at this the available copying of department publication of in partial fulfilment University of British Columbia, for this or thesis reference thesis by this for his thesis and scholarly or for her Department The University of British Columbia Vancouver, Canada DE-6 (2/88) I I further purposes gain the shall requirements agree that agree may representatives. financial permission. Date study. of be It not that the Library by understood be an advanced shall permission for granted is for allowed the make extensive head that without it of copying my my or written Abstract Continuous forest inventory (CFI) with partial replacement of sampling units (SPR) has been established as an efficient inventory technique. Previous C F I with S P R theory was based on equal probability sampling which might not be as efficient as the unequal probability sampling that is now widely used in forest inventories. In this study, the general theory of multistage sampling combined with C F I using unequal probabilities in a SPR structure on two occasions is discussed. Estimators of the current means and the change in means for both one-stage and two-stage cases are given, along with the estimators for the overall variances. The estimators derived for estimating the current means take the same form as those developed by other researchers. However, the estimators proposed for predicting the change in means take simpler forms than those presented in the literature. This might result in greater simplicity when used in data processing. The application of C F I with S P R using multistage unequal probability sampling to a Chinese national forest inventory is presented for illustrative purposes. Some particular points are addressed regarding specific situations in China. Although no precise theoretical demonstration of the gains from using this method is given, the combination of highly efficient sampling with unequal probabilities and the very practical multistage sampling in C F I with SPR inventories can certainly provide an efficient alternative to traditional Chinese inventory systems. 11 Table of Contents Abstract ii List of Tables v Acknowledgements vi 1 Introduction 1 2 Literature R e v i e w 8 3 One-stage S P R w i t h U n e q u a l Probabilities 18 3.1 Definition of the Sample 18 3.2 Estimator of the Current Mean 24 3.3 Estimator of the Change 28 4 5 Two-stage S P R w i t h U n e q u a l Probabilities 34 4.1 Definition of the sample 34 4.2 Estimator of the Current Mean 38 4.3 Estimator of the Change 42 A p p l i c a t i o n to Chinese Forest Inventory 48 5.1 Background on Chinese Forests and Inventory Problems 48 5.2 Design of a Two-stage Sample . 53 5.2.1 General Considerations 53 5.2.2 The Shape and Size of Sample Plots iii : 55 5.2.3 Sample Size 57 5.2.4 The Layout and Survey of Sample Plots 61 6 Conclusions 63 7 Literature Cited 66 A Derivation of Equation 3.6 73 B Derivation of Equation 3.11 76 C Derivation of Equation 4.19 79 iv List of Tables 3.1 Various observations and sample plots 19 4.2 Various observations and sample plots 35 5.3 Forest types and their amounts in China 49 5.4 Forest land utilization in China 49 v Acknowledgements I am most grateful to Dr. Peter L. Marshall, my supervisor, for his guidance, encouragement, enthusiasm, patience, and understanding throughout the period of my gratuate studies at the University of British Columbia. Special thanks to the members of my thesis committee, Dr. A . Kozak, and Dr. J . H. G. Smith for reviewing the thesis draft and making very useful comments. I am very grateful to the State Educational Committee, Government of the People's Republic of China, for financial assistance in form of the National Overseas Graduate Scholarship for two years, and to the Department of Forest Resource Management at the University of British Columbia for providing a teaching assistantship. I thank Nanjing Forestry University and the Ministry of Forestry of the People's Republic of China for giving me the opportunity to study abroad. I extend my sincere appreciation to the Faculty of Forestry, especially the Department of Forest Resource Management at the University of British Columbia for providing such an excellent learning environment, and to my teachers and my fellow biometrics students at the University of British Columbia. I thank Dr. Val Lemay of the University of British Columbia, and Professor D'eyin Zhou of Nanjing Forestry University for reviewing the thesis draft and for their comments, as well as Mr. Basivi Reddy and Mr. James Thrower for assisting during the early stages of writing this thesis. M y appreciation also goes to my fellow Chinese students in Forestry, and to the General Consulate of the People's Republic of China in Vancouver. Finally, I am very grateful to my mother, my grandmother, and my friends in China. vi I shall always feel indebted for their continuing encouragement and moral support. vn To my mother, I dedicate this thesis. Chapter 1 Introduction Originally, inventory was a commercial term meaning the preparation of a detailed descriptive list of articles with number, quantity and value of each item (Loetsch and Haller, 1964). In forest resource management, the term inventory has a wider scope. Forest inventories require specialized and often complicated information and consequently expensive techniques. In recent years, the increase in intensive forest management has stimulated development of integrated or multiresource inventories (Beers, 1978; Lund, 1979; McClure, 1979). Forest inventories have expanded to cover a broader range of topics. A complete forest inventory not only provides information about the estimates of timber quantity and quality, but also provides estimates of growth and drain. When needed, additional information may be collected according to the purpose of the inventory. Many statistical methods have been used in forest inventories to summarize data and estimate resource characteristics. Relying on sampling for collecting important series of data that are published at regular intervals has become common (Cochran, 1977). Forests are dynamic populations and a. survey at infrequent intervals is of limited use. Highly precise information about the characteristics of a population in 1970 and 1980 may not help much in planning that demands a knowledge of the population in 1990. A series of small samples at shorter intervals may be more valuable. A wide variety of sampling techniques and associated statistical procedures is applied to forest inventory designs for both single-occasion and continuous forest inventories. The 1 Chapter 1. 2 Introduction forester is often faced with the question of which sampling technique to use. The answer to this question depends in part on the sampling objectives. Often, the objective is to minimize cost subject to precision limits, or maximize precision subject to cost limits (Fontaine, 1973; Cochran, 1977). Cost is measured by the resources required to conduct the survey. The precision of a sampling procedure is usually measured as the variance of the sample estimator. The best of the several alternative procedures is the one that most nearly meets survey objectives. Stott and Semmens (1962) defined the term "Continuous Forest Inventory" (CFI) as "a precise frequently repeated, and directly comparable measurement of all commercial trees in systematically placed sample plots. These plots have fixed radii and are permanently located in the forest. Their treatment, and the treatment of surrounding forest must be analoguous." CFI's use repeated sampling to estimate quantities and characteristics of the forest present at different occasions and the change in the forest during the intervening period. They provide information on forest populations at various points of time. Single-occasion forest inventories can only provide information on the current values of the forest population at a given time. The limitations of single-occasion inventories in terms of estimating resource changes are obvious (Gillis, 1988). There are three different strategies that CFI's can follow. These are: • The sample plots at later occasions are all different from those at previous occasions. A new sample is drawn at each occasion, and the means and totals, as well as standard errors are calculated separately. Estimates of the change or growth would be the difference in the means for the two inventories. Early CFI's usually used this method. It can provide the current forest information in a fast and relatively cheap way. • The same set of sample plots is used on all occasions. This involves the concept Chapter 1. Introduction 3 of permanent sample plots. Estimates of the means, totals, and standard errors for each inventory would again be found as in the case of two separate inventories. Similarly, the differences between the means for each inventory would indicate the change in the forest. • Partial replacement of sample plots from occasion to occasion. This means that a certain portion of the sample plots observed on the previous occasions will be remeasured at a later occasion as well as a new set of sample plots. This technique has been accepted as a valid CFI technique for estimating current forest resource values and changes in these values over time. Detailed accounts of the development and application of this procedure to forest inventory are given by many researchers and many articles have been published on the general statistical theory of sampling with partial replacement. These will be discussed in the literature review. Partial replacement of sample units requires selection of sample fractions on each occasion (i.e., the proportion of permanent and temporary sample plots). These fractions can be selected either for the purpose of maxmizing precision or minimizing cost. Cochran (1977) made the following statements about replacement policy for the purpose of maximizing precision: 1. For estimating change, it is best to retain the same sample throughout all occasions; 2. For estimating the average over all occasions, it is best to draw a new sample on each occasion; 3. For current estimates, equal precision is obtained either by keeping the sample or by changing it on every occasion. Replacement of part of the sample on each occasion may be better than these alternatives. Chapter 1. Introduction 4 Since there is nearly always a positive correlation between the measurements on the same unit on successive occasions, the first two replacement policies are obvious. The third replacement policy is less obvious and it can be very flexible. Much of the research has been directed towards this aspect and the partial replacement policy is often decided by some kind of cost function because of the difficulty of optimizing the precision. In a typical C F I there are three quantities which we may wish to estimate. These are: 1 . The change in values of a parameter from one occasion to the next (e.g., the difference between the total volumes on different occasions). One of the main purposes of a CFI is to obtain accurate information about the changes in values from occasion to occasion. 2. The average value of a parameter over all occasions (e.g., the average volume for the whole population at a certain inventory period). This value may be adequate for a population in which changes over time are slow, but this will not often be the case in forest populations so the overall average value may not be as useful as the current average. 3. The average value of a parameter for the most recent occasion (e.g., the average volume for the whole population at the time of inventory). This is the most important item in most CFI's. Most inventories concentrate on the current average and an accurate estimate of this value will ensure not only an accurate estimate of the changes in values, but also an accurate average value over all occasions. In national and other large forest inventories covering broad forest areas, each unit in the population can be divided into a number of smaller units, or subunits. For instance, the forest areas in China can be divided into TV different units according to forest region. These units can be subdivided into subunits according to forest types or ages. Chapter 1. 5 Introduction Suppose a sample of n units has been selected. If subunits within a selected unit give similar results, it seems uneconomical to measure them all. A common practice is to select and measure a sample of the subunits in any chosen unit (Cochran, 1977). This technique is called subsampling or multistage sampling since the unit is not measured completely but is itself sampled. The principal advantage of multistage sampling is that it is more flexible than single stage sampling (Hansen et al. 1953; Freese, 1962). It allows less intensive sampling, and can take the advantage of prior knowledge of a population which can be used to increase the precision or usefulness of the sampling. Multistage sampling is also a cost effective technique (Deming, 1950; Titus, 1979). In many forest sampling situations, locating and travelling to a sampling unit is expensive, while measurement of the unit is relatively cheap (this is partically true in Chinese forest inventories). In these circumstances, it seems logical to make measurements on two or three units at or near each location. It can yield estimators of a given precision at a cost lower than a completely random sample. This eventually reduces to a balance between statistical precision and cost. Substantial work has been done using multistage sampling in CFI's with SPR. Singh (1968), Chakrabarty and Rana (1974), and Rana (1978) used two-stage or multistage successive sampling. Omule (1981) provided an extensive review and derived a complete set of formulas using multistage sampling in CFI's with SPR. One practical problem of multistage sampling in an extensive population is that primary units that vary in size are encountered frequently. If the primary unit sizes do not vary greatly, one method of overcoming this problem is to stratify by the size of the primary unit so that the units within a stratum become equal, or nearly equal (Cochran, 1977). The multistage sample may then be combined with any type of sampling of the primary unit. Results are given for stratified sampling of the primary units in a two-stage Chapter 1. Introduction 6 sample by Hansen et al. (1953) and by Cochran (1977). Often, however, substantial differences in size remain within some strata, and sometimes it is advisable to base the stratification on other variables (Cochran, 1977). Concentrated effort is required to obtain a good working knowledge of multistage sampling when the units vary in size, because the technique is flexible. The units may be chosen either with equal probabilities or with probabilities proportional to size (PPS) or to a prediction of size (3P). Various rules can be devised to determine the sampling and subsampling fractions, and various methods of estimation are available (Kish, 1965; Cochran, 1977; Brewer and Hanif, 1983). The problems of multistage sampling with units of unequal size are of common occurrence. The sampling procedures in CFI's are more complicated. Since the same population is sampled repeatedly, the sampler is in an ideal position to make realistic estimates both of costs and of variances and to apply the techniques that lead to optimum sampling efficiency. One important question is how frequently and in what manner the sample should be changed as time progresses. Many considerations affect the decision. In forest inventories, loss of permanent sample plots and other unexpected occurences may make the inventory very complicated. The question of replacement of part of the sample and the related question of making estimates from a series of repeated samples should be carefully considered. Much of the work in CFI's with SPR has been done using different equal probability sampling techniques. Relatively little research has been done involving unequal prob- ability sampling in CFI's. Sampling with unequal probabilities on single occasions has been studied in detail by Grosenbaugh (1964, 1965, 1971, 1976), Schreuder et al. (1968, 1971), Schreuder (1984). Furnival et al. (1987) and many other researchers. Althrough Van Hooser (1972, 1973) and Yandle and White (1977) presented some general theories about sampling with unequal probabilities in two-stage sampling, and Scott (1977, 1979, 7 Chapter 1. Introduction 1984) discussed midcycle updating using PPS and 3P sampling, some questions remain unsolved. This is particularly true for large scale applications, optimum sample designs, sample size calculations and sample allocations. A general theory of multistage sampling with unequal probabilities in CFI's with SPR has not been developed for forestry. Many researchers have shown that sampling with unequal probabilities can be much more efficient than sampling with equal probability, particularly when the sample units have large variations which is often the case in forest inventories. Consideration of cost often dictates the use of multistage sampling. Hence the combination of multistage sampling with unequal probabilities is particularly meaningful in CFI's. The main objectives of this study are: 1. To summarize the use of CFI's with SPR and discuss multistage sampling using unequal probabilities of selecting sampling plots in CFI's with SPR, and combine the general theory of multistage sampling with unequal probabilities and SPR in CFI's. 2. To illustrate the application of multistage sampling with unequal probabilities in CFI's with specific reference to Chinese forest situations. In the following chapters, a literature review is given first, and then, the combination of multistage sampling with unequal probabilities and partial replacement is discussed. The estimators are given for current means and change in means for both one-stage and two-stage SPR with unequal probabilities are given. A n application is discussed with special reference to Chinese forest situations and some suggestions and recommendations are given. Finally, the conclusions are presented. Chapter Literature 2 Review The concept of C F I was originated by Stott (1947). It involves the concept of permanent sample plots (PSP's) which means that the established sample plots are measured more than once. Based on a PSP system, periodic measurements are made to determine primarily tree growth and yield with time or treatment. The original purpose of these CFI plots was to collect volume, growth and mortality information and provide a base for better forest management (Stott and Semmens, 1962; Stott, 1968; Solomon, 1979). In the early stages of CFI development, the plots established at the first occasion were all remeasured at later occasions (Spurr, 1952; Cutter, 1955; Hall, 1959). Through these repetitive measurements, estimates of present volume and change in volume, basal area, ingrowth, and diameter distribution were obtained. Although remeasurement of all plots could provide more precise information, and many researchers argued for its advantages, the general theory of C F I in the early years was not very clear, and the high cost limited its use. With the rapid development of sampling techniques, many more efficient procedures were presented. There is a considerable volume of literature which addressed the statistical theories used in CFI's. Contributions have been made by many researchers. Repeated sampling with partial replacement (SPR) was introduced by Jessen (1942) and its application resulted in a great gain in sampling precision for estimating the current mean. Patterson (1950), Cochran (1953), Yates (1960) and many other researchers continued to work on the general theory of repetitive sampling. A detailed discussion of the statistical 8 Chapter 2. Literature Review 9 development was given by Omule (1981). Most of the papers reviewed here address the development of C F I with SPR in forestry situations. Bickford (1956) first introduced the concept of S P R into CFI's. In his proposed design for a continuous inventory of the Northeastern United States, he presented the idea of a partial replacement sampling scheme for forest inventories and pointed out some of its advantages. Three years later, Bickford (1959) reported an application of C F I with SPR to the Allegheny National Forest in Pennsylvania. He concluded that inventories need repeated sampling and good designs. He further noted that for an area as large as a national forest, there are usually sizable differences in forest cover and stratified sampling is, therefore, a more efficient procedure. Several sampling strategies deserving consideration were presented including SPR. His proposed method of continuous inventory was almost the same as the procedure later fully developed for C F I using SPR theory by many other researchers. The method consisted of remeasurement of enough of the plots established at the first occasion to obtain the needed accuracy in periodic net growth and establishment and measurement of enough new plots to obtain the desired accuracy in current volume. He then estimated the current volume by combining the independent estimates obtained from the different types of sample plots. A test was made in the Allegheny National Forest. The test made use of: aerial photographs, double sampling, point sampling, remeasured plots, new plots, and adjustment through regression^ To be able to use stratified sampling, area by stratum was obtained independently from aerial photographs, classifying a large number of plots, or by some other adaptions of double sampling. This may be the earliest example of an application of C F I with SPR in forestry. One of the most influential papers on C F I with SPR was by Ware and Cunia (1962). They discussed the sampling and statistical aspects of the use of remeasured permanent plots and partial replacement of the initial sample for forest inventory. A thorough 10 Chapter 2. Literature Review conceptual presentation and theoretical derivation made the efficient sampling design for CFI with SPR initially proposed by Bickford (1956, 1959) appear to be a very attractive scheme. Simple random sampling was combined with partial replacement of sampling units on two successive occasions. The most efficient estimator for current volume at the second occasion was given, along with its variance. Ware and Cunia (1962) concluded that this weighted estimate was the most efficient linear estimate possible given the sample information under their assumptions of known variance. This means in statistical terms, it is the best linear unbiased estimate (BLUE) obtainable. In estimating the change in volume, Ware and Cunia (1962) presented several estimators and their ancillary variances according to various remeasurement procedures. The most efficient of these was the difference between the most efficient estimator of current mean volume at the second occasion and the most efficient estimator available at occasion two for the mean volume at occasion one. When using S P R it is necessary to determine how many plots to replace. This involves calculating the optimum sampling fraction. One set of optimum solutions was given for estimating current mean volume, and one for simultaneously estimating both change and current mean volume by Ware and Cunia (1962). They applied a graphical solution which minimized a linear function of the variances. A n example of the application of the derived theory was presented for a large scale forest inventory and it showed an increased efficiency and cost saving when using CFI with SPR. Bickford (1963) clarified the application of SPR to CFI by explaining what was involved in successive inventories and providing a basis for identifying the most efficient sampling procedure in relation to particular needs. He noted that sampling can be a major problem when facing the various situations, and that, although several theoretical studies had been completed, much was still unknown. The formulae given in his Chapter 2. Literature Review 11 paper were similar to those of Ware and Cunia (1962). However, he gave some rationale for choosing among the many types of successive inventories in order to obtain the information desired in useable forms, and at approximately least cost. Bickford et al. (1963) presented an efficient sampling design for CFI which was developed to provide data to meet a pre-set standard of accuracy at least cost for the forest survey of the Northeastern United States as the second survey was made. The essentials of the method were: an initial stratified inventory with some remeasurable plots, remeasurement of some of the initial plots, and establishment of some new plots in a stratified sample upon the second occasion. It was essentially a stratified double sample at the first occasion^ followed by S P R at the second occasion. Formulae were given for estimating total volumes and areas and their sampling error(variances of the estimates), and for calculating number of plots to establish and remeasure in order to meet an accuracy goal. The design was readilj amendable to local intensification, and the principle of SPR was r useful even if the first sample was unstratified. They stated that this design should be of particular interest to those concerned with repeated inventories of extensive forested area who would like to obtain the required data at least cost. Cunia (1964) continued work on SPR. He urged the application of CFI with SPR and pursuaded the user to ignore the complicated statistical aspects of the theory. At this stage of development, although C F I with SPR theory had been advocated by many researchers, very few applications of SPR to forest inventory had been reported. The literarure also differed as to which estimators should be used, and many researchers still worked on this aspect, trying to improve the theory. Cunia (1965) extended the theory of sampling populations on two occasions with S P R to multiple regression estimates. He showed that multiple regressions usually are better than simple linear regression estimates. Formulae were derived for the best estimators of current averages and changes in these averages from the first to second occasion and the Chapter 2. Literature Review 12 standard errors of these estimates were also given. An example was given of its application for the Canadian International Paper Company's limits in Quebec. The sampling method used was basically cluster sampling with systematic subsampling. Cunia (1965) also suggested that the theory of S P R on two occasions with multiple regression could be extended to SPR on multiple occasions. Letourneau (1966) reported a forest inventory by continuous forest control. Pleines and Letourneau (1970) further discussed C F I control systems, and later, Letourneau (1979) used simulation techniques to update C F I results. Frayer (1966) conducted a rigorous analysis of available data to test the validity of updating timber volume estimates by means of least-squares regression. Since CFI with SPR involves remeasurement of a portion of the initial sample, certain assumptions involving the distribution of regression residuals must be made to justify the use of regression techniques. These assumptions include independence, normality, and homogeneous variance of the residuals as well as freedom from error in the independent variables. Data from West Virginia and New Hampshire were tested to determine if the CFI's with SPR model developed by Ware and Cunia (1962) violates any of the assumptions of least-squares regression. Frayer (1966) found that all assumptions were justified except homogeneity of the variance. He employed a function of the expected current volume as a weight, and used weighted regression to stabilize the variance. His work suggested that the weighting procedure for estimating variances was more efficient. Frayer and Furnival (1967) presented a method of calculating changes in area from sampling with partial replacement. Previously S P R had been concerned primarily with volume estimates. They showed how the SPR formulae could be used to estimate area and they gave an example based on data taken from the forest survey of Vermont to demonstrate the procedures and efficiency of the estimate. Cunia. and Chevrou (1969) extended theory from two "successive occasions to three or Chapter 2. Literature Review 13 more remeasurements of the same forest population. This was mainly built on the work of Ware and Cunia (1962) and Cunia (1965). Formulae were given for the best estimates of the current values and change in these values from one to any of the subsequent measurements, and the variances of these estimates. However, C F I with SPR on three or more occasions seemed to be very complicated statistically as well as in its practical applications. For consecutive remeasurement on more than two occasions, the correlation structures may be complex and the regression assumptions may not be met. This limits the application of this theory. As we will see later, the C F I with S P R theory has many alternatives, but the applications remained mostly for two occasion cases. Frayer et al. (1971) gave an example of the application of S P R to timber inventories in the Central Rocky Mountains. The results indicated a substantial savings in time and money with C F I using SPR, especially when the correlation between the remeasurement plots is high. See (1974) examined SPR theory on successive occasions to establish the validity of the claims of increased efficiency in comparison with conventional systems. Several examples were discussed and he concluded that for estimation of a mean on more than one occasion, SPR was superior to conventional independent surveys in cost, and for estimation of both current mean and changes in means on successive occasions, SPR was again a superior system. Newton et al. (1974) extended the work on CFI theory from the univariate case to the multivariate case. General formulae were given and used to calculate stand tables for a particular forest area measured on two occasions. Problems related to practical application were discussed. Simple random sampling was the basic sampling design used in gathering the data. Hazard and Promnitz (1974) proposed convex mathematical programming as a method Chapter 2. Literature Review 14 of optimally allocating forest inventory sampling resources under different sampling plans to meet specified precision requirements for several variables.' Prior to this, optimum design was only discussed in a graphical sense by Ware and Cunia (1962). The solutions presented in their paper illustrate that the optimum replacement fraction in C F I with SPR can vary from complete remeasurement to large replacement fractions depending upon the specified precision levels, the population parameters, and the relative costs of obtaining information. Results of a sensitivity analysis were also given, along with an example. Dixon and Howitt (1979) used Kalman linear filters to estimate the current means. General formulas are given and they showed that the partial replacement estimator of Ware and Cunia (1962) was a special case of the Kalman estimator. They also proved that the best variance estimate of Ware and Cunia (1962) was always greater than or equal to the variance of the corresponding Kalman estimator. The theory was developed solely for measurement on two occasions, and only current value estimators were discussed. Peng and Zhu (1985) used Kalman linear filters to estimate the change in means. Two estimated values of the linear filter were used. In addition, the linear filtering formula, which was used in various SPR projects for C F I on two occasions, was given and some associated problems discussed based.on the results of forest inventories made in a forestry bureau of Jilin Province, China. Omule (1981) extended the work of Ware and Cunia (1962) from one-stage SPR to multistage SPR. He noted that theoretical developments in multistage S P R had been restricted to situations in which the sample size and variance are constant over time, and no application of multistage SPR had been reported in forestry. He derived the theory of multistage SPR without the restrictive assumptions of constant sample size or variance on the successive occasions. For practical purposes, Omule (1981) considered only the case in which partial replacement occurs at the primary stage of the multistage design. Chapter 2. Literature Review 15 The multistage design was based on simple random sampling within each stage. For estimators of the current population mean and change in means between two successive occasions, Omule (1981) derived estimators which were B L U E , together with their respective variances for a multistage structure. He also gave the biased estimators of the ratio form of the current population mean and of the change in means between two successive occasions, and made some comparisons among the different estimation methods. A n example of two-stage design was given based on data obtained from the Cranbrook Public Sustained Yield Unit in British Columbia. Results shown were convincing. He also suggested that this multistage SPR theory maybe particularly useful for inventories of large and diverse forest areas. Omule and Williams (1982) employed dynamic programing to determine replacement fractions. Most of this work was derived from Omule (1981). One example was given and the solutions showed that the optimum results were almost identical to those of other optimization techniques such as the graphical procedure. However, they explained that dynamic programming could be more appropriate since it could handle optimization problems in which the sampling extended over several occasions. Omule and Kozak (1982) discussed different estimators for successive forest sampling with SPR. They assumed that the coefficient of variation of timber volume remains constant over successive occasions and analysed previous estimators for the current mean volume. A n efficient estimator of the current mean volume was developed; it was a modification of the estimator by Ware and Cunia (1962). Results showed that the new estimator could provide gains in precision up to 200 percent when used for even-aged and normally distributed forest populations. The estimator was obtained using a simple random design. Omule (1984) presented the multistage SPR theory in a simpler and more understandable form. A two-stage sampling design on two successive occasions was used as Chapter 2. Literature Review 16 an example. He also noted the possibility of improving the precision by using variable probabilities to select the sampling units at the various stages of the multistage design, especially when the sampling units within each stage are of unequal size. Cunia (1987a) presented some theory on the error of CFI estimates and gave an illustrative numerical example. A n approach was shown that introduced the error of volume tables into the total error of the C F I estimates of average volume per tree, average volume per acre, and growth components. Cunia (1987b) extended the approach to C F I systems using SPR. He obtained S P R estimators which were slightly different from the classical estimators. The formulae are derived for the case of SPR on two occasions when (1) volume tables are constructed from linear regressions for which an estimate of the covariance matrix of the regression coefficients is known, and (2) the sample plots or points are selected by simple random sampling independently of the given volume regression function. Recently, Gillis (1988) discussed some problems in estimating change from successive static forest inventories. He noted that the estimated results were frequently masked by artificial differences and CFI should be done with detailed background knowledge. It can be seen from the literature reviewed above that there are many ways that CFI's can be conducted and that many estimation formulae exist. However, the basic sampling design in CFI's with SPR is the same. Equal probability sampling was always used whether the actual design was simple random sampling, stratified random sampling, or multistage sampling. In the following sections, multistage unequal probability sampling with SPR will be described. The S P R estimate procedures used here originated from Ware and Cunia (1962), Cunia (1965), and Omule (1981). The sampling structures for multistage S P R are assumed to be the same as those developed by Omule (1981). The unequal probability sampling theory, which was first developed by Hansen and Hurvitz (1943), is summarized Chapter 2. Literature Review mainly from Cochran (1977). 17 A detailed discussion of unequal probability sampling theory can be found in papers by Hendricks (1956), Kish (1965), and Brewer and Hanif (1983). Chapter 3 One-stage SPR with Unequal Probabilities 3.1 Definition of the Sample Suppose a population consists of N units of unequal size. A random sample of size n 1 of the ./V units will be selected with probability proportional to a measure or an estimate of size with replacement (PPSWR) on the first occasion. The probability of selection of the zth population unit is p xi = Ma/Mi, where Mu is the measure of size of that unit and Mi — Y^iLi Mu is the total measure of size of the N units in the population on the first occasion. Further, suppose a random sample (selected by P P S W R ) of size n of the n sample p units is retained for remeasurement on the second occasion. 1 The remaining part of the rii sample unit is n , and is measured only on the first occasion. In addition, a tl random sample of size n t2 of the TV — rii other sample units is selected by ppswr on the second occasion. The probability of selecting the zth unit on the second occasion is p i — M /M , 2 2i 2 where M i is the measure of size of that unit and M — S i L i " Mu. Thus, 1 2 there are rii — n + n p tl 2 sample units on the first occasion, and n — n + n 2 p t2 on the second occasion. Assume the variable of interest on the first occasion is A" and the same variable of interest on the second occasion is Y. The various observations and sample units can be defined in Table 3.1. In CFI's, interest centers on estimating the current population mean fly a n d the change in means A = fiy — fix- Suppose the observations on both permanent sample 18 Chapter 3. One-stage SPR with Unequal Probabilities 19 Table 3.1: Various observations and sample plots name occasion 1 Permanent sample units Temporary sample units Observations on permanent sample units Observations on temporary sample plots n occasion 2 n p p Tit, *i (i = 1, . .. vi Tip) (i = i , . . . ti) (» = Tl ) p Vi n ! , . . . , nt2) units and temporary sample units on the first occasion are described by the linear model: Si = Px + eii for i = l,...,N where X{ is the observation on the ith. unit, fix is overall mean of the observations, and en is the difference of the ith observation from the mean. Since the observations are made for different units, all the e^'s are assumed to be independent random variables with expected value 0 and variance <r\ . On the second occasion, observations on both permanent sample units and temporary sample units are described by a similar linear model: Vi - Mr + e 2i for i = 1,... , N- where yi is the observation on the ith unit, fly is overall mean of the observations, and £2i is the difference of the ith observation from the mean. Again, all the e i's are independent random variables with expected value 0 and variance 2 Chapter 3. One-stage SPR with Unequal Probabilities 20 <r\ . Further, it will be assumed that Corr(en, e ;) = p. The covariance between Xj and 2 yi therefore is Cov(xi, ) yi = E{[xi-E(xi)][yi-E(y )]} = E(e ,ex) { u = P<r <r . ci t3 From the sample observations after both occasions, the following four sample estimators are obtained: 1 " p n x1 p v i = 1 u Xi u i n n n i= i pu ' n nj 1 p i = 1 1 n t 2 t p i 2 2 i = l V2i where is an unbiased estimate of the population total based on the observations X p matched on occasion 1 and 2, t, is an unbiased estimate of the population total based on the observations unmatched on occasion 1. Y is an unbiased estimate of the population total based on the observations p matched on occasion 2 and 1, and Y t2 is an unbiased estimate of the population total based on the observations unmatched on occasion 2. The above unbiased estimates of the population total can be converted into unbiased estimates of the population mean by dividing the population size N. Chapter 3. One-stage SPR with Unequal Probabilities * n 1 21 l - = J _ f 1L ^ Nn fr[ p -I P 2 i "to i = l P2i t2 Nn To derive the variances of these estimators,- let ti be the number of times that the i t h unit appears i n a specific sample n . For sampling with replacement, the i t h unit may p appear 0, 1, 2,.. . ,n p times in the sample. The probability of the zth unit being selected is pu in each draw. This is a multinomial problem. Consequently, distribution of the the joint frequency for all N units in the population is the multinomial expression: . , , , • .P11P12 • - -PIN- For a multinomial distribution, the following properties are well known: E(U) = V(U) = Cov(ti,tj) np, p u n pu(l p = p ),and u -Uppupxj. It follows that, 1 Nn . Xi -(tipn p 1 N n p A2 h tP12 h • • • PIN + ijv" N Y "1 PU 2 The only difference between the formulas X ^-Nn~ ti^ P and 1 Xif ^ . X i Mnp^l 'pu Chapter 3. 22 One-stage SPR with Unequal Probabilities is that the sum extends over all units in the population in the later formula. In repeated sampling the £ 's are the random variables, while the Xi and pu are a set of fixed numbers. t Since E(ti) = n pu, hence, p E(x ) = £ ( ^ _ f > p p P l j ) ^ ) = (X) = fix • E Therefore, the variance of x is: p V(x- ) = P A 1 vi—Zu Xi —) pf^l Pli Nn \£(-) V(U) 1 -jH- ^ + 2 2J:^Cov(t ,t )] i E(-) p«(i - P * ) a 2 - E - ^ p i ^ i , ] £2 P % £i P i i (E —-* ) Ji P!i 2 -^ p i l l P l i 2n 1 Y N £ P I . ( - - * ) N Tlp itl 2 S Pli where 1 s ^ i = 1H52PK{— A'- »=i P i i i V - X) ,and 2 The unbiased estimate of V(x ) is: p Similarly. i E(x ) tl = 1 E(—— N n u N , pu)—) , X; — E(X) = ^2(n i=i PIT ti j p x Chapter 3. One-stage SPR with Unequal Probabilities E(y ) = E{^-Y,{n p )^~) = E(Y) = fi E(Vt ) = E(-^J2(n p )^) = E(Y) = P 2 p 2i t2 Y 2l i. l Y and The unbiased estimators of these variances are: 1 iV n (n 2 t l -1) t l 7 = -. A/2, N n (n 2 p p \ ^Mi-gpf nt .1 ^ n (n . -1) 2 2 t 2 t 2 and the covariance between mean of the sample plots observed on both occasions Cov(x ,y ) p p = — N Tl 2 ~[fr[Plr^P2i 2 p »2 where n 1 p "r 1 The CFI estimators can now be derived using these terms and formulas. Chapter 3. One-stage SPR with Unequal Probabilities 3.2 24 Estimator of the Current M e a n The most commonly used (Ware and Cunia, 1962; Omule, 1981) estimate of the current mean py on the second occasion is obtained by a linear estimator of the form: y = ax b where a ,b ,c , x x x x + bx + cy + dy u x p x p x t2 and d are constants. x To make this estimator an unbiased estimator of u-y, it follows that: E(y ) = py. b Since, E(x ) tl = E(x ) = p p x E(y ) = E(y ) = fly t2 p therefore, E(y~b) = = a-ipx + b p x (<*i + b )fi t x X + dp ci/xy + x Y + ( c i + d )p x Y then: x = 0 x = a +b x c -\- d x 1. Hence. Vb = a (x x tl - x ) + c y -f (1 - c )y . p x p x (3.1) t2 The variance of y;, is: V(y ) b = al[V(x )-r = g 2 ( i L + f l ) + 2 j j + (1 -Cl n n rip . V(x )) + c\V%) + (1 - c ) V(y ) 2 tl p C tl p x )2ii2a n t2 t2 e l C i - 2a c Cov(x ,y ) x P ^ . n a l p x p • p (3.2) (3.3) Chapter 3. One-stage SPR with Unequal Probabilities 25 In order to calculate ay and Cy and minimize V(y ), the equation (3.3) is differentiated b with respect to ay and Cy and the equations set equal to zero. That is: d[V(y )] b = d[V(y )} d(cy) 0 = 0. b Therefore, 2 f l l (il n i l ) - 2 + n tl n Tii p C l p ^ ^ = 0 np Tip 2 Solving the above two formulas simultaneously, after simplification, yields: _ T)s n\n 7 n ssnn p d = 1 2 1 2 l £l tl (3.4) 2 i2 1 (3.5) n Ti T) - 2 8pa n n r] L2_LP — np s- S2n n where r/ = ix - 2 2 1 tl i2 a . £2 The relationship of a and Cy can be written as: a TjTl Ci tl a = y SyTlyTlp . Substituting the ay and cy into (3.1) yields the best estimate of the current mean: y b = ay(x tl - x) + cy .= x p {x tl p + (1 - - Xp) + cyy p cy)y t2 4- (1 - cy)y t2 SyTlyTlp = (x - x ) +y ] cy{ VTltl synyn tl p p + (1 - cy)y . t2 p Notice that the overall mean of the population on the first occasion, x— n x tl tl + n Xp p Tly n ti - — x Tly tl H, P-x . n Tly p Chapter 3. 26 One-stage SPR with Unequal Probabilities Since n = n — n , the above formulas can be written as: p x tl x — xv (x ni - x ). tl p Thus, Vb = c i [ {x - x ) + y ] + (1 - c )y . p SlTl p x t2 p Note that: r, sn x _ Cov(x ,y )n _ — Pi 2 p p p V(x )n n p p p p where j3i is similar to a regression coefficient. Therefore, C ! [ ^ i ( i - x ) +y ] + (1 - ci)y Vb = p where y = /3 (x - x ) x x p Ci-Vi + (1 - c )y x t3 t2 +y . p p Substituting the calculated values of (3.4) and (3.5) into (3.3) yields the general formula for calculating the minmized variance of 'V(i/ ): 0 n p ^ i f i z n x ^ - n n rj' tl n sis n n - n n r] 2 p + 2 x n t2 n 2 2 tl t2 ^ n tl p p ^ S 7 i p s ^ ^ i " ' 2 — n n rj 1 tl 2 p 1 2 t? t2 - n n T]^ n ^s s nin 2 n ssnn A t2 tl 2 n 2 2 p x 2 x Notice s s n nln x 1- d - 2 : x - tn TOp s 5 n Ti 2 1 2 1 2 " - n n rj •-- • tl t2 n n rj tl t2 2 2 - n n T]* tl t2 n* v Chapter 3. One-stage SPR with Unequal Probabilities 27 After simplification (See Appendix A ) , V(yi) can be expressed as: y(fc) = — ( 3 . 6 ) If n = n — n and n — n = n (the number of sample plots on both occasions is the tl t2 t x 2 same), then: T]s n nl 2 nnss 2 t — 2 x 2 nnss 2 2 2 2 t - 2 1 n r] n 2 t n ' 2 Substituting these values into ( 3 . 3 ) yields: n n£s s — nil] z 1 i n 1 2 i 2n ns s f 2 p n^n^Si - nfrj ns n n li 2 t n n SiS s 2 ns s nl p x 2 2 — n r} 2 p x 2 2 t rj 2 nn ss 2 p 2 n p — n\j] 2 n 2 2 n 2 2 2 — nfr z s 3 n n£s- s t 2 n 2 p Notice _ nn n sis - 2 t p 2 n n£siS ^bi) c a n 2 2 — n^rj^ z After simphfication, the variance V(y~b) n r) 2 be expressed as: = -(l-c ). n n t Another assumption often made in CFI's states that the variances on both occasions are the same. This means: = s n, = n n = n •Si t 1 = s 2 2 h - n — n. t Chapter 3. One-stage SPR with Unequal Probabilities 28 The ai and C\ values in this case become: •qsn n 2 v t 0-12 = 9 9 9 n n s' z p i 2 9 — nirj- ns nl 2 Cl2 = nns 2 — n rj 2 2 p 2 Substituting these values into (3.3) yields: V(y ) b2 = ( ^ ^ +( w 2 2 . Vsn n Noti ce: n n p ns n 2 t z — 2 2 p 9 9 9 n ripS' 2 fL) -+(i-a- n^n^s* — rtfTj^ V \ 2 9 'V — niT]* 9 •> 9 n n~s z nn n s 2 t nns 1 2 2 2 — 2 2 z — t 3 p v 2 9 / — nirj z 9 ' n p nn 2 n rj t 2 2 2 After simplification, the variance V(yb) can be expressed as: V(y ) b2 3.3 = - ( 1 - c ). 12 Estimator of the Change The most commonly used (Ware and Cunia, 1962; Omule, 1981) estimation of the change in means between two occasions, A = py — px, is by a linear estimator of the form: A b = a[x tl + b' x + c[y + d\y x p p t7 where a\, b[, c[, d[ are constants. To make this estimator be the unbiased estimator of A , it follows that: E(A ) = u. - px = A . b Y Chapter 3. One-stage SPR with Unequal Probabilities 29 Since, E(x ) = E(x ) = p u p x E(y ) = E(y ) = p t2 p Y therefore, E(A ) b + Kfix = a[nx + W + = [a\ + b[)px + (c[ + d'^py C d W then. a[ + b[ = -1 'i + = 1- c Hence, A (3.7) = b[x + c[y + (1 - c p y , - (1 + b[)x . p b p tl t The variance of A(, is: V(A ) = b b^^ + cfVtiJ + (1 + b[) V{x ) + 2b' c' Cov(x ,y ). + 2 tl il-tfVfa) 1 1 p p Substituting the different variances and covariance into the formula yields: V'(A ) 6 = b[ fl n 2 + cftl + (1 - e p z i l + (1 + 6 p 2 . f i -r 2 6 ' ^ ^ ' ^ . (3.8) n n n np1 p p t 2 ( l In order to calculate b\ and c[ and minimize the variance V(A ), b the above formula is differentiated with respect to b\ and c[ and the equations set equal to zero. That is d[V(A )} b d{b[) d(c[) = 0 = 0. Chapter 3. One-stage SPR with Unequal Probabilities 30 Therefore, 26' ^+2(l+&' )-^ + 2 c ; ^ ^ = 0 1 1 Tip Tit i Tip 2 ^ - 2 ( 1 - ^ +2 6 1 ^ ^ n Tip n h - 0. 2 p Solving the above two formulas simultaneously and simplifying yields: , _ -nls ( n Tip 2 s s n n nl 1 , _ ei 1 n n Tj 2 t2 tl + n r]) t2 S j s ^ n ^ l (3.10) n n ri ' 2 - tl and the relationship between e2 tl - 2 nls^s^rip 1 where r/ = 8pa a , 2 n T)) + 2 Sl i2 and c[ can be seen to be: y _ _ r i p _ n rjc\ tl Tii 1 <SlTiiTip Substituting the calculated values of (3.9) and (3.10) into (3.7) yields: A b = b'^Xp - x ) + c[y + (1 - c'^yt, - = (— tl ){ P ~ t ) X 'i{y c P x tl + c-iVv + (i - i)yt - x c 1 Tli = p 2 tl -^^-{xp ~ * t j ] - — (S - s ) + (1 - ci)i/t - S P 5lTiiTi tl Tii p Notice that the overall mean on the first occasion is: x n = x tl + ripip tl ^ x t l ^ ^ l x p + Tii = n (x Xp + tl Til Hence, x S\TliTlp l - Xp). 2 t l Chapter 3. One-stage SPR with Unequal Probabilities 31 and x - p 1- = —(x -x \ ). 7 1 1 x tl tl Also, V Cov(x ,y )nl SjUp V(xp)n n p p p "•- Pi • p Therefore, &b = c[[y - Pi(x - x)] - (x - x ) + (1 - c[)y - x = C'IV'I + = y'b-x p tl p tl t2 (i - i)y<2 - * c where: A ( P - x), and V'l =y - s P y' = c' y' + b 1 (l-c' )y . l 1 t2 This result means that the best estimate of the change on both occasions is the best estimate of the current mean on second occasion minus the overall mean on the first occasion. Now substituting (3.9) and (3.10) into (3.8), we get the general formula for calculating the minimized variance of V (Ab)- n V(A ) B = 2 p s ( 2 S n l I n 2 + p n n ) u ; s s n n n^ 1 2 1 - n n T) 2 tl +l ~ 2 t? 2 p s 2 2 n ( s 1 1 n n 2 n ^ 2 n p 1 2 s 2 - n n t i l t n t 2 Tj T ] ) ^ i n i t 2 r} 2 tl 1 n^ (5 n n •+ 1 2 1 2 7 n t2 p 1 2 + p n t l rj) i 2 s i p 2 — -n^ (ain n 2 1 2 - n n rj p h ssnnn l ; l 1- 1 H n z s h P 2 ; - \ - n - SiS-iU^nl s s nin n p + n^ri) 2 s i s ^ ^ - n n ^ i ( s n i » + n v) ; l n z n l s ^ s ^ r i p l S 2 p s n 2 n t 2 - n n r} tl t2 2 1 n 2 n * - n t i n t 2 rj 2 S l — n t 2 rj . v) n 2 p This can be simplified (See Appendix B) to give: .V(A ) = — ( l - c i ) + - ^ . ( l - ^ c i . ) . n TI a k t2 (3.11) 32 Chapter 3. One-stage SPR with Unequal Probabilities The best estimate of variance is a. weighted estimate of the variances divided by the number of temporary sample plots on the second occasion and the total number of sample plots on the first occasion. When n t l =n f 2 = n and Tij = n = n, b\ and c\ become: t 2 -s nl(s nn 2 b n , x = sis n n _ n s (s nn 2 2 2 11 p x n rj 2 + p 2_9 SiS n n z t - 2 2 2 + n rf) p p 2 n n) t 2 9 — ' ni7] i Substituting the above values into (3.8) yields: \/(A V{A ) bl \ = -s nl(s nn — 2 1 + n ri) p 2 SiS^Tlp t _ S l — 2 — TifT/^ _ nl (s nn Sl 2 S\s n n l 2 v -s^Ksynrip I 2 9 - nfri l s n z 2 2 Tip p x p „ , + n n) t 2 9 s^n^n* - nfrj t2 x 2 H 2 S J l J n tl + n v) y p t 9 9 — niT)* SiS n n~ Z 2 -s n (s nn n s (s nn p 2 n{TI [1 T 2 t 29 — 2 T 2 t P r~2 2 J +n T/) 99 2 SiS^Tl* tV 99 n rj) s 2~1~$ ~ n*ai(* wra + Tip +n ) p '[ H 2 9 J 9' — niv/ SiS n n^ z 2 n p After simplification, the variance V"(A ) can be expressed as: V(A ) = ^ ( l - c ' ) + ^ ( l - / 3 c i ) . Dl D l 1 1 Ti 1 n t 1 In CFI's, the assumption is often made that the population variances are the same on both occasions in addition to equivalent sample size. This means: S1 S2 — s -— fit, — n = n, tn n-i n = n. — 2 The b[ and c[ values become: , _ -5Tip(^TiTip + snn 12 2 2 2 n rj) t — nfrj n s(snn +n rj) 2 p snn 2 2 2 t — nfrj 2 2 33 Chapter 3. One-stage SPR with Unequal Probabilities Substituting the above values of b' and c' into (3.8) yields: 12 _ ' \*-*b2 ) -sn (snn + n n) 2 p [ p n 12 l 5 7 l n^(a7in snn 2 2 s 2 p p J 2 9 s n n^ — nfr) n s(snn n z + n 7/) t — n 7/ 2 p ' I s*n*n* — nil} j 2 h + p 5 7Z n 2 snn 2 t p — 71 7? 2 + -^(j57t7i t 2 2 2 2 2 2 2 M 1 2 1 - ^ ) - p + ^7?) p 2 — t 2 This variance can again be simplified to: v(A ) = i - ( i - c ' ) + r ( 2 + n 7/) 77 5 7 Z n — 7l 77 2 n 2 sn (sran + n rj) n s(snn 2 s t 2 9 J 9 9 9 p n 2 2 +2[- + n n) 2 t 7Z 2 7i r/ 2 2 s n t Chapter 4 Two-stage S P R with Unequal Probabilities 4.1 Definition of the sample As in the one-stage case, suppose the population consists of TV primary units (pu). W i t h i n each of the TV sample units, there are M secondary units (su). Both the TV pu and the M su within each pu are of unequal size. O n the first occasion, a random sample of size rii of the TV pu is selected by ppswr. The probability of selection of the i t h pu at any of these draws is pu — Mu/Mi, where Mu is the measure of size of that unit and M i = Y^iLi Mu is the total measure of the size of the TV units. W i t h i n each selected rii of the TV pu, a random sample of size m of the M units is assumed to be drawn with probability proportional sample to a measure or an estimate of size without replacement (ppswor). If a particular su is drawn more than once, the whole sample is replaced, and a new independent drawing of m subunits, again without replacement, is made. Suppose a random sample of size n p of the n pu is retained for 1 remeasurement on the second occasion together with its respective su drawn from the first occasion. A n additional random sample of size n t2 of the TV — rii other sample units is selected with ppswr on the second occasion. W i t h i n each selected pu, m su are selected again by ppswor. The probability of selecting the i t h pu on the second occasion is p i = M i/M , 2 The 2 2 where M ; is the measure of size of that unit and 2 M 2 = Hili" 1 M i. 2 various observations and sample units are summarized i n Table 4.2. The total number of sample units on the first occasion is riim = (n + n )m p 34 tl and Chapter 4. Two-stage SPR with Unequal Probabilities 35 Table 4.2: Various observations and sample plots name Number of matched pu Number of unmatched pu Number of matched su Number of unmatched su Observations on permanent sample plots Observations on temporarj sample plots r occasion 1 n p u nm n m n occasion 2 n n nm n m p h p p tl t2 ™' (i = 1,... , n ) (j = 1, . . . ,m) p (i = (j = (i = 1, ... ,n ) (i = (j = 1, .. . ,m) (;' = tl V'ij l,...,n ) 1, . . . ,m) Vij 1, . . . , n ) 1,... ,m) p i2 the total on the second occasion is n m = {n -f n )m. From the observations on both 2 p t2 occasions, the following four unbiased estimators of the population means are obtained: 1 ripMi i= l Pii 1 v MuXi n t , M i Zpu i=l 1 M y'i Tij, 2i 7l M p P2i 2 1 Mnyi n M 2- Pu t7 2 where x is an unbiased estimate of the population mean based on the observations p matched on occasion 1 and 2, x is an unbiased estimate of the population mean based on the observations tl unmatched on occasion 1, y is an unbiased estimate of the population mean based on the observations p matched on occasion 2 and 1, and Chapter 4. Two-stage SPR with Unequal ProbabiHties y t2 36 is an unbiased estimate of the population mean based on the observations unmatched on occasion 2. To find the variances of these estimators, let ti be the number of times that the zth unit appears in the specific sample n . Since ppswr was assumed i n selecting of pu, the p ith pu may appear 0,1,2,... ,n times i n the sample. The probability of the i t h unit p being selected is pu in each draw. This is a multinomial problem as i n the one-stage case. The between-unit component can be written as: 1 Y N Its variance, V j , can be obtained by the same method as i n the one-stage situation, that is: 1 Y N pMi Pi* ~i n The contribution of the su within a pu to the variance for a unit that is drawn once is 1 M (l-f )Sl u 2 where f = ^ is the sampling fraction in second stage, and 2 1 M '* l j=i 1 Its variance, V , is calculated as: 2 n M\ £ p^m p • A d d i n g the variances obtained from pu's and su's to get the total variance of the twostage sampling yields: ivi \p n i i = 1 V^i^ = 1+4 P n y Pit Ml(l-f )Sli 2 p m. 2 u Chapter 4. Two-stage SPR with Unequal Probabilities 37 where 1 £ M^^V- " 1 ,Xi and 1 ^M {l-h)Sl E u The unbiased estimate of V(x ) is obtained as: p where / = jjj is the samphng fraction in first stage, and X{ is the sample total in the ith x pu and x is the average of the X{. Similarly, M = i hfri Pu n n tiMifr{ 2 P U -fLlfi VQ \ f 1 M i n n P f^i ( Yl P2i Y) + 1 2 nM f -'>* 2 fr[ 2 p pliiri ' 4- s " 2 _ m ' 2 A 2 ^ 2 s S and the unbiased estimators of these variances are: H " p) = (n -l)l M2) ^~* 2 P h )2 5 38 Chapter 4. Two-stage SPR with Unequal Probabilities The co variance between the remeasured sample plots, using the correlation structure given by Omule (1981), is expressed as: Cov(x' ,y' ) p = (3 /n p 2 p where (3 = {m-picr^a*, + p a a )/m. 2 2 ei £2 Pi is the true population correlation coefficient between the primary units on the two occasions and p is the true population correlation coefficient between the secondary 2 units within the primarj' units on the two occasions. The CFI estimators can now be derived using these terms and formulas. E s t i m a t o r of the C u r r e n t M e a n 4.2 The most commonly used (Ware and Cuina, 1962; Omule, 1981) estimate of the current mean py on the second occasion is obtained by a. linear estimator of the form: Vb = 2^t, + b &p + c y + d y a 2 where a , b ,c ,and 2 2 2 2 p 2 i2 d are constants. 2 To make this estimator an unbiased estimator of uy, it follows that: E{y )=p b Y Since, E(x ) = E(x ) = p E(y ) = E{y ) = tl h p p x Chapter 4. Two-stage SPR with Unequal Probabilities 39 therefore, = ap E{y ) b 2 = + b px + c p x 2 2 Y + dp 2 Y (<l2 + & 2 W ' + (c + d )p 2 2 Y Thus, a +6 2 = 0 2 C2+^2 = 1. Hence, Vb = 2i h ~ P) + 2 2 / x a X + (1 ~ 2 ) y - (4.12) c C t2 P The variance of y is b ffo) = al{V(x ) + V(x. )} + clV(§ ). tl p p +(l-c ) V(y: )-2a c Cov(x ,§ ) 2 (4.13) = t2 2 2 2 p p 2/4+4' , 4+4'\ , 4 + 4' . „ \2 2 + 4' S 2 = a{ 1 2 J+ c h(.l —c ) 2 2 0 „ „ ^2 7l .. ., 2a c — -(4.14) 2 r 2 In order to calculate a and c and minimize V(2/f,), equation (4.14) is differentiated 2 2 with respect to a and c and the equations set equal to zero. That is 2 2 d[V(y )] d(a ) = 0 d[V(y )] d(c ) = 0. b 2 b 2 Therefore, 2 o ! ( f i ± n 2c 4 + 4' 2 2 2 Tip i I + £ i ± i n tl , - 2(1 - c ) 2 I ) 2 c f t n p ,4+4' „ 2 2 n = - 2a — 2 i 2 o p 02 ra p = 0. Chapter 4. Two-stage SPR with Unequal ProbabUities 40 Solving the above two formulas simultaneously, after simplification, yields: K + s i)( 2 4>in - n n ^ + s 2 tl { t2 ' ; The relationship of a and c can be written as: 2 2 _ f3 n c 2 tl 2 ° " (*i + * i > i ' 2 Substituting (4.15) and (4.16) into (4.12) yields the best estimate of the current mean: Vb = 2{x - i ) + c y + (1 - c )y a 2 p h /3 nt,c 2 p 2 t2 2 Notice that the overall mean of the population on the first occasion, = —aj rii t l + — V ri! Since n = Tlx — n , the above formulas can be written as: p tl = _ = - ^ i i ^ X Xp y^t\ — _ = ^ *^p)' Hi Thus, yb ~ C2L „Ax ~ x ) + y ) + (1 - ci)y . P \ 1 -T l ) S S Suppose, ft ft p t2 Chapter 4. Two-stage SPR with Unequal Probabilities 41 where j3 is similar to a regression coefficient. Therefore, r Vb = c [/3 (l - x ) + y ] + (1 - c )y 2 = p r p c y + (1 2 2 t2 c )y , 2 2 h where y = f3 (x - x ) + y . 2 r p p Substituting the calculated values of (4.15) and (4.16) into (4.14) yields the general formula for calculating the minimized variance of V(y ): b V(v) - ( /3 ^' S 2 '^ ^ + S W-> | -»V U (*i + *"){*2 + 4)nin M (s\ + s")(s + 4 ' ) i 2 - n n f3i n 2 + _ 2 /j _ (*i + s")(*2 / ^2(4 + (*i + s t2 tl ")( 2 s , 2 (5 2 w n ) p + agViTip + s' ') 2 (*1 + *l)(*2 + •S2) l -p n ~ n n (3% (si + sJ'Xs'rj + s ).*in - + 2) i 2 3 1 + s'i) n p n | 1 n i3 tl 2 2 jfo T i n n 0l t2 tl n' p Notice (s[ + s'{(s' + s'Drixnt, - n n /3 2 1- c 2 tl (s[ + s")(s + s ' ' ) i 2 - n n f3 n n 2 2 2 t2 t2 tl : After simplification, V(y ) can be expressed as: b V(y )= -^^(l-c ): { b If n tl 2 — n = n and n = n — n, then: t2 t x 2 a /3 (s' + s )n n 2 21 n ^ 2 2 + '/(5' S 2 2 p + 4') t - V& 2 71(5; + 4)(4 + 5'')n 2 C21 p n W+«I)('2+4)-'H% ' 2 2 Substituting these values into (4.14), after simplification, the variance V(y ) can be b expressed as: m 1 ) = ^ ± ^ ( i 71, C 2 1 ) . Chapter 4. 42 Two-stage SPR with Unequal Probabilities Another assumption often made in CFI's states that the variances on both occasions are the same. This means: i i a ' 711 = i " n. 71 = 2 The a and c values in this case become, 2 2 (3 sn n 2 &22 p n V t n /3 - 2 2 ns n 2 c p ns 2 22 - 2 n (3l' 2 t Substituting these values into (4.14) yields: v -fasripnt x r 2 l n 5 2 2 s 2 - n /? 2 0 sn n 2 n 2 ™ 5 2 t p t 2 - n /3 2 2 5 n 2 2 P 2 , ,+ ( ns n 2 p n , -n /3 2 2 2 (3 sn n 2 p n s -n {3 [ 2 2 2 2 ) ( - 1 C 2 2 ns n t 2 n s 2 ^ 2 - n (3 2 ) J3 p ) [ 5 2 + 2 2 t s ^ 2 } 2 2 n p Notice nn s - n /3 ns — nt /3 2 1—c 22 t = 2 2 2 2 t 2 2 After simplification, the variance V(y ) can be expressed as: b V($ ) = b2 4.3 -(l-c 71 2 2 ). 1 Estimator of the Change The most commonly used (Ware and Cunia, 1962; Omule, 1981) estimation of the change in means between two occasions, A = py — px > is estimated by a linear estimator Chapter 4. Two-stage SPR with Unequal Probabihties of the form: A = a' x + b' x + c' y + d' y 2 b tl 2 p 2 p 2 t2 where a' ,fc ,c ,d are constants. 2 2 2 2 To make this estimator be the unbiased estimator of A, it follows that E(A ) = A. = ny - p b x Since, = E(x ) = fi E(x ) tl p x Therefore, — a fi E(Ab) 2 = + b' p + c' p + d' p 2 x x 2 Y 2 Y + (4' + d' )p . (2 + a 2 Y Thus, a = -1 +K 2 c' + d 2 2 ~ 1. Hence. A b = b' x + c' y + (1 - c' )y - (1 + b' )x . 2 p 2 v 2 t2 2 tl The variance of Ab is V(Ab) •= b' V(x ) + c' V(y ) + 2 2 (l-c' YV(y ) 2 p 2 p 2 t2 +(1 + o ' ) V ( I ) i 26 c C (x ,z7 ). 2 2 = tl 2 2 OT p p Chapter 4. 44 Two-stage SPR with Unequal Probabilities Substituting the different variances and the co variance into the above formula yields: ' f(A±fll 2b' c' ^. {1+b 2 + 2 (4.18) 2 In order to calculate b' and c and minimize the variance V(A ), 2 2 b the above formula is differentiated with respect to b and c' and the equations set equal to zero. That is 2 2 d[V(A )} b d{b> ) = 0 = 0. 2 d(c> ) 2 Therefore, 26 ;W±ii) + 2 ( 1+ ^W±ii) rip 24' ' a n ^-2(l-4)^ + + Tip ft + 24 n tl p ''+2^ 5 Tlt „ = = 0. Up 2 Solving the above two formulas simultaneously and simplifying yields: b' -n (s' p = + s' '){{s\ + s'{)n 2 2 (s[ + s'{)(s n (4 (4 + 4')(4 + s ')n i 2 + s ) in tlf 2 tl t2 +n (3 ) 1 t2 - n 2 n 3) n n f3l n 2 +s'l){(s' p = + + 5- ') i 2 - 2 c, 2 n 2 2 n n f3l tl t2 The relationship of b' and c' can be seen to be: 2 2 rii (4+4>i 2 n (3 c' TJ T Til 2 2 n : Substituting the calculated values of b' and c' into (4.17) yields: 2 A fc 2 = b' (x - i t , ) + c y + (1 - c )j/ - x 2 p Tip n r 2 n (3 c tl 2 p 2 2 •)(lp - l , = ( 5 ; + 8'{)n t2 t ] tl ) + c' y + (1 - c )y = 2 p 2 t2 -i t x n uPi ,= • — \ 1 n , p = = .• , H ,, = l 45 Chapter 4. Two-stage SPR with Unequal Probabilities Notice that the overall mean on the first occasion is: nx tl + nx tl p T> p , , -•• n , (ri! - n ) tl t l = — i ^ti Tii 7i! = _L hL(= Xp i v""'i — — = I f =^ V)' T X Hence, and Xp - x = t l x J. —-(i - f Also, Substituting the above results into A yields: b A fc = c' [y - 0 {x - x)\ - (1 - l 2 p T = c y' + = y'b-x 2 v t l ) + (1 - c' )y - a 2 t2 £l (l-c' )y -x 2 2 t2 where y ' = y - f3 (x - x). 2 Now p T p substituting the calculated values of b' and c' into (4.18), we get the general 2 2 formula for calculating the minmized variance of: = + s")n +n /3 )M ^-Ms'z-r tl 2 + s'{) 2 ( s i - f s " ) ( s + s ) r i i 7 i 2 - n n (3 2 tl 2 t2 n 2 p n (s[ + s';)(( ' + s'Dm + n /3 ) f Y + s ') P + (si S + s'{)(s' + 2 2 la s )7ii7i 2 2 -n n Bl tx i2 2 l 2 2 2 n p 46 Chapter 4. Two-stage SPR with Unequal Probabilities + [1 _ ntM M + s' ') (s[ + s'{)(s' + a2 )ra1n2 - n n /3 n -n (s' + s' ')((si + s'{)n + n /3 ). {s[ + s'() % ( - i+ s 4>i + + 2 / 2 tl 2 p fl [ 2 2 ( i + 5 ,-n (5 p 2 2 • i')( 2 s 5 t2 2 2 tl 2 + •s' ')n n - n n / ? 2 1 + s ')((sj + s j > 2 (s[ + s'{)(s' + s'i)nin 2 2 2 tl tl 2 2 n 2 t2 + 7i /3 ) 2 2 t ] 7i (sj + sj')((s + s > ! + n f o ) 2 p (s[ + s'^s!, + - n n /3^ tl tj t2 2 s ')ni7i 2 ta 2 - n n f3l tl t2 n p This can be simplified (See Appendix C) to give: \/(A ) = ^ ± ^ ( l - c fc 2 ) + M+iil(i-^ ). Ca (4.19) The best estimate of the change in means can be seen to be a weighted linear combination of a regression estimate, the mean obtained from the temporary sample units on the second occasion, and the overall mean on the first occasion. The best estimate of variance is a weighted estimate of the variances divided by the number of sample units on the different occasions. When n tl =n t2 = n and n — n — n, b' and c' become: t x 2 1 _ ~ 2 2 ~(s' + 4 ) p ( ( * l + 2) + tP2) {s[ + s'{)(s' + s > - 7 i / 3 n S n n 2 2 2 n c 2 p( 'i + (s> +s )(s s 1 1 2 + s' ')n + n (3 ) + s'i)n*-n?ft ' 2 2 2 2 t 2 Substituting the above values into (4.18), after simphfication, the variance V(Abi) can be expressed as: F(A !) = ^ ± ^ ( l - c b Tif 2 1 ) + ^±il(l-/5 c 71 r 2 1 ). (4.20) In CFI's, the assumption is often made that the population variances are the same on both occasions in addition to equivalent sample size. This means: (si + s") = (s' -f s ) = s 2 n, = •n t2 t Tii — 2 = n n — n. 2 t Chapter 4. Two-stage SPR with Unequal Probabilities The b' and c' values become: 2 2 _ -sn (sn 2 2 ~ sn 2 2 ~ + n (3 ) p 2 t 2 - n f3 2 2 2 s n -n (3 2 2 2 ' 2 Substituting the above values of b' and c' into (4.18) yields: 22 , . ^ > b 2 -sn (sn + n f3 ) p L 2 S n 2 22 t _ n 2 / 3 2 J 2 s 2 U n s{sn +n /3 ) p + i p _ n s(sn +n /3 ) t p + l 2 5 n 2 - n /3 2 2 2 2 n n 2 i _ 2p% 2 n a 2 J g t n s 2 p -ara (aTi p + [ 2 t 2 2 2 - , s n ( 5 n 4- n /3 ) n s(sn + n (3 ) (3 t p [ s n -n /3 2 2 2 2 2 p l [ t s n -n (3 2 2 2 2 2 2 l n p This variance can be simplified to: y(A ) = - ( l - c fc2 n t 2 2 )4--(l-&c n 4- n ^ ) ^ - , s n - n (3 + 2 2 ). 2 1 ' Chapter 5 A p p l i c a t i o n to Chinese Forest Inventory 5.1 B a c k g r o u n d on Chinese Forests and Inventory P r o b l e m s China is the world's third largest country covering an area of 9.6 million square km. Its land extends 5,500 km from north to south, crossing the tropics, subtropics, warm temperate zone and cold temperate zone, and 5200 kilometers from east to west, stretching upwards from the Pacific all the way to top of the world-Mount Everest situated in the hinterland of the Asian continent. China's major landscapes are mountains, high plateaus, deserts and extensive river plains. Its varied geographical and climatic condi-. tions give it perhaps the greast variety of plants in the world. The forests in China are highly variable. In the south (Hainandao, Guangdong, Guangxi, Yunnan and Taiwan), there are large areas of tropical monsoon rain forest and small stretches of tropical rain forest. Going progressively northwards, one finds subtropical humid evergreen broadleaf forests, warm temperate zone deciduous broad leaf forests, temperate zone mixed broadleaf-conifer forests and finally, temperate and frigid zone conifer forests. In addition, there are large mountain conifer forests on the Tibet-Qinghai Plateau and the Inner Mongolia-Xinjing Plateau and large diversiformleaved poplar forests in the Talimu River Valley in Xinjing. There are also desert forests on the fringes of deserts and dryland areas. According to the latest forest inventory results released by the Ministry of Forestry in 1984 (Ministry of Forestry, 1984), China's stocked forests cover 115,277,400 ha which 48 49 Chapter 5. Application to Chinese Forest Inventory is 12 percent of the total land area. Thirteen provinces and autonomous regions are over 20 percent forested. Table 5.3 illustrates the forest types and their amounts. Table 5.3: Forest types and their amounts in China Type of forest Area (10,000 ha) Percent Protective forests Timber forests Fuel forests Special forests Economic forests Bamboo forests Total 1000.24 8062.84 369.09 129.90 1128.04 319.96 11527.74 9.1 73.2 3.4 1.2 10.2 2.9 100.0 It is estimated that China has about 250 million ha (26 percent of its total land area) suitable for forests (Hsiung, 1980), but one-third of that still remains unforested. Table 5.4 illustrates the problems. Table 5.4: Forest land utilization in China Utilization Million ha Percent Nurseries New plantations Thinly-stocked forests Shrub lands Nonforested lands Stocked forests Total 0.12 4.15 15.63 29.75 85.82 115.27 251.01 0.08 1.8 6.2 11.9 34.2 45.9 100.0 There are six forest regions in China. These forest regions are divided according to the administrative authorities and forest types. A brief description of each follows. Chapter 5. Application to Chinese Forest Inventory 50 1. The northern coniferous forest region. Most of China's remaining forest is in the north. Located in the extreme north of Manchuria, this region includes the northern part of Daxinganling. Adminis- tratively, the area covers some 200,000 square kilometers, and embraces part of Heilungjian province and the Inner Mongolian autonomous region. The vegetation type in this region is strongly reminiscent of the northern coniferous forests of Europe and North America, with many common or closely related species. The region is almost entirely forested. The forest is mainly overmature with corresponding deficiencies in young-age classes. The distribution of these forests is irregular both in volume and area. 2. The mixed coniferous and deciduous broadleaved forest region. This type forms a transition between the coniferous forest in the north of Manchuria and the deciduous broadleaved forest in the south. It covers much of the Xiaoxdnanling and the Changba mountains. China's richest timber resources are located here. Vast forests remain which provide some 30 percent of China's present timber requirements. Like the northern coniferous forest, the forests in this area are mainly overmature with corresponding deficiencies in younger age classes. However, the distribution of these forests is more regular than that of the northern coniferous forest region. 3. The deciduous broadleaved forest region. This zone skirts the southwestern end of the Manchurian plain to take in the peninsulas of Shantung and Liaoning, the northern China plain, the northern slopes of the Tsingling mountains and the western half of the Shansi highlands. It covers almost 1,000,000 square km in nine provinces. Most of the accessible forest in the deciduous broadleaved region has been logged and the present potential for timber Chapter 5. Application to Chinese Forest Inventory 51 production is not high. This area has a long history of drought and flood. There is an obvious need for afforestation of protection forests in conjunction with a water conservancy scheme, and considerable scope for the development of farm forestry and shelter planting on the plains. At the present time, China has an ambitious reforestation program in the region; 4. The mixed deciduous and evergreen broadleaf forest region. This region is a transition between the deciduous and the evergreen broadleaved forest zone and forms a belt of 300,000 square kilometers along the Hanjing river in the west and the lower reaches of the Yangtze in the east. It runs through parts of eight provinces. The existing forest, as in most accessible parts of China, has been cut over a long period and most of the natural forest has disappeared. Only the protected areas around the secluded temples in the mountains and hills give an indication of the natural forest. This area is densely populated and forms part of the traditional rice bowl of China. 5. The evergreen broadleaved forest region. This forest region in China extends in a broad belt some 1,750,000 square km in area across the southern provinces and include parts of Zejing, Anhuai, Jiangsu, and Guangxi autonomous regions. Much of the natural forest of this area has been cleared for crop cultivation. Most of the evergreen species are utilizable timber trees and, where readily accessible, they have been extensively logged. Excessive tree cutting and poor methods of cultivation and reclaimation have made the forests in this region into lean secondary forest types. 6. The tropical monsoon rain forest region. The tropical monsoon rain forest is restricted to southern China and forms a narrow belt 250,000 square km in area along the sourthern edge of Fujian and Guangdong Chapter 5. Application to Chinese Forest Inventory 52 provinces. The tropical regions in general, and the rain forest in particular, have been a great reservoir of plant stocks. Although the tropical forests are different in many respects from those of the temperate forests, they have important similarities. While irreversibilities are more common in the tropics, in the vast majority of cases the tropical forest has the resiliency to restore itself. However, often the new forest will be different than the original in the composition of species, etc. This is particularly true when the forest areas harvested or destroyed are large. Existing natural forest in the tropical zone in China is limited to the hills of some provinces. In the southern part of this region, there is still a large amount of unused forest. However, it is difficult to access these forests. There was no planned inventory of China's forests for a long time. The first inventory was conducted between 1953 and 1961. This survey was actually conducted in different forest regions, using point sampling and visual estimation. No permanent sample plots were established and the data obtained at this period were only used to provide basic information such as the approximate, total forest areas. From 1973 to 1977, a more detailed forest inventory was conducted by provinces. An average volume of 78.26 m / h a (Ministry of Forestry, 1980) was found. In several forest 3 provinces such as Guanxi and Heilongjiang, a basic CFI system was established in this period. Some detailed checks and remeasurements have been made on a small scale in recent years. However, forest inventory is still a weak link in managing Chinese forests. Due to a lack of adequately trained personnel and equipment, the data obtained through the early forest inventories have seen very limited use and are not precise (Huei, 1980; Zhou, 1980). Many problems still exist. Lack of suitable inventory information has made forest management decisions very difficult to make. Chapter 5. Application to Chinese Forest Inventory 53 Forestry is playing an increasingly important role in recent years in China's social and economic development. Demands for timely and accurate forest inventory data are becoming more urgent. Forest inventory systems are now being established on a national basis to provide the required data. New and improved sampling methods are being tried in different forest types and regions. The following section will outline the steps necessary to apply the methodology developed in the previous chapters to China and discuss aspects of the design of a national forest inventory for China based on CFI principles. The purpose is to illustrate points which should be addressed, rather than providing a design which should be implemented. 5.2 Design of a Two-stage Sample 5.2.1 G e n e r a l Considerations It is very important to define clearly the various uses of CFI's in forests. In order to design and implement an inventory which best solves the problems, a clear idea, of the objectives must be decided before starting the work. In China, objectives center on current average volume, forest areas, stand structures, and growth patterns. These objectives should not be defined by an inventory specialist alone; they must be defined jointly by the people who will make use of inventory results and the inventory specialist. The inventory specialist should design an inventory which will provide the users with the information they need in a suitable form and with the required precision (Fontaine, 1981). Because CFI's can produce information on a number of variables besides volume, the priority of C F I objectives to be met has to be clearly assessed before designing an inventory. Additional requirements also have to be known before starting the inventory. Fontaine (1973, 1981), and Loetsch et al. (1973) gave a detailed discussion of these Chapter 5. Application to Chinese Forest Inventory aspects. 54 Some major points are summarized here to illustrate how CFI's should be applied. 1. Acquire general information. For example, determine the authority responsible for the inventory, available information and data on the area to be inventoried from past surveys, reports, maps, general descriptions of the forest, and conditions of terrain, transportation facilities, accessibility, and other necessary information. 2. Determine the exact purpose of the inventory. Detailed specifications of the objectives should be listed, for example, nature of the information required, exact limits and size of the area to be inventoried, and precision required. A clearly defined purpose provides a focus for the inventory foresters in designing and executing the details of the inventory. 3. Decide on which type of design to use. Many sampling designs can be used with CFI's. The selection of a sampling design is not limited only to the theoretical optimum. It is the combination of many theoretical and practical considerations. Some design possibilities include systematic or random locations, clustered or single locations, points or plots, one-stage or multistage. 4. Establish correct measurement procedures. For example, knowing the instruments that will be used in the inventory, provide detailed instructions on all techniques. Field test the instruments and the measurement procedures that will be used and provide forms for recording of observations. Outline a procedure for locating and establishing sample plots, etc. 5. Review compilation procedures and production of the final report. For example, provide detailed formulae for estimating means, totals, and their sampling errors, skeleton tables, and detailed descriptions of all phases of calculation from raw data Chapter 5. Application to Chinese Forest Inventory 55 on original forms to final results. It is also important to consider the analyses which will be made on the results in the final report. 5.2.2 T h e S h a p e a n d Size of S a m p l e Plots It appears that no single kind of sample plot shape is the best for Chinese forest inventories. The establishment of CFI sample plots should differ according to the different natural and economic situations in the various regions in China. It is recommended that in the eastern part of China (the deciduous broadleaved forest region and the mixed deciduous and evergreen broadleaf forest region), which largely consists of plains, the use of clusters of circular sample plots is better. Cluster sample plots are popular in northern Europe, especially in Sweden, where they are used in national forest inventories. This approach can profit from the advantages of circular plots while having at the same time sufficiently large sampling units. Hagberg (1957) discussed this method in detail. Janz (1975) made some improvements on Hagberg's method. Although cluster sample plots have many advantages, they also have some disadvantages, especially in mountainous regions. Also, for the same size of plot, the cluster of circular plots may have a longer total perimeter and consequently more borderline trees than equivalent rectangular sampling plots. However, these deficiencies will not influence the adoption of cluster sample plots in the eastern part of China which largely consists of plains. A cluster of 2 to 5 circular plots can be used. In southern China (the evergreen broadleaved forest region and the tropical monsoon rain forest region) which largely consist of mountains and rain forests, while the transportation systems are very limited, variable area single plots should be used. Variable area plots are mostly used in North America based on point sampling theory. In the southern provinces of China where the rain forests are common, these sample plots are used widely at the present time. It has been proven that this type of plot is very efficient. 56 Chapter 5. Application to Chinese Forest Inventory It can gain the required precision at lower cost because in each sample plot, fewer trees are measured. Plot centers can be relocated very easily when the sample plots become permanent plots. On the negative side, this kind of sample plot can cause some trees in the sample location to be missed unless caution is taken. Small mistakes can cause a big difference in the estimate of the current volume and especially in estimates of the change in volume. It should also be mentioned that the larger the sample plot, the more precise the results. However, the cost of the plots will also increase. It is therefore necessary to stipulate basal area factors (BAF's) for different regions. In Guangxi and Guangdong where the tropical forests are more intensively managed and need more accurete inventory data, smaller BAF's (more trees per plot) are desirable. In Yunnan, Hainandao, and part of Gueizhou where the forests are mainly unexplored and the accessability is not good, larger B A F ' s (less trees per plot) can better serve inventory purposes. In the northern part of China (the northern coniferous forest region and the mixed coniferous and deciduous broadleaved forest region) where forests are regulated well and the topography is simple, the use of rectangular sample plots would appear best. Rectangular sample plots are used extensively in most Chinese provinces. They usually cover an area of 0.06 to 1.0 ha. Rectangular plots have some advantages. They have clear boundaries and each tree has a smaller probability of being ignored or double measured than with other types of plots. The disadvantage of this kind of plot is that sometimes there are too many trees in a sample plot which may increase survey costs. Since the forests in China have a great diversity even in a very small area (especially in the southern and eastern parts), any small sample plot will not provide representative data. It is recommended that a relatively large plot size with a relatively small number of total sample plots should be used for Chinese forest inventories. A large plot size may increase the measurement cost,.but the relatively small number of total sample plots will Chapter 5. Application to Chinese Forest Inventory 57 greatly reduce the cost of transportation, a major cost in Chinese forest inventories. The plot sizes suggested for different regions are as follow. They should be jointly considered when used in designing sample plots. 1. In the evergreen broadleaved forest region and the tropical monsoon rain forest region, a large plot size of 1 ha to 1.5 ha should be used. Such a large plot may increase surveying cost by a large amount in many other countries. However, the low labor cost in China (approximately 2 to 3 Canadian dollars per day) will keep the costs reasonable. The large plot size will reduce the total number of sample plots in these regions. This will greatly reduce the cost of transportation and minimize the total inventory costs. 2. In the northern coniferous forest region and the mixed coniferous and deciduous broadleaved forest region, a sample size of 0.8 ha to 1 ha should meet the inventory purposes. Northern forests are well regulated and the topography is simple. Inventories are much easer than in the southern tropical regions. The transportation facilities in this area are also readily accessble. 3. In the deciduous broadleaved forest region and the mixed deciduous and evergreen broadleaved forest region, a sample size of 0.2 ha to 1 ha should be used. The better economic situation and transportation capacities in these areas make an inventory relatively easy to conduct. Past inventory results showed that a relatively small sample plot in this region not only met the inventory purposes, but also were very practical. 5.2.3 Sample Size Althrough the recommended plot size varied from region to region, an individual CFI sample plot of 1 ha will be assumed for sake of simplicity in subsequent discussions. Chapter 5. Application to Chinese Forest Inventory 58 As seen previously, the total forest area is approximately 115,277,400 ha. Suppose a sampling intensity of 0.02 percent is used on a national level (previous sampling intensities have ranged from 0.005 percent in the south to 0.1 percent in the north for large scale inventories). Then, on the national level, a total area of approximately 23,600 ha should be sampled, requiring 23,600 sample plots. The average cost for establishing and measuring a sample plot in China varies from region to region. In the tropical south where the forest conditions are very complex, a sample plot usually costs 500 Yuan ( 3 Yuan= 1 Canadian dollar). The average cost in the northern regions is 400 Yuan while in the eastern regions it is 350 Yuan. Suppose an average cost of 400 Yuan per plot is used for the national inventory system. For a total sample size of 23,600, a total cost of 9.44 million Yuan is needed for the first stage of the national inventory system. This cost will include: 1. Transportation and locating sample plots. 2. Measurement costs and compensation for the workers. 3. Auxiliary materials (maps, recording papers and tables, pens, etc.) and equipment (raincoats, bags, etc.). The cost for measurement equipment is not included. Of the 23,600 sample plots, a certain portion will remain as permanent sample plots and will be remeasured on a second occasion. This portion is not only decided by the inventory results from the first occasion, but also by the total budget available at the second occasion. For the sampling design presented in Chapter 3 and 4, the number of sample plots that should be remeasured on the second occasion is also decided by the individual values of each sample plot and the total values of all sample plots, as well as additional requirements such as the required precision, expected population variance, coefficient of variation, correlation coefficient, etc. Some rough data for these requirements 59 Chapter 5. Application to Chinese Forest Inventory in China are provided by Guang et al. (1983) from past inventory results. The coefficient of variation usually ranges from 20 percent to 50 percent, the correlation coefficient ranges from 0.6 to 0.9, and the required precision for a large scale inventory is 5 to 10 percent. This aspect will not be discussed in detail and only the basic sample size calculation formula used by many researchers in a simple random design will be given as a rough guide for calculating the portion that should be left as permanent sample plots: where n is the number of permanent sample plots, n is the total number of sample plots, p a p is the correlation coefficient between the sample plots measured on the first occasion and remeasured on the second occasion, C is the cost for establishing and measuring p permanent sample plot, and C is the cost for establishing and measuring a temporary t sample plot. For example, suppose the cost of relocating and remeasuring the established sample plots is 450 Yuan. The average cost of establishing and measuring the new plots on the second occasion is the same as on the first occasion, 400 Yuan. A correlation coefficient between the sample plots measured on the first occasion and remeasured on the second occasion is approximately 0.8. Hence, C = 400, C = 450, p = 0.8, and n = 23, 600. t p a Substituting the values into the formula yields n =7585, which means that within the p 23,600 established sample plots, a total number of 7,585 sample plots should remain as permanent sample plots and be remeasured on the second occasion. It is therefore recommended that for a national CFI system in China, a total number of 8,000 to 10,000 permanent sample plots should meet all the purposes of the inventories. In order to retain a total number of 23,600 sample plots on the second occasion, 13,600 to 15,600 temporary sample plots would need to be established. 60 Chapter 5. Application to Chinese Forest Inventory It must be emphasized that the above results can only serve as a rough guide for a national CFI system in China. Determining the "best" sample size involves optimizing the sample design. This means finding the design that for a specific cost gives the smallest error for the parameter to be estimated or for a specified error yields the least cost. This is very complicated for a two-stage design using unequal probabilities with S P R on more than one occasion, since good information on the sampling and measurement errors, as well as biases is never simple to obtain. In many textbooks and manuals the calculation of the optimum is restricted to the estimation of the optimum number of sample units at each stage based on a simple cost formulation. Loetsch and Haller (1964) gave an example of optimization of a two-stage sample based on equal probability for the one occasion case using inventory data collected in Thailand. Other researchers have used different cost functions in calculating the optimum sample size in a two-stage sample design. These include: 1. C = C\n + C nm. 2 where C is the total cost for a inventory, Cjii is the total cost for the number of primar3 units in the sample, and C2nm is the total cost for the total number of r second-stage units. This cost function is discussed in detail by Cochran (1977) and is often used when the travel costs between units are unimportant. 2. C = C-^n + C Jn + C nm. ty 2 This cost function is presented by Hansen et al. (1953) and used when the cost of travel (C ) between primary units is substantial. t 3 where C w is the cost of walking a unit distance, C is the cost of establishing a plot or point, C e m is the cost of measuring the characteristic of interest, yJnA is the average minimum distance among randomly drawn plots in the area A , and 61 Chapter 5. Application to Chinese Forest Inventory T is the average number of elements observed and measured at a location. This function was presented by O'Regan and Arvanitis (1966). A general discussion concerning the optimum sample size was given by Fontaine (1973). He assumed that the standard error of the estimate (SE) of a given sample, and the total cost (C) are functions of various design characteristics (x , x , ... , Xk). x 2 For a given cost of C , to minimize the standard error SE, it must follow: 0 SE{xi, x ,Xk) = 2 C(x ,x , 1 2 ...,x ) k minimum = C 0 Optimum sample size is obtained by taking the derivative with regards to every characteristic involved. 5.2.4 T h e Layout and Survey of Sample Plots The layout of a two-stage sample design should be carefullj r planned. There are man3 layout forms that could be used in a Chinese national C F I system. In order to r provide more precise estimates, it seems reasonable that different forest regions should be treated as distinct populations and different layouts should be employed in each region. However, the forests and topography in China are so complicated that this approach is seldom possible in practice. Given the present situation in China, systematic layout of plots has the most promise for a national C F I system. The vast majority of previous inventories have used some form of systematic selection. Systematic sampling provides good coverage and good representation. It is also good for mapping. The following steps are necessary in designing the layout of the sample plots: 1. Divide the whole nation into uniform portions. The total forest area in China is 115,277,400 ha. This is comprised of forest production farms of approximately Chapter 5. Application to Chinese Forest Inventory 62 50,000 ha in size. Assuming a uniform size of by 50,000 ha means that there will be 115,277,400/50,000=2306 farm size regions in the country. 2. From the 2,306 production farms, systematically select 200 from which to be sampled further. 3. Classify each of the selected farms according to forest types and ages. Then divide the farm into different units according to the classification. These units will be used as primary sampling units and they usually cover an area of 20 ha to 1,000 ha. 4. Choose the primary units according to probability proportional to the area of the units. Fifteen primary units should be selected in each farm using this method. 5. Select eight plots (secondary units) w ithin each primary units by probability proT portional to predicted volume. The above procedures will ensure: 1. A total sample size of 8x15x200=24,000 will be selected. This number is approximately equal to the required total number of sample plots which is 23,600. 2. A two-stage unequal probability sampling estimate procedure presented in Chapter 4 can be used in each of the selected 200 farms. 3. A final national inventory result can be obtained by using the results from the 200 production farms. Chapter 6 Conclusions The previous chapters discussed the general theory of multistage sampling combined with CFI using unequal probabilities in a S P R structure on two occasions. Only one-stage and two-stage S P R with unequal probabilities were discussed. Estimators of the current means and the change in means for both one-stage and two-stage cases were given. Although no precise theoretical demonstration of the gains from using the presented method was given, the combination of highly efficient sampling with unequal probabilities and the very practical multistage sampling in CFI with SPR in Chinese forest inventories can certainly provide an efficient alternative to traditional systems. General theories of CFI with SPR, unequal probability sampling, and multistage sampling have been developed by manj' researchers. Their applications in forestry have been examined in depth. CFI with S P R estimation structure has existed for a long time and the main difference among the various CFI methods is the way to obtain the least variance or to obtain the least combination of the variances from each component. In this study, multistage unequal probability sampling was the method used to obtain these variances. The results in Chapter 3 showed that for one-stage SPR with unequal probabilities, the estimator of the current mean takes the same form as presented by many researchers. However, after using unequal probability sampling, the variances obtained and the constants involved in estimating the current means differ from those reported for equal probability sampling. The estimator of the current mean in .two-stage SPR with unequal 63 64 Chapter 6. Conclusions probabilities also showed a similar form as the one-stage case, with the variances and constants differing from the one-stage case. This resulted from a similar linear form and estimation structure. It is difficult to tell whether the estimators derived in this study are better or worse than the estimators obtained using other methods. Further analysis and comparison might be done in this area. In estimating the change in means, the conventional linear estimate form is used again, and the results from one-stage SPR with unequal probabilities showed that the best estimator of the changes on both occasions is the difference between the best estimate of the current mean on second occasion and the overall mean on the first occasion. A similar result is also obtained for two-stage SPR with unequal probabilities. The overall variances in estimating the change in means for both one-stage and twostage SPR with unequal probabilities on two occasions are obtained by substituting the different variances from independently observed sample plots (temporary sample plots) and from dependently observed sample plots (permanent sample plots). Results showed that in the one-stage case, the best estimate of the total variance is a weighted estimate of the different variances. For two-stage sampling, a similar result is obtained with the variances and constants differing from the one-stage case. The estimators of overall variances developed take simpler forms than those presented in the literature which might result in greater simplicity when used in data processing. The application of C F I with S P R using multistage unequal probability sampling in Chinese national forest inventories is discussed for illustrative purposes. Inventory systems are a weak link in Chinese forest resource management and the nation is now undergoing an ambitious program trying to establish a good inventory system. Some researchers have suggested that large scale remote sensing images combined with ground sampling design is a better way to achieve accurate results. Others have suggested that current social and economic situations might not allow the wide use of remote sensing Chapter 6. Conclusions 65 images in forestry, and thus, more emphasis should put on developing and using more efficient and practical sampling designs. It will be benefical if a CFI with S P R using multistage unequal probability sampling could be combined with the use of available remote sensing results, especially in forest region classifications, forest type identifications, and area calculations. The establishment of a national C F I system in China should use all kinds of useful auxiliary facilities and information that might reduce the complexity of the inventories. China's inventory problems can never be solved by a uniform sampling design. Some suggestions made in Chapter 5 may be used as a rough guide in actually designing the inventories. More research is needed regarding how the design can be implemented more efficiently and how the optimum sample size could be calculated. With the improving understanding of current inventory methods and with the practice of new inventory designs, a better inventory system will eventually be established in China, and forest resource management decision making will benefit from improved information on the current status of the resource. Chapter 7 Literature Cited Beers, T. W . 1978. Developing efficient estimation techniques for integrated inventories. USDA For. Serv. Gen. Tech. Rep. RM-55. pp. 270-275. Bickford, C. A . 1956. Proposed design for continuous inventory: a system of perpetual forest survey for Northeast. USDA For. Sev. Eastern Technical Meeting, Forest Survey. Cumberland Falls, K Y . 37pp. Bickford, C. A . 1959. A test of continuous inventory for national forest management based upon aerial photographs, double sampling and remeasured plots. SAF proceedings, 1959:143-148. Bickford, C. A . 1963. On successive forest inventories. Proc. Soc. Amer. For., 25-30pp. Bickford, C. A . , C. E. Mayer, and K . D. Ware. 1963. An efficient sampling design for forest inventory: the Northwest forest survey. J. For. 61:826-833. Brewer, K . R. W., and M . Hanif. 1983. Sampling with unequal probabilities. SpringerVerlag. New York. 164pp. Chakrabarty, R. P., and D. S. Rana. 1974. Multi-stage sampling with partial replacement of the sample on successive occasions. Amer. Stat. Assoc., Proc. Social Stat. Sec.:289-291. Cochran, W. G. 1953. Sampling techniques.. John Wiley and Sons, Inc., New York. 330pp. Cochran, W. G. 1977. Sampling techniques. 3rd Edition. John Wiley and Sons, Inc., NewYork. 428pp. Cunia, T. 1964. What is sampling with paxtial replacement and why use it in continuous 66 Chapter 7. Literature Cited 67 forest inventory. Proc. Soc. Amer. For. 207-21 lpp. Cunia, T. 1965. Continuous forest inventory, partial replacement of samples and multiple regression. For. Sci. 11:480-502. Cunia, T. 1987a. On the error of continuous forest inventory estimates. Can. J . For. Res. 17:436-441. Cunia, T. 1987b. On the error of estimates from successive forest inventories. Can. J . For. Res. 17:446-447. Cunia, T. and R. B. Chevrou 1969. Sampling with partial replacement on three or more occasions. For. Sci. 15:204-224. Cutter, D. D. 1955. A permanent plot system of survey for the continuous inventory on ponderosa pine stands in the southwest. J . For. 53: 186-189. Deming, W . E . 1950. Some theory of sampling. John Wiley and Sons, Inc. New York. 602pp. Dixon, B . L., and R. E . Howitt. 1979. Continuous forest inventory using a linear filter. For. Sci. 25:675-689. Fontaine, R. G . 1973. Manual of forest inventory with special reference to mixed tropical forests. FAO manual, Rome. 200pp. Fontaine, R. G . 1981. Manual of forest inventory with special reference to mixed tropical forests. FAO Forestry paper, Rome. 200pp. Frayer, W. E . 1966. Weighted regression in successive forest inventories. For. Sci. 12:464472. Frayer, W . E . 1967. A systematic bias in the interpretation of CFI results. USDA For. Serv. Res. Note NE-60. 4pp. Frayer, W . E . , and G. M . Furnival. 1967. Area change estimates from sampling with partial replacement. For. Sci. 13:72-77. Frayer, W . E., R. C. Van Aken and R. D. Sullivan. 1971. Application of sampling with Chapter 7. Literature Cited 68 partial replacement to timber inventories, central Rocky Mountains. For. Sci. 17: 160-162. Freese, F. 1962. Elementary forest sampling. USDA For. Serv. Southern For. Exp. Stat. 91pp. Furnival, G. M . , T. G . Gregoire, L. R. Grosenbaugh. 1987. Adjusted inclusion probabilities with 3P sampling. For. Sci. Vol. 33, No. 3. pp:617-631. Gillis, M . D. 1988. Estimating change from successive static forest inventories. For. Chron. Vol. 64:352-354. Grosenbaugh, L. R. 1964. STX-Fortran-4 program for estimates of tree populations from 3P sample-tree measurements. USDA Fer. Serv. Res. Pap. PSW-B, Berkeley, C A 49pp. —— 1965. Three-P sampling theory and program T H R P for computer generation of selection criteria. USDA For. Serv. Res. Pap. PSW-21. 1971. STX1-11-71 for dendrometry or multistage 3P samples. USDA For. Serv. Pub. FS-277. 63p 1976. Approximate sampling variance of adjusted 3P estimates. For. Sci. 22:172-176. Guang, Y . X . et al.. 1983. Forest mensuration. Vol. 2. 191-345pp. Hagberg, E . 1957. The Swedish national forest survey. Unasylva. 11:3-8. Hall, 0 . F. 1959. The contribution of remeasured sample plots to the precision of growth estimates. J . of For. 57:807-811. Hansen, M . H., and W . N . Hurwitz. 1943. On the theory of sampling from finite population. Ann. Math. Stat., 14:333-362. Hansen, M . H . , W. N . Hurwitz, and W. G. Madow. 1953. Sampling survey methods and theory. Volume LMethods and applications. Volume ILTheory. John Wiley and Sons., Inc., 638pp. and 332pp. Chapter 7. Literature Cited 69 Hazard, J . W . and L. C. Promnitz 1974. Design of successive forest inventories: optimization by convex mathematical programing. For. Sci. 20:117-127. Hendricks, W . A . 1956. The mathematical theory of sampling. Scarecrow Press. 364pp. Hsiung, W . Y . 1980. Forests and forestry in China. The H . R. MacMillan lectureship in forestry. Faculty of Forestry, The University of British Columbia. Vancouver, B.C. 18pp. Huei, Z. 1980. Analysis of the data obtained through the national forest inventories in the past years in China. J . of Forest Management and Surveying (Chinese). 1980(2):32-38. Janz, K . 1975. The Swedish national forest survey. Second FAO/SIDA training course on forest inventory. Jessen, R. J . 1942. Statistical investigation of a sample survey for obtaining farm facts. Iowa Agric. Exp. Stn. Res. Bull 304, 104p. Kish, L.T965. Survey sampling. John Wiley and Sons., Inc. New York. 643pp. Letourneau, L . 1966. A forest inventory by continuous forest control. For. Chron. 42: 414-419. Letourneau, L . 1979. Updating continuous forest inventory by simulation. Forest Resources Inventories: Workshop Proceedings. Colorado State Univ., Fort Collins. Vol. I. pp:327-332. Loetsch, F. and K . E . Haller. 1964. Forest inventory. Vol. I. Munchen. B E R N WEEN, 436pp. Loetsch, F., F. Zohrer, and K . E. Haller. 1973. Forest inventory. Vol. II. Munchen. B E R N W I E N . 469pp. Lund, H. G. 1979. Linking inventories. Forest Resources Inventories: Workshop Proceedings. Colorado State Univ., Fort Collins. Vol. I. pp:2-9. McClure, J . P. 1979. Multi-resource inventories-a new concept for forest survey in the Chapter 7. Literature Cited 70 southeast. Forest Resources Inventories: Workshop Proceedings. Colorado State Univ., Fort Collins. Vol. I. pp:23-28. Ministry of Forestry (China). 1980. China's forestry and its role in social development. Forestry Press. 43pp. Ministry of Forestry (China). 1984. A brief account of China's forestry. Forestry Press. 20pp. Newton, C. W . , T. Cunia, C. A . Bickford. 1974. Multivariate estimators for sampling v/ith partial replacement on two occasions. For. Sci. 20:106-116. Omule, S. A. Y . 1981. Successive forest sampling using multistage with partial replacement of unit. Ph.D. thesis. Faculty of For., The Univ. of B. C. Vancouver, B . C. 117pp. Omule, S. A. Y . 1984. Multistage sampling with partial replacement. Can. J . For. Res. 14:869-873. Omule, S. A . Y . , and A . Kozak. 1982. Estimators for successive forest sampling with partial replacement. Can. J . For. Res. 12:753-760. Omule, S. A. Y . , and D. H . Williams. 1982. Optimun allocation by dynamic programming for sampling on successive occasions with partial replacement of units. Can. J . For. Res. 12:264-269. O'Regan, W . G., and L . G . Arvanitis. 1966. Cost effectiveness in forest sampling. For. Sci. 12:406-414. Patterson, H . D . 1950. Sampling on successive occasions with partial replacement of units. J . Royal Stat. Soc, Series B , 12:241-255. Peng, S., and M . Zhu. 1985. Kalman's linear filtering and its application to the estimation of forest resource dynamics. Scientia Silvae Sinicae. Vol. 21, No. 2:122-131. Pleines, W . and L . Letourneau. 1970. The continuous forest control inventory system. For. Chron. 46: 39-43. Chapter 7. Literature Cited 71 Rana, D. S. 1978. Ratio method of estimation in multi-stage successive sampling on two occasions. Amer. Stat. Assoc., Proc. Social Stat. Sec.:289-291. Schreuder, H . T., J . Sedransk, and K . D. Ware. 1968. 3-p sampling and some alternatives, I. For. Sci. 14:429-453. Schreuder, H . T., J . Sedransk, K . D. Ware, D. A . Hamilton. 1971. 3-p sampling and some alternatives, II. For. Sci. 17:103-118. Schreuder, H . T. 1984. Alternative estimators for point-Poisson sampling. For. Sci. 30:803812. Scott, C. T. 1977. Unequal probability sampling for updating inventory estimates. M.S. thesis. Univ. of Ga., Athens. 62pp. Scott, C. T. 1979. Midcycle updating: some practical suggestions. Forest Resources Inventories: Workshop Proceedings. Colorado State Univ., Fort Collins. Vol. I. pp:327332. Scott, C. T. 1984. A new look at sampling with partial replacement. For. Sci. 30:157-166. See, T. 1974. Forest sampling on two occasions with partical replacement of sample units. M.Sc. thesis. Faculty of Forestry, The Univ. of B.C., Vancouver, B . C . 63pp. Singh, D. 1968. Estimates in successive sampling using a multi-stage design. American Stat. Assoc. J . 63(321):99-112. Solomon, D . S. 1979. Permanent plots in forest research. Forest Resources Inventories:Workshop Proceedings. Colorado State Univ. Fort Collins. Vol. I. pp:327-332. Spurr, S. H. 1952. Forest inventory. The Ronald Press Co., New York. 476pp. Stott, C. B. 1947. Permanent growth and mortality plots in half the time. J . For. 37:669673. Stott, C. B. 1968. A short history of continuous forest inventory east of the Mississippi. J. For. 66:834-837. Stott, C. B., and G. Semmens. 1962. Our changing inventory methods and the C F I system Chapter 7. Literature Cited 72 in North America. Proc. 5th. World For. Congr., Seattle, 1:451-454. Titus, S. J . 1979. Multi-stage sampling: what's it all about? Univ. of Alberta, Agri-For. Bui. 2(3):16-18. Van Hooser, D . D . 1972. Evaluation of two-stage 3P sampling for forest surveys. USDA For. Serv. Res. Pap. No. SO-77. Van Hooser, D. D. 1973. Field evaluation of two-stage 3P sampling. USDA for. Serv. Res. Pap. SO-86. 15pp. Ware, K . D. and T. Cunia 1962. Continuous forest inventory with partial replacement of samples. For. Sci. Monog 3. 40p. Yandle, D. 0. and F. M . White 1977. A n application of two-stage forest sampling. South J. Appl. For. (3):27-32. Yates, F. 1960. Sampling methods for censuses and surveys. Griffin and Co., London. 440pp. Zhou, C. Q. 1980. Some comments about current forest management situations in China. J . of Forest Management and Surveying (Chinese). 1980(2):28-31. Appendix A Derivation of E q u a t i o n 3.6 The variance of y is b V(y ) = b a [ y ( x J + y ( x ) ] + c?y(y ) + ( l - ) 1 / ( y ) 2 2 f p p C l t 2 -2aic Cov(x ,y ) 1 = p (A.22) p + f l ) + > f i + (1 n n c n tl p C f i l - 2 a n l p l C l t2 ^ ^ n (A.23) l p and, „ . = n sis nin p = C l 2 2 -. - n n rj l 2 . 2 4 ) 2 — n sis nin p ( A - n^n^rj z (A.25 / 2 tl t2 where rj = 8pa a . ei £2 The relationship of a\ and c\ can be written as: Vu n ai = c sinin i • p Substituting the calculated values of (A.24) and (A.25) into (A.23) yields the general formula for calculating the minimized variance of V (y ): b „ > V{Vb) r\s n\n , , . 1 1 . ' 2 + — n s s n n2 — n n rj n n . t = 2 — ti — i p z 1 2 1 tl 7i s s "-in 2 p 1 n 2 - tj n n r] tl p . . . _ T>. n. •-ie. *s csnn- T n - n n r)' 1 n ssnn 2 2 1 2 1 2 1 2 tl t2 73 n t2 - n n r] tl p n 2 t2 , p p 2 t? 2 h n ssnn 2 p l 2 1 - n n rj 2 2 u t2 n 2 p 74 Appendix A. Derivation of Equation 3.6 Therefore, y s ny vr- ^ V{y ) = 2 2 1 tiSl -j— {n s s n n -n n T} ) s s\n ni b ^ 2 2 2 p l 2 ,N — + —) n n 2 l 2 tl 2 t2 2 1 ti p 2 _| 12 1 p (n SiS Tiin2 - n n r/ ) / aig2wm; 2 2 p + t l 2 2 t 2 2 x n 2 p s s ra n 1 2 1 n n r] 2 2 tl h s 2 n t2 2r] sln n n s 2 i p (n s s n n p 1 2 1 tl 1 1 - 2 2 n n v)' 2 tl 2 t2 So, 7/ 5 7Z 7 l 5 l T// - N 2 (n s s 7i n - n n 7)*y 2 p 2 1 (n 2 p 1 2 + (1- 1 - 2 n n r/ ) 2 f l (n 2 p 2 p p 5 i5 7i n 2 2 s s n n 1 2 1 1 - 2 - 2 2 + (1- 1 t l 2 S l S 2 7 l i n n s 5 ri 7i 1 p 1 2 tl 7/ ) 2 r — [7— 2 > 2 " ) - n n n n z 2 tl t2 h 1 I 2 s 5 n n 1 2 2 t 2 S lil 2 t2 t2 2 n n S p n t2 2 2 s 2 n n T}* tl 2 ri p n n n) s s 7i 7i -n ra 1 n tl j ) ' ^ , ! T) slnlrit^nx (n 1 f 2 2 n n t2 tl sis n ri 2 1 tl 2 1 - n n 7]' 2 h N 2 n 1 t2 t2 •s s n n^ 1 2 1 2 2 r -(ni«is n n 2 2 t Notice, •SiS27ian 7i 2 1- c = a ra s s n n 2 p 1 2 1 2 t2 - n n T} - n n T} tl 2 t2 2 tl t2 - 7/ n n t l t 2 ) + (1 - c j 75 Appendix A. Derivation of Equation 3.6 Thus, V(yb) = — —i ~ = i l ^ - d X d +l - d ) j ( l - c i ) + (1-ci) Appendix B Derivation of E q u a t i o n 3.11 The general formula for calculating the minmized variance of V(Ab): - n 5 ( 5 7 i 7 i p + n r]) 2 2 1 tl 2 S l s^n^nl - n n rj n ^ nls (s n n + n rj) s 2 tl ! l 2 1 p p s^sin^nl s 2 t2 2 - n n j] n 2 tl t2 [ i _ lA 2 i P n s + t2 n n n + p , t v) 2 s-is^^nl - n n T) —n s (s\n n + n Tj) s 2 2 n 2 tl p + [1 + 2 p tl 5 5 n n 7 i 2 - n n r] ^-%s (sin n + n r]) 2 2 ssnnn 1 2 1 1 2 tl 2 p tl u t2 + n v) t2 v s-^s^n^n - n n rj 2 2 1 n ,2 , n^As^rip t2 - n n rj 2 t2 s 2 2 1 ( 2 t2 2 t2 n 2 tl t2 Notice this variance can be split into two parts: V(A ) = V + V . b A B where , nlsi(s nin 2 p n Tj) + 5 s n ra n| - n n v 1 + 2 1 2 t2 t l 2 s, 2 n 2 t 2 p ^ _ nlsAs^rip + n ri) ^ i2 s s n n n 2 - n n rj 1 2 1 ^-n s (s n n tl 2 2 p 2 A S-LS^I^TI 2 tl 1 t2 and 1 2 n 2 2 2 tl t2 76 t2 ssnnn 2 5i5 7i n n - n n v 2 n rj) 71^1(^2^1% + tl - n n rj 2 t2 + n 77) p s 2 n 2 p 2 1 2 2 t2 - n n T} tl T2 2 7 '? n 2 p Appendix B. Derivation of Equation 3.11 j + -s nl(s n n 1+ 2 1 2 + n T]) p tl - n n rf 5i5 nin n^ 2 ]2 S l 2 tl 2 +[ 5i5 7i Ti n 2 a 2 p n t2 -n s (sin n + n r)) - n n rj 2 7i^i(a Tt TO 2 tl 2 s^s^^n 2 2 tl tl t2 1 + rt^T/) p - n n Tj 2 2 t2 tl T/ n 2 v Hence T . 5!S 7l 1 2 2 ^ = _ F>2 iy2(7lt Wp5iS Tli + 2 - n n^77 )l -5 n n 7i r/ 1 A = where Q 2 'l'52 p 1 1 Tl p r , Q2("t Tt 2 2 $2 = •5 7iin 2 1 2 2 t2 tl - c'i) ~ 2 2\1 i f i I \2 72.^71^7? ) + (1 - C j ) 2 5^271! - T i 5 i 7 l 7/ t 2 and n n T) , p — SiUpT) t 2 (1 + tl 2 s 5 n n 7i = 2 t2 U 2 + Ti ,7/. t p Note that: , _ n nls s n 1 t2 1 2 - n^s^y 1 n ny - 2 t2 tl So the variance V 4 can be written as: V™t ™t 2 2 = ^-(l-cj^-.+a-cj] —(1-4). = Ti t j Similarly, the variance VB can be written as: T , -sassn 1 2 ' = —[yiKV * ™ + B 1 n 2 = where Q•= x ^ 7 ^ 5 5 n 7i 7i 1 p 2 1 2 2 2 n t i S 2 n p v - n n ,y )] + (1 + & i ) — -s n n n ri 2 2 2 2 tl ti t —[<?iK n 5 5 7i 2 1 — nny, 2 tl u 1 2 and 2 - n s n r] 2 tl 2 n n T} )} 2 tl t2 + (1 + fc'J 2 Appendix B. Derivation of Equation 3.11 Qi = Sin n 2 p + n T]. tl Note that: _ n n ss n 2 tl p l 2 - n 5 n^T/ - 2 tl 2 n n ri tl 2 t2 Q So the variance VB can be written as: v B = ^ [ < ? i ( i Qn _ i l 7l - _ i l ! ^ Tl t l _ & Til 2 5 * c i 'HR\-^liifi Tij tl _£I!^L _ ^ Tl Tii Tii i 5 i C Tlx i = -(l-M). Tii Hence, the variance V(Ab) can be expressed as: ^0 = ^(1-0 Tl t 2 tl Til i l _ (i _ t l K ) n tl f l 7l - + K)] + ( i + + ^(1-^). Tij Appendix C Derivation of E q u a t i o n 4.19 From the general formula for calculating the minmized variance: -n (s' + 4 X ( 4 + s'{)n + n MM + s'l) (4 + 4 X 4 + 4 ) M ™ 2 - n n f3l n ( p 2 2 t r tl | p t2 rx (4 + 4)((4 + 4 ) ^ + n M p (4 + s'{)(s'2 + 4 ' ) i 2 n t {s' + s'l) 2 2 n n Pl ra tl n t2 p M 4 + 4X(4 + 4 > i + ™t&) (4 + 4) + [1- 2 l2 (4 + 4 X 4 + 4 ) i 2 - n n f3 n -n (4 + 4 X ( 4 + 4 > 2 + ritMM + 4) + [1 + (4 + 4 X 4 + 4 ' ) i 2 - ntxnt2fil ntl ,2 r - n p ( 4 + 4 X ( 4 + 4 > 2 w /3 ) np(4 + 4)((4 + s'^n, 2 i[ ^ (-si + 4 X 4 + 4 ) i 2 - ntlnl2(3 (s[ + s'j(s2 -f 4 ' ) i 2 7 l n 2 tl tt t2 p n n + n tl 2 n n n + nM p n n Mn t 2 tl t Notice this variance can be split into two parts: V{A ) = V + y b 1 (C.30) 2 where ^ ( 4 + 4) ((4 + 4 > i + " ,j9 ) 2 V 1 = t (4 + 4 ) 2 2 - n n 0i) n +[i _ V 4 + 4 X ( 4 + 4'K + n M M + 4') (•si + 4 X 4 + 4') M 2 - n n Pl n -n (4 + 4 X 4 + 4)&((4 + 4 > + n M(i4 + 4 > i + n,/3) 2 ((4 + 4 X 4 + 4 ' K ™ 2 tl t2 p t r n tl p t2 t2 2 t t ((4 + 4')(4+4>m2-nt nt /?I) l and n (4 + 4) ((4 + 4 ' K + %,/3 ) (,i + s'l) ((4 + 4 X 4 + 4 > i ^ . -n n f3 y 2 2 2 2 2 tl l2 79 a 2 2 p Appendix C. Derivation of Equation 4.19 80 , , ~ K + * 2 > P ( ( * I + *1>2 + ntlP2),M + -si) [2 ( i + *i)(*2 + a ' ) i 2 5 n tl tj tl + *"){*2 + ^ M ^ ' i + a?)ra + n /3 )((s' + s'Dm + n 3 ) -rh^'i { n ra ra /3f n 2 2 {{s[ + s'{)(s + s'l) n 2 ni tl 2 nt^Pt) - 2 2 taf 2 2 Hence, \r Vi ( i + i)( 2 + "SD ? 1 ri u i i " W ' i jr + i)(( 2+ s = s QsK^i 71 s 2 V , / i "-42 „ / \2 fl re 2 1 2 2 tl 2 ( 2 + 2) 5 + 1 ~ 2> 1 ii\ \ a \f i , "\ a a'i 2) i +«-t P2j(-s + s J n / 3 - ri /3 s f l C 5 n t2 {s[+s'{)(s' + s'i)n 1 2 p [QaKK + s'l){s' 2 <? ra , 2 + s'^n, - n {s\ + s'{)n (3 - n n (S )} 2 ti p tl 2 t2 2 t ( l _ ' ) 2 ( 2 + 4) t 5 + c i n 2 where Q = (s[ + s")(s' + 5 ')n. n — n n fil, 2 g 3 2 1 2 ti and t2 = (a + s > + n /3 . 2 2 1 <2 2 Note that: , _ ra (^ + t2 ° 2 Q ( 4 + 4')n x ~ : 71^(^1 + 4)n /3 - p 2 - n n {3 2 tl Q — t2 ~ ~ - So the variance Vi can be written as: • _ (s[+s';)(s' + s'^n »'i = 2 j , M + *Z) Qs(l - c2)J + (1 - c2) p n V™<2 - —— Tl (!- 2) 7j c f 2 r- (1 — c )J 2 = W±i2(i-4). Similarly, the variance V can be written as: 2 T/ V 2 = (^l + )(' 2 + 2) P — S S n -(s> + s';) n A 2 ni t 4[ tl 1 2 2 tl t2 2 tl 2 ^ ^ 1 - n n /3 )] + (1 + nt, [<?4(ra (si + <? + ( 1 4 - O ^ n Unr, ( J A ' n"\( J A o"V~ 1 „ 2/ / , ll\Q U (n (s + s )(s + s )n + n (s + s )/3 •' 2 2 2 u 2 2 2 ^iflliS «i)( 2 + • 2 ) 2 - ra (s' + s )n 0 s s n tl 2 2 p 2 ra n /?|)] tl t2 Appendix C. Derivation of Equation 4.19 where Q - (s[ + 4 X 4 + s )n-in n - n n f3 , and 2 2 Q 4 81 2 tl t2 2 = (4 + 4 > + n / 3 . 2 tl 2 Note that: 1 + 6; = n (4 + 4 X 4 + 4 ) w tl 2 - n (fl t l 2 + s'l)n p p 2 - n n P ti 2 t2 So the variance V can be written as: 2 (4 + 4)(4 + 4 > p [g Q(l+6 )] 4 li+ I) i (1 + t 4+4')/ , (4±i|K Q )| n n 4+4') 2 P : n tiA-4 (4 + 4 K (l + 6' ) + 1 + 1+ i 2 , 2 ( s ' 1 + S ' ' 1 ) ] ) #.(4 + 4 ) 4 4 + 4) _ _ np _ # ( 4 + 4 ) 4 ( 1 4 +4 W A-(4 + 4 ) 4 'i+<) (i-#4). Hence, the variance V(A{,) can be expressed as: ™t 2 (C.31)
- Library Home /
- Search Collections /
- Open Collections /
- Browse Collections /
- UBC Theses and Dissertations /
- Continuous forest inventory using multistage unequal...
Open Collections
UBC Theses and Dissertations
Featured Collection
UBC Theses and Dissertations
Continuous forest inventory using multistage unequal probability sampling with partial replacement Huang, Shongming 1988
pdf
Page Metadata
Item Metadata
Title | Continuous forest inventory using multistage unequal probability sampling with partial replacement |
Creator |
Huang, Shongming |
Publisher | University of British Columbia |
Date Issued | 1988 |
Description | Continuous forest inventory (CFI) with partial replacement of sampling units (SPR) has been established as an efficient inventory technique. Previous CFI with SPR theory was based on equal probability sampling which might not be as efficient as the unequal probability sampling that is now widely used in forest inventories. In this study, the general theory of multistage sampling combined with CFI using unequal probabilities in a SPR structure on two occasions is discussed. Estimators of the current means and the change in means for both one-stage and two-stage cases are given, along with the estimators for the overall variances. The estimators derived for estimating the current means take the same form as those developed by other researchers. However, the estimators proposed for predicting the change in means take simpler forms than those presented in the literature. This might result in greater simplicity when used in data processing. The application of CFI with SPR using multistage unequal probability sampling to a Chinese national forest inventory is presented for illustrative purposes. Some particular points are addressed regarding specific situations in China. Although no precise theoretical demonstration of the gains from using this method is given, the combination of highly efficient sampling with unequal probabilities and the very practical multistage sampling in CFI with SPR inventories can certainly provide an efficient alternative to traditional Chinese inventory systems. |
Genre |
Thesis/Dissertation |
Type |
Text |
Language | eng |
Date Available | 2010-08-13 |
Provider | Vancouver : University of British Columbia Library |
Rights | For non-commercial purposes only, such as research, private study and education. Additional conditions apply, see Terms of Use https://open.library.ubc.ca/terms_of_use. |
DOI | 10.14288/1.0075388 |
URI | http://hdl.handle.net/2429/27348 |
Degree |
Master of Forestry - MF |
Program |
Forestry |
Affiliation |
Forestry, Faculty of |
Degree Grantor | University of British Columbia |
Campus |
UBCV |
Scholarly Level | Graduate |
Aggregated Source Repository | DSpace |
Download
- Media
- 831-UBC_1989_A6 H82.pdf [ 4.58MB ]
- Metadata
- JSON: 831-1.0075388.json
- JSON-LD: 831-1.0075388-ld.json
- RDF/XML (Pretty): 831-1.0075388-rdf.xml
- RDF/JSON: 831-1.0075388-rdf.json
- Turtle: 831-1.0075388-turtle.txt
- N-Triples: 831-1.0075388-rdf-ntriples.txt
- Original Record: 831-1.0075388-source.json
- Full Text
- 831-1.0075388-fulltext.txt
- Citation
- 831-1.0075388.ris
Full Text
Cite
Citation Scheme:
Usage Statistics
Share
Embed
Customize your widget with the following options, then copy and paste the code below into the HTML
of your page to embed this item in your website.
<div id="ubcOpenCollectionsWidgetDisplay">
<script id="ubcOpenCollectionsWidget"
src="{[{embed.src}]}"
data-item="{[{embed.item}]}"
data-collection="{[{embed.collection}]}"
data-metadata="{[{embed.showMetadata}]}"
data-width="{[{embed.width}]}"
async >
</script>
</div>
Our image viewer uses the IIIF 2.0 standard.
To load this item in other compatible viewers, use this url:
http://iiif.library.ubc.ca/presentation/dsp.831.1-0075388/manifest