Computing the Stochastic Complexity of Simple Probabilistic Graphical Models

Minimum Description Length (MDL) is an information-theoretic principle that can be used for model selection and other statistical inference tasks. There are various ways to use the principle in practice. One theoretically valid way is to use the normalized maximum likelihood (NML) criterion. Due to computational difficulties, this approach has not been used very often. This thesis presents efficient floating-point algorithms that make it possible to compute the NML for multinomial, Naive Bayes and Bayesian forest models. None of the presented algorithms rely on asymptotic analysis and with the first two model classes we also discuss how to compute exact rational number solutions.

Production process of Scots pine.

The accompanying collective research report is the result of the research project in 1986­90 between The Finnish Academy and the former Soviet Academy of Sciences. The project was organized around common field work in Finland and in the former Soviet Union and theoretical analyses of tree growth determining processes. Based on theoretical analyses, dynamic stand growth models were made and their parameters were determined utilizing the field results. Annual cycle affects the tree growth. Our theoretical approach was based on adaptation to local climate conditions from Lapland to South Russia. The initiation of growth was described as a simple low and high temperature accumulation driven model. Linking the theoretical model with long term temperature data allowed us to analyze what type of temperature response produced favorable outcome in different climates. Initiation of growth consumes the carbohydrate reserves in plants. We measured the dynamics of insoluble and soluble sugars in the very northern and Karelian conditions. Clear cyclical pattern was observed but the differences between locations were surprisingly small. Analysis of field measurements of CO2 exchange showed that irradiance is the dominating factor causing variation in photosynthetic rate in natural conditions during summer. The effect of other factors is so small that they can be omitted without any considerable loss of accuracy. A special experiment carried out in Hyytiälä showed that the needle living space, defined as the ratio between the shoot cylindric volume and needle surface area, correlates with the shoot photosynthesis. The penetration of irradiance into Scots pine canopy is a complicated phenomenon because of the movement of the sun on the sky and the complicated structure of branches and needles. A moderately simple but balanced forest radiation regime submodel was constructed. It consists of the tree crown and forest structure, the gap probability calculation and the consideration of spatial and temporal variation of radiation inside the forest. The common field excursions in different geographical regions resulted in a lot of experimental data of regularities of woody structures. The water transport seems to be a good common factor to analyse these properties of tree structure. There are evident regressions between cross-sectional areas measured at different locations along the water pathway from fine roots to needles. The observed regressions have clear geographical trends. For example, the same cross-sectional area can support three times higher needle mass in South Russia than in Lapland. Geographical trends can also be seen in shoot and needle structure. Analysis of data published by several Russian authors show, that one ton of needles transpire 42 ton of water a year. This annual amount of transpiration seems to be independent of geographical location, year and site conditions. The produced theoretical and experimental material is utilised in the development of stand growth model that describes the growth and development of Scots pine stands in Finland and the former Soviet Union. The core of the model is carbon and nutrient balances. This means that carbon obtained in photosynthesis is consumed for growth and maintenance and nutrients are taken according to the metabolic needs. The annual photosynthetic production by trees in the stand is determined as a function of irradiance and shading during the active period. The utilisation of the annual photosynthetic production to the growth of different components of trees is based on structural regularities. Since the fundamental metabolic processes are the same in all locations the same growth model structure can be applied in the large range of Scots pine. The annual photosynthetic production and structural regularities determining the allocation of resources have geographical features. The common field measurements enable the application of the model to the analysis of growth and development of stands growing on the five locations of experiments. The model enables the analysis of geographical differences in the growth of Scots pine. For example, the annual photosynthetic production of a 100-year-old stand at Voronez is 3.5 times higher than in Lapland. The share consumed to needle growth (30 %) and to growth of branches (5 %) seems to be the same in all locations. In contrast, the share of fine roots is decreasing when moving from north to south. It is 20 % in Lapland, 15 % in Hyytiälä Central Finland and Kentjärvi Karelia and 15 % in Voronez South Russia. The stem masses (115­113 ton/ha) are rather similar in Hyytiälä, Kentjärvi and Voronez, but rather low (50 ton/ha) in Lapland. In Voronez the height of the trees reach 29 m being in Hyytiälä and Kentjärvi 22 m and in Lapland only 14 m. The present approach enables utilization of structural and functional knowledge, gained in places of intensive research, in the analysis of growth and development of any stand. This opens new possibilities for growth research and also for applications in forestry practice.

