Estimating Two-Stage Models for Genetic Influences on Alcohol, Tobacco or Drug Use Initiation and Dependence Vulnerability in Twin and Family Data

Genetic research on risk of alcohol, tobacco or drug dependence must make allowance for the partial overlap of risk-factors for initiation of use, and risk-factors for dependence or other outcomes in users. Except in the extreme cases where genetic and environmental risk-factors for initiation and dependence overlap completely or are uncorrelated, there is no consensus about how best to estimate the magnitude of genetic or environmental correlations between Initiation and Dependence in twin and family data. We explore by computer simulation the biases to estimates of genetic and environmental parameters caused by model misspecification when Initiation can only be defined as a binary variable. For plausible simulated parameter values, the two-stage genetic models that we consider yield estimates of genetic and environmental variances for Dependence that, although biased, are not very discrepant from the true values. However, estimates of genetic (or environmental) correlations between Initiation and Dependence may be seriously biased, and may differ markedly under different two-stage models. Such estimates may have little credibility unless external data favor selection of one particular model. These problems can be avoided if Initiation can be assessed as a multiple-category variable (e.g. never versus early-onset versus later onset user), with at least two categories measurable in users at risk for dependence. Under these conditions, under certain distributional assumptions., recovery of simulated genetic and environmental correlations becomes possible, Illustrative application of the model to Australian twin data on smoking confirmed substantial heritability of smoking persistence (42%) with minimal overlap with genetic influences on initiation.

Genetic consequences of sequential founder events by an island-colonizing bird

The importance of founder events in promoting evolutionary changes on islands has been a subject of long-running controversy. Resolution of this debate has been hindered by a lack of empirical evidence from naturally founded island populations. Here we undertake a genetic analysis of a series of historically documented, natural colonization events by the silvereye species-complex (Zosterops lateralis), a group used to illustrate the process of island colonization in the original founder effect model. Our results indicate that single founder events do not affect levels of heterozygosity or allelic diversity, nor do they result in immediate genetic differentiation between populations. Instead, four to five successive founder events are required before indices of diversity and divergence approach that seen in evolutionarily old forms. A Bayesian analysis based on computer simulation allows inferences to be made on the number of effective founders and indicates that founder effects are weak because island populations are established from relatively large flocks. Indeed, statistical support for a founder event model was not significantly higher than for a gradual-drift model for all recently colonized islands. Taken together, these results suggest that single colonization events in this species complex are rarely accompanied by severe founder effects, and multiple founder events and/or long-term genetic drift have been of greater consequence for neutral genetic diversity.

A methodology for building energy modelling and calibration in warm climates

In the last 7 years, a method has been developed to analyse building energy performance using computer simulation, in Brazil. The method combines analysis of building design plans and documentation, walk-through visits, electric and thermal measurements and the use of an energy simulation tool (DOE-2.1E code), The method was used to model more than 15 office buildings (more than 200 000 m(2)), located between 12.5degrees and 27.5degrees South latitude. The paper describes the basic methodology, with data for one building and presents additional results for other six cases. (C) 2002 Elsevier Science Ltd. All rights reserved.

Subsídio ao aprimoramento do manejo da irrigação de consórcios agroflorestais em situação de escassez hídrica

A demanda por água do setor agrícola vem crescendo a cada ano, bem como o aperfeiçoamento do manejo da irrigação. Ainda assim, em locais onde a água encontra-se escassa, os conflitos entre usos e usuários devido a disputas para acesso à água tendem a se intensificar. Na maioria das vezes esse problema está relacionado à ausência de cobertura florestal, levando, por vezes, a problemas socioambientais. Nesse contexto, a pesquisa buscou avaliar o impacto da implantação de sistemas agroflorestais (SAF’s) irrigados na condição hídrica local (de escassez), estabelecendo-se estratégias que os torne atrativos ambientalsocial- e economicamente. Assim, por meio de simulação computacional, alternativas de consórcios agroflorestais foram avaliadas, tomando-se como referência uma região piloto, típica da agricultura espírito-santense, constituída por pequenas propriedades agrícolas de base familiar, inseridas parcialmente em APPs. As espécies banana, pupunha e goiaba foram selecionadas para compor os cenários agroflorestais. Os resultados da pesquisa mostram que os SAF’s irrigados são uma alternativa no sentido de minimizar os conflitos por demandas de água em regiões de escassez hídrica, com a redução de tais demandas o comprometimento de rendimentos financeiros. Além disso, os SAF’s são alternativas para diversificação da renda e para tentar controlar a sazonalidade dos preços de mercado.

