Bases para o estudo da genética de populações dos Hymenoptera em geral e dos Apinae sociais em particular

This paper deals with problems on population genetics in Hymenoptera and particularly in social Apidae. 1) The studies on populations of Hymenoptera were made according to the two basic types of reproduction: endogamy and panmixia. The populations of social Apinae have a mixed method of reproduction with higher percentage of panmixia and a lower of endogamy. This is shown by the following a) males can enter any hive in swarming time; b) males of Meliponini are expelled from hives which does not need them, and thus, are forced to look for some other place; c) Meliponini males were seen powdering themselves with pollen, thus becoming more acceptable in any other hive. The panmixia is not complete owing to the fact that the density of the breeding population as very low, even in the more frequent species as low as about 2 females and 160 males per reproductive area. We adopted as selection values (or survival indices) the expressions according to Brieger (1948,1950) which may be summarised as follows; a population: p2AA + ²pq Aa + q2aa became after selection: x p2AA + 2pq Aa + z q²aa. For alge-braics facilities Brieger divided the three selective values by y giving thus: x/y p2 AA + y/y 2 pq Aa + z/y q²aa. He called x/y of RA and z/y of Ra, that are survival or selective index, calculated in relation to the heterozygote. In our case all index were calculated in relation to the heterozygote, including the ones for haploid males; thus we have: RA surveval index of genotype AA Ra surveval index of genotype aa R'A surveval index of genotype A R'a surveval index of genotype a 1 surveval index of genotype Aa The index R'A ande R'a were equalized to RA and Ra, respectively, for facilities in the conclusions. 2) Panmitic populations of Hymenoptera, barring mutations, migrations and selection, should follow the Hardy-Weinberg law, thus all gens will be present in the population in the inicial frequency (see Graphifc 1). 3) Heterotic genes: If mutation for heterotic gene ( 1 > RA > Ra) occurs, an equilibrium will be reached in a population when: P = R A + Ra - 2R²a _____________ (9) 2(R A + Ra - R²A - R²a q = R A + Ra - 2R²A _____________ (10) 2(R A + Ra - R²A - R²a A heterotic gene in an hymenopteran population may be maintained without the aid of new mutation only if the survival index of the most viable mutant (RA) does not exced the limiting value given by the formula: R A = 1 + √1+Ra _________ 4 If RA has a value higher thah the one permitted by the formula, then only the more viable gene will remain present in the population (see Graphic 10). The only direct proof for heterotic genes in Hymenoptera was given by Mackensen and Roberts, who obtained offspring from Apis mellefera L. queens fertilized by their own sons. Such inbreeding resulted in a rapid loss of vigor the colony; inbred lines intercrossed gave a high hybrid vigor. Other fats correlated with the "heterosis" problem are; a) In a colony M. quadrifasciata Lep., which suffered severely from heat, the percentage of deths omong males was greater .than among females; b) Casteel and Phillips had shown that in their samples (Apis melifera L). the males had 7 times more abnormalities tian the workers (see Quadros IV to VIII); c) just after emerging the males have great variation, but the older ones show a variation equal to that of workers; d) The tongue lenght of males of Apis mellifera L., of Bombus rubicundus Smith (Quadro X), of Melipona marginata Lep. (Quadro XI), and of Melipona quadrifasciata Lep. Quadro IX, show greater variationthan that of workers of the respective species. If such variation were only caused by subviables genes a rapid increasse of homozigoty for the most viable alleles should be expected; then, these .wild populations, supposed to be in equilibrium, could .not show such variability among males. Thus we conclude that heterotic genes have a grat importance in these cases. 4) By means of mathematical models, we came to the conclusion tht isolating genes (Ra ^ Ra > 1), even in the case of mutations with more adaptability, have only the opor-tunity of survival when the population number is very low (thus the frequency of the gene in the breeding population will be large just after its appearence). A pair of such alleles can only remain present in a population when in border regions of two races or subspecies. For more details see Graphics 5 to 8. 5) Sex-limited genes affecting only females, are of great importance toHymenoptera, being subject to the same limits and formulas as diploid panmitic populations (see formulas 12 and 13). The following examples of these genes were given: a) caste-determining genes in the genus Melipona; b) genes permiting an easy response of females to differences in feeding in almost all social Hymenoptera; c) two genes, found in wild populations, one in Trigona (Plebéia) mosquito F. SMITH (quadro XII) and other in Melipona marginata marginata LEP. (Quadro XIII, colonies 76 and 56) showing sex-limited effects. Sex-limited genes affecting only males do not contribute to the plasticity or genie reserve in hymenopteran populations (see formula 14). 6) The factor time (life span) in Hymenoptera has a particular importance for heterotic genes. Supposing one year to be the time unit and a pair of heterotic genes with respective survival indice equal to RA = 0, 90 and Ra = 0,70 to be present; then if the life time of a population is either one or two years, only the more viable gene will remain present (see formula 11). If the species has a life time of three years, then both alleles will be maintained. Thus we conclude that in specis with long lif-time, the heterotic genes have more importance, and should be found more easily. 7) The colonies of social Hymenoptera behave as units in competition, thus in the studies of populations one must determine the survival index, of these units which may be subdivided in indice for egg-laying, for adaptive value of the queen, for working capacity of workers, etc. 8) A study of endogamic hymenopteran populations, reproduced by sister x brother mating (fig. 2), lead us to the following conclusions: a) without selection, a population, heterozygous for one pair of alleles, will consist after some generations (theoretically after an infinite number of generation) of females AA fecundated with males A and females aa fecundated with males a (see Quadro I). b) Even in endogamic population there is the theoretical possibility of the presence of heterotic genes, at equilibrium without the aid of new mutations (see Graphics 11 and 12), but the following! conditions must be satisfied: I - surveval index of both homozygotes (RA e Ra) should be below 0,75 (see Graphic 13); II - The most viable allele must riot exced the less viable one by more than is permited by the following formula (Pimentel Gomes 1950) (see Gra-fic 14) : 4 R5A + 8 Ra R4A - 4 Ra R³A (Ra - 1) R²A - - R²a (4 R²a + 4 Ra - 1) R A + 2 R³a < o Considering these two conditions, the existance of heterotic genes in endogamic populations of Hymenoptera \>ecames very improbable though not - impossible. 9) Genie mutation offects more hymenopteran than diploid populations. Thus we have for lethal genes in diploid populations: u = q2, and in Hymenoptera: u = s, being u the mutation ratio and s the frequency of the mutant in the male population. 10) Three factors, important to competition among species of Meliponini were analysed: flying capacity of workers, food gathering capacity of workers, egg-laying of the queen. In this connection we refer to the variability of the tongue lenght observed in colonies from several localites, to the method of transporting the pollen in the stomach, from some pots (Melliponi-ni storage alveolus) to others (e. g. in cases of pillage), and to the observation that the species with the most populous hives are almost always the most frequent ones also. 11) Several defensive ways used for Meliponini to avoid predation are cited, but special references are made upon the camouflage of both hive (fig. 5) and hive entrance (fig. 4) and on the mimetism (see list in page ). Also under the same heading we described the method of Lestrimelitta for pillage. 12) As mechanisms important for promoting genetic plasticity of hymenopteran species we cited: a) cytological variations and b) genie reserve. As to the former, duplications and numerical variations of chromosomes were studied. Diprion simile ATC was cited as example for polyploidy. Apis mellife-ra L. (n •= 16) also sugests polyploid origen since: a) The genus Melipona, which belongs to a" related tribe, presents in all species so far studied n = 9 chromosomes and b) there occurs formation of dyads in the firt spermatocyte division. It is su-gested that the origin of the sex-chromosome of Apis mellifera It. may be related to the possible origin of diplo-tetraploidy in this species. With regards to the genie reserve, several possible types of mutants were discussed. They were classified according to their survival indices; the heterotic and neutral mutants must be considered as more important for the genie reserve. 13) The mean radius from a mother to a daghter colony was estimated as 100 meters. Since the Meliponini hives swarm only once a year we may take 100 meters a year as the average dispersion of female Meliponini in ocordance to data obtained from Trigona (tetragonisca) jaty F. SMITH and Melipona marginata LEP., while other species may give different values. For males the flying distance was roughly estimated to be 10 times that for females. A review of the bibliography on Meliponini swarm was made (pg. 43 to 47) and new facts added. The population desity (breeding population) corresponds in may species of Meliponini to one male and one female per 10.000 square meters. Apparently the males are more frequent than the females, because there are sometimes many thousands, of males in a swarm; but for the genie frequency the individuals which have descendants are the ones computed. In the case of Apini and Meliponini, only one queen per hive and the males represented by. the spermatozoos in its spermateca are computed. In Meliponini only one male mate with the queen, while queens of Apis mellijera L. are fecundated by an average of about 1, 5 males. (Roberts, 1944). From the date cited, one clearly sees that, on the whole, populations of wild social bees (Meliponini) are so small that the Sewall Wright effect may become of great importance. In fact applying the Wright's formula: f = ( 1/aN♂ + 1/aN♀) (1 - 1/aN♂ + 1/aN♀) which measures the fixation and loss of genes per generation, we see that the fixation or loss of genes is of about 7% in the more frequent species, and rarer species about 11%. The variation in size, tergite color, background color, etc, of Melipona marginata Lep. is atributed to this genetic drift. A detail, important to the survival of Meliponini species, is the Constance of their breeding population. This Constance is due to the social organization, i. e., to the care given to the reproductive individuals (the queen with its sperm pack), to the way of swarming, to the food storage intended to control variations of feeding supply, etc. 14) Some species of the Meliponini are adapted to various ecological conditions and inhabit large geographical areas (e. g. T. (Tetragonisca jaty F. SMITH), and Trigona (Nanno-trigona testaceicornis LEP.) while others are limited to narrow regions with special ecological conditions (e. g. M. fuscata me-lanoventer SCHWARZ). Other species still, within the same geographical region, profit different ecological conditions, as do M. marginata LEP. and M. quadrifasciata LEP. The geographical distribution of Melipona quadrifasciata LEP. is different according to the subspecies: a) subsp anthidio-des LEP. (represented in Fig. 7 by black squares) inhabits a region fron the North of the S. Paulo State to Northeastern Brazil, ,b) subspecies quadrifasciata LEP., (marked in Fig. 7 with black triangles) accurs from the South of S. Paulo State to the middle of the State of Rio Grande do Sul (South Brazil). In the margined region between these two areas of distribution, hi-brid colonies were found (Fig. 7, white circles); they are shown with more details in fig. 8, while the zone of hybridization is roughly indicated in fig. 9 (gray zone). The subspecies quadrifasciata LEP., has 4 complete yellow bands on the abdominal tergites while anthidioides LEP. has interrupted ones. This character is determined by one or two genes and gives different adaptative properties to the subspecies. Figs. 10 shows certains meteorological isoclines which have aproximately the same configuration as the limits of the hybrid zone, suggesting different climatic adaptabilities for both genotypes. The exis-tance of a border zone between the areas of both subspecies, where were found a high frequency of hybrids, is explained as follows: being each subspecies adapted to a special climatic zone, we may suppose a poor adaptation of either one in the border region, which is also a region of intermediate climatic conditions. Thus, the hybrids, having a combination of the parent qualities, will be best adapted to the transition zone. Thus, the hybrids will become heterotic and an equilibrium will be reached with all genotypes present in the population in the border region.