The structure of Norway spruce (Picea abies [L.] Karst.) stems in relation to wood properties of sawn timber

An important challenge in forest industry is to get the appropriate raw material out from the forests to the wood processing industry. Growth and stem reconstruction simulators are therefore increasingly integrated in industrial conversion simulators, for linking the properties of wooden products to the three-dimensional structure of stems and their growing conditions. Static simulators predict the wood properties from stem dimensions at the end of a growth simulation period, whereas in dynamic approaches, the structural components, e.g. branches, are incremented along with the growth processes. The dynamic approach can be applied to stem reconstruction by predicting the three-dimensional stem structure from external tree variables (i.e. age, height) as a result of growth to the current state. In this study, a dynamic growth simulator, PipeQual, and a stem reconstruction simulator, RetroSTEM, are adapted to Norway spruce (Picea abies [L.] Karst.) to predict the three-dimensional structure of stems (tapers, branchiness, wood basic density) over time such that both simulators can be integrated in a sawing simulator. The parameterisation of the PipeQual and RetroSTEM simulators for Norway spruce relied on the theoretically based description of tree structure developing in the growth process and following certain conservative structural regularities while allowing for plasticity in the crown development. The crown expressed both regularity and plasticity in its development, as the vertical foliage density peaked regularly at about 5 m from the stem apex, varying below that with tree age and dominance position (Study I). Conservative stem structure was characterized in terms of (1) the pipe ratios between foliage mass and branch and stem cross-sectional areas at crown base, (2) the allometric relationship between foliage mass and crown length, (3) mean branch length relative to crown length and (4) form coefficients in branches and stem (Study II). The pipe ratio between branch and stem cross-sectional area at crown base, and mean branch length relative to the crown length may differ in trees before and after canopy closure, but the variation should be further analysed in stands of different ages and densities with varying site fertilities and climates. The predictions of the PipeQual and RetroSTEM simulators were evaluated by comparing the simulated values to measured ones (Study III, IV). Both simulators predicted stem taper and branch diameter at the individual tree level with a small bias. RetroSTEM predictions of wood density were accurate. For focusing on even more accurate predictions of stem diameters and branchiness along the stem, both simulators should be further improved by revising the following aspects in the simulators: the relationship between foliage and stem sapwood area in the upper stem, the error source in branch sizes, the crown base development and the height growth models in RetroSTEM. In Study V, the RetroSTEM simulator was integrated in the InnoSIM sawing simulator, and according to the pilot simulations, this turned out to be an efficient tool for readily producing stand scale information about stem sizes and structure when approximating the available assortments of wood products.

Plant secondary compounds and soil microbial processes in carbon and nitrogen cycling in relation to tree species

The aim of this study was to explore soil microbial activities related to C and N cycling and the occurrence and concentrations of two important groups of plant secondary compounds, terpenes and phenolic compounds, under silver birch (Betula pendula Roth), Norway spruce (Picea abies (L.) Karst) and Scots pine (Pinus sylvestris L.) as well as to study the effects of volatile monoterpenes and tannins on soil microbial activities. The study site, located in Kivalo, northern Finland, included ca. 70-year-old adjacent stands dominated by silver birch, Norway spruce and Scots pine. Originally the soil was very probably similar in all three stands. All forest floor layers (litter (L), fermentation layer (F) and humified layer (H)) under birch and spruce showed higher rates of CO2 production, greater net mineralisation of nitrogen and higher amounts of carbon and nitrogen in microbial biomass than did the forest floor layers under pine. Concentrations of mono-, sesqui-, di- and triterpenes were higher under both conifers than under birch, while the concentration of total water-soluble phenolic compounds as well as the concentration of condensed tannins tended to be higher or at least as high under spruce as under birch or pine. In general, differences between tree species in soil microbial activities and in concentrations of secondary compounds were smaller in the H layer than in the upper layers. The rate of CO2 production and the amount of carbon in the microbial biomass correlated highly positively with the concentration of total water-soluble phenolic compounds and positively with the concentration of condensed tannins. Exposure of soil to volatile monoterpenes and tannins extracted and fractionated from spruce and pine needles affected carbon and nitrogen transformations in soil, but the effects were dependent on the compound and its molecular structure. Monoterpenes decreased net mineralisation of nitrogen and probably had a toxic effect on part of the microbial population in soil, while another part of the microbes seemed to be able to use monoterpenes as a carbon source. With tannins, low-molecular-weight compounds (also compounds other than tannins) increased soil CO2 production and nitrogen immobilisation by soil microbes while the higher-molecular-weight condensed tannins had inhibitory effects. In conclusion, plant secondary compounds may have a great potential in regulation of C and N transformations in forest soils, but the real magnitude of their significance in soil processes is impossible to estimate.

