Statistical analysis of the associations of hereditary factors with the risk of Type 1 diabetes and Diabetic Nephropathy based on the familial data collected through ascertaiment from population-based registers

In genetic epidemiology, population-based disease registries are commonly used to collect genotype or other risk factor information concerning affected subjects and their relatives. This work presents two new approaches for the statistical inference of ascertained data: a conditional and full likelihood approaches for the disease with variable age at onset phenotype using familial data obtained from population-based registry of incident cases. The aim is to obtain statistically reliable estimates of the general population parameters. The statistical analysis of familial data with variable age at onset becomes more complicated when some of the study subjects are non-susceptible, that is to say these subjects never get the disease. A statistical model for a variable age at onset with long-term survivors is proposed for studies of familial aggregation, using latent variable approach, as well as for prospective studies of genetic association studies with candidate genes. In addition, we explore the possibility of a genetic explanation of the observed increase in the incidence of Type 1 diabetes (T1D) in Finland in recent decades and the hypothesis of non-Mendelian transmission of T1D associated genes. Both classical and Bayesian statistical inference were used in the modelling and estimation. Despite the fact that this work contains five studies with different statistical models, they all concern data obtained from nationwide registries of T1D and genetics of T1D. In the analyses of T1D data, non-Mendelian transmission of T1D susceptibility alleles was not observed. In addition, non-Mendelian transmission of T1D susceptibility genes did not make a plausible explanation for the increase in T1D incidence in Finland. Instead, the Human Leucocyte Antigen associations with T1D were confirmed in the population-based analysis, which combines T1D registry information, reference sample of healthy subjects and birth cohort information of the Finnish population. Finally, a substantial familial variation in the susceptibility of T1D nephropathy was observed. The presented studies show the benefits of sophisticated statistical modelling to explore risk factors for complex diseases.

Study of Friction and Transfer Layer Formation in Copper-Steel Tribo-System: Role of Surface Texture and Roughness Parameters

In the present investigation, experiments were conducted on a tribological couple-copper pin against steel plate-using an inclined pin-on-plate sliding tester to understand the role of surface texture and roughness parameters of the plate on the coefficient friction and transfer layer formation. Two surface characteristics of the steel plates-roughness and texture-were varied in the tests. It was observed that the transfer layer formation and the coefficient of friction along with its two components, namely, the adhesion and plowing, are controlled by the surface texture of the plate. The plowing component of friction was highest for the surface texture that promotes plane strain conditions while it was lowest for the texture that favors plane stress conditions at the interface. Dimensionless quantifiable roughness parameters were formulated to describe the degree of plowing and hence the plane strain/stress type deformations taking place at the asperity level.

Compressible 2nd-order boundary-layers for 3-dimensional stagnation point flows with mass-transfer

All the second-order boundary-layer effects have been studied for the steady laminar compressible 3-dimensional stagnation-point flows with variable properties and mass transfer for both saddle and nodal point regions. The governing equations have been solved numerically using an implicit finite-difference scheme. Results for the heat transfer and skin friction have been obtained for several values of the mass-transfer rate, wall temperature, and also for several values of parameters characterizing the nature of stagnation point and variable gas properties. The second-order effects on the heat transfer and skin friction at the wall are found to be significant and at large injection rates, they dominate over the results of the first-order boundary layer, but the effect of large suction is just the opposite. In general, the second-order effects are more pronounced in the saddle-point region than in the nodal-point region. The overall heat-transfer rate for the 3-dimensional flows is found to be more than that of the 2-dimensional flows.

Carbon grazing: Managing carbon flows for resilient pasture in variable climate

Work with consultants to help graziers understand implications of carbon in their pastures.

Simultaneous prediction of acidity parameters (pH and titratable acidity) in Kefir using near infrared reflectance spectroscopy.

Acidity in terms of pH and titratable acids influences the texture and flavour of fermented dairy products, such as Kefir. However, the methods for determining pH and titratable acidity (TA) are time consuming. Near infrared (NIR) spectroscopy is a non-destructive method, which simultaneously predicts multiple traits from a single scan and can be used to predict pH and TA. The best pH NIR calibration model was obtained with no spectral pre-treatment applied, whereas smoothing was found to be the best pre-treatment to develop the TA calibration model. Using cross-validation, the prediction results were found acceptable for both pH and TA. With external validation, similar results were found for pH and TA, and both models were found to be acceptable for screening purposes.