Comportamento dinâmico da rede eléctrica da ilha de São Vicente em Cabo Verde: estabilidade de tensão

Neste trabalho é efectuado, não só o diagnóstico em regime permanente, mas também o estudo, simulação e análise do comportamento dinâmico da rede eléctrica da ilha de São Vicente em Cabo Verde. Os estudos de estabilidade transitória desempenham um importante papel, tanto no planeamento como na operação dos sistemas de potência. Tais estudos são realizados, em grande parte, através de simulação digital no domínio do tempo, utilizando integração numérica para resolver as equações não-lineares que modelam a dinâmica do sistema e dependem da existência de registos reais de perturbação (ex: osciloperturbografia). O objectivo do trabalho será também verificar a aplicabilidade dos requisitos técnicos que as unidades geradoras devem ter, no que concerne ao controlo de tensão, estabelecidos na futura regulamentação europeia desenvolvida pela ENTSO-E (European Network Transmission System Operator for Electricity). De entre os requisitos analisou-se a capacidade das máquinas existentes suportarem cavas de tensão decorrentes de curto-circuitos trifásicos simétricos, Fault Ride Through, no ponto de ligação à rede. Identificaram-se para o efeito os factores que influenciam a estabilidade desta rede, em regime perturbado nomeadamente: (i) duração do defeito, (ii) caracterização da carga, com e sem a presença do sistema de controlo de tensão (AVR) em unidades de geração síncronas. Na ausência de registos reais sobre o comportamento do sistema, conclui-se que este é sensível à elasticidade das cargas em particular do tipo potência constante, existindo risco de perda de estabilidade, neste caso, para defeitos superiores a 5ms sem AVR. A existência de AVR nesta rede afigura-se como indispensável para garantir estabilidade de tensão sendo contudo necessário proceder a uma correcta parametrização.

Monte Carlo modeling and simulations of the High Definition (HD120) micro MLC and validation against measurements for a 6 MV beam

Purpose: The most recent Varian® micro multileaf collimator(MLC), the High Definition (HD120) MLC, was modeled using the BEAMNRCMonte Carlo code. This model was incorporated into a Varian medical linear accelerator, for a 6 MV beam, in static and dynamic mode. The model was validated by comparing simulated profiles with measurements. Methods: The Varian® Trilogy® (2300C/D) accelerator model was accurately implemented using the state-of-the-art Monte Carlo simulation program BEAMNRC and validated against off-axis and depth dose profiles measured using ionization chambers, by adjusting the energy and the full width at half maximum (FWHM) of the initial electron beam. The HD120 MLC was modeled by developing a new BEAMNRC component module (CM), designated HDMLC, adapting the available DYNVMLC CM and incorporating the specific characteristics of this new micro MLC. The leaf dimensions were provided by the manufacturer. The geometry was visualized by tracing particles through the CM and recording their position when a leaf boundary is crossed. The leaf material density and abutting air gap between leaves were adjusted in order to obtain a good agreement between the simulated leakage profiles and EBT2 film measurements performed in a solid water phantom. To validate the HDMLC implementation, additional MLC static patterns were also simulated and compared to additional measurements. Furthermore, the ability to simulate dynamic MLC fields was implemented in the HDMLC CM. The simulation results of these fields were compared with EBT2 film measurements performed in a solid water phantom. Results: Overall, the discrepancies, with and without MLC, between the opened field simulations and the measurements using ionization chambers in a water phantom, for the off-axis profiles are below 2% and in depth-dose profiles are below 2% after the maximum dose depth and below 4% in the build-up region. On the conditions of these simulations, this tungsten-based MLC has a density of 18.7 g cm− 3 and an overall leakage of about 1.1 ± 0.03%. The discrepancies between the film measured and simulated closed and blocked fields are below 2% and 8%, respectively. Other measurements were performed for alternated leaf patterns and the agreement is satisfactory (to within 4%). The dynamic mode for this MLC was implemented and the discrepancies between film measurements and simulations are within 4%. Conclusions: The Varian® Trilogy® (2300 C/D) linear accelerator including the HD120 MLC was successfully modeled and simulated using the Monte CarloBEAMNRC code by developing an independent CM, the HDMLC CM, either in static and dynamic modes.

