Azathioprine and 6-Mercaptopurine use in the Swiss IBD cohort : adverse effects, causes of discontinuation and risk of "flares" according to 6-TG levels

Background: To characterize and analyze in the Swiss IBD Cohort: a) reported Azathioprine (AZA) and 6-Mercaptopurine (6-MP) adverse effects (AE), b) causes of discontinuation and c) response to therapy according to gastroenterologists' clinical judgment, d) whether level of 6-TGN < 235pmol/8 x108 red blood cells (RBC) is associated with a higher risk of "flare" occurrence. Methods: Retrospective statistical description, Cox model and Kaplan-Meier survival estimation. Results: 1499 patients with Crohn's Disease (CD) and 1066 with Ulcerative colitis (UC).

Effect of Cumulating Exposure to Abacavir on the Risk of Cardiovascular Disease Events in Patients From the Swiss HIV Cohort Study.

BACKGROUND: Patients with HIV exposed to the antiretroviral drug abacavir may have an increased risk of cardiovascular disease (CVD). There is concern that this association arises because of a channeling bias. Even if exposure is a risk, it is not clear how that risk changes as exposure cumulates. METHODS: We assess the effect of exposure to abacavir on the risk of CVD events in the Swiss HIV Cohort Study. We use a new marginal structural Cox model to estimate the effect of abacavir as a flexible function of past exposures while accounting for risk factors that potentially lie on a causal pathway between exposure to abacavir and CVD. RESULTS: A total of 11,856 patients were followed for a median of 6.6 years; 365 patients had a CVD event (4.6 events per 1000 patient-years). In a conventional Cox model, recent--but not cumulative--exposure to abacavir increased the risk of a CVD event. In the new marginal structural Cox model, continued exposure to abacavir during the past 4 years increased the risk of a CVD event (hazard ratio = 2.06; 95% confidence interval: 1.43 to 2.98). The estimated function for the effect of past exposures suggests that exposure during the past 6-36 months caused the greatest increase in risk. CONCLUSIONS: Abacavir increases the risk of a CVD event: the effect of exposure is not immediate, rather the risk increases as exposure cumulates over the past few years. This gradual increase in risk is not consistent with a rapidly acting mechanism, such as acute inflammation.

Prognostic factors for progression-free and overall survival in advanced biliary tract cancer.

BACKGROUND: Biliary tract cancer is an uncommon cancer with a poor outcome. We assembled data from the National Cancer Research Institute (UK) ABC-02 study and 10 international studies to determine prognostic outcome characteristics for patients with advanced disease. METHODS: Multivariable analyses of the final dataset from the ABC-02 study were carried out. All variables were simultaneously included in a Cox proportional hazards model, and backward elimination was used to produce the final model (using a significance level of 10%), in which the selected variables were associated independently with outcome. This score was validated externally by receiver operating curve (ROC) analysis using the independent international dataset. RESULTS: A total of 410 patients were included from the ABC-02 study and 753 from the international dataset. An overall survival (OS) and progression-free survival (PFS) Cox model was derived from the ABC-02 study. White blood cells, haemoglobin, disease status, bilirubin, neutrophils, gender, and performance status were considered prognostic for survival (all with P < 0.10). Patients with metastatic disease {hazard ratio (HR) 1.56 [95% confidence interval (CI) 1.20-2.02]} and Eastern Cooperative Oncology Group performance status (ECOG PS) 2 had worse survival [HR 2.24 (95% CI 1.53-3.28)]. In a dataset restricted to patients who received cisplatin and gemcitabine with ECOG PS 0 and 1, only haemoglobin, disease status, bilirubin, and neutrophils were associated with PFS and OS. ROC analysis suggested the models generated from the ABC-02 study had a limited prognostic value [6-month PFS: area under the curve (AUC) 62% (95% CI 57-68); 1-year OS: AUC 64% (95% CI 58-69)]. CONCLUSION: These data propose a set of prognostic criteria for outcome in advanced biliary tract cancer derived from the ABC-02 study that are validated in an international dataset. Although these findings establish the benchmark for the prognostic evaluation of patients with ABC and confirm the value of longheld clinical observations, the ability of the model to correctly predict prognosis is limited and needs to be improved through identification of additional clinical and molecular markers.