Influence of plant traits on production costs and profitability of strawberry in southeast Queensland.

In Queensland the subtropical strawberry ( Fragaria * ananassa) breeding program aims to combine traits into novel genotypes that increase production efficiency. The contribution of individual plant traits to cost and income under subtropical Queensland conditions was investigated, with the overall goal of improving the profitability of the industry through the release of new strawberry cultivars. The study involved specifying the production and marketing system using three cultivars of strawberry that are currently widely grown annually in southeast Queensland, developing methods to assess the economic impact of changes to the system, and identifying plant traits that influence outcomes from the system. From May through September P (price; $ punnet -1), V (monthly mass; tonne of fruit on the market) and M (calendar month; i.e. May=5) were found to be related ( r2=0.92) by the function (SE) P=4.741(0.469)-0.001630(0.0005) V-0.226(0.102) M using data from 2006 to 2010 for the Brisbane central market. Both income and cost elements in the gross margin were subject to sensitivity analysis. 'Harvesting' and 'Handling/Packing' 'Groups' of 'Activities' were the major contributors to variable costs (each >20%) in the gross margin analysis. Within the 'Harvesting Group', the 'Picking Activity' contributed most (>80%) with the trait 'display of fruit' having the greatest (33%) influence on the cost of the 'Picking Activity'. Within the 'Handling/Packing Group', the 'Packing Activity' contributed 50% of costs with the traits 'fruit shape', 'fruit size variation' and 'resistance to bruising' having the greatest (12-62%) influence on the cost of the 'Packing Activity'. Non-plant items (e.g. carton purchases) made up the other 50% of the costs within the 'Handling/Packing Group'. When any of the individual traits in the 'Harvesting' and 'Handling/Packing' groups were changed by one unit (on a 1-9 scale) the gross margin changed by up to 1%. Increasing yield increased the gross margin to a maximum (15% above present) at 1320 g plant -1 (94% above present). A 10% redistribution of total yield from September to May increased the gross margin by 23%. Increasing fruit size increased gross margin: a 75% increase in fruit size (to ~30 g) produced a 22% increase in the gross margin. The modified gross margin analysis developed in this study allowed simultaneous estimation of the gross margin for the producer and gross value of the industry. These parameters sometimes move in opposite directions.

A four-variable model for the pattern-forming mechanism in Hydra

A generalized Gierer-Meinhardt model has been used to account for the transplantation experiments in Hydra. In this model, a cross inhibition between the two organizing centres (namely, head and foot) are assumed to be the only mode of interaction in setting up a stable morphogen distribution for the pattern formation in Hydra.

Soil carbon modelling as a tool for carbon balance studies in forestry

Soils represent a remarkable stock of carbon, and forest soils are estimated to hold half of the global stock of soil carbon. Topical concern about the effects of climate change and forest management on soil carbon as well as practical reporting requirements set by climate conventions have created a need to assess soil carbon stock changes reliably and transparently. The large spatial variability of soil carbon commensurate with relatively slow changes in stocks hinders the assessment of soil carbon stocks and their changes by direct measurements. Models therefore widely serve to estimate carbon stocks and stock changes in soils. This dissertation aimed to develop the soil carbon model YASSO for upland forest soils. The model was aimed to take into account the most important processes controlling the decomposition in soils, yet remain simple enough to ensure its practical applicability in different applications. The model structure and assumptions were presented and the model parameters were defined with empirical measurements. The model was evaluated by studying the sensitivities of the model results to parameter values, by estimating the precision of the results with an uncertainty analysis, and by assessing the accuracy of the model by comparing the predictions against measured data and to the results of an alternative model. The model was applied to study the effects of intensified biomass extraction on the forest carbon balance and to estimate the effects of soil carbon deficit on net greenhouse gas emissions of energy use of forest residues. The model was also applied in an inventory based method to assess the national scale forest carbon balance for Finland’s forests from 1922 to 2004. YASSO managed to describe sufficiently the effects of both the variable litter and climatic conditions on decomposition. When combined with the stand models or other systems providing litter information, the dynamic approach of the model proved to be powerful for estimating changes in soil carbon stocks on different scales. The climate dependency of the model, the effects of nitrogen on decomposition and forest growth as well as the effects of soil texture on soil carbon stock dynamics are areas for development when considering the applicability of the model to different research questions, different land use types and wider geographic regions. Intensified biomass extraction affects soil carbon stocks, and these changes in stocks should be taken into account when considering the net effects of forest residue utilisation as energy. On a national scale, soil carbon stocks play an important role in forest carbon balances.