Study and implementation of some quantitative trading models

Quantitative or algorithmic trading is the automatization of investments decisions obeying a fixed or dynamic sets of rules to determine trading orders. It has increasingly made its way up to 70% of the trading volume of one of the biggest financial markets such as the New York Stock Exchange (NYSE). However, there is not a signi cant amount of academic literature devoted to it due to the private nature of investment banks and hedge funds. This projects aims to review the literature and discuss the models available in a subject that publications are scarce and infrequently. We review the basic and fundamental mathematical concepts needed for modeling financial markets such as: stochastic processes, stochastic integration and basic models for prices and spreads dynamics necessary for building quantitative strategies. We also contrast these models with real market data with minutely sampling frequency from the Dow Jones Industrial Average (DJIA). Quantitative strategies try to exploit two types of behavior: trend following or mean reversion. The former is grouped in the so-called technical models and the later in the so-called pairs trading. Technical models have been discarded by financial theoreticians but we show that they can be properly cast into a well defined scientific predictor if the signal generated by them pass the test of being a Markov time. That is, we can tell if the signal has occurred or not by examining the information up to the current time; or more technically, if the event is F_t-measurable. On the other hand the concept of pairs trading or market neutral strategy is fairly simple. However it can be cast in a variety of mathematical models ranging from a method based on a simple euclidean distance, in a co-integration framework or involving stochastic differential equations such as the well-known Ornstein-Uhlenbeck mean reversal ODE and its variations. A model for forecasting any economic or financial magnitude could be properly defined with scientific rigor but it could also lack of any economical value and be considered useless from a practical point of view. This is why this project could not be complete without a backtesting of the mentioned strategies. Conducting a useful and realistic backtesting is by no means a trivial exercise since the \laws" that govern financial markets are constantly evolving in time. This is the reason because we make emphasis in the calibration process of the strategies' parameters to adapt the given market conditions. We find out that the parameters from technical models are more volatile than their counterpart form market neutral strategies and calibration must be done in a high-frequency sampling manner to constantly track the currently market situation. As a whole, the goal of this project is to provide an overview of a quantitative approach to investment reviewing basic strategies and illustrating them by means of a back-testing with real financial market data. The sources of the data used in this project are Bloomberg for intraday time series and Yahoo! for daily prices. All numeric computations and graphics used and shown in this project were implemented in MATLAB^R scratch from scratch as a part of this thesis. No other mathematical or statistical software was used.