FMRI 3D registration based on Fourier space subsets using neural networks

In this work, we present a neural network (NN) based method designed for 3D rigid-body registration of FMRI time series, which relies on a limited number of Fourier coefficients of the images to be aligned. These coefficients, which are comprised in a small cubic neighborhood located at the first octant of a 3D Fourier space (including the DC component), are then fed into six NN during the learning stage. Each NN yields the estimates of a registration parameter. The proposed method was assessed for 3D rigid-body transformations, using DC neighborhoods of different sizes. The mean absolute registration errors are of approximately 0.030 mm in translations and 0.030 deg in rotations, for the typical motion amplitudes encountered in FMRI studies. The construction of the training set and the learning stage are fast requiring, respectively, 90 s and 1 to 12 s, depending on the number of input and hidden units of the NN. We believe that NN-based approaches to the problem of FMRI registration can be of great interest in the future. For instance, NN relying on limited K-space data (possibly in navigation echoes) can be a valid solution to the problem of prospective (in frame) FMRI registration.

Mutation accumulation in Tetrahymena

Background - The rate and fitness effects of mutations are key in understanding the evolution of every species. Traditionally, these parameters are estimated in mutation accumulation experiments where replicate lines are propagated in conditions that allow mutations to randomly accumulate without the purging effect of natural selection. These experiments have been performed with many model organisms but we still lack empirical estimates of the rate and effects of mutation in the protists. Results - We performed a mutation accumulation (MA) experiment in Tetrahymena thermophila, a species that can reproduce sexually and asexually in nature, and measured both the mean decline and variance increase in fitness of 20 lines. The results obtained with T. thermophila were compared with T. pyriformis that is an obligate asexual species. We show that MA lines of T. thermophila go to extinction at a rate of 1.25 clonal extinctions per bottleneck. In contrast, populations of T. pyriformis show a much higher resistance to extinction. Variation in gene copy number is likely to be a key factor in explaining these results, and indeed we show that T. pyriformis has a higher mean copy number per cell than T. thermophila. From fitness measurements during the MA experiment, we infer a rate of mutation to copy number variation of 0.0333 per haploid MAC genome of T. thermophila and a mean effect against copy number variation of 0.16. A strong effect of population size in the rate of fitness decline was also found, consistent with the increased power of natural selection. Conclusions - The rate of clonal extinction measured for T. thermophila is characteristic of a mutational degradation and suggests that this species must undergo sexual reproduction to avoid the deleterious effects detected in the laboratory experiments. We also suggest that an increase in chromosomal copy number associated with the phenotypic assortment of amitotic divisions can provide an alternative mechanism to escape the deleterious effect of random chromosomal copy number variation in species like T. pyriformis that lack the resetting mechanism of sexual reproduction. Our results are relevant to the understanding of cell line longevity and senescence in ciliates.

Besnoitia besnoiti protein disulfide isomerase (BbPDI): molecular characterization, expression and in silico modelling

Besnoitia besnoiti is an apicomplexan parasite responsible for bovine besnoitiosis, a disease with a high prevalence in tropical and subtropical regions and re-emerging in Europe. Despite the great economical losses associated with besnoitiosis, this disease has been underestimated and poorly studied, and neither an effective therapy nor an efficacious vaccine is available. Protein disulfide isomerase (PDI) is an essential enzyme for the acquisition of the correct three-dimensional structure of proteins. Current evidence suggests that in Neosporacaninum and Toxoplasmagondii, which are closely related to B. besnoiti, PDI play an important role in host cell invasion, is a relevant target for the host immune response, and represents a promising drug target and/or vaccine candidate. In this work, we present the nucleotide sequence of the B. besnoiti PDI gene. BbPDI belongs to the thioredoxin-like superfamily (cluster 00388) and is included in the PDI_a family (cluster defined cd02961) and the PDI_a_PDI_a'_c subfamily (cd02995). A 3D theoretical model was built by comparative homology using Swiss-Model server, using as a template the crystallographic deduced model of Tapasin-ERp57 (PDB code 3F8U chain C). Analysis of the phylogenetic tree for PDI within the phylum apicomplexa reinforces the close relationship among B. besnoiti, N. caninum and T. gondii. When subjected to a PDI-assay based on the polymerisation of reduced insulin, recombinant BbPDI expressed in E. coli exhibited enzymatic activity, which was inhibited by bacitracin. Antiserum directed against recombinant BbPDI reacted with PDI in Western blots and by immunofluorescence with B. besnoiti tachyzoites and bradyzoites.