Suunnittelunäkemys metsäorganisaatioissa ja sen vaikutus tietojärjestelmälle asetettaviin vaatimuksiin

Department of Forest Resource Management in the University of Helsinki has in years 2004?2007 carried out so-called SIMO -project to develop a new generation planning system for forest management. Project parties are organisations doing most of Finnish forest planning in government, industry and private owned forests. Aim of this study was to find out the needs and requirements for new forest planning system and to clarify how parties see targets and processes in today's forest planning. Representatives responsible for forest planning in each organisation were interviewed one by one. According to study the stand-based system for managing and treating forests continues in the future. Because of variable data acquisition methods with different accuracy and sources, and development of single tree interpretation, more and more forest data is collected without field work. The benefits of using more specific forest data also calls for use of information units smaller than tree stand. In Finland the traditional way to arrange forest planning computation is divided in two elements. After updating the forest data to present situation every stand unit's growth is simulated with different alternative treatment schedule. After simulation, optimisation selects for every stand one treatment schedule so that the management program satisfies the owner's goals in the best possible way. This arrangement will be maintained in the future system. The parties' requirements to add multi-criteria problem solving, group decision support methods as well as heuristic and spatial optimisation into system make the programming work more challenging. Generally the new system is expected to be adjustable and transparent. Strict documentation and free source code helps to bring these expectations into effect. Variable growing models and treatment schedules with different source information, accuracy, methods and the speed of processing are supposed to work easily in system. Also possibilities to calibrate models regionally and to set local parameters changing in time are required. In future the forest planning system will be integrated in comprehensive data management systems together with geographic, economic and work supervision information. This requires a modular method of implementing the system and the use of a simple data transmission interface between modules and together with other systems. No major differences in parties' view of the systems requirements were noticed in this study. Rather the interviews completed the full picture from slightly different angles. In organisation the forest management is considered quite inflexible and it only draws the strategic lines. It does not yet have a role in operative activity, although the need and benefits of team level forest planning are admitted. Demands and opportunities of variable forest data, new planning goals and development of information technology are known. Party organisations want to keep on track with development. One example is the engagement in extensive SIMO-project which connects the whole field of forest planning in Finland.

Genetic diversity of Atlantic salmon (Salmo salar) and European whitefish (Coregonus lavaretus) in the Baltic Sea

We described the patterns and extent of microsatellite DNA variation in historical and present-day Atlantic salmon (Salmo salar L.) stocks in the Baltic Sea and neighbouring areas, and in European whitefish (Coregonus lavaretus) ecotypes, populations and run-timing types in Finland. Moreover, the amount and pattern of genetic diversity in historical salmon populations before human impact were described, and the proportion of diversity maintained in the present hatchery stocks evaluated. Salmon populations in the Baltic Sea were, on average, significantly less variable than eastern Atlantic populations, and the diversity of landlocked populations (Lakes Vänern, Saimaa, Onega and Ladoga) was in turn significantly lower than that of anadromous salmon populations in the Baltic Sea populations. Within the Baltic Sea, the anadromous populations of Atlantic salmon formed three clear groups, corresponding to the northern (Gulf of Bothnia), eastern (Gulf of Finland and eastern Baltic Main Basin) and southern (western Baltic Main Basin) regions. Based on microsatellite data, three salmon population groups in the Baltic Sea were considered potentially different colonization lineages. In short- and long-term breeding programmes of Atlantic salmon, the average observed rate of loss of alleles was 4.9% and 2.0% per generation and the average rate of loss of heterozygosity was 1.4% and 1% per generation, respectively. When comparing the genetic parameters of stocks before and after hatchery breeding of several successive generations (Rivers Iijoki and Oulujoki), statistically significant changes in allele frequencies were common, while large wild stock in the Teno River has remained temporally very stable over 56 years. Despite the observed losses of genetic diversity in broodstock breeding, a large proportion of the genetic resources of the extirpated stocks are still conserved in the broodstocks. Genetic differentiation among European whitefish ecotypes was generally low, thus giving support to the hypothesis of one native European whitefish species in Fennoscandia. Among the ecotypes, the northern, large sparsely rakered, bottom-dwelling whitefish was the most unique. The known genetic differences in quantitative traits have thus either developed independently of potential phylogenetic lineages, or the lineages have mixed and the quantitative traits of the ecotypes, like gill-raker number, have later changed according to environment and selection pressures. Overall, genetic distances between the anadromous whitefish populations along the Finnish coast, especially in the Bothnian Bay area, were small. Wild whitefish populations studied had slightly higher allelic diversity than hatchery-reared populations in corresponding rivers.