Lignin biosynthesis in Norway spruce: from a model system to the tree

Lignin is a hydrophobic polymer that is synthesised in the secondary cell walls of all vascular plants. It enables water conduction through the stem, supports the upright growth habit and protects against invading pathogens. In addition, lignin hinders the utilisation of the cellulosic cell walls of plants in pulp and paper industry and as forage. Lignin precursors are synthesised in the cytoplasm through the phenylpropanoid pathway, transported into the cell wall and oxidised by peroxidases or laccases to phenoxy radicals that couple to form the lignin polymer. This study was conducted to characterise the lignin biosynthetic pathway in Norway spruce (Picea abies (L.) Karst.). We focused on the less well-known polymerisation stage, to identify the enzymes and the regulatory mechanisms that are involved. Available data for lignin biosynthesis in gymnosperms is scarce and, for example, the latest improvements in precursor biosynthesis have only been verified in herbaceous plants. Therefore, we also wanted to study in detail the roles of individual gene family members during developmental and stress-induced lignification, using EST sequencing and real-time RT-PCR. We used, as a model, a Norway spruce tissue culture line that produces extracellular lignin into the culture medium, and showed that lignin polymerisation in the tissue culture depends on peroxidase activity. We identified in the culture medium a significant NADH oxidase activity that could generate H2O2 for peroxidases. Two basic culture medium peroxidases were shown to have high affinity to coniferyl alcohol. Conservation of the putative substrate-binding amino acids was observed when the spruce peroxidase sequences were compared with other peroxidases with high affinity to coniferyl alcohol. We also used different peroxidase fractions to produce synthetic in vitro lignins from coniferyl alcohol; however, the linkage pattern of the suspension culture lignin could not be reproduced in vitro with the purified peroxidases, nor with the full complement of culture medium proteins. This emphasised the importance of the precursor radical concentration in the reaction zone, which is controlled by the cells through the secretion of both the lignin precursors and the oxidative enzymes to the apoplast. In addition, we identified basic peroxidases that were reversibly bound to the lignin precipitate. They could be involved, for example, in the oxidation of polymeric lignin, which is required for polymer growth. The dibenzodioxocin substructure was used as a marker for polymer oxidation in the in vitro polymerisation studies, as it is a typical substructure in wood lignin and in the suspension culture lignin. Using immunolocalisation, we found the structure mainly in the S2+S3 layers of the secondary cell walls of Norway spruce tracheids. The structure was primarily formed during the late phases of lignification. Contrary to the earlier assumptions, it appears to be a terminal structure in the lignin macromolecule. Most lignin biosynthetic enzymes are encoded for by several genes, all of which may not participate in lignin biosynthesis. In order to identify the gene family members that are responsible for developmental lignification, ESTs were sequenced from the lignin-forming tissue culture and developing xylem of spruce. Expression of the identified lignin biosynthetic genes was studied using real-time RT-PCR. Candidate genes for developmental lignification were identified by a coordinated, high expression of certain genes within the gene families in all lignin-forming tissues. However, such coordinated expression was not found for peroxidase genes. We also studied stress-induced lignification either during compression wood formation by bending the stems or after Heterobasidion annosum infection. Based on gene expression profiles, stress-induced monolignol biosynthesis appeared similar to the developmental process, and only single PAL and C3H genes were specifically up-regulated by stress. On the contrary, the up-regulated peroxidase genes differed between developmental and stress-induced lignification, indicating specific responses.

