Local power of some tests in exponential family nonlinear models

In this paper we obtain asymptotic expansions up to order n(-1/2) for the nonnull distribution functions of the likelihood ratio, Wald, score and gradient test statistics in exponential family nonlinear models (Cordeiro and Paula, 1989), under a sequence of Pitman alternatives. The asymptotic distributions of all four statistics are obtained for testing a subset of regression parameters and for testing the dispersion parameter, thus generalising the results given in Cordeiro et al. (1994) and Ferrari et al. (1997). We also present Monte Carlo simulations in order to compare the finite-sample performance of these tests. (C) 2010 Elsevier B.V. All rights reserved.

Genetic heterogeneity of within-family variance of body weight in Atlantic salmon (Salmo salar)

BACKGROUND: Canalization is defined as the stability of a genotype against minor variations in both environment and genetics. Genetic variation in degree of canalization causes heterogeneity of within-family variance. The aims of this study are twofold: (1) quantify genetic heterogeneity of (within-family) residual variance in Atlantic salmon and (2) test whether the observed heterogeneity of (within-family) residual variance can be explained by simple scaling effects. RESULTS: Analysis of body weight in Atlantic salmon using a double hierarchical generalized linear model (DHGLM) revealed substantial heterogeneity of within-family variance. The 95% prediction interval for within-family variance ranged from ~0.4 to 1.2 kg2, implying that the within-family variance of the most extreme high families is expected to be approximately three times larger than the extreme low families. For cross-sectional data, DHGLM with an animal mean sub-model resulted in severe bias, while a corresponding sire-dam model was appropriate. Heterogeneity of variance was not sensitive to Box-Cox transformations of phenotypes, which implies that heterogeneity of variance exists beyond what would be expected from simple scaling effects. CONCLUSIONS: Substantial heterogeneity of within-family variance was found for body weight in Atlantic salmon. A tendency towards higher variance with higher means (scaling effects) was observed, but heterogeneity of within-family variance existed beyond what could be explained by simple scaling effects. For cross-sectional data, using the animal mean sub-model in the DHGLM resulted in biased estimates of variance components, which differed substantially both from a standard linear mean animal model and a sire-dam DHGLM model. Although genetic differences in canalization were observed, selection for increased canalization is difficult, because there is limited individual information for the variance sub-model, especially when based on cross-sectional data. Furthermore, potential macro-environmental changes (diet, climatic region, etc.) may make genetic heterogeneity of variance a less stable trait over time and space.

Fitting spatial models in the R package: hglm

We present a new version of the hglm package for fittinghierarchical generalized linear models (HGLM) with spatially correlated random effects. A CAR family for conditional autoregressive random effects was implemented. Eigen decomposition of the matrix describing the spatial structure (e.g. the neighborhood matrix) was used to transform the CAR random effectsinto an independent, but heteroscedastic, gaussian random effect. A linear predictor is fitted for the random effect variance to estimate the parameters in the CAR model.This gives a computationally efficient algorithm for moderately sized problems (e.g. n<5000).

Comparison of treatment of incomplete abortion with misoprostol by physicians and midwives at district level in Uganda : a randomised controlled equivalence trial

BACKGROUND: Misoprostol is established for the treatment of incomplete abortion but has not been systematically assessed when provided by midwives at district level in a low-resource setting. We investigated the effectiveness and safety of midwives diagnosing and treating incomplete abortion with misoprostol, compared with physicians. METHODS: We did a multicentre randomised controlled equivalence trial at district level at six facilities in Uganda. Eligibility criteria were women with signs of incomplete abortion. We randomly allocated women with first-trimester incomplete abortion to clinical assessment and treatment with misoprostol either by a physician or a midwife. The randomisation (1:1) was done in blocks of 12 and was stratified for study site. Primary outcome was complete abortion not needing surgical intervention within 14-28 days after initial treatment. The study was not masked. Analysis of the primary outcome was done on the per-protocol population with a generalised linear-mixed effects model. The predefined equivalence range was -4% to 4%. The trial was registered at ClinicalTrials.gov, number NCT01844024. FINDINGS: From April 30, 2013, to July 21, 2014, 1108 women were assessed for eligibility. 1010 women were randomly assigned to each group (506 to midwife group and 504 to physician group). 955 women (472 in the midwife group and 483 in the physician group) were included in the per-protocol analysis. 452 (95·8%) of women in the midwife group had complete abortion and 467 (96·7%) in the physician group. The model-based risk difference for midwife versus physician group was -0·8% (95% CI -2·9 to 1·4), falling within the predefined equivalence range (-4% to 4%). The overall proportion of women with incomplete abortion was 3·8% (36/955), similarly distributed between the two groups (4·2% [20/472] in the midwife group, 3·3% [16/483] in the physician group). No serious adverse events were recorded. INTERPRETATION: Diagnosis and treatment of incomplete abortion with misoprostol by midwives is equally safe and effective as when provided by physicians, in a low-resource setting. Scaling up midwives' involvement in treatment of incomplete abortion with misoprostol at district level would increase access to safe post-abortion care. FUNDING: The Swedish Research Council, Karolinska Institutet, and Dalarna University.