The MicroArray Quality Control (MAQC)-II study of common practices for the development and validation of microarray-based predictive models.

Gene expression data from microarrays are being applied to predict preclinical and clinical endpoints, but the reliability of these predictions has not been established. In the MAQC-II project, 36 independent teams analyzed six microarray data sets to generate predictive models for classifying a sample with respect to one of 13 endpoints indicative of lung or liver toxicity in rodents, or of breast cancer, multiple myeloma or neuroblastoma in humans. In total, >30,000 models were built using many combinations of analytical methods. The teams generated predictive models without knowing the biological meaning of some of the endpoints and, to mimic clinical reality, tested the models on data that had not been used for training. We found that model performance depended largely on the endpoint and team proficiency and that different approaches generated models of similar performance. The conclusions and recommendations from MAQC-II should be useful for regulatory agencies, study committees and independent investigators that evaluate methods for global gene expression analysis.

Epidemic models with an infected-infectious period

The introduction of an infective-infectious period on the geographic spread of epidemics is considered in two different models. The classical evolution equations arising in the literature are generalized and the existence of epidemic wave fronts is revised. The asymptotic speed is obtained and improves previous results for the Black Death plague

Formation of the long range Dpp morphogen gradient.

The TGF-β homolog Decapentaplegic (Dpp) acts as a secreted morphogen in the Drosophila wing disc, and spreads through the target tissue in order to form a long range concentration gradient. Despite extensive studies, the mechanism by which the Dpp gradient is formed remains controversial. Two opposing mechanisms have been proposed: receptor-mediated transcytosis (RMT) and restricted extracellular diffusion (RED). In these scenarios the receptor for Dpp plays different roles. In the RMT model it is essential for endocytosis, re-secretion, and thus transport of Dpp, whereas in the RED model it merely modulates Dpp distribution by binding it at the cell surface for internalization and subsequent degradation. Here we analyzed the effect of receptor mutant clones on the Dpp profile in quantitative mathematical models representing transport by either RMT or RED. We then, using novel genetic tools, experimentally monitored the actual Dpp gradient in wing discs containing receptor gain-of-function and loss-of-function clones. Gain-of-function clones reveal that Dpp binds in vivo strongly to the type I receptor Thick veins, but not to the type II receptor Punt. Importantly, results with the loss-of-function clones then refute the RMT model for Dpp gradient formation, while supporting the RED model in which the majority of Dpp is not bound to Thick veins. Together our results show that receptor-mediated transcytosis cannot account for Dpp gradient formation, and support restricted extracellular diffusion as the main mechanism for Dpp dispersal. The properties of this mechanism, in which only a minority of Dpp is receptor-bound, may facilitate long-range distribution.

