Massey Documents by Type

Permanent URI for this communityhttps://mro.massey.ac.nz/handle/10179/294

Browse

Search Results

Now showing 1 - 2 of 2
  • Item
    Applications of linear modelling in energy analysis : a thesis presented in partial fulfilment of the requirements for the degree of Doctor of Philosophy in Technology at Massey University
    (Massey University, 1984) Patterson, Murray Graham
    The primary objective of this study was to explore the use of classical linear models in Energy Analysis; so as to resolve some of the methodological problems associated with Energy Analysis, and to extend the scope and potential of Energy Analysis as a scientific discipline. This was undertaken in the form of two separate yet related discourses. The first discourse provided a basis for resolving the energy quality problem encountered in Energy Analysis. A general equation was hypothesized and tested: [Equation] The matrix of simultaneous linear equations represented by this equation is usually overdetermined. Therefore, an appropriate solution method is a fitting procedure, such as regression. Further, in order to solve this general equation, one coefficient must arbitrarily be given a value equal to unity. Hence, all estimated coefficients are expressed in terms of multiples of that coefficient (termed quality equivalents). The general equation was first tes ted for the 1 97 6 New Zealand economy, so as to estimate 'actual' quality coefficients. Subsequently, the general equation was tested for a notional 'energy efficient' New Zealand economy, so as to estimate 'long run' quality coefficients, which reflected thermodynamic limits. Generally very accurate estimates of the coefficients were obtained. The solutions to the equations indicated that hydroelectricity was the highest quality primary energy source, followed by natural gas, oil, coal and then wood. The second discourse examined the 'optimal' use of primary energy resources in the New Zealand food system, using the formalism of Linear Programming. A preliminary discussion concluded that the concept of 'optimality' had greater potency than the concept of 'efficiency' in evaluating the use of energy resources, particularly in food systems. For each food sector (Production, Processing, Export-Import, Distribution, Catering and Household), coefficient matrices were assembled, drawing on literature data. Various combinations of constraints and objective functions were applied, in different Runs. The main objective functions used were minimising energy inputs for providing a nutritionally adequate diet , or maximising net energy gain from exporting agro-food products. The most critical constraints were found to be land area and market demand constraints. Detailed results of the Linear Programming runs are presented and discussed. An energetically 'optimal' diet was found to consist of large amounts of cereals, significant amounts of fresh fruit and vegetables and dairy products, and a very small amount of meat. Such a diet provided the Recommended Daily Allowances, for all nutritional elements, for the New Zealand population. Meat, Fish and Dairy products were found to have a particularly important function as commodities to be traded for imports of oil and energy intensive goods. A final discussion reviewed the use of Linear Models in Energy Analysis, and future directions for growth and development in Energy Analysis.
  • Item
    On the geometry of generalized linear models : submitted in partial fulfillment of the requirements for the degree of Doctor of Philosophy at Massey University, Palmerston North, New Zealand
    (Massey University, 2003) Luo, Dongwen
    The perspective afforded by Euclidean geometry led to the rapid development of linear models in the early stages of the twentieth century: Fisher saw the data as a point in finite-dimensional Euclidean space, the model as a subspace and least squares fitting as projection of the observation vector onto the model space. From the late 1960s to early 1970s, Fienberg revealed geometry underlying loglinear models for two-way tables, while Haberman discussed geometry for the log-transformed case. Generalized linear models, however, have largely eluded geometers until recently. In 1997 an extension of Fisher's view to generalized linear models was given by Kass and Vos, using the language of differential geometry. The aim of this work is to develop a simple, general geometric framework for generalized linear models, closely related to the thinking of Fienberg and Haberman. Whereas Kass and Vos developed a geometric view which leads to the usual scoring method, we develop geometry which leads to a new algorithm. A linearization of this new algorithm yields the scoring method. The geometry discussed by Kass and Vos is based on the log-likelihood function whereas the geometry developed here depends on sufficiency. In the geometry of generalized linear models, developed through chapters 1 to 3, an observation with n values is viewed as a vector in Euclidian space Rn. This Euclidian space Rn is partitioned into two orthogonal spaces, the sufficiency space S and the auxiliary space A, with respect to a new basis. We focus on two mean sets relating to generalized linear models, one for the untransformed model space and another for the link-transformed model space. There are two critical properties of the maximum likelihood estimate of the parameters of a generalized linear model with canonical link. The first property is that the coefficients of the basis of the sufficiency space, the sufficient statistics, are preserved in the untransformed model space in the fitting process. The second property is that the coefficients of the basis of the auxiliary space are zeroed in the link-transformed model space in the fitting process. Linear models and loglinear models serve as special cases of generalized linear models with identity and log link respectively. Based on the geometric framework discussed in the thesis, a new algorithm is constructed for fitting generalized linear models with canonical link in Chapter 4. This algorithm, which relies on sufficient statistics for the parameters in the model rather than the likelihood function, takes two projections alternately, orthogonal projection onto a sufficiency affine plane and non-orthogonal projection onto the transformed model space. In the process, we match the model space and sufficient statistics iteratively until convergence. Linearization of the new algorithm induces the scoring method. In Chapter 5 we pay special attention to a subset of loglinear models, graphical loglinear models, those which are the intersection of a finite set of conditional independence statements. The model space of one conditional independence statement is described through the notions of "corresponding point convex hull" and "set convex hull". The fitting of one conditional independence statement is considered geometrically using a direct fitting method and the familiar iterative proportional fitting method.