Reproduction and population structure in European aspen

The European aspen (Populus tremula) is a keystone species for biodiversity in boreal forests. However, the future of aspen may be threatened, because large aspens have mostly been removed from managed forests, whereas regeneration and the long-term persistence of mature trees are subjects of concern in protected areas. Aspen is a pioneer tree, and it can reproduce both sexually by seed and asexually by root suckers. Through asexual reproduction aspen forms clones, groups of genetically identical trees (ramets). In my thesis, I have studied the structure of aspen populations in terms of number, size, clonal and demographic properties. Additionally, I have investigated the emergence and survival of seedlings as well as the seed quantity and quality in crosses between the European and hybrid aspen. To study the regeneration and population structure, mature aspens were recorded in old-growth and managed forests in eastern Finland based on a large-scale inventory (11 400 ha). In addition, small aspen trees were surveyed on sample plots. Clonal structure was investigated both by morphological characters and by DNA-based markers (microsatellites). Seedling emergence and survival was studied with two sowing experiments. With crosses between European and hybrid aspens we wanted to study whether elevated temperatures due to climate change would benefit the different crosses of European and hybrid aspen unequally and thus affect the gene flow between the two species. The average volumes of mature aspen were 5.3 m3/ha in continuous old-growth, and 0.8 m3/ha in managed forests. Results indicate also that large aspen trees in managed forests are a legacy of the past less intensively managed forest landscapes. Long-term persistence of aspen in protected areas can only be secured by restoration measures creating sufficiently large gaps for regeneration. More emphasis should be given to sparing aspens in thinnings and to retaining of mature aspens in regeneration cutting in managed forests. Aspen was found to be spatially aggregated in the landscape. This could be explained by site type, disturbance history and / or limitations in seed dispersal. Clonal structure does not explain the spatial aggregation, since average size of the clones was only 2.3 ramets, and most clones (70 %) consisted of just one ramet. The small size of the clones suggests that most of them are relatively young. Therefore, sexual reproduction may be more common than has previously been thought. Seedling emergence was most successful in mineral soil especially, when the site had been burned. Only few seedlings occurred on humus. Survival of the seedlings was low, and strongly dependent on moisture, but also on seedbed conditions. The seeds were found to maintain their germinability longer than has earlier been thought to be possible. Interspecific crosses produced more seeds with higher quality than intraspecific crosses. When temperature was elevated, germination of hybrid aspen seeds increased more than seeds from P. tremula x P. tremula crosses. These results suggest that hybrid aspen may have a significant genetic impact on the European aspen, and this effect may become strengthened by climate warming.

Microsite occupancy and spatial structure of regeneration in three late-successional Norway spruce forests in northern Europe

A comparison of microsite occupancy and the spatial structure of regeneration in three areas of late-successional Norway spruce dominated forest. Pallas-Ylläs is understood to have been influenced only by small-scale disturbance; Dvina-Pinega has had sporadic larger-scale disturbances; Kazkim has been affected by fire. All spruce and birch trees with diameter at breast height (DBH) ?10 cm were mapped in five stands on 40 m x 400 m transects, and those with DBH < 10 cm on 2 or 4 m x 400 m subplots. Microsite type was inventoried at 1m intervals along the centre line and for each tree with DBH < 10 cm. At all study areas small seedlings (h < 0.3 m, DBH < 10 cm) preferentially occupied disturbed microsites. In contrast, spruce saplings (h ? 1.3 m, DBH <10 cm) at all study areas showed less, or no, preference. At Pallas-Ylläs spruce seedlings (h < 1.3 m, DBH < 10 cm) and saplings (h ? 1.3 m, DBH < 10 cm) exhibited spatial correlation at scales from 32-52 m. At Dvina-Pinega saplings of both spruce and birch exhibited spatial correlation at scales from 32-81 m. At Kazkim spatial correlation of seedlings and saplings of both species was exhibited over variable distances. No spatial cross-correlation was found between overstorey basal area (DBH ? 10 cm) and regeneration (h ? 1.3 m, DBH < 10 cm) at any study area. The results confirm the importance of disturbed microsites for seedling establishment, but suggest that undisturbed microsites may sometimes be more advantageous for long-term tree survival. The regeneration gap concept may not be useful in describing the regeneration dynamics of late-successional boreal forests.