Mortalidade infantil e condições socioeconômicas nas microrregiões do Nordeste Brasileiro

The study aims to answer the following question: what are the different profiles of infant mortality, according to demographic, socioeconomic, infrastructure and health care, for the micro-regions at the Northeast of Brazil? Thus, the main objective is to analyze the profiles or typologies associated mortality levels sociodemographic conditions of the micro-regions, in the year 2010. To this end, the databases of birth and death certificates of SIM and SINASC (DATASUS/MS), were taken from the 2010 population Census microdata and from SIDRA/IBGE. As a methodology, a weighted multiple linear regression model was used in the analysis in order to find the most significant variables in the explanation child mortality for the year 2010. Also a cluster analysis was performed, seeking evidence, initially, of homogeneous groups of micro-regions, from of the significant variables. The logit of the infant mortality rate was used as dependent variable, while variables such as demographic, socioeconomic, infrastructure and health care in the micro-regions were taken as the independent variables of the model. The Bayesian estimation technique was applied to the database of births and deaths, due to the inconvenient fact of underreporting and random fluctuations of small quantities in small areas. The techniques of Spatial Statistics were used to determine the spatial behavior of the distribution of rates from thematic maps. In conclusion, we used the method GoM (Grade of Membership), to find typologies of mortality, associated with the selected variables by micro-regions, in order to respond the main question of the study. The results points out to the formation of three profiles: Profile 1, high infant mortality and unfavorable social conditions; Profile 2, low infant mortality, with a median social conditions of life; and Profile 3, median and high infant mortality social conditions. With this classification, it was found that, out of 188 micro-regions, 20 (10%) fits the extreme profile 1, 59 (31.4%) was characterized in the extreme profile 2, 34 (18.1%) was characterized in the extreme profile 3 and only 9 (4.8%) was classified as amorphous profile. The other micro-regions framed up in the profiles mixed. Such profiles suggest the need for different interventions in terms of public policies aimed to reducing child mortality in the region

Estimativas de correlações genéticas entre escores visuais e características de carcaça medidas por ultrassonografia em bovinos Nelore utilizando modelos bayesianos linear-limiar

O objetivo neste estudo foi estimar as correlações genéticas entre escores visuais e características de carcaça medidas por ultrassonografia em bovinos da raça Nelore utilizando a estatística bayesiana por meio da Amostragem de Gibbs, sob modelo animal linear-limiar. Foram estudadas as características categóricas morfológicas de musculosidade, estrutura física, conformação e sacro, avaliadas aos 15 e 22 meses de idade. Para as características de carcaça, foram avaliadas as características área de olho-de-lombo, espessura de gordura subcutânea, espessura de gordura subcutânea na garupa e altura na garupa. Os escores visuais devem ser empregados como critérios de seleção para aumentar o progresso genético para a característica área de olhode-lombo e, consequentemente, melhorar o rendimento de carcaça. As estimativas de correlação genética obtidas para musculosidade com espessura de gordura subcutânea e espessura de gordura subcutânea na garupa indicaram que a seleção para musculosidade pode levar a animais com melhor acabamento de carcaça. A seleção para a estrutura física e conformação aos 15 e 22 meses de idade pode promover resposta correlacionada para o aumento da altura na garupa.