Acute leukemias in Piauí: comparison with features observed in other regions of Brazil

Differences in age and sex distribution as well as FAB (French-American-British classification) types have been reported for acute leukemias in several countries. We studied the demographics and response to treatment of patients with acute myeloid leukemia (AML) and acute lymphoblastic leukemia (ALL) between 1989 and 2000 in Teresina, Piauí, and compared these results with reports from Brazil and other countries. Complete data concerning 345 patients (230 ALL, 115 AML) were reviewed. AML occurred predominantly in adults (77%), with a median age of 34 years, similar to that found in the southeast of Brazil but lower than the median age in the United States and Europe (52 years). FAB distribution was similar in children and adults and FAB-M2 was the most common type, as also found in Japan. The high frequency of FAB-M3 described in most Brazilian studies and for Hispanics in the United States was not observed. Overall survival for adults was 40%, similar to other studies in Brazil. A high mortality rate was observed during induction. No clinical or hematological parameter influenced survival in the Cox model. ALL presented the characteristic peak of incidence between 2-8 years. Most of the cases were CD10+ pre-B ALL. In 25%, abnormal expression of myeloid antigens was observed. Only 10% of the patients were older than 30 years. Overall survival was better for children. Age and leukocyte count were independent prognostic factors. These data demonstrate that, although there are regional peculiarities, the application of standardized treatments and good supportive care make it possible to achieve results observed in other countries for the same chemotherapy protocols.

An improved strategy for evaluating the extent of chronic arterial baroreceptor denervation in conscious rats

There is no index or criterion of aortic barodenervation, nor can we differentiate among rats that have suffered chronic sham, aortic or sino-aortic denervation. The objective of this study was to develop a procedure to generate at least one quantitative, reproducible and validated index that precisely evaluates the extent of chronic arterial barodenervation performed in conscious rats. Data from 79 conscious male Wistar rats of about 65-70 days of age with diverse extents of chronic arterial barodenervation and used in previous experiments were reanalyzed. The mean arterial pressure (MAP) and the heart rate (HR) of all rats were measured systematically before (over 1 h) and after three consecutive iv bolus injections of phenylephrine (PHE) and sodium nitroprusside (SNP). Four expressions of the effectiveness of barodenervation (MAP lability, PHE ratio, SNP ratio, and SNP-PHE slope) were assessed with linear fixed models, three-level average variance, average separation among levels, outlier box plot analysis, and overlapping graphic analysis. The analysis indicated that a) neither MAP lability nor SNP-PHE slope was affected by the level of chronic sodium intake; b) even though the Box-Cox transformations of both MAP lability [transformed lability index (TLI)] and SNP-PHE slope [transformed general sensitivity index (TGSI), {((3-(ΔHRSNP-ΔHRPHE/ΔMAPSNP-ΔMAPPHE))-0.4-1)/-0.04597}] could be two promising indexes, TGSI proved to be the best index; c) TLI and TGSI were not freely interchangeable indexes for this purpose. TGSI ranges that permit differentiation between sham (10.09 to 11.46), aortic (8.40 to 9.94) and sino-aortic (7.68 to 8.24) barodenervated conscious rats were defined.

Morphometric evaluation of nitric oxide synthase isoforms and their cytokine regulators predict pulmonary dysfunction and survival in systemic sclerosis

Because histopathological changes in the lungs of patients with systemic sclerosis (SSc) are consistent with alveolar and vessel cell damage, we presume that this interaction can be characterized by analyzing the expression of proteins regulating nitric oxide (NO) and plasminogen activator inhibitor-1 (PAI-1) synthesis. To validate the importance of alveolar-vascular interactions and to explore the quantitative relationship between these factors and other clinical data, we studied these markers in 23 cases of SSc nonspecific interstitial pneumonia (SSc-NSIP). We used immunohistochemistry and morphometry to evaluate the amount of cells in alveolar septa and vessels staining for NO synthase (NOS) and PAI-1, and the outcomes of our study were cellular and fibrotic NSIP, pulmonary function tests, and survival time until death. General linear model analysis demonstrated that staining for septal inducible NOS (iNOS) related significantly to staining of septal cells for interleukin (IL)-4 and to septal IL-13. In univariate analysis, higher levels of septal and vascular cells staining for iNOS were associated with a smaller percentage of septal and vascular cells expressing fibroblast growth factor and myofibroblast proliferation, respectively. Multivariate Cox model analysis demonstrated that, after controlling for SSc-NSIP histological patterns, just three variables were significantly associated with survival time: septal iNOS (P=0.04), septal IL-13 (P=0.03), and septal basic fibroblast growth factor (bFGF; P=0.02). Augmented NOS, IL-13, and bFGF in SSc-NSIP histological patterns suggest a possible functional role for iNOS in SSc. In addition, the extent of iNOS, PAI-1, and IL-4 staining in alveolar septa and vessels provides a possible independent diagnostic measure for the degree of pulmonary dysfunction and fibrosis with an impact on the survival of patients with SSc.