Exactly solvable phase oscillator models with syncrhonization dynamics

Populations of phase oscillators interacting globally through a general coupling function f(x) have been considered. We analyze the conditions required to ensure the existence of a Lyapunov functional giving close expressions for it in terms of a generating function. We have also proposed a family of exactly solvable models with singular couplings showing that it is possible to map the synchronization phenomenon into other physical problems. In particular, the stationary solutions of the least singular coupling considered, f(x) = sgn(x), have been found analytically in terms of elliptic functions. This last case is one of the few nontrivial models for synchronization dynamics which can be analytically solved.

Dynamic models of viral replication and latency.

PURPOSE OF REVIEW: HIV targets primary CD4(+) T cells. The virus depends on the physiological state of its target cells for efficient replication, and, in turn, viral infection perturbs the cellular state significantly. Identifying the virus-host interactions that drive these dynamic changes is important for a better understanding of viral pathogenesis and persistence. The present review focuses on experimental and computational approaches to study the dynamics of viral replication and latency. RECENT FINDINGS: It was recently shown that only a fraction of the inducible latently infected reservoirs are successfully induced upon stimulation in ex-vivo models while additional rounds of stimulation make allowance for reactivation of more latently infected cells. This highlights the potential role of treatment duration and timing as important factors for successful reactivation of latently infected cells. The dynamics of HIV productive infection and latency have been investigated using transcriptome and proteome data. The cellular activation state has shown to be a major determinant of viral reactivation success. Mathematical models of latency have been used to explore the dynamics of the latent viral reservoir decay. SUMMARY: Timing is an important component of biological interactions. Temporal analyses covering aspects of viral life cycle are essential for gathering a comprehensive picture of HIV interaction with the host cell and untangling the complexity of latency. Understanding the dynamic changes tipping the balance between success and failure of HIV particle production might be key to eradicate the viral reservoir.