Análise das relações da curva de crescimento e eficiência produtiva de vacas da raça Holandesa

O objetivo neste trabalho foi analisar as relações entre parâmetros da curva de crescimento (peso à maturidade e taxa de maturação) e medidas da eficiência produtiva de 333 vacas da raça Holandesa nascidas de 1992 a 2002. Dados de idade ao primeiro parto, produção de leite e duração da primeira lactação, primeiro intervalo de partos e produção de leite por dia de intervalo de partos foram analisados por meio de um modelo linear generalizado com o efeito fixo de grupo contemporâneo (ano-estação de nascimento), os efeitos aleatórios de pai da vaca e erro e os efeitos lineares e quadráticos do peso à maturidade e da taxa de maturação. Nas análises de longevidade e de duração da vida útil, incluiu-se o efeito fixo de motivo de descarte da vaca. Houve efeitos linear e quadrático da taxa de maturação sobre a maioria das características, exceto o primeiro intervalo de partos e a produção de leite por dia de intervalo de partos. As estimativas das taxas de maturação ótimas variaram de 0,0896 a 0,1187 kg/kg de peso vivo/mês. O peso à maturidade influenciou a longevidade de forma linear e quadrática. A combinação ótima foi de 701 kg de peso à maturidade e 0,0934 kg/kg de peso vivo/mês de taxa de maturação. O peso à maturidade correlacionou-se de forma desfavorável com a duração da primeira lactação (-0,11). Coeficientes de correlação significativos foram obtidos entre taxa de maturação e idade ao primeiro parto (-0,22), produção de leite na primeira lactação (0,18), duração da primeira lactação (0,18) e longevidade (0,14), o que confirma a hipótese de que vacas com taxas de maturação mais rápidas apresentam maior eficiência produtiva.

Resposta à Seleção para Características de Desempenho em um Rebanho de Seleção de Suínos

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Modelo de programação linear para seleção de pulverizadores agrícolas de barras

A seleção de pulverizadores agrícolas que se adaptem às necessidades da propriedade, é um processo trabalhoso, sendo uma das etapas mais importantes dentro do processo produtivo. O objetivo do presente trabalho foi o de desenvolver e utilizar um modelo de programação linear para auxiliar na seleção de pulverizadores agrícolas de barras, baseado no menor custo horário do equipamento. Foram utilizadas as informações técnicas referentes a 20 modelos de pulverizadores disponíveis no mercado, sendo quatro autopropelidos, oito de arrasto e oito do tipo montado. A análise de sensibilidade dos componentes dos custos operacionais mostrou que as taxas de reparo e depreciação foram os fatores que mais interferiram na variação do custo horário do conjunto trator-pulverizador. O modelo matemático desenvolvido facilitou a realização da análise de sensibilidade que foi processada em um tempo muito pequeno.

Caracterização e modelagem dos sistemas de produção de caprinos leiteiros

The objectives of this research were characterizing the dairy goat production systems and model it using linear program. On the first step of this research, the model was developed using data from farms that was affiliated in the ACCOMIG/Caprileite, used a similar dairy goat production systems and have a partnership program with Universidade Federal de Minas Gerais . The data of research were from a structured questionnaire applied with farmers and monitoring of production systems during a guided visit on their farms. The results permitted identify that all farms were classified as a small and have a intensive production system. The average herd size had 63.75 dairy goats on lactation; it permits a production of 153, 38 kg of goat milk per day. It was observed that existing more than one channel of commercialization for the goat milk and their derivative products. The data obtained, on the first step of this research, was used to develop a linear program model. It was evaluated in two goat production systems, called P1 and P2. The results showed that the P1 system, with an annual birth and lactation during approximately 300 days was the best alternative for business. These results were compared with a mixed (beef and dairy) goat system in the semiarid region, which indicated merged with both systems. Therefore, to achieve profits and sustainability of the system, in all simulations it was necessary a minimum limit of funding of U.S. $ 10,000.00; this value permit earning of U.S. $ 792.00 per month and pay the investment within 5 years