New statistical methods to assess the effect of time-dependent exposures in case-control studies

Contexte. Les études cas-témoins sont très fréquemment utilisées par les épidémiologistes pour évaluer l’impact de certaines expositions sur une maladie particulière. Ces expositions peuvent être représentées par plusieurs variables dépendant du temps, et de nouvelles méthodes sont nécessaires pour estimer de manière précise leurs effets. En effet, la régression logistique qui est la méthode conventionnelle pour analyser les données cas-témoins ne tient pas directement compte des changements de valeurs des covariables au cours du temps. Par opposition, les méthodes d’analyse des données de survie telles que le modèle de Cox à risques instantanés proportionnels peuvent directement incorporer des covariables dépendant du temps représentant les histoires individuelles d’exposition. Cependant, cela nécessite de manipuler les ensembles de sujets à risque avec précaution à cause du sur-échantillonnage des cas, en comparaison avec les témoins, dans les études cas-témoins. Comme montré dans une étude de simulation précédente, la définition optimale des ensembles de sujets à risque pour l’analyse des données cas-témoins reste encore à être élucidée, et à être étudiée dans le cas des variables dépendant du temps. Objectif: L’objectif général est de proposer et d’étudier de nouvelles versions du modèle de Cox pour estimer l’impact d’expositions variant dans le temps dans les études cas-témoins, et de les appliquer à des données réelles cas-témoins sur le cancer du poumon et le tabac. Méthodes. J’ai identifié de nouvelles définitions d’ensemble de sujets à risque, potentiellement optimales (le Weighted Cox model and le Simple weighted Cox model), dans lesquelles différentes pondérations ont été affectées aux cas et aux témoins, afin de refléter les proportions de cas et de non cas dans la population source. Les propriétés des estimateurs des effets d’exposition ont été étudiées par simulation. Différents aspects d’exposition ont été générés (intensité, durée, valeur cumulée d’exposition). Les données cas-témoins générées ont été ensuite analysées avec différentes versions du modèle de Cox, incluant les définitions anciennes et nouvelles des ensembles de sujets à risque, ainsi qu’avec la régression logistique conventionnelle, à des fins de comparaison. Les différents modèles de régression ont ensuite été appliqués sur des données réelles cas-témoins sur le cancer du poumon. Les estimations des effets de différentes variables de tabac, obtenues avec les différentes méthodes, ont été comparées entre elles, et comparées aux résultats des simulations. Résultats. Les résultats des simulations montrent que les estimations des nouveaux modèles de Cox pondérés proposés, surtout celles du Weighted Cox model, sont bien moins biaisées que les estimations des modèles de Cox existants qui incluent ou excluent simplement les futurs cas de chaque ensemble de sujets à risque. De plus, les estimations du Weighted Cox model étaient légèrement, mais systématiquement, moins biaisées que celles de la régression logistique. L’application aux données réelles montre de plus grandes différences entre les estimations de la régression logistique et des modèles de Cox pondérés, pour quelques variables de tabac dépendant du temps. Conclusions. Les résultats suggèrent que le nouveau modèle de Cox pondéré propose pourrait être une alternative intéressante au modèle de régression logistique, pour estimer les effets d’expositions dépendant du temps dans les études cas-témoins

Facteurs de risque de chutes chez les aînés vivant dans la communauté et ayant recours aux services de soutien à domicile : covariables dépendantes du temps et événements récurrents

