Survival analysis: a tool in the study of post-harvest diseases in peaches

Survival analysis is applied when the time until the occurrence of an event is of interest. Such data are routinely collected in plant diseases, although applications of the method are uncommon. The objective of this study was to use two studies on post-harvest diseases of peaches, considering two harvests together and the existence of random effect shared by fruits of a same tree, in order to describe the main techniques in survival analysis. The nonparametric Kaplan-Meier method, the log-rank test and the semi-parametric Cox's proportional hazards model were used to estimate the effect of cultivars and the number of days after full bloom on the survival to the brown rot symptom and the instantaneous risk of expressing it in two consecutive harvests. The joint analysis with baseline effect, varying between harvests, and the confirmation of the tree effect as a grouping factor with random effect were appropriate to interpret the phenomenon (disease) evaluated and can be important tools to replace or complement the conventional analysis, respecting the nature of the variable and the phenomenon.

Liver morphology with emphasis on bile ducts changes and survival analysis in mice submitted to multiple Schistosoma mansoni infections and chemotherapy

In an attempt to be as close as possible to the infected and treated patients of the endemic areas of schistosomiasis (S. mansoni) and in order to achieve a long period of follow-up, mice were repeatedly infected with a low number of cercariae. Survival data and histological variables such as schistosomal granuloma, portal changes, hepatocellular necrosis, hepatocellular regeneration, schistosomotic pigment, periductal fibrosis and chiefly bile ducts changes were analysed in the infected treated and non treated mice. Oxamniquine chemotherapy in repeatedly infected mice prolonged survival significantly when compared to non-treated animals (chi-square 9.24, p = 0.0024), thus confirming previous results with a similar experimental model but with a shorter term follow-up. Furthermore, mortality decreased rapidly after treatment suggesting an abrupt reduction in the severity of hepatic lesions. A morphological and immunohistochemical study of the liver was carried out. Portal fibrosis, with a pattern resembling human Symmers fibrosis was present at a late phase in the infected animals. Bile duct lesions were quite close to those described in human Mansonian schistosomiasis. Schistosomal antigen was observed in one isolated altered bile duct cell. The pathogenesis of the bile duct changes and its relation to the parasite infection and/or their antigens are discussed.

Mathematical models to describe the volumetric shrinkage rate of red beans during drying

The present study aimed to determine the volumetric shrinkage rate of bean (Phaseolus vulgaris L.) seeds during air-drying under different conditions of air, temperature and relative humidity, and to adjust several mathematical models to the empiric values observed, and select the one that best represents the phenomenon. Six mathematical models were adjusted to the experimental values to represent the phenomenon. It was determined the degree of adjustment of each model from the value of the coefficient of determination, the behavior of the distribution of the residuals, and the magnitude of the average relative and estimated errors. The rate of volumetric shrinkage that occurred in bean seeds during drying is between 25 and 37%. It basically depends on the final moisture content, regardless of the air conditions during drying. The Modified Bala & Woods' model best represented the process.

Survival analysis of female dogs with mammary tumors after mastectomy: epidemiological, clinical and morphological aspects

Abstract: Mammary gland tumors are the most common type of tumors in bitches but research on survival time after diagnosis is scarce. The purpose of this study was to investigate the relationship between survival time after mastectomy and a number of clinical and morphological variables. Data was collected retrospectively on bitches with mammary tumors seen at the Small Animal Surgery Clinic Service at the University of Brasília. All subjects had undergone mastectomy. Survival analysis was conducted using Cox's proportional hazard method. Of the 139 subjects analyzed, 68 died and 71 survived until the end of the study (64 months). Mean age was 11.76 years (SD=2.71), 53.84% were small dogs. 76.92% of the tumors were malignant, and 65.73% had both thoracic and inguinal glands affected. Survival time in months was associated with age (hazard rate ratios [HRR] =1.23, p-value =1.4x10-4), animal size (HRR between giant and small animals =2.61, p-value =0.02), nodule size (HRR =1.09, p-value =0.03), histological type (HRR between solid carcinoma and carcinoma in a mixed tumor =2.40, p-value =0.02), time between diagnosis and surgery (TDS, with HRR =1.21, p-value =2.7x10-15), and the interaction TDS*follow-up time (HRR =0.98, p-value =1.6x10-11). The present study is one of the few on the subject matter. Several important covariates were evaluated and age, animal size, nodule size, histological type, TDS and TDS*follow up time were identified as significantly associated to survival time.

Application of mathematical models for the prediction of adsorption isotherms in solid mixture for mango powder refreshment