Application of theoretical models to financial innovation

The present paper is aimed at providing a general strategic overview of the existing theoretical models that have applications in the field of financial innovation. Whereas most financialdevelopments have relied upon traditional economic tools, a new stream of research is defining a novel paradigm in which mathematical models from diverse scientific disciplines are being applied to conceptualize and explain economic and financial behavior. Indeed, terms such as ‘econophysics’ or ‘quantum finance’ have recently appeared to embrace efforts in this direction. As a first contact with such research, the project will present a brief description of some of the main theoretical models that have applications in finance and economics, and will try to present, if possible, potential new applications to particular areas in financial analysis, or new applicable models. As a result, emphasiswill be put on the implications of this research for the financial sector and its future dynamics.

Monte Carlo Methods in Parameter Estimation of Nonlinear Models

Yksi keskeisimmistä tehtävistä matemaattisten mallien tilastollisessa analyysissä on mallien tuntemattomien parametrien estimointi. Tässä diplomityössä ollaan kiinnostuneita tuntemattomien parametrien jakaumista ja niiden muodostamiseen sopivista numeerisista menetelmistä, etenkin tapauksissa, joissa malli on epälineaarinen parametrien suhteen. Erilaisten numeeristen menetelmien osalta pääpaino on Markovin ketju Monte Carlo -menetelmissä (MCMC). Nämä laskentaintensiiviset menetelmät ovat viime aikoina kasvattaneet suosiotaan lähinnä kasvaneen laskentatehon vuoksi. Sekä Markovin ketjujen että Monte Carlo -simuloinnin teoriaa on esitelty työssä siinä määrin, että menetelmien toimivuus saadaan perusteltua. Viime aikoina kehitetyistä menetelmistä tarkastellaan etenkin adaptiivisia MCMC menetelmiä. Työn lähestymistapa on käytännönläheinen ja erilaisia MCMC -menetelmien toteutukseen liittyviä asioita korostetaan. Työn empiirisessä osuudessa tarkastellaan viiden esimerkkimallin tuntemattomien parametrien jakaumaa käyttäen hyväksi teoriaosassa esitettyjä menetelmiä. Mallit kuvaavat kemiallisia reaktioita ja kuvataan tavallisina differentiaaliyhtälöryhminä. Mallit on kerätty kemisteiltä Lappeenrannan teknillisestä yliopistosta ja Åbo Akademista, Turusta.

Steady-state global optimization of metabolic non-linear dynamic models through recasting into power-law canonical models

Background: Design of newly engineered microbial strains for biotechnological purposes would greatly benefit from the development of realistic mathematical models for the processes to be optimized. Such models can then be analyzed and, with the development and application of appropriate optimization techniques, one could identify the modifications that need to be made to the organism in order to achieve the desired biotechnological goal. As appropriate models to perform such an analysis are necessarily non-linear and typically non-convex, finding their global optimum is a challenging task. Canonical modeling techniques, such as Generalized Mass Action (GMA) models based on the power-law formalism, offer a possible solution to this problem because they have a mathematical structure that enables the development of specific algorithms for global optimization. Results: Based on the GMA canonical representation, we have developed in previous works a highly efficient optimization algorithm and a set of related strategies for understanding the evolution of adaptive responses in cellular metabolism. Here, we explore the possibility of recasting kinetic non-linear models into an equivalent GMA model, so that global optimization on the recast GMA model can be performed. With this technique, optimization is greatly facilitated and the results are transposable to the original non-linear problem. This procedure is straightforward for a particular class of non-linear models known as Saturable and Cooperative (SC) models that extend the power-law formalism to deal with saturation and cooperativity. Conclusions: Our results show that recasting non-linear kinetic models into GMA models is indeed an appropriate strategy that helps overcoming some of the numerical difficulties that arise during the global optimization task.