Les chutes chez les personnes âgées représentent un problème majeur. Il n’est donc pas étonnant que l’identification des facteurs qui en accroissent le risque ait mobilisé autant d’attention. Les aînés plus fragiles ayant besoin de soutien pour vivre dans la communauté sont néanmoins demeurés le parent pauvre de la recherche, bien que, plus récemment, les autorités québécoises en aient fait une cible d’intervention prioritaire. Les études d’observation prospectives sont particulièrement indiquées pour étudier les facteurs de risque de chutes chez les personnes âgées. Leur identification optimale est cependant compliquée par le fait que l’exposition aux facteurs de risque peut varier au cours du suivi et qu’un même individu peut subir plus d’un événement. Il y a 20 ans, des chercheurs ont tenté de sensibiliser leurs homologues à cet égard, mais leurs efforts sont demeurés vains. On continue aujourd’hui à faire peu de cas de ces considérations, se concentrant sur la proportion des personnes ayant fait une chute ou sur le temps écoulé jusqu’à la première chute. On écarte du coup une quantité importante d’information pertinente. Dans cette thèse, nous examinons les méthodes en usage et nous proposons une extension du modèle de risques de Cox. Nous illustrons cette méthode par une étude des facteurs de risque susceptibles d’être associés à des chutes parmi un groupe de 959 personnes âgées ayant eu recours aux services publics de soutien à domicile. Nous comparons les résultats obtenus avec la méthode de Wei, Lin et Weissfeld à ceux obtenus avec d’autres méthodes, dont la régression logistique conventionnelle, la régression logistique groupée, la régression binomiale négative et la régression d’Andersen et Gill. L’investigation est caractérisée par des prises de mesures répétées des facteurs de risque au domicile des participants et par des relances téléphoniques mensuelles visant à documenter la survenue des chutes. Les facteurs d’exposition étudiés, qu’ils soient fixes ou variables dans le temps, comprennent les caractéristiques sociodémographiques, l’indice de masse corporelle, le risque nutritionnel, la consommation d’alcool, les dangers de l’environnement domiciliaire, la démarche et l’équilibre, et la consommation de médicaments. La quasi-totalité (99,6 %) des usagers présentaient au moins un facteur à haut risque. L’exposition à des risques multiples était répandue, avec une moyenne de 2,7 facteurs à haut risque distincts par participant. Les facteurs statistiquement associés au risque de chutes incluent le sexe masculin, les tranches d’âge inférieures, l’histoire de chutes antérieures, un bas score à l’échelle d’équilibre de Berg, un faible indice de masse corporelle, la consommation de médicaments de type benzodiazépine, le nombre de dangers présents au domicile et le fait de vivre dans une résidence privée pour personnes âgées. Nos résultats révèlent cependant que les méthodes courantes d’analyse des facteurs de risque de chutes – et, dans certains cas, de chutes nécessitant un recours médical – créent des biais appréciables. Les biais pour les mesures d’association considérées proviennent de la manière dont l’exposition et le résultat sont mesurés et définis de même que de la manière dont les méthodes statistiques d’analyse en tiennent compte. Une dernière partie, tout aussi innovante que distincte de par la nature des outils statistiques utilisés, complète l’ouvrage. Nous y identifions des profils d’aînés à risque de devenir des chuteurs récurrents, soit ceux chez qui au moins deux chutes sont survenues dans les six mois suivant leur évaluation initiale. Une analyse par arbre de régression et de classification couplée à une analyse de survie a révélé l’existence de cinq profils distinctifs, dont le risque relatif varie de 0,7 à 5,1. Vivre dans une résidence pour aînés, avoir des antécédents de chutes multiples ou des troubles de l’équilibre et consommer de l’alcool sont les principaux facteurs associés à une probabilité accrue de chuter précocement et de devenir un chuteur récurrent. Qu’il s’agisse d’activité de dépistage des facteurs de risque de chutes ou de la population ciblée, cette thèse s’inscrit dans une perspective de gain de connaissances sur un thème hautement d’actualité en santé publique. Nous encourageons les chercheurs intéressés par l’identification des facteurs de risque de chutes chez les personnes âgées à recourir à la méthode statistique de Wei, Lin et Weissfeld car elle tient compte des expositions variables dans le temps et des événements récurrents. Davantage de recherches seront par ailleurs nécessaires pour déterminer le choix du meilleur test de dépistage pour un facteur de risque donné chez cette clientèle.