Solid mixtures for refreshment are already totally integrated to the Brazilian consumers' daily routine, because of their quick preparation method, yield and reasonable price - quite lower if compared to 'ready-to-drink' products or products for prompt consumption, what makes them economically more accessible to low-income populations. Within such a context, the aim of this work was to evaluate the physicochemical and mineral composition, as well as the hygroscopic behavior of four different brands of solid mixture for mango refreshment. The BET, GAB, Oswim and Henderson mathematical models were built through the adjustment of experimental data to the isotherms of adsorption. Results from the physiochemical evaluation showed that the solid mixtures for refreshments are considerable sources of ascorbic acid and reductor sugar; and regarding mineral compounds, they are significant sources of calcium, sodium and potassium. It was also verified that the solid mixtures for refreshments of the four studied brands are considered highly hygroscopic.

A mathematical model for malaria transmission relating global warming and local socioeconomic conditions

OBJECTIVE: Sensitivity analysis was applied to a mathematical model describing malaria transmission relating global warming and local socioeconomic conditions. METHODS: A previous compartment model was proposed to describe the overall transmission of malaria. This model was built up on several parameters and the prevalence of malaria in a community was characterized by the values assigned to them. To assess the control efforts, the model parameters can vary on broad intervals. RESULTS: By performing the sensitivity analysis on equilibrium points, which represent the level of malaria infection in a community, the different possible scenarios are obtained when the parameters are changed. CONCLUSIONS: Depending on malaria risk, the efforts to control its transmission can be guided by a subset of parameters used in the mathematical model.

Draft prediction models for soil engaging tines in two soils of Rio Grande do Sul, Brazil

The draft forces of soil engaging tines and theoretical analysis compared to existing mathematical models, have yet not been studied in Rio Grande do Sul soils. From the existing models, those which can get the closest fitting draft forces to real measure on field have been established for two of Rio Grande do Sul soils. An Albaqualf and a Paleudult were evaluated. From the studied models, those suggested by Reece, so called "Universal Earthmoving Equation", Hettiaratchi and Reece, and Godwin and Spoor were the best fitting ones, comparing the calculated results with those measured "in situ". Allowing for the less complexity of Reece's model, it is suggested that this model should be used for modeling draft forces prediction for narrow tines in Albaqualf and Paleudut.

Survival of adult AIDS patients in a reference hospital of a metropolitan area in Brazil

OBJECTIVE: To evaluate the influence of sociodemographic, clinical, and epidemiological factors in AIDS patients survival in a reference hospital. METHODS: A sample of 502 adult AIDS patients out of 1,494 AIDS cases registered in a hospital in Fortaleza, Brazil, was investigated between 1986 and 1998. Sixteen cases were excluded due to death at the moment of the AIDS diagnosis and 486 were analyzed in the study. Socioeconomic and clinical epidemiological were the variables studied. Statistical analysis was conducted using the Kaplan-Meier survival analysis and the Cox proportional hazards model. RESULTS: Three hundred and sixty two out of the 486 patients studied took at least one antiretroviral drug and their survival was ten times longer than those who did not take any drug (746 and 79 days, respectively, p <0.001). Patients who took two nucleoside reverse transcriptase inhibitors (NRTI) plus protease inhibitor were found to have higher survival rates (p <0.001). The risk of dying in the first year was significantly lower for patients who took NRTI and a protease inhibitor compared to those who took only NRTI. In addition, this risk was much lower from the second year on (0.10; 95%CI: 0.42-0.23). The risk of dying in the first year was significantly higher for less educated patients (15.58; 95%CI: 6.64-36.58) and those who had two or more systemic diseases (3.03; 95%CI: 1.74-5.25). After the first year post-diagnosis, there was no risk difference for these factors. CONCLUSIONS: Higher education revealed to exert a significant influence in the first-year survival. Antiretroviral drugs had a greater impact in the survival from the second year on. A more aggressive antiretroviral therapy started earlier could benefit those patients.

Dengue and the risk of urban yellow fever reintroduction in São Paulo State, Brazil

OBJECTIVE: To propose a mathematical method for the estimation of the Basic Reproduction Number, R0, of urban yellow fever in a dengue-infested area. METHODS: The method is based on the assumption that, as the same vector (Aedes aegypti) causes both infections, all the quantities related to the mosquito, estimated from the initial phase of dengue epidemic, could be applied to yellow fever dynamics. It is demonstrated that R0 for yellow fever is, on average, 43% lower than that for dengue. This difference is due to the longer dengue viremia and its shorter extrinsic incubation period. RESULTS: In this study the analysis was expanded to the epidemiological situation of dengue in São Paulo in the year 2001. The total number of dengue cases increased from 3,582 in 2000 to 51,348 in 2001. It was then calculated R0 for yellow fever for every city which have shown R0 of dengue greater than 1. It was also estimated the total number of unprotected people living in highly risky areas for urban yellow fever. CONCLUSIONS: Currently there is a great number of non-vaccinated people living in Aedes aegypti infested area in the state of São Paulo.