Primary power distribution systems planning taking into account reliability, operation and expansion costs

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Estudo de simulação de modelos lineares mistos com distribuição normal contaminada no melhoramento genético animal

Objetivou-se com esse trabalho comparar estimativas de componentes de variâncias obtidas por meio de modelos lineares mistos Gaussianos e Robustos, via Amostrador de Gibbs, em dados simulados. Foram simulados 50 arquivos de dados com 1.000 animais cada um, distribuídos em cinco gerações, em dois níveis de efeito fixo e três valores fenotípicos distintos para uma característica hipotética, com diferentes níveis de contaminação. Exceto para os dados sem contaminação, quando os modelos foram iguais, o modelo Robusto apresentou melhores estimativas da variância residual. As estimativas de herdabilidade foram semelhantes em todos os modelos, mas as análises de regressão mostraram que os valores genéticos preditos com uso do modelo Robusto foram mais próximos dos valores genéticos verdadeiros. Esses resultados sugerem que o modelo linear normal contaminado oferece uma alternativa flexível para estimação robusta em melhoramento genético animal.

The effect of implant design on insertion torque and immediate micromotion

Objectives: To evaluate the effect of insertion torque on micromotion to a lateral force in three different implant designs. Material and methods: Thirty-six implants with identical thread design, but different cutting groove design were divided in three groups: (1) non-fluted (no cutting groove, solid screw-form); (2) fluted (901 cut at the apex, tap design); and (3) Blossomt (Patent pending) (non-fluted with engineered trimmed thread design). The implants were screwed into polyurethane foam blocks and the insertion torque was recorded after each turn of 901 by a digital torque gauge. Controlled lateral loads of 10N followed by increments of 5 up to 100N were sequentially applied by a digital force gauge on a titanium abutment. Statistical comparison was performed with two-way mixed model ANOVA that evaluated implant design group, linear effects of turns and displacement loads, and their interaction. Results: While insertion torque increased as a function of number of turns for each design, the slope and final values increased (Po0.001) progressively from the Blossomt to the fluted to the non-fluted design (M +/- standard deviation [SD] = 64.1 +/- 26.8, 139.4 +/- 17.2, and 205.23 +/- 24.3 Ncm, respectively). While a linear relationship between horizontal displacement and lateral force was observed for each design, the slope and maximal displacement increased (Po0.001) progressively from the Blossomt to the fluted to the non-fluted design (M +/- SD 530 +/- 57.7, 585.9 +/- 82.4, and 782.33 +/- 269.4 mm, respectively). There was negligible to moderate levels of association between insertion torque and lateral displacement in the Blossomt, fluted and non-fluted design groups, respectively. Conclusion: Insertion torque was reduced in implant macrodesigns that incorporated cutting edges, and lesser insertion torque was generally associated with decreased micromovement. However, insertion torque and micromotion were unrelated within implant designs, particularly for those designs showing the least insertion torque.

Choosing between Cox proportional hazards and logistic models for interval-censored data via bootstrap

This work develops a new methodology in order to discriminate models for interval-censored data based on bootstrap residual simulation by observing the deviance difference from one model in relation to another, according to Hinde (1992). Generally, this sort of data can generate a large number of tied observations and, in this case, survival time can be regarded as discrete. Therefore, the Cox proportional hazards model for grouped data (Prentice & Gloeckler, 1978) and the logistic model (Lawless, 1982) can befitted by means of generalized linear models. Whitehead (1989) considered censoring to be an indicative variable with a binomial distribution and fitted the Cox proportional hazards model using complementary log-log as a link function. In addition, a logistic model can be fitted using logit as a link function. The proposed methodology arises as an alternative to the score tests developed by Colosimo et al. (2000), where such models can be obtained for discrete binary data as particular cases from the Aranda-Ordaz distribution asymmetric family. These tests are thus developed with a basis on link functions to generate such a fit. The example that motivates this study was the dataset from an experiment carried out on a flax cultivar planted on four substrata susceptible to the pathogen Fusarium oxysoprum. The response variable, which is the time until blighting, was observed in intervals during 52 days. The results were compared with the model fit and the AIC values.

Body size and flight distance in stingless bees (Hymenoptera: Meliponini): inference of flight range and possible ecological implications

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)