Validation des modèles statistiques tenant compte des variables dépendantes du temps en prévention primaire des maladies cérébrovasculaires

L’intérêt principal de cette recherche porte sur la validation d’une méthode statistique en pharmaco-épidémiologie. Plus précisément, nous allons comparer les résultats d’une étude précédente réalisée avec un devis cas-témoins niché dans la cohorte utilisé pour tenir compte de l’exposition moyenne au traitement : – aux résultats obtenus dans un devis cohorte, en utilisant la variable exposition variant dans le temps, sans faire d’ajustement pour le temps passé depuis l’exposition ; – aux résultats obtenus en utilisant l’exposition cumulative pondérée par le passé récent ; – aux résultats obtenus selon la méthode bayésienne. Les covariables seront estimées par l’approche classique ainsi qu’en utilisant l’approche non paramétrique bayésienne. Pour la deuxième le moyennage bayésien des modèles sera utilisé pour modéliser l’incertitude face au choix des modèles. La technique utilisée dans l’approche bayésienne a été proposée en 1997 mais selon notre connaissance elle n’a pas été utilisée avec une variable dépendante du temps. Afin de modéliser l’effet cumulatif de l’exposition variant dans le temps, dans l’approche classique la fonction assignant les poids selon le passé récent sera estimée en utilisant des splines de régression. Afin de pouvoir comparer les résultats avec une étude précédemment réalisée, une cohorte de personnes ayant un diagnostique d’hypertension sera construite en utilisant les bases des données de la RAMQ et de Med-Echo. Le modèle de Cox incluant deux variables qui varient dans le temps sera utilisé. Les variables qui varient dans le temps considérées dans ce mémoire sont iv la variable dépendante (premier évènement cérébrovasculaire) et une des variables indépendantes, notamment l’exposition

Risque d’acquisition du virus de l’immunodéficience humaine (VIH) chez les femmes utilisant des hormones contraceptives orales et injectables

Objectif : Étudier l'association entre l’utilisation de contraceptifs hormonaux et le risque d'acquisition du VIH-1 chez les femmes au Malawi, en Afrique du Sud, en Zambie et au Zimbabwe. Devis : Analyses secondaires de 2887 femmes âgées de 17-55 ans ayant participé à l’étude HPTN 035, une étude de phase II/IIb sur l’efficacité de deux gels microbicides pour prévenir la transmission du VIH chez les femmes à risque. Méthodes : L'association entre l'utilisation de contraceptifs hormonaux et le risque d'acquisition du VIH-1 a été évaluée en utilisant des modèles de Cox. Des risques relatifs sont estimés où le groupe de référence est celui des femmes qui n’utilisent pas de contraceptifs hormonaux. De plus, un modèle multivarié de Cox est utilisé afin de contrôler pour les facteurs potentiellement confondants. Résultats : Les contraceptifs injectables ont été utilisés par 52,1% des femmes, alors que les contraceptifs oraux ont été utilisés par 20,7% de celles-ci. Pendant l'étude, il y a eu 192 séroconversions. L'incidence observée du VIH était de 2,28; 4,19 et 4,69 pour 100 personne-années pour les contraceptifs oraux, injectables et non hormonaux, respectivement. Lors de l’analyse multivariée, nous n'avons trouvé aucune association significative entre l’usage des contraceptifs hormonaux et l’acquisition du VIH-1. Le risque relatif ajusté (RRa) pour les contraceptifs oraux est de 0,573 (IC de 95% : [0,31-1,06]) et 0,981 (IC de 95% : [0,69 ; 1,39]) pour les contraceptifs injectables. Conclusions : Bien que cette étude ne démontre pas d’association entre l’usage des contraceptifs hormonaux et le VIH-1, nous concluons toutefois que ces méthodes de contraception ne protègent pas contre le VIH-1, et il est ainsi recommandé aux femmes utilisant des hormones contraceptives de toujours utiliser le condom pour prévenir l'infection au VIH-1.

Marine04 marine radiocarbon age calibration, 0-26 cal kyr BP