Algorithms for Exact Structure Discovery in Bayesian Networks

Bayesian networks are compact, flexible, and interpretable representations of a joint distribution. When the network structure is unknown but there are observational data at hand, one can try to learn the network structure. This is called structure discovery. This thesis contributes to two areas of structure discovery in Bayesian networks: space--time tradeoffs and learning ancestor relations. The fastest exact algorithms for structure discovery in Bayesian networks are based on dynamic programming and use excessive amounts of space. Motivated by the space usage, several schemes for trading space against time are presented. These schemes are presented in a general setting for a class of computational problems called permutation problems; structure discovery in Bayesian networks is seen as a challenging variant of the permutation problems. The main contribution in the area of the space--time tradeoffs is the partial order approach, in which the standard dynamic programming algorithm is extended to run over partial orders. In particular, a certain family of partial orders called parallel bucket orders is considered. A partial order scheme that provably yields an optimal space--time tradeoff within parallel bucket orders is presented. Also practical issues concerning parallel bucket orders are discussed. Learning ancestor relations, that is, directed paths between nodes, is motivated by the need for robust summaries of the network structures when there are unobserved nodes at work. Ancestor relations are nonmodular features and hence learning them is more difficult than modular features. A dynamic programming algorithm is presented for computing posterior probabilities of ancestor relations exactly. Empirical tests suggest that ancestor relations can be learned from observational data almost as accurately as arcs even in the presence of unobserved nodes.

The Phraseology and Structure of Latin Building Inscriptions in Roman North Africa

This study analyses the diction of Latin building inscriptions. Despite its importance, this topic has rarely been discussed before: the most substantial contribution on the subject is a short dissertation by Klaus Gast (1965) that focuses on 100 inscriptions dating mostly from the Republican period. Marietta Horster (2001) also touched upon this theme in her thesis on imperial building inscriptions. I have collected my source material in North Africa because more Latin building inscriptions dating from the Imperial period have survived there than in any other area of the Roman Empire. By means of a thorough and independent survey, I have assembled all relevant African Latin building inscriptions datable to the Roman period (between 146 BC and AD 425), 1002 texts, into a corpus. These inscriptions are all fully edited in Appendix 1; Appendix 2 contains references to earlier editions. To facilitate search operations, both are also available in electronic form. They are downloadable from the address http://www.helsinki.fi/hum/kla/htm/jatkoopinnot.htm. Chapter one is an introduction dealing with the nature of building inscriptions as source material. Chapter two offers a statistical overview of the material. The following main section of the work falls into five chapters, each of which analyses one main part of a building inscription. An average building inscription can be divided into five parts: the starting phrase opens the inscription (a dedication to gods, for example), the subject part identifies the builder, the object part describes the constructed or repaired building, the predicate part records the building activity and the supplement part offers additional information on the project (it can specify the funding, for instance). These chapters are systematic and chronological and their purpose is to register and interpret the phrases used, to analyse reasons for their use and for their popularity among the different groups of builders. Chapter eight, which follows the main section of the work, creates a typology of building inscriptions based on their structure. It also presents the most frequently attested types of building inscriptions. The conclusion describes, on a general level, how the diction of building inscriptions developed during the period of study and how this striking development resulted from socio-economic changes that took place in Romano-African society during Antiquity. This study shows that the phraseology of building inscriptions had a clear correlation both with the type of builder and with the date of carving. Private builders tended to accentuate their participation (especially its financial side) in the project; honouring the emperor received more emphasis in the building inscriptions set up by communities; the texts produced by the army were concise. The chronological development is so clear that it enables stylistic dating. At the beginning of the imperial period the phrases were clear, concrete, formal and stereotyped but by Late Antiquity they have become vague, subjective, flexible, varied and even rhetorically or poetically coloured.

Structure informationnelle et constructions du kabyle : Etude de trois types de phrase dans le cadre de la grammaire constructionnelle