A measure of stability as a criterion for the verification and analysis of simulation models

The aim of this study is to define a new statistic, PVL, based on the relative distance between the likelihood associated with the simulation replications and the likelihood of the conceptual model. Our results coming from several simulation experiments of a clinical trial show that the PVL statistic range can be a good measure of stability to establish when a computational model verifies the underlying conceptual model. PVL improves also the analysis of simulation replications because only one statistic is associated with all the simulation replications. As well it presents several verification scenarios, obtained by altering the simulation model, that show the usefulness of PVL. Further simulation experiments suggest that a 0 to 20 % range may define adequate limits for the verification problem, if considered from the viewpoint of an equivalence test.

A comparison of electricity spot prices simulation using ARMA-GARCH and mean-reverting models

The aim of this work is to compare two families of mathematical models for their respective capability to capture the statistical properties of real electricity spot market time series. The first model family is ARMA-GARCH models and the second model family is mean-reverting Ornstein-Uhlenbeck models. These two models have been applied to two price series of Nordic Nord Pool spot market for electricity namely to the System prices and to the DenmarkW prices. The parameters of both models were calibrated from the real time series. After carrying out simulation with optimal models from both families we conclude that neither ARMA-GARCH models, nor conventional mean-reverting Ornstein-Uhlenbeck models, even when calibrated optimally with real electricity spot market price or return series, capture the statistical characteristics of the real series. But in the case of less spiky behavior (System prices), the mean-reverting Ornstein-Uhlenbeck model could be seen to partially succeeded in this task.

Models and applications for risk assessment and prediction of Asian soybean rust epidemics

Asian rust of soybean [Glycine max (L.) Merril] is one of the most important fungal diseases of this crop worldwide. The recent introduction of Phakopsora pachyrhizi Syd. & P. Syd in the Americas represents a major threat to soybean production in the main growing regions, and significant losses have already been reported. P. pachyrhizi is extremely aggressive under favorable weather conditions, causing rapid plant defoliation. Epidemiological studies, under both controlled and natural environmental conditions, have been done for several decades with the aim of elucidating factors that affect the disease cycle as a basis for disease modeling. The recent spread of Asian soybean rust to major production regions in the world has promoted new development, testing and application of mathematical models to assess the risk and predict the disease. These efforts have included the integration of new data, epidemiological knowledge, statistical methods, and advances in computer simulation to develop models and systems with different spatial and temporal scales, objectives and audience. In this review, we present a comprehensive discussion on the models and systems that have been tested to predict and assess the risk of Asian soybean rust. Limitations, uncertainties and challenges for modelers are also discussed.

Mathematical model to estimate of the deterioration of wooden poles in contact with soil used in rural areas

In São Paulo State, mainly in rural areas, the utilization of wooden poles is observed for different purposes. In this context, wood in contact with the ground presents faster deterioration, which is generally associated to environmental factors and, especially to the presence of fungi and insects. With the use of mathematical models, the useful life of wooden structures can be predicted by obtaining "climatic indexes" to indicate, comparatively among the areas studied, which have more or less tendency to fungi and insects attacks. In this work, by using climatological data of several cities at São Paulo State, a simplified mathematical model was obtained to measure the aggressiveness of the wood in contact with the soil.

Adjust of regression models to estimate the rectal temperature of broilers for the first 14 days of life

The broiler rectal temperature (t rectal) is one of the most important physiological responses to classify the animal thermal comfort. Therefore, the aim of this study was to adjust regression models in order to predict the rectal temperature (t rectal) of broiler chickens under different thermal conditions based on age (A) and a meteorological variable (air temperature - t air) or a thermal comfort index (temperature and humidity index -THI or black globe humidity index - BGHI) or a physical quantity enthalpy (H). In addition, through the inversion of these models and the expected t rectal intervals for each age, the comfort limits of t air, THI, BGHI and H for the chicks in the heating phase were determined, aiding in the validation of the equations and the preliminary limits for H. The experimental data used to adjust the mathematical models were collected in two commercial poultry farms, with Cobb chicks, from 1 to 14 days of age. It was possible to predict the t rectal of conditions from the expected t rectal and determine the lower and superior comfort thresholds of broilers satisfactorily by applying the four models adjusted; as well as to invert the models for prediction of the environmental H for the chicks first 14 days of life.