New radiocarbon calibration curves, IntCal04 and Marine04, have been constructed and internationally ratified to replace the terrestrial and marine components of IntCal98. The new calibration data sets extend an additional 2000 yr, from 0-26 cal kyr BP (Before Present, 0 cal. BP = AD 1950), and provide much higher resolution, greater precision, and more detailed structure than IntCal98. For the Marine04 curve, dendrochronologically-dated tree-ring samples, converted with a box diffusion model to marine mixed-layer ages, cover the period from 0-10.5 call kyr BR Beyond 10.5 cal kyr BP, high-resolution marine data become available from foraminifera in varved sediments and U/Th-dated corals. The marine records are corrected with site-specific C-14 reservoir age information to provide a single global marine mixed-layer calibration from 10.5-26.0 cal kyr BR A substantial enhancement relative to IntCal98 is the introduction of a random walk model, which takes into account the uncertainty in both the calendar age and the C-14 age to calculate the underlying calibration curve (Buck and Blackwell, this issue). The marine data sets and calibration curve for marine samples from the surface mixed layer (Marine04) are discussed here. The tree-ring data sets, sources of uncertainty, and regional offsets are presented in detail in a companion paper by Reimer et al. (this issue).

Irreversible compressible work and APE dissipation in turbulent stratified fluid

Although it plays a key role in the theory of stratified turbulence, the concept of available potential energy (APE) dissipation has remained until now a rather mysterious quantity, owing to the lack of rigorous result about its irreversible character or energy conversion type. Here, we show by using rigorous energetics considerations rooted in the analysis of the Navier-Stokes for a fully compressible fluid with a nonlinear equation of state that the APE dissipation is an irreversible energy conversion that dissipates kinetic energy into internal energy, exactly as viscous dissipation. These results are established by showing that APE dissipation contributes to the irreversible production of entropy, and by showing that it is a part of the work of expansion/contraction. Our results provide a new interpretation of the entropy budget, that leads to a new exact definition of turbulent effective diffusivity, which generalizes the Osborn-Cox model, as well as a rigorous decomposition of the work of expansion/contraction into reversible and irreversible components. In the context of turbulent mixing associated with parallel shear flow instability, our results suggests that there is no irreversible transfer of horizontal momentum into vertical momentum, as seems to be required when compressible effects are neglected, with potential consequences for the parameterisations of momentum dissipation in the coarse-grained Navier-Stokes equations.

Coronary reserve impairment prevents the improvement of left ventricular dysfunction and adversely affects the long-term outcome of patients with hypertensive dilated cardiomyopathy

In hypertension, left ventricular (LV) hypertrophy develops as an adaptive mechanism to compensate for increased afterload and thus preserve systolic function. Associated structural changes such as microvascular disease might potentially interfere with this mechanism, producing pathological hypertrophy. A poorer outcome is expected to occur when LV function is put in jeopardy by impaired coronary reserve. The aim of this study was to evaluate the role of coronary reserve in the long-term outcome of patients with hypertensive dilated cardiomyopathy. Between 1996 and 2000, 45 patients, 30 of them male, with 52 +/- 11 years and LV fractional shortening <30% were enrolled and followed until 2006. Coronary flow velocity reserve was assessed by transesophageal Doppler of the left anterior descending coronary artery. Sixteen patients showed >= 10% improvement in LV fractional shortening after 17 +/- 6 months. Coronary reserve was the only variable independently related to this improvement. Total mortality was 38% in 10 years. The Cox model identified coronary reserve (hazard ratio = 0.814; 95% CI = 0.72-0.92), LV mass, low diastolic blood pressure, and male gender as independent predictors of mortality. In hypertensive dilated cardiomyopathy, coronary reserve impairment adversely affects survival, possibly by interfering with the improvement of LV dysfunction. J Am Soc Hypertens 2010;4(1):14-21. (C) 2010 American Society of Hypertension. All rights reserved.

A family of autoregressive conditional duration models

This paper develops a family of autoregressive conditional duration (ACD) models that encompasses most specifications in the literature. The nesting relies on a Box-Cox transformation with shape parameter λ to the conditional duration process and a possibly asymmetric shocks impact curve. We establish conditions for the existence of higher-order moments, strict stationarity, geometric ergodicity and β-mixing property with exponential decay. We next derive moment recursion relations and the autocovariance function of the power λ of the duration process. Finally, we assess the practical usefulness of our family of ACD models using NYSE transactions data, with special attention to IBM price durations. The results warrant the extra flexibility provided either by the Box-Cox transformation or by the asymmetric response to shocks.