Male traits and herd reproductive capability in tropical beef cattle. 2. Genetic parameters of bull traits

A total of 4063 young bulls of two tropical genotypes (1639 Brahman and 2424 Tropical Composite) raised in northern Australia were evaluated for a comprehensive range of production and reproduction traits up to 24 months of age. Prior to weaning, peripheral blood concentrations of luteinising hormone (LH) and inhibin were measured at 4 months of age. At weaning (6 months) blood insulin-like growth factor-1 (IGF-I) and flight time were recorded. Body composition traits of fat depth and eye-muscle area were determined by ultrasonography at 15 months of age when additional measurements of liveweight, hip height and body condition score were recorded. Bull breeding soundness was evaluated at similar to 12, 18 and 24 months of age when measurements of scrotal circumference, sheath score, semen mass activity, progressive motility of individual sperm and percent morphologically normal sperm were recorded. Magnitude of heritability and genetic correlations changed across time for some traits. Heritability of LH, inhibin, IGF-I and of 18-month scrotal circumference, mass activity, progressive motility and percent normal sperm was 0.31, 0.74, 0.44, 0.75, 0.24, 0.15 and 0.25, respectively, for Brahmans and 0.48, 0.72, 0.36, 0.43, 0.13, 0.15 and 0.20, respectively, for Tropical Composites. Inhibin and IGF-I had moderate genetic association with percent normal sperm at 24 months in Brahmans but low to negligible associations in Tropical Composites. Body condition score in Brahmans and sperm motility (mass and individual) traits in both genotypes had moderate to strong genetic correlation with percent normal sperm and may prove useful candidates for indirect selection. There is scope to increase scrotal circumference by selection and this will be associated with favourable correlated responses of improved semen quality in both genotypes. The lack of genetic antagonism among bull traits indicates that selection for improved semen quality will not adversely affect other production traits.

Marker selection by Akaike information criterion and Bayesian information criterion

We carried out a discriminant analysis with identity by descent (IBD) at each marker as inputs, and the sib pair type (affected-affected versus affected-unaffected) as the output. Using simple logistic regression for this discriminant analysis, we illustrate the importance of comparing models with different number of parameters. Such model comparisons are best carried out using either the Akaike information criterion (AIC) or the Bayesian information criterion (BIC). When AIC (or BIC) stepwise variable selection was applied to the German Asthma data set, a group of markers were selected which provide the best fit to the data (assuming an additive effect). Interestingly, these 25-26 markers were not identical to those with the highest (in magnitude) single-locus lod scores.

Optimal simultaneous maximum a posterioriestimation of states, noise statistics and parameters I. Algorithm

The simultaneous state and parameter estimation problem for a linear discrete-time system with unknown noise statistics is treated as a large-scale optimization problem. The a posterioriprobability density function is maximized directly with respect to the states and parameters subject to the constraint of the system dynamics. The resulting optimization problem is too large for any of the standard non-linear programming techniques and hence an hierarchical optimization approach is proposed. It turns out that the states can be computed at the first levelfor given noise and system parameters. These, in turn, are to be modified at the second level.The states are to be computed from a large system of linear equations and two solution methods are considered for solving these equations, limiting the horizon to a suitable length. The resulting algorithm is a filter-smoother, suitable for off-line as well as on-line state estimation for given noise and system parameters. The second level problem is split up into two, one for modifying the noise statistics and the other for modifying the system parameters. An adaptive relaxation technique is proposed for modifying the noise statistics and a modified Gauss-Newton technique is used to adjust the system parameters.

Optimal simultaneous maximum a posterioriestimation of states, noise statistics and parameters II. Numerical performance

A very general and numerically quite robust algorithm has been proposed by Sastry and Gauvrit (1980) for system identification. The present paper takes it up and examines its performance on a real test example. The example considered is the lateral dynamics of an aircraft. This is used as a vehicle for demonstrating the performance of various aspects of the algorithm in several possible modes.