Efeito da exposição às temperaturas extremas na população portuguesa - um contributo da epidemiologia

Backgroud: O International Panel on Climate Change prevê que o aumento da temperatura média global, até ao ano de 2100, varie entre 1,4 e 5,8ºC desconhecendo-se a evolução da adaptação da população a esta subida da temperatura. Em Portugal morre-se mais no Inverno que no Verão. Mas existem evidências de repercussões na mortalidade atribuíveis ao calor extremo. Este estudo procura conhecer os grupos etários e/ou populacionais que parecem revelar vulnerabilidade acrescida à exposição a temperaturas extremas e identificar indicadores de saúde apropriados para revelar esses mesmos efeitos. Métodos: Foram analisados dados de internamentos hospitalar e mortalidade por doenças cardiovasculares, respiratórias, renais, efeitos directos do frio e do calor, na população com 75 e mais anos de idade, nos distritos de Beja, Bragança e Faro, nos meses de Janeiro e Junho. Para os dados de morbilidade o período de análise foi 2002 a 2005 e para os de mortalidade de 2002 a 2004. Os dados meteorológicos analisados corresponderam aos valores da temperatura máxima e percentis da temperatura máxima, nos meses de Janeiro (P10) e Junho (P90). Os excessos de internamentos hospitalares, definidos como os dias em que ocorreram internamentos acima do valor da média mais 2 desvio padrão, foram relacionados com a distribuição das temperaturas extremas (frias abaixo do P10, quentes acima do P90.Os dias com óbitos acima do valor da média foram relacionados com a distribuição das temperaturas extremas (frias abaixo do P10, quentes acima do P90). Os indicadores propostos foram baseados em Odds Ratios e intervalos de confiança que sugeriam as estimativas mais precisas. Resultados: O grupo que revelou maior vulnerabilidade às temperaturas extremas foi o grupo dos 75 e mais anos, com doenças cardiovasculares quando exposto a temperaturas extremas, nos 3 distritos observados.O nº de dias de excesso de óbitos por doenças cardiovasculares relacionados com temperaturas extremas foi o mais elevado comparado com as restantes causas de morte. O grupo etário dos 75 e mais anos com de doenças respiratórias também é vulnerável, às temperaturas extremas frias, nos 3 distritos. Verificaram-se dias de excessos de internamentos hospitalares e óbitos por esta causa de morte, relacionados com a exposição às temperaturas extremas frias. Em Junho, não se verificou excesso de mortalidade associado à exposição a temperaturas extremas por esta causa, em qualquer dos distritos analisados. Apenas se verificou a associação entre os dias de ocorrência de internamentos hospitalares por doenças renais e o calor extremo, em Bragança. Conclusões: Foram encontradas associações estatísticas significativas entre dias de excesso de ocorrência de internamentos hospitalares ou óbitos por causa e exposição a temperaturas extremas frias e quentes possibilitando a identificação de um conjunto de indicadores de saúde ambiental apropriados para monitorizar a evolução dos padrões de morbilidade, mortalidade e susceptibilidade das populações ao longo do tempo.-------------------- Backgroud: International Panel on Climate Change estimates that the rise of mean global temperature varies between 1,4 e 5,8ºC until 2100, with unknowing evolution adaptation of populations. In Portugal we die more in Winter than in Summer time. But there are several evidences of mortality attributable to extreme eat. The proposal of this study is to know the age and/or populations groups that reveal more vulnerability to exposure to extreme temperature and identifying proper health indicators to reveal those effects. Methods: Data from hospital admissions and mortality caused by cardiovascular, respiratory, renal diseases and direct effects from direct exposure to extreme cold and heat, in population with 75 and more years, in Beja, Bragança and Faro districts, during January and June, were analysed. Analysis period for morbidity data was from 2002 to 2005 and form mortality was 2002 to 2004. Meteorological data analysed were maximum temperature and percentile of maximum temperature, from January (P10) and June (P90. Relationship between excess of hospital admission, defined as the days that occurred hospital admissions above mean value more 2 standards desviation and distribution of extreme temperatures were established (cold under P10 and heat above P90. Proposal indicators were based on Odds Ratios and confidence intervals, suggesting the most precises estimatives. Results: The most vulnerable group to extreme temperature were people with 75 or more years older with cardiovascular diseases, observed in the 3 districts. Number of days caused by excess cardiovascular mortality and extreme temperature were the most number of days between the other causes. The group with 75 or more years old with respiratory diseases is vulnerable too, especially to cold extreme temperature, in all the 3 districts. There were excess of days of hospital admissions and days with deaths, for this cause relating to extreme cold temperature. In June, does not funded excess of mortality associated to extreme temperature by this cause in any district of the in observation. Just was found relationship between days of hospital admissions caused by renal diseases in Bragança in days with extreme heat. Conclusions: Were found statistically significant associations between days of excess of hospital admissions or deaths and exposure to extreme cold and heat temperatures giving the possibility of identifying a core of environmental indicators proper to monitoring patterns and trends evolutions on morbidity, mortality and susceptibly of populations for a long time.

Can red cell distribution width be used as a marker of Crohn's disease activity?

Introduction: Recently, it has been suggested an association between red cell distribution width (RDW) and Crohn’s disease activity index (CDAI), but its use is not yet performed in daily clinical practice. Objectives: To determine whether RDW can be used as a marker of Crohn’s disease (CD) activity. Methods: This was a cross-sectional study including patients with CD, observed consecutively in an outpatient setting between January 1st and September 30th 2013. Blood cell indices, erythrocyte sedimentation rate (ESR), and C-reactive protein were measured. CD activity was determined by CDAI (active disease if CDAI ≥ 150). Associations were analyzed using logistic regression (SPSS version 20). Results: 119 patients (56% female) were included in the study with a mean age of 47 years (SD 15.2). Twenty patients (17%) had active disease. The median RDW was 14.0 (13---15). There was an association between RDW and disease activity (p = 0.044). After adjustment for age and gender, this association remained consistent (OR 1.20, 95% CI 1.03---1.39, p = 0.016). It was also found that the association between RDW and disease activity was independent of hemoglobin and ESR (OR 1.36, 95% CI 1.08---1.72, p = 0.01) and of biologic therapy (OR 1.19, 95% CI 1.03---1.37, p = 0.017). A RDW cutoff of 16% had a specificity and negative predictive value for CDAI ≥ 150 of 88% and 86%, respectively. Conclusion: In this study, RDW proved to be an independent and relatively specific marker of CD activity. These results may contribute to the implementation of this simple parameter, in clinical practice, aiming to help therapeutic decisions.

Factorization of elliptic boundary value problems by invariant embedding and application to overdetermined problems

Dissertação para obtenção do Grau de Doutor em Matemática

Optimal Cut-Off Value for Homeostasis Model Assessment (HOMA) Index of Insulin-Resistance in a Population of Patients Admitted Electively in a Portuguese Cardiology Ward

INTRODUCTION: Insulin resistance is the pathophysiological key to explain metabolic syndrome. Although clearly useful, the Homeostasis Model Assessment index (an insulin resistance measurement) hasn't been systematically applied in clinical practice. One of the main reasons is the discrepancy in cut-off values reported in different populations. We sought to evaluate in a Portuguese population the ideal cut-off for Homeostasis Model Assessment index and assess its relationship with metabolic syndrome. MATERIAL AND METHODS: We selected a cohort of individuals admitted electively in a Cardiology ward with a BMI < 25 Kg/m2 and no abnormalities in glucose metabolism (fasting plasma glucose < 100 mg/dL and no diabetes). The 90th percentile of the Homeostasis Model Assessment index distribution was used to obtain the ideal cut-off for insulin resistance. We also selected a validation cohort of 300 individuals (no exclusion criteria applied). RESULTS: From 7 000 individuals, and after the exclusion criteria, there were left 1 784 individuals. The 90th percentile for Homeostasis Model Assessment index was 2.33. In the validation cohort, applying that cut-off, we have 49.3% of individuals with insulin resistance. However, only 69.9% of the metabolic syndrome patients had insulin resistance according to that cut-off. By ROC curve analysis, the ideal cut-off for metabolic syndrome is 2.41. Homeostasis Model Assessment index correlated with BMI (r = 0.371, p < 0.001) and is an independent predictor of the presence of metabolic syndrome (OR 19.4, 95% CI 6.6 - 57.2, p < 0.001). DISCUSSION: Our study showed that in a Portuguese population of patients admitted electively in a Cardiology ward, 2.33 is the Homeostasis Model Assessment index cut-off for insulin resistance and 2.41 for metabolic syndrome. CONCLUSION: Homeostasis Model Assessment index is directly correlated with BMI and is an independent predictor of metabolic syndrome.

Spatial distribution and esterase activity in populations of Aedes (Stegomyia) aegypti (Linnaeus) (Diptera: Culicidae) resistant to temephos

INTRODUCTION:The need for studies that describe the resistance patterns in populations of Aedes aegypti (Linnaeus) in function of their region of origin justified this research, which aimed to characterize the resistance to temephos and to obtain information on esterase activity in populations of Aedes aegypticollected in municipalities of the State of Paraíba.METHODS:Resistance to temephos was evaluated and characterized from the diagnostic dose of 0.352mg i.a./L and multiple concentrations that caused mortalities between 5% and 99%. Electrophoresis of isoenzymes was used to verify the patterns of esterase activity among populations of the vector.RESULTS:All populations of Aedes aegypti were resistant to temephos, presenting a resistance rate (RR) greater than 20. The greatest lethal dose 50% of the sample (CL50) was found for the municipality of Lagoa Seca, approximately forty-one times the value of CL50 for the Rockefeller population. The populations characterized as resistant showed two to six regions of α and β-esterase, called EST-1 to EST-6, while the susceptible population was only seen in one region of activity.CONCLUSIONS:Aedes aegyptiis widely distributed and shows a high degree of resistance to temephos in all municipalities studied. In all cases, esterases are involved in the metabolism and, consequently, in the resistance to temephos.

Urinary tract infection in full-term newborn infants: value of urine culture by bag specimen collection

OBJETIVE: to evaluate the efficacy of urine culture by bag specimen for the detection of neonatal urinary tract infection in full-term newborn infants. Retrospective study (1997) including full-term newborn infants having a positive urine culture (>100,000 CFU/ml) by bag specimen collection. The urinary tract infection diagnosis was confirmed by positive urine culture (suprapubic bladder aspiration method). The select cases were divided into three groups, according to newborn infant age at the bag specimen collection: GI (< 48 h, n = 17), GII (48 h to 7 d, n = 35) and GIII (> 7 d, n = 9). Sixty one full-term newborn infants were studied (5.1 % of total infants). The diagnosis was confirmed on 19/61 (31.1 %) of full-term infants born alive. Distribution among the groups was: GI = 2/17 (11.8 %), GII = 10//35 (28.6 %), and GIII = 7/9 (77.7 %). The most relevant clinical symptoms were: fever (GI - 100 %, GII - 91.4 %) and weight loss (GI - 35.3 %, GII - 45.7 %). Urine culture results for specimens collected by suprapubic aspiration were: E. coli GI (100 %), GII (40 %) and GIII (28.6 %), E. faecalis GI (30%), Staphylococcus coagulase-negative GII (20 %) and GIII (42.8 %), and Staphylococcus aureus GII (10 %). Correlation between positive urine culture collection (bag specimen method) and urinary tract infection diagnosis, using relative risk analysis, produced the following results: GI=0.30 (CI95% 0.08-1.15), GII=0.51 (CI 95% 0.25-1.06) and GIII=3.31 (CI95% 1.8-6.06) The most frequent urinary tract infection clinical signs in the first week were fever and weight loss, while non-specific symptomatology occurred later. E. coli was most frequent infectious agent, although from the 7th day of life, staphylococcus was noted. The urine culture (bag specimen method) was effective in detecting urinary tract infection only after the 7th day of life.

Cross-scale determinants of palm species distribution

We investigate palm species distribution, richness and abundance along the Mokoti, a seasonally-dry river of southeastern Amazon and compare it to the patterns observed at a large scale, comprising the entire Brazilian territory. A total of 694 palms belonging to 10 species were sampled at the Mokoti River basin. Although the species showed diverse distribution patterns, we found that local palm abundance, richness and tree basal area were significantly higher from the hills to the bottomlands of the study region, revealing a positive association of these measures with moisture. The analyses at the larger spatial scale also showed a strong influence of vapor pressure (a measure of moisture content of the air, in turn modulated by temperature) and seasonality in temperature: the richest regions were those where temperature and humidity were simultaneously high, and which also presented a lower degree of seasonality in temperature. These results indicate that the distribution of palms seems to be strongly associated with climatic variables, supporting the idea that, by 'putting all the eggs in one basket' (a consequence of survival depending on the preservation of a single irreplaceable bud), palms have become vulnerable to extreme environmental conditions. Hence, their distribution is concentrated in those tropical and sub-tropical regions with constant conditions of (mild to high) temperature and moisture all year round.

Systematics, functional morphology and distribution of a bivalve (Apachecorbula muriatica gen. et sp. nov.) from the rim of the ‘Valdivia Deep’ brine pool in the Red Sea

The deep brine pools of the Red Sea comprise extreme, inhospitable habitats yet house microbial communities that potentially may fuel adjacent fauna. We here describe a novel bivalve from a deep-sea (1525 m) brine pool in the Red Sea, where conditions of high salinity, lowered pH, partial anoxia and high temperatures are prevalent. Remotely operated vehicle (ROV) footage showed that the bivalves were present in a narrow (20 cm) band along the rim of the brine pool, suggesting that it is not only tolerant of such extreme conditions but is also limited to them. The bivalve is a member of the Corbulidae and named Apachecorbula muriatica gen. et sp. nov. The shell is atypical of the family in being modioliform and thin. The semi-infaunal habit is seen in ROV images and reflected in the anatomy by the lack of siphons. The ctenidia are large and typical of a suspension feeding bivalve, but the absence of guard cilia and the greatly reduced labial palps suggest that it is non-selective as a response to low food availability. It is proposed that the low body mass observed is a consequence of the extreme habitat and low food availability. It is postulated that the observed morphology of Apachecorbula is a result of paedomorphosis driven by the effects of the extreme environment on growth but is in part mitigated by the absence of high predation pressures.

Pulse wave velocity distribution in a cohort study: from arterial stiffness to early vascular aging

BACKGROUND: By contrast with other southern European people, north Portuguese population registers an especially high prevalence of hypertension and stroke incidence. We designed a cohort study to identify individuals presenting accelerated and premature arterial aging in the Portuguese population. METHOD: Pulse wave velocity (PWV) was measured in randomly sampled population dwellers aged 18-96 years from northern Portugal, and used as a marker of early vascular aging (EVA). Of the 3038 individuals enrolled, 2542 completed the evaluation. RESULTS: Mean PWV value for the entire population was 8.4?m/s (men: 8.6?m/s; women: 8.2?m/s; P?values adjusted for age (using normal European reference values as comparators). The overall prevalence of EVA was 12.5%; 26.1% of individuals below 30 years presented this feature and 40.2% of individuals in that same age strata were placed above the 90th percentile of PWV; and 18.7% of the population exhibited PWV values above 10?m/s, with male predominance (17.2% of men aged 40-49 years had PWV?>?10?m/s). Logistic regression models indicated gender differences concerning the risk of developing large artery damage, with women having the same odds of PWV above 10?m/s 10 years later than men. CONCLUSION: The population PWV values were higher than expected in a low cardiovascular risk area (Portugal). High prevalence rates of EVA and noteworthy large artery damage in young ages were found.

Influence of emphysema distribution on pulmonary function parameters in COPD patients

Objective: To evaluate the impact that the distribution of emphysema has on clinical and functional severity in patients with COPD. Methods: The distribution of the emphysema was analyzed in COPD patients, who were classified according to a 5-point visual classification system of lung CT findings. We assessed the influence of emphysema distribution type on the clinical and functional presentation of COPD. We also evaluated hypoxemia after the six-minute walk test (6MWT) and determined the six-minute walk distance (6MWD). Results: Eighty-six patients were included. The mean age was 65.2 ± 12.2 years, 91.9% were male, and all but one were smokers (mean smoking history, 62.7 ± 38.4 pack-years). The emphysema distribution was categorized as obviously upper lung-predominant (type 1), in 36.0% of the patients; slightly upper lung-predominant (type 2), in 25.6%; homogeneous between the upper and lower lung (type 3), in 16.3%; and slightly lower lung-predominant (type 4), in 22.1%. Type 2 emphysema distribution was associated with lower FEV1 , FVC, FEV1 /FVC ratio, and DLCO. In comparison with the type 1 patients, the type 4 patients were more likely to have an FEV1 < 65% of the predicted value (OR = 6.91, 95% CI: 1.43-33.45; p = 0.016), a 6MWD < 350 m (OR = 6.36, 95% CI: 1.26-32.18; p = 0.025), and post-6MWT hypoxemia (OR = 32.66, 95% CI: 3.26-326.84; p = 0.003). The type 3 patients had a higher RV/TLC ratio, although the difference was not significant. Conclusions: The severity of COPD appears to be greater in type 4 patients, and type 3 patients tend to have greater hyperinflation. The distribution of emphysema could have a major impact on functional parameters and should be considered in the evaluation of COPD patients.

Defining habitat characteristics influencing the distribution, density and growth of Plaice (Pleuronectes Platessa) and Dab (Limanda Limanda) on west of Ireland nursery grounds

Plaice (Pleuronectes platessa, L.) and dab (Limanda limanda, L.) are among the most abundant flatfishes in the north-eastern Atlantic region and the dominant species in shallow coastal nursery grounds. With increasing pressures on commercial flatfish stocks in combination with changing coastal environments, better knowledge of population dynamics during all life stages is needed to evaluate variability in year-class strength and recruitment to the fishery. The aim of this research was to investigate the complex interplay of biotic and abiotic habitat components influencing the distribution, density and growth of plaice and dab during the vulnerable juvenile life stage and to gain insight in spatial and temporal differences in nursery habitat quality along the west coast of Ireland. Intraspecific variability in plaice diet was observed at different spatial scales and showed a link with condition, recent growth and morphology. This highlights the effect of food availability on habitat quality and the need to consider small scale variation when attempting to link habitat quality to feeding, growth and condition of juvenile flatfish. There was evidence of trophic, spatial and temporal resource partitioning between juvenile plaice and dab allowing the co-existence of morphologically similar species in nursery grounds. In the limited survey years there was no evidence that the carrying capacity of the studied nursery grounds was reached but spatial and interannual variations in fish growth indicated fluctuating environments in terms of food availability, predator densities, sediment features and physico-chemical conditions. Predation was the most important factor affecting habitat quality for juvenile plaice and dab with crab densities negatively correlated to fish condition whereas shrimp densities were negatively associated with densities of small-sized juveniles in spring. A comparison of proxies for fish growth showed the advantage of Fulton’s K for routine use whereas RNA:DNA ratios proved less powerful when short-term environmental fluctuations are lacking. This study illustrated how distinct sets of habitat features can drive spatial variation in density and condition of juvenile flatfish highlighting the value of studying both variables when modeling habitat requirements. The habitat models generated in this study also provide a powerful tool to predict potential climate and anthropogenic impacts on the distribution and condition of juveniles in flatfish nurseries. The need for effective coastal zone management was emphasized to ensure a sustainable use of coastal resources and successful flatfish recruitment to the fishery.

Population dynamics of 0-group flounder (Platichthys flesus L.) in Galway Bay, West of Ireland and the value of flounder from shore angling tourism in Ireland

The economic value of flounder from shore angling around Ireland was assessed. Flounder catches from shore angling tournaments around Ireland were related to domestic and overseas shore angling expenditure in order to determine an economic value for the species. Temporal trends in flounder angling catches, and specimen (trophy) flounder reports were also investigated. Flounder was found to be the most caught shore angling species in competitions around Ireland constituting roughly one third of the shore angling competition catch although this did vary by area. The total value of flounder from shore angling tourism was estimated to be of the order of €8.4 million. No significant temporal trends in flounder angling catches and specimen reports were found. Thus there is no evidence from the current study for any decline in flounder stocks. The population dynamics of 0-group flounder during the early benthic stage was investigated at estuarine sites in Galway Bay, west of Ireland. Information was analysed from the March to June sampling period over five years (2002 to 2006). Spatial and temporal variations in settlement and population length structure were analysed between beach and river habitats and sites. Settlement of flounder began from late March to early May of each year, most commonly in April. Peak settlement was usually in April or early May. Settlement was recorded earlier than elsewhere, although most commonly was similar to the southern part of the UK and northern France. Settlement was generally later in tidal rivers than on sandy beaches. Abundance of 0-group flounder in Galway Bay did not exhibit significant inter -annual variability. 0-group flounder were observed in dense aggregations of up to 105 m'2, which were patchy in distribution. Highest densities of 0-group flounder were recorded in limnetic and oligohaline areas as compared with the lower densities in polyhaline and to a lesser extent mesohaline areas. Measurements to of salinity allowed the classification of beaches, and tidal river sections near the mouth, into a salinity based scheme for length comparisons. Beaches were classified as polyhaline,the lower section of rivers as mesohaline, and the middle and upper sections as oligohaline. Over the March to June sampling period 0-group flounder utilised different sections at different length ranges and were significantly larger in more upstream sections. During initial settlement in April, 0-group flounder of 8-10 mm (standard length, SL) were present in abundance on polyhaline sandy beaches. By about 10mm (SL), flounder were present in all polyhaline, mesohaline and (oligohaline) sections. 0-group flounder became absent or in insignificant numbers in polyhaline and mesohaline sections in a matter of weeks after first appearance. From April to June, 0-group flounder of 12-30mm (SL) were found in more upstream locations in the oligohaline sections. About one month (May or June) after initial settlement, 0-group flounder became absent from the oligohaline sections. Concurrently, flounder start to reappear in mesohaline and polyhaline areas at approximately 30mm (SL) in June. The results indicate 0-group flounder in the early benthic stage are associated with low salinity areas, but as they grow, this association diminishes. Results strongly suggest that migration of 0-group flounder between habitats takes place during the early benthic phase.

Prognostic Value of TIMI Score versus GRACE Score in ST-segment Elevation Myocardial Infarction

Background: The TIMI Score for ST-segment elevation myocardial infarction (STEMI) was created and validated specifically for this clinical scenario, while the GRACE score is generic to any type of acute coronary syndrome. Objective: Between TIMI and GRACE scores, identify the one of better prognostic performance in patients with STEMI. Methods: We included 152 individuals consecutively admitted for STEMI. The TIMI and GRACE scores were tested for their discriminatory ability (C-statistics) and calibration (Hosmer-Lemeshow) in relation to hospital death. Results: The TIMI score showed equal distribution of patients in the ranges of low, intermediate and high risk (39 %, 27 % and 34 %, respectively), as opposed to the GRACE Score that showed predominant distribution at low risk (80 %, 13 % and 7%, respectively). Case-fatality was 11%. The C-statistics of the TIMI score was 0.87 (95%CI = 0.76 to 0.98), similar to GRACE (0.87, 95%CI = 0.75 to 0.99) - p = 0.71. The TIMI score showed satisfactory calibration represented by χ2 = 1.4 (p = 0.92), well above the calibration of the GRACE score, which showed χ2 = 14 (p = 0.08). This calibration is reflected in the expected incidence ranges for low, intermediate and high risk, according to the TIMI score (0 %, 4.9 % and 25 %, respectively), differently to GRACE (2.4%, 25% and 73%), which featured middle range incidence inappropriately. Conclusion: Although the scores show similar discriminatory capacity for hospital death, the TIMI score had better calibration than GRACE. These findings need to be validated populations of different risk profiles.

A determinação dos números de indivíduos mínimos necessários na experimentação genética

The main object of the present paper consists in giving formulas and methods which enable us to determine the minimum number of repetitions or of individuals necessary to garantee some extent the success of an experiment. The theoretical basis of all processes consists essentially in the following. Knowing the frequency of the desired p and of the non desired ovents q we may calculate the frequency of all possi- ble combinations, to be expected in n repetitions, by expanding the binomium (p-+q)n. Determining which of these combinations we want to avoid we calculate their total frequency, selecting the value of the exponent n of the binomium in such a way that this total frequency is equal or smaller than the accepted limit of precision n/pª{ 1/n1 (q/p)n + 1/(n-1)| (q/p)n-1 + 1/ 2!(n-2)| (q/p)n-2 + 1/3(n-3) (q/p)n-3... < Plim - -(1b) There does not exist an absolute limit of precision since its value depends not only upon psychological factors in our judgement, but is at the same sime a function of the number of repetitions For this reasen y have proposed (1,56) two relative values, one equal to 1-5n as the lowest value of probability and the other equal to 1-10n as the highest value of improbability, leaving between them what may be called the "region of doubt However these formulas cannot be applied in our case since this number n is just the unknown quantity. Thus we have to use, instead of the more exact values of these two formulas, the conventional limits of P.lim equal to 0,05 (Precision 5%), equal to 0,01 (Precision 1%, and to 0,001 (Precision P, 1%). The binominal formula as explained above (cf. formula 1, pg. 85), however is of rather limited applicability owing to the excessive calculus necessary, and we have thus to procure approximations as substitutes. We may use, without loss of precision, the following approximations: a) The normal or Gaussean distribution when the expected frequency p has any value between 0,1 and 0,9, and when n is at least superior to ten. b) The Poisson distribution when the expected frequecy p is smaller than 0,1. Tables V to VII show for some special cases that these approximations are very satisfactory. The praticai solution of the following problems, stated in the introduction can now be given: A) What is the minimum number of repititions necessary in order to avoid that any one of a treatments, varieties etc. may be accidentally always the best, on the best and second best, or the first, second, and third best or finally one of the n beat treatments, varieties etc. Using the first term of the binomium, we have the following equation for n: n = log Riim / log (m:) = log Riim / log.m - log a --------------(5) B) What is the minimun number of individuals necessary in 01der that a ceratin type, expected with the frequency p, may appaer at least in one, two, three or a=m+1 individuals. 1) For p between 0,1 and 0,9 and using the Gaussean approximation we have: on - ó. p (1-p) n - a -1.m b= δ. 1-p /p e c = m/p } -------------------(7) n = b + b² + 4 c/ 2 n´ = 1/p n cor = n + n' ---------- (8) We have to use the correction n' when p has a value between 0,25 and 0,75. The greek letters delta represents in the present esse the unilateral limits of the Gaussean distribution for the three conventional limits of precision : 1,64; 2,33; and 3,09 respectively. h we are only interested in having at least one individual, and m becomes equal to zero, the formula reduces to : c= m/p o para a = 1 a = { b + b²}² = b² = δ2 1- p /p }-----------------(9) n = 1/p n (cor) = n + n´ 2) If p is smaller than 0,1 we may use table 1 in order to find the mean m of a Poisson distribution and determine. n = m: p C) Which is the minimun number of individuals necessary for distinguishing two frequencies p1 and p2? 1) When pl and p2 are values between 0,1 and 0,9 we have: n = { δ p1 ( 1-pi) + p2) / p2 (1 - p2) n= 1/p1-p2 }------------ (13) n (cor) We have again to use the unilateral limits of the Gaussean distribution. The correction n' should be used if at least one of the valors pl or p2 has a value between 0,25 and 0,75. A more complicated formula may be used in cases where whe want to increase the precision : n (p1 - p2) δ { p1 (1- p2 ) / n= m δ = δ p1 ( 1 - p1) + p2 ( 1 - p2) c= m / p1 - p2 n = { b2 + 4 4 c }2 }--------- (14) n = 1/ p1 - p2 2) When both pl and p2 are smaller than 0,1 we determine the quocient (pl-r-p2) and procure the corresponding number m2 of a Poisson distribution in table 2. The value n is found by the equation : n = mg /p2 ------------- (15) D) What is the minimun number necessary for distinguishing three or more frequencies, p2 p1 p3. If the frequecies pl p2 p3 are values between 0,1 e 0,9 we have to solve the individual equations and sue the higest value of n thus determined : n 1.2 = {δ p1 (1 - p1) / p1 - p2 }² = Fiim n 1.2 = { δ p1 ( 1 - p1) + p1 ( 1 - p1) }² } -- (16) Delta represents now the bilateral limits of the : Gaussean distrioution : 1,96-2,58-3,29. 2) No table was prepared for the relatively rare cases of a comparison of threes or more frequencies below 0,1 and in such cases extremely high numbers would be required. E) A process is given which serves to solve two problemr of informatory nature : a) if a special type appears in n individuals with a frequency p(obs), what may be the corresponding ideal value of p(esp), or; b) if we study samples of n in diviuals and expect a certain type with a frequency p(esp) what may be the extreme limits of p(obs) in individual farmlies ? I.) If we are dealing with values between 0,1 and 0,9 we may use table 3. To solve the first question we select the respective horizontal line for p(obs) and determine which column corresponds to our value of n and find the respective value of p(esp) by interpolating between columns. In order to solve the second problem we start with the respective column for p(esp) and find the horizontal line for the given value of n either diretly or by approximation and by interpolation. 2) For frequencies smaller than 0,1 we have to use table 4 and transform the fractions p(esp) and p(obs) in numbers of Poisson series by multiplication with n. Tn order to solve the first broblem, we verify in which line the lower Poisson limit is equal to m(obs) and transform the corresponding value of m into frequecy p(esp) by dividing through n. The observed frequency may thus be a chance deviate of any value between 0,0... and the values given by dividing the value of m in the table by n. In the second case we transform first the expectation p(esp) into a value of m and procure in the horizontal line, corresponding to m(esp) the extreme values om m which than must be transformed, by dividing through n into values of p(obs). F) Partial and progressive tests may be recomended in all cases where there is lack of material or where the loss of time is less importent than the cost of large scale experiments since in many cases the minimun number necessary to garantee the results within the limits of precision is rather large. One should not forget that the minimun number really represents at the same time a maximun number, necessary only if one takes into consideration essentially the disfavorable variations, but smaller numbers may frequently already satisfactory results. For instance, by definition, we know that a frequecy of p means that we expect one individual in every total o(f1-p). If there were no chance variations, this number (1- p) will be suficient. and if there were favorable variations a smaller number still may yield one individual of the desired type. r.nus trusting to luck, one may start the experiment with numbers, smaller than the minimun calculated according to the formulas given above, and increase the total untill the desired result is obtained and this may well b ebefore the "minimum number" is reached. Some concrete examples of this partial or progressive procedure are given from our genetical experiments with maize.

The ghost of past species occurrence: improving species distribution models for presence-only data

1. Model-based approaches have been used increasingly in conservation biology over recent years. Species presence data used for predictive species distribution modelling are abundant in natural history collections, whereas reliable absence data are sparse, most notably for vagrant species such as butterflies and snakes. As predictive methods such as generalized linear models (GLM) require absence data, various strategies have been proposed to select pseudo-absence data. However, only a few studies exist that compare different approaches to generating these pseudo-absence data. 2. Natural history collection data are usually available for long periods of time (decades or even centuries), thus allowing historical considerations. However, this historical dimension has rarely been assessed in studies of species distribution, although there is great potential for understanding current patterns, i.e. the past is the key to the present. 3. We used GLM to model the distributions of three 'target' butterfly species, Melitaea didyma, Coenonympha tullia and Maculinea teleius, in Switzerland. We developed and compared four strategies for defining pools of pseudo-absence data and applied them to natural history collection data from the last 10, 30 and 100 years. Pools included: (i) sites without target species records; (ii) sites where butterfly species other than the target species were present; (iii) sites without butterfly species but with habitat characteristics similar to those required by the target species; and (iv) a combination of the second and third strategies. Models were evaluated and compared by the total deviance explained, the maximized Kappa and the area under the curve (AUC). 4. Among the four strategies, model performance was best for strategy 3. Contrary to expectations, strategy 2 resulted in even lower model performance compared with models with pseudo-absence data simulated totally at random (strategy 1). 5. Independent of the strategy model, performance was enhanced when sites with historical species presence data were not considered as pseudo-absence data. Therefore, the combination of strategy 3 with species records from the last 100 years achieved the highest model performance. 6. Synthesis and applications. The protection of suitable habitat for species survival or reintroduction in rapidly changing landscapes is a high priority among conservationists. Model-based approaches offer planning authorities the possibility of delimiting priority areas for species detection or habitat protection. The performance of these models can be enhanced by fitting them with pseudo-absence data relying on large archives of natural history collection species presence data rather than using randomly sampled pseudo-absence data.