Projecto suspensão formula student

Trabalho Final de Mestrado para obtenção do grau de Mestre em Engenharia Mecânica

Controlo clássico vs. preditivo na regulação de frequência em sistemas de potência

Dissertação para obtenção do grau de Mestre em Engenharia Electrotécnica Ramo Automação e Electrónica Industrial

Avaliação e optimização operacional: estudo do terminal de contentores de Alcântara

Trabalho Final de Mestrado para obtenção do grau de Mestre em Engenharia Mecânica - Ramo Manutenção e Produção

Leilão de direitos de transmissão financeiros de energia eléctrica em ambiente de mercado

Os direitos de transmissão (Transmission Rigths – TRs) correspondem, na sua essência, a contratos que conferem aos seus proprietários o direito de transmitir energia eléctrica, por um determinado caminho, a um preço fixo. Na actualidade vive-se uma era de liberalização dos mercados de energia eléctrica nos quais, no caso concreto dos mercados de TRs, os diversos agentes podem ceder os TRs que possuem a outros agentes desde que cumpram certos requisitos impostos pelo sistema na figura do seu operador de sistema (Independent System Operator – ISO). Neste sentido, o ISO oferece, aos diversos agentes do mercado, algumas ferramentas que lhes permite transaccionar, sob sua orientação, os seus respectivos TRs fazendo-se cumprir todos os requisitos indispensáveis para o efeito. A mais popular dessas ferramentas, nos principais mercados energético da actualidade, é o leilão. Com o presente trabalho de dissertação é pretendido apresentar-se um modelo para a resolução do problema inerente a um leilão de TRs em ambiente de mercado, neste caso concreto de direitos de transmissão financeiros de energia eléctrica (Financial Transmission Rigths – FTRs). Neste sentido foi desenvolvido um simulador informático (SIM_AuctFTR) que implementa um modelo para este tipo de problemas. Este trabalho foi estruturado essencialmente em três etapas com objectivos inerentes. Assim, numa primeira fase da realização deste trabalho, foi realizado um estudo de diversos conceitos e metodologias inerentes ao problema de leilão de FTRs em ambiente de mercado, suportado por uma sólida base bibliográfica. A segunda fase tratou-se do desenvolvimento do algoritmo da aplicação computacional que solucione um problema de um leilão de FTRs. O SIM_AuctFTR foi desenvolvido sobre o pressuposto da maximização dos rendimentos financeiros provenientes da atribuição dos FTRs propostos a leilão, tendo estes últimos de coabitar no sistema respeitando as restrições técnicas a que este são inerentes, mesmo na eventualidade da ocorrência de contingências de nível ‘n-1’. Por último, numa terceira fase, partindo da aplicação desenvolvida, foram realizados 4 casos de estudos com outras tantas redes eléctricas de forma a testar a robustez da ferramenta desenvolvida.

Wind-CSP short-term coordination by MILP approach

This paper is on the maximization of total profit in a day-ahead market for a price-taker producer needing a short-term scheduling for wind power plants coordination with concentrated solar power plants, having thermal energy storage systems. The optimization approach proposed for the maximization of profit is a mixed-integer linear programming problem. The approach considers not only transmission grid constraints, but also technical operating constraints on both wind and concentrated solar power plants. Then, an improved short-term scheduling coordination is provided due to the more accurate modelling presented in this paper. Computer simulation results based on data for the Iberian wind and concentrated solar power plants illustrate the coordination benefits and show the effectiveness of the approach.