Aprimoramento estatístico da regionalização de vazões máximas e médias. Aplicação a bacias hidrográficas do Rio Grande do Sul e Santa Catarina

Os registros de vazões líquidas obtidos das estações fluviométricas são válidos, unicamente, para o local de onde foram coletados ou muito próximo a ele. Na maioria das vezes, a região de influência deste não inclui o local de interesse para o desenvolvimento de projetos de recursos hídricos. Este inconveniente, geralmente, pode ser resolvido através do uso de métodos de regionalização hidrológica. Para determinar os coeficientes da equação de regionalização, o procedimento mais usado consiste em: i) estabelecer uma relação do tipo exponencial multivariada entre a variável dependente (vazão média de longo prazo ou média anual de cheia) e as covariáveis (variáveis climáticas e fisiográficas da bacia hidrográfica); ii) linearizar a equação anterior mediante a transformação logarítmica de ambos os membros; iii) utilizar modelos lineares de regressão para estimar os coeficientes, geralmente, o método dos mínimos quadrados ordinários; e iv) aplicar a transformação inversa para definir a equação. A aplicação deste procedimento implica assumir certas propriedades dos dados (assimetria positiva, registros da mesma extensão e que os mesmos possuem o mesmo período de início e fim, entre outros) que dificilmente podem ser atendidas, prejudicando a verificação das hipóteses nas quais estão baseados os métodos lineares de regressão e, em conseqüência, seu adequado uso, bem como a confiabilidade dos resultados obtidos. Esta pesquisa apresenta um aprimoramento dos métodos de regionalização de vazões geralmente empregados, incluindo-se técnicas que levam em consideração as limitações anteriores. Estas técnicas foram: i) uso da transformada de Box-Cox na linearização da equação exponencial multivariada; ii) determinação dos coeficientes da equação de regionalização usando mínimos quadrados ponderados; e iii) verificação se os resíduos da regressão estão correlacionados ou não. Para o desenvolvimento e verificação da metodologia proposta foram usados somente registros fluviométricos de Bacias Hidrográficas Brasileiras, que drenam suas águas para o Rio Grande do Sul e/ou que estejam localizadas dentro dele. Geograficamente, a área de estudo inclui a totalidade do estado do Rio Grande do Sul e parte de Santa Catarina. As equações de regionalização foram definidas usando dados de vazões médias de longo prazo e média de cheia, para tempo de retorno de 2,33 e 50 anos. Neste último caso, as freqüências foram estimadas através do método dos momentos-L.Comparando os resultados obtidos utilizando o modelo de regionalização hidrológica proposto neste trabalho (transformada de Box-Cox / mínimos quadrados ponderados) junto a seus similares gerados usando a metodologia convencional (transformada logarítmica / mínimos quadrados ordinários) e de modelos intermediários (transformada logarítmica / mínimos quadrados ponderados e transformada de Box-Cox / mínimos quadrados ordinários), os mesmos podem ser considerados satisfatórios, visto que em todas as simulações realizadas o modelo proposto forneceu melhores resultados que aqueles obtidos com os outros modelos, sendo utilizado como padrão de comparação: 1) a qualidade do ajuste, 2) o grau de verificação das hipóteses dos métodos lineares de regressão e 3) os erros na estimativa das descargas, em termos de vazão específica. Nas simulações realizadas usando os modelos intermediários, observou-se que: i) na regionalização de vazões médias, o ganho de considerar a heterogeneidade temporal dos dados é maior do que corrigir a assimetria dos mesmos; ii) quando são usadas séries de descargas máximas, ocorre o efeito contrário, visto que o ganho de corrigir a assimetria das séries é maior do que o efeito da heterogeneidade temporal dos dados. Com relação aos resíduos da regressão, contrariamente ao esperado, os mesmos não sugerem estar correlacionados; isto pode ser conseqüência de utilizar como variável dependente um único registro por estação (vazão média de longo prazo ou média anual de cheia).