Information structure and Kabyle constructions Three sentence types in the Construction Grammar framework The study examines three Kabyle sentence types and their variants. These sentence types have been chosen because they code the same state of affairs but have different syntactic structures. The sentence types are Dislocated sentence, Cleft sentence, and Canonical sentence. I argue first that a proper description of these sentence types should include information structure and, second, that a description which takes into account information structure is possible in the Construction Grammar framework. The study thus constitutes a testing ground for Construction Grammar for its applicability to a less known language. It constitutes a testing ground notably because the differentiation between the three types of sentences cannot be done without information structure categories and, consequently, these categories must be integrated also in the grammatical description. The information structure analysis is based on the model outlined by Knud Lambrecht. In that model, information structure is considered as a component of sentence grammar that assures the pragmatically correct sentence forms. The work starts by an examination of the three sentence types and the analyses that have been done in André Martinet s functional grammar framework. This introduces the sentence types chosen as the object of study and discusses the difficulties related to their analysis. After a presentation of the state of the art, including earlier and more recent models, the principles and notions of Construction Grammar and of Lambrecht s model are introduced and explicated. The information structure analysis is presented in three chapters, each treating one of the three sentence types. The analyses are based on spoken language data and elicitation. Prosody is included in the study when a syntactic structure seems to code two different focus structures. In such cases, it is pertinent to investigate whether these are coded by prosody. The final chapter presents the constructions that have been established and the problems encountered in analysing them. It also discusses the impact of the study on the theories used and on the theory of syntax in general.

Biophysical Characterization of Oxidized Phospholipids : Implications for Altered Membrane Structure and Function

The present study aims to elucidate the modifications in the structure and functionality of the phospholipid matrix of biological membranes brought about by free radical-mediated oxidative damage of its molecular constituents. To this end, the surface properties of two oxidatively modified phospholipids bearing an aldehyde or carboxyl function at the end of truncated sn-2 acyl chain were studied using a Langmuir balance. The results obtained reveal both oxidized species to have a significant impact on the structural dynamics of phospholipid monolayers, as illustrated by the progressive changes in force-area isotherms with increasing mole fraction of the oxidized lipid component. Moreover, surface potential measurements revealed considerable modifications in the electric properties of oxidized phospholipid containing monolayers during film compression, suggesting a packing state-controlled reorientation of the intramolecular electric dipoles of the lipid headgroups and acyl chains. Based on the above findings, a model describing the conformational state of oxidized phospholipid molecules in biological membranes is proposed, involving the protrusion of the acyl chains bearing the polar functional groups out from the hydrocarbon phase to the surrounding aqueous medium. Oxidative modifications alter profoundly the physicochemical properties of unsaturated phospholipids and are therefore readily anticipated to have important implications for their interactions with membrane-associating molecules. Along these lines, the carboxyl group bearing lipid was observed to bind avidly the peripheral membrane protein cytochrome c. The binding was reversed following increase in ionic strength or addition of polyanionic ATP, thus suggesting it to be driven by electrostatic interactions between cationic residues of the protein and the deprotonated lipid carboxyl exposed to the aqueous phase. The presence of aldehyde function bearing oxidized phospholipid was observed to enhance the intercalation of four antimicrobial peptides into phospholipid monolayers and liposomal bilayers. Partitioning of the peptides to monolayers was markedly attenuated by the aldehyde scavenger methoxyamine, revealing it to be mediated by the carbonyl moiety possibly through efficient hydrogen bonding or, alternatively, formation of covalent adduct in form of a Schiff base between the lipid aldehydes and primary amine groups of the peptide molecules. Lastly, both oxidized phospholipid species were observed to bind with high affinity three small membrane-partitioning therapeutic agents, viz. chlorpromazine, haloperidol, and doxorubicin. In conclusion, the results of studies conducted using biomimetic model systems support the notion that oxidative damage influences the molecular architecture as well as the bulk physicochemical properties of phospholipid membranes. Further, common polar functional groups carried by phospholipids subjected to oxidation were observed to act as molecular binding sites at the lipid-water interface. It is thus plausible that oxidized phospholipid species may elicit cellular level effects by modulating integration of various membrane-embedded and surface-associated proteins and peptides, whose conformational state, oligomerization, and functionality is known to be controlled by highly specific lipid-protein interactions and proper physical state of the membrane environment.