Survival of a cohort of women with cervical cancer diagnosed in a Brazilian cancer center

OBJECTIVE: To assess overall survival of women with cervical cancer and describe prognostic factors associated. METHODS: A total of 3,341 cases of invasive cervical cancer diagnosed at the Brazilian Cancer Institute, Rio de Janeiro, southeastern Brazil, between 1999 and 2004 were selected. Clinical and pathological characteristics and follow-up data were collected. There were performed a survival analysis using Kaplan-Meier curves and a multivariate analysis through Cox model. RESULTS: Of all cases analyzed, 68.3% had locally advanced disease at the time of diagnosis. The 5-year overall survival was 48%. After multivariate analysis, tumor staging at diagnosis was the single variable significantly associated with prognosis (p<0.001). There was seen a dose-response relationship between mortality and clinical staging, ranging from 27.8 to 749.6 per 1,000 cases-year in women stage I and IV, respectively. CONCLUSIONS: The study showed that early detection through prevention programs is crucial to increase cervical cancer survival.

Some predictors of cardiovascular mortality among the elderly population of Botucatu (SP)

OBJECTIVE: To detect factors associated with cardiovascular mortality in the elderly of Botucatu. METHODS: We evaluated 29 variables of interest in a cohort of patients aged ³60 using data from a survey conducted between 1983/84. The elderly cohort was analyzed in 1992 to detect the occurrence of cardiovascular deaths. Survival analysis was performed using the Kaplan-Meier method, the log-rank test, and Cox regression analysis. Three models were adapted for each group of variables, and a final model was chosen from those variables selected from each group. RESULTS: We identified predictor for cardiovascular death according to age for elderly males not supporting the family, not possessing a vehicle, and previous cardiovascular disease. In elderly females, the predictor variables were previous cardiovascular disease and diabetes mellitus. CONCLUSION: Socioeconomic indicators (family heading and vehicle ownerrship) may be added to well stabilished medical factors (diabete mellitus and hypertension to select target groups for programs intended to reduce deaths due to cardiovascular diseases in elderly people.

Survival and Predictive Factors of Lethality in Hemodyalisis: D/I Polymorphism of The Angiotensin I-Converting Enzyme and of the Angiotensinogen M235T Genes

Background: End-stage kidney disease patients continue to have markedly increased cardiovascular disease morbidity and mortality. Analysis of genetic factors connected with the renin-angiotensin system that influences the survival of the patients with end-stage kidney disease supports the ongoing search for improved outcomes. Objective: To assess survival and its association with the polymorphism of renin-angiotensin system genes: angiotensin I-converting enzyme insertion/deletion and angiotensinogen M235T in patients undergoing hemodialysis. Methods: Our study was designed to examine the role of renin-angiotensin system genes. It was an observational study. We analyzed 473 chronic hemodialysis patients in four dialysis units in the state of Rio de Janeiro. Survival rates were calculated by the Kaplan-Meier method and the differences between the curves were evaluated by Tarone-Ware, Peto-Prentice, and log rank tests. We also used logistic regression analysis and the multinomial model. A p value ≤ 0.05 was considered to be statistically significant. The local medical ethics committee gave their approval to this study. Results: The mean age of patients was 45.8 years old. The overall survival rate was 48% at 11 years. The major causes of death were cardiovascular diseases (34%) and infections (15%). Logistic regression analysis found statistical significance for the following variables: age (p = 0.000038), TT angiotensinogen (p = 0.08261), and family income greater than five times the minimum wage (p = 0.03089), the latter being a protective factor. Conclusions: The survival of hemodialysis patients is likely to be influenced by the TT of the angiotensinogen M235T gene.

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.

New tools for epidemiology: a space odyssey

Geographical Information Systems (GIS) facilitate access to epidemiological data through visualization and may be consulted for the development of mathematical models and analysis by spatial statistics. Variables such as land-cover, land-use, elevations, surface temperatures, rainfall etc. emanating from earth-observing satellites, complement GIS as this information allows the analysis of disease distribution based on environmental characteristics. The strength of this approach issues from the specific environmental requirements of those causative infectious agents, which depend on intermediate hosts for their transmission. The distribution of these diseases is restricted, both by the environmental requirements of their intermediate hosts/vectors and by the ambient temperature inside these hosts, which effectively govern the speed of maturation of the parasite. This paper discusses the current capabilities with regard to satellite data collection in terms of resolution (spatial, temporal and spectral) of the sensor instruments on board drawing attention to the utility of computer-based models of the Earth for epidemiological research. Virtual globes, available from Google and other commercial firms, are superior to conventional maps as they do not only show geographical and man-made features, but also allow instant import of data-sets of specific interest, e.g. environmental parameters, demographic information etc., from the Internet.

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.