Avaliação económica em transplantação renal : abordagem estratégica do processo de doação e análise de custo-utilidade do programa de transplantação renal versus hemodiálise

SUMÁRIO - O desafio atual da Saúde Pública é assegurar a sustentabilidade financeira do sistema de saúde. Em ambiente de recursos escassos, as análises económicas aplicadas à prestação dos cuidados de saúde são um contributo para a tomada de decisão que visa a maximização do bem-estar social sujeita a restrição orçamental. Portugal é um país com 10,6 milhões de habitantes (2011) com uma incidência e prevalência elevadas de doença renal crónica estadio 5 (DRC5), respetivamente, 234 doentes por milhão de habitantes (pmh) e 1.600 doentes/pmh. O crescimento de doenças associadas às causas de DRC, nomeadamente, diabetes Mellitus e hipertensão arterial, antecipam uma tendência para o aumento do número de doentes. Em 2011, dos 17.553 doentes em tratamento substitutivo renal, 59% encontrava-se em programa de hemodiálise (Hd) em centros de diálise extra-hospitalares, 37% viviam com um enxerto renal funcionante e 4% estavam em diálise peritoneal (SPN, 2011). A lista ativa para transplante (Tx) renal registava 2.500 doentes (SPN 2009). O Tx renal é a melhor modalidade terapêutica pela melhoria da sobrevida, qualidade de vida e relação custo-efetividade, mas a elegibilidade para Tx e a oferta de órgãos condicionam esta opção. Esta investigação desenvolveu-se em duas vertentes: i) determinar o rácio custo-utilidade incremental do Tx renal comparado com a Hd; ii) avaliar a capacidade máxima de dadores de cadáver em Portugal, as características e as causas de morte dos dadores potenciais a nível nacional, por hospital e por Gabinete Coordenador de Colheita e Transplantação (GCCT), e analisar o desempenho da rede de colheita de órgãos para Tx. Realizou-se um estudo observacional/não interventivo, prospetivo e analítico que incidiu sobre uma coorte de doentes em Hd que foi submetida a Tx renal. O tempo de seguimento mínimo foi de um ano e máximo de três anos. No início do estudo, colheram-se dados sociodemográficos e clínicos em 386 doentes em Hd, elegíveis para Tx renal. A qualidade de vida relacionada com a saúde (QVRS) foi avaliada nos doentes em Hd (tempo 0) e nos transplantados, aos três, seis, 12 meses, e depois, anualmente. Incluíram-se os doentes que por falência do enxerto renal transitaram para Hd. Na sua medição, utilizou-se um instrumento baseado em preferências da população, o EuroQol-5D, que permite o posterior cálculo dos QALY. Num grupo de 82 doentes, a QVRS em Hd foi avaliada em dois tempos de resposta o que permitiu a análise da sua evolução. Realizou-se uma análise custo-utilidade do Tx renal comparado com a Hd na perspetiva da sociedade. Identificaram-se os custos diretos, médicos e não médicos, e as alterações de produtividade em Hd e Tx renal. Incluíram-se os custos da colheita de órgãos, seleção dos candidatos a Tx renal e follow-up dos dadores vivos. Cada doente transplantado foi utilizado como controle de si próprio em diálise. Avaliou-se o custo médio anual em programa de Hd crónica relativo ao ano anterior à Tx renal. Os custos do Tx foram avaliados prospetivamente. Considerou-se como horizonte temporal o ciclo de vida nas duas modalidades. Usaram-se taxas de atualização de 0%, 3% e 5% na atualização dos custos e QALY e efetuaram-se análises de sensibilidade one way. Entre 2008 e 2010, 65 doentes foram submetidos a Tx renal. Registaram-se, prospetivamente, os resultados em saúde incluíndo os internamentos e os efeitos adversos da imunossupressão, e o consumo dos recursos em saúde. Utilizaram-se modelos de medidas repetidas na avaliação da evolução da QVRS e modelos de regressão múltipla na análise da associação da QVRS e dos custos do transplante com as características basais dos doentes e os eventos clínicos. Comparativamente à Hd, observou-se melhoria da utilidade ao 3º mês de Tx e a qualidade de vida aferida pela escala EQ-VAS melhorou em todos os tempos de observação após o Tx renal. O custo médio da Hd foi de 32.567,57€, considerado uniforme ao longo do tempo. O custo médio do Tx renal foi de 60.210,09€ no 1º ano e 12.956,77€ nos anos seguintes. O rácio custo-utilidade do Tx renal vs Hd crónica foi de 2.004,75€/QALY. A partir de uma sobrevivência do enxerto de dois anos e cinco meses, o Tx associou-se a poupança dos custos. Utilizaram-se os dados nacionais dos Grupos de Diagnóstico Homogéneos e realizou-se um estudo retrospectivo que abrangeu as mortes ocorridas em 34 hospitais com colheita de órgãos, em 2006. Considerou-se como dador potencial o indivíduo com idade entre 1-70 anos cuja morte ocorrera a nível hospitalar, e que apresentasse critérios de adequação à doação de rim. Analisou-se a associação dos dadores potenciais com características populacionais e hospitalares. O desempenho das organizações de colheita de órgãos foi avaliado pela taxa de conversão (rácio entre os dadores potenciais e efetivos) e pelo número de dadores potenciais por milhão de habitantes a nível nacional, regional e por Gabinete Coordenador de Colheita e Transplantação (GCCT). Identificaram-se 3.838 dadores potenciais dos quais 608 apresentaram códigos da Classificação Internacional de Doenças, 9.ª Revisão, Modificações Clínicas (CID- 9-MC) que, com maior frequência, evoluem para a morte cerebral. O modelo logit para dados agrupados identificou a idade, o rácio da lotação em Unidades de Cuidados Intensivos e lotação de agudos, existência de GCCT e de Unidade de Transplantação, e mortalidade por acidente de trabalho como fatores preditivos da conversão dum dador potencial em efetivo e através das estimativas do modelo logit quantificou-se a probabilidade dessa conversão. A doação de órgãos deve ser assumida como uma prioridade e as autoridades em saúde devem assegurar o financiamento dos hospitais com programas de doação, evitando o desperdício de órgãos para transplantação, enquanto um bem público e escasso. A colheita de órgãos deve ser considerada uma opção estratégica da atividade hospitalar orientada para a organização e planeamento de serviços que maximizem a conversão de dadores potenciais em efetivos incluindo esse critério como medida de qualidade e efetividade do desempenho hospitalar. Os resultados deste estudo demonstram que: 1) o Tx renal proporciona ganhos em saúde, aumento da sobrevida e qualidade de vida, e poupança de custos; 2) em Portugal, a taxa máxima de eficácia da conversão dos dadores cadavéricos em dadores potenciais está longe de ser atingida. O investimento na rede de colheita de órgãos para Tx é essencial para assegurar a sustentabilidade financeira e promover a qualidade, eficiência e equidade dos cuidados em saúde prestados na DRC5.

Modelação dinâmica do embate de embarcações sobre pilares de pontes

Após o embate de uma embarcação contra o maciço do pilar da Ponte 25 de Abril em Lisboa e tendo em conta a importância desta estrutura, levantaram-se algumas preocupações referentes à segurança da ponte quanto a este tipo de solicitações. Como tal, no seguimento desta problemática, o principal objetivo do presente trabalho é avaliar a resposta da Ponte 25 de Abril sujeita ao embate de embarcações. Sabe-se que a Ponte 25 de Abril foi dimensionada tendo como base o sismo de El Centro de 1940, não tendo sido feito um estudo pormenorizado na temática do embate de uma embarcação uma vez que os momentos na base do maciço, resultantes do sismo estudado no projeto, eram largamente superiores aos mesmos momentos provocados por uma força no topo do maciço que pretendia simular o embate de uma embarcação. No entanto, as expressões que pretendem simular o efeito de um embate de uma embarcação em pilares de pontes evoluíram desde a data do projeto. Sendo os maciços dos pilares os elementos estruturais diretamente sujeitos ao evento do embate, transmitindo posteriormente os deslocamentos e forças aos respetivos pilares, o presente estudo focou-se nestes elementos. Para se estudar o comportamento do maciço de um pilar da ponte sujeito a estas ações dinâmicas, foram necessários dados para se definir a estrutura em causa e a respetiva solicitação, como também foi necessário um programa de cálculo não-linear de estruturas indicado para o presente caso. Obteve-se o projeto da Ponte 25 de Abril assim como dados relativos a todas as embarcações que cruzaram a Ponte no ano de 2014. A modelação da estrutura e respetivas solicitações foram feitas com recurso ao programa de cálculo não-linear de estruturas Extreme Loading for Structures, que utiliza o Método dos Elementos Aplicados. A partir do modelo do maciço do pilar da ponte traçaram-se as curvas de capacidade dinâmica e estática da estrutura e simularam-se tanto o sismo de El Centro de 1940 como o embate da embarcação mais condicionante que poderia ter ocorrido no ano de 2014. Com base nos resultados obtidos nos ensaios efetuados no presente trabalho verificou-se que, embora o embate de uma embarcação possa ser mais condicionante do que o sismo para o qual foi dimensionada, a estrutura do maciço resiste a ambas as solicitações não apresentando danos significativos.

Análise da demora média pré-operatória no Serviço Nacional de Saúde

RESUMO - INTRODUÇÃO A escassez de recursos e os crescentes custos em saúde conduzem à adoção de medidas que visam um melhor desempenho das organizações de saúde. O facto da despesa corrente na saúde ter aumentado nos últimos anos e até 2010 resultou, entre outras, das seguintes condicionantes: (i) aumento da esperança média de vida, pressupondo-se uma maior canalização e utilização de recursos para a prestação de cuidados de saúde à população idosa; (ii) inovação na medicina, ao nível de diagnóstico, de tratamento e de tecnologias da saúde, repercutindo-se no eventual aumento de expectativas e exigências da população em relação aos seus níveis de saúde, acessibilidade e qualidade dos cuidados de saúde. Atualmente, os gestores hospitalares deparam-se com o imperativo da utilização racional e eficiente de recursos na saúde. O presente trabalho resultou de dois fatores, designadamente: (i) a tendência crescente da despesa pública relativa à saúde e (ii) os níveis de eficiência dos hospitais do SNS. O objetivo geral consistiu em identificar o potencial de melhoria na eficiência hospitalar em relação à DMPO observada no universo do SNS, medida através de dias de internamento desnecessários. Para tal, pretendeu-se: (i) descrever e analisar a variação na demora média pré-operatória de doentes submetidos a cirurgia eletiva nos hospitais que integram o Serviço Nacional de Saúde, identificando os fatores inerentes aos doentes e de gestão que influenciam a demora média pré-operatória; (ii) estimar os dias de internamento pré-operatórios desnecessários por hospital, por Grande Categoria de Diagnóstico, por Grupo de Diagnóstico Homogéneo e por dia da semana da admissão do doente; e (iii) calcular a valorização, em termos monetários, dos dias de internamento pré-operatório desnecessários ao nível dos Grupos de Diagnósticos Homogéneos. METODOLOGIA Pretendeu-se analisar a demora média pré-operatória para o triénio 2009-2011, segundo as variáveis independentes sexo, idade, número de diagnósticos secundários, Grupo de Diagnóstico Homogéneo, hospital, dia da semana da admissão do doente, recorrendo a testes de hipóteses envolvendo até duas variáveis e ao modelo de regressão linear múltipla. Também se pretendeu estimar os dias de internamento pré-operatório desnecessários por hospital, por Grande Categoria de Diagnóstico, por Grupo de Diagnóstico Homogéneo e por dia da semana da admissão do doente. Por último, procedeu-se à valorização dos dias de internamento pré-operatório desnecessários em termos monetários. RESULTADOS Os principais resultados apontam para a relação entre as variáveis independentes consideradas e a demora média pré-operatória, sendo que o poder explicativo das mesmas em relação à demora média pré-operatória reduzido (7,3%), o que pressupõe que haja outros fatores com maior influência na demora média pré-operatória. Em termos globais a demora média pré-operatória rondou os 1,08 dias, sendo que existe margem para melhorar. CONCLUSÃO Atendendo à demora média pré-operatória, parece haver margem para incrementar a produção hospitalar com os recursos de que o Serviço Nacional de Saúde dispõe atualmente (e.g. minimizando os dias de internamento pré-operatórios desnecessários), sem prejuízo dos níveis de qualidade; os hospitais do Serviço Nacional de Saúde podem incrementar os níveis de desempenho no que respeita à eficiência, mormente em relação à demora média pré-operatória; e aparentemente, de um modo um pouco grosseiro, a valorização dos dias de internamento pré-operatório desnecessários para os Grupos de Diagnósticos Homogéneos analisados repercutiu-se em cerca de 85 milhões de €.

A PESCA ARTESANAL NO BAIXO AMAZONAS: ANÁLISE MULTIVARIADA DA CAPTURA POR ESPÉCIE.

Dados sobre esforço e captura por espécie correspondentes à produção pesqueira desembarcada no porto de Santarém, no ano de 1993, foram submetidos a duas técnicas de análise multivariada: uma análise de fatores segundo o método de componentes principais e outra de co-variância conforme o modelo linear geral (GLM). Os resultados indicam que a atividade pesqueira na região está influenciada pelas características dos ciclos de vida das espécies-alvo, pelo ciclo hidrológico e condição climáticas de sistema, e ainda pelas preferências culturais e interesses econômicos do mercado consumidor. Os resultados indicam o direcionamento da atividade pesqueira para determinados grupos de espécies, e as diferentes variáveis incluídas no modelo permitiram uma explicação aproximada dos padrões desse direcionamento. Esses padrões são: (a) pesca de grandes bagres, alvo de pesca para a exportação (FATOR 1); (b) pesca de entre-safra do mapará (Hypophthalmusspp.) e da pescada (Plagioscionspp.), realizada nos lagos com um componente para exportação e outro para o consumo local (FATOR4); (c) pesca de peixes de escamas de hábitos sedentários e/ou migratórios, alvo da pesca comercial de pequena escala e de importância no mercado local (FATOR2, FATOR3 e FATOR5), que inclui pescarias importantes como a do tambaqui (Colossoma macropomum)e do pirarucu(Arapaima gigas)nos lagos, ou a do jaraqui (Semaprochilodusspp.) e pacu (Mylossomaspp. e Metynnisspp.) nos rios.

Detecção de cicatrizes de áreas queimadas baseada no modelo linear de mistura espectral e imagens índice de vegetação utilizando dados multitemporais do sensor MODIS/TERRA no estado do Mato Grosso, Amazônia brasileira

O objetivo desta pesquisa foi avaliar os dados do sensor MODIS para detectar e monitorar cicatrizes de áreas recém queimadas. Utilizamos imagens da reflectância de superfície do sensor MODIS: produto MOD09 (dia 5 de outubro) e produto MOD13A1 (meses de outubro e novembro). Foi avaliada também uma série temporal de um ano dos índices de vegetação (IV) EVI e NDVI (produto MOD13A1). Uma imagem do sensor ETM+ (dia 5 de outubro) foi utilizada como base para a delimitação dos polígonos amostrais e avaliação dos dados MODIS devido a sua melhor resolução espacial. A metodologia focou na aplicação do modelo linear de mistura espectral nas imagens reflectância para a geração das imagens fração sombra. Análises de regressão foram efetuadas para comparação entre o percentual de sombra derivado da imagem ETM+ e das imagens MODIS. As alterações multitemporais nas imagens IV foram avaliadas com base no teste de Tukey. Os resultados mostraram que a imagem fração sombra gerada a partir do produto MOD09 apresentou um R² = 0,66 (p < 0,01) em relação aos dados ETM+. Para as imagens do produto MOD13A1 não foram identificadas relações significativas. Os IV dentro dos mesmos polígonos apresentaram uma variação sazonal durante o ano. No entanto, não houve uma diminuição significativa dos valores destes índices nos meses onde foram observadas as cicatrizes de áreas recém queimadas. Portanto, o produto MOD09 mostrou-se mais eficiente que o produto MOD13A1 para a detecção de cicatrizes de áreas recém queimadas. A análise multitemporal dos IV sugeriu que não foi possível detectar este mesmo padrão na área de estudo.

Instrumentos de políticas públicas do turismo: uma análise empírica dos municípios portugueses

Tese de Doutoramento em Ciências da Administração

Desenvolvimento e otimização do projeto mecatrónico para um Tribómetro Linear

Dissertação de mestrado integrado em Engenharia Mecânica

Modelação de processos de negócio: análise comparativa de linguagens

Dissertação de mestrado integrado em Engenharia e Gestão de Sistemas de Informação

Seismic behavior of concrete block masonry buildings

Tese de Doutoramento em Engenharia Civil

A estética audiovisual: o caso do Canal 180

Relatório de estágio de mestrado em Ciências da Comunicação (área de especialização em Audiovisual e Multimédia)

Análise da função ventricular de pacientes com doença pulmonar obstrutiva crônica avançada pela ressonância magnética

OBJETIVO: Avaliar os possíveis mecanismos envolvidos na evolução da disfunção ventricular direita e esquerda em pacientes com doença pulmonar obstrutiva crônica (DPOC). MÉTODOS: A ressonância magnética foi aplicada em 27 pacientes com DPOC divididos em grupos, DPOCc e DPOCs, de acordo com a presença ou ausência de disfunção ventricular direita, respectivamente, e 11 controles (grupo C). O exame Doppler ecocardiográfico foi empregado para análise do grau de hipertensão pulmonar. RESULTADOS: O diâmetro do ventrículo direito foi similar nos 3 grupos, DPOCs, DPOCc e C, 29±8 mm vs 31±7 mm vs 30±6 mm; p NS, respectivamente. Foram observados hipertrofia ventricular direita somente nos grupos DPOC (8±2 mm e 9±3 mm vs 5±1 mm; p<0,01), percentual de espessamento sistólico da parede lateral do ventrículo direito similar (%RVLWT) (86±82% vs 41±35% vs 86±89%; p NS) e diferentes frações de ejeção ventricular esquerda entre os 3 grupos (69±9% vs 55±16% vs 76±6%; p<0,01); correlação linear positiva e significante entre o diâmetro diastólico do ventrículo esquerdo (VE) e o débito sistólico do VE (r=0,72, p < 0,01). Não houve correlação entre os volumes pulmonares e a gasometria arterial com a função ventricular. CONCLUSÃO: Não houve correlação entre a gravidade da função pulmonar e o grau de comprometimento da função ventricular. Fica a ser elucidado se um preservado %RVLWT significa a possibilidade de reversibilidade da função ventricular direita. No entanto, confirmamos a presença do fenômeno da interdependência ventricular.

Análise combinada de fatores genéticos e ambientais na hipertensão essencial em um município da região Amazônica

OBJETIVO: Este estudo avaliou a contribuição de seis polimorfismos genéticos presentes em genes do sistema renina-angiotensina-aldosterona (SRAA) e fatores de risco clínicos para o desenvolvimento da hipertensão arterial essencial em um município da região Amazônica. MÉTODOS: Oitenta e dois indivíduos hipertensos e setenta e oito indivíduos normotensos foram genotipados quanto à presença de polimorfismos REN-G1051A (renina), AGT-M235T (angiotensinogênio), ECA-Alu I/D (enzima conversora de angiotensina I), AGTR1-A1166C (receptor tipo 1 da angiotensina II) e CYP11B2-C344T (aldosterona sintetase) pela técnica de reação em cadeia da polimerase (PCR), com análise de restrição quando necessário. A influência de polimorfismos genéticos e fatores de risco clínicos na variação da pressão arterial foi avaliada por meio de regressão linear stepwise. RESULTADOS: Relatamos a co-ocorrência de fatores de risco clínicos e polimorfismo do gene da enzima conversora de angiotensina (ECA) na população de um município da região amazônica. Nossos resultados mostram que a elevação da pressão arterial sistólica é favorecida pelo alelo D do polimorfismo de inserção/deleção do gene da ECA e pelo aumento da idade, enquanto consumo de bebida alcoólica e envelhecimento estão associados ao aumento da pressão arterial diastólica (PAD). CONCLUSÃO: Esses achados indicam que os moradores de Santa Isabel do Rio Negro que possuem o alelo D da ECA ou têm o hábito de beber apresentam valores mais elevados de PAS e PAD, respectivamente, com o passar dos anos.

Efeitos agudos do tabagismo sobre a modulação autonômica: análise por meio do plot de poincaré

FUNDAMENTO: O tabagismo altera a função autonômica. OBJETIVO: Investigar os efeitos agudos do tabagismo sobre a modulação autonômica e a recuperação dos índices de variabilidade de frequência cardíaca (VFC) pós-fumo, por meio do plot de Poincaré e índices lineares. MÉTODOS: Foram avaliados 25 fumantes jovens, os quais tiveram a frequência cardíaca analisada, batimento a batimento, na posição sentada, após 8 horas de abstinência, por 30 minutos em repouso, 20 minutos durante o fumo e 30 minutos pós-fumo. Análise de variância para medidas repetidas, seguido do teste de Tukey, ou teste de Friedman seguido do teste de Dunn foram aplicados dependendo da normalidade dos dados, com p < 0,05. RESULTADOS: Durante o fumo, houve redução dos índices SD1 (23,4 ± 9,2 vs 13,8 ± 4,8), razão SD1/SD2 (0,31 ± 0,08 vs 0,2 ± 0,04), RMSSD (32,7 ± 13 vs 19,1 ± 6,8), SDNN (47,6 ± 14,8 vs 35,5 ± 8,4), HFnu (32,5 ± 11,6 vs 19 ± 8,1) e do intervalo RR (816,8 ± 89 vs 696,5 ± 76,3) em relação ao repouso, enquanto que aumentos do índice LFnu (67,5 ± 11,6 vs 81 ± 8,1) e da razão LF/HF (2,6 ± 1,7 vs 5,4 ± 3,1) foram observados. A análise visual do plot mostrou menor dispersão dos intervalos RR durante o fumo. Com exceção da razão SD1/SD2, os demais índices apresentaram recuperação dos valores, 30 minutos após o tabagismo. CONCLUSÃO: O tabagismo produziu agudamente modificações no controle autonômico, caracterizadas por ativação simpática e retirada vagal, com recuperação 30 minutos após o fumo.

Análise da tendência da mortalidade por acidente vascular cerebral no Brasil no século XXI

FUNDAMENTO: Embora seja mundialmente a segunda principal causa de óbitos, o Acidente Vascular Cerebral (AVC) vem apresentando uma importante redução das taxas de mortalidade nas últimas décadas. OBJETIVO: Avaliar a tendência da taxa de mortalidade por acidente vascular cerebral no Brasil, em ambos os sexos, a partir dos 30 anos de idade, entre 2000 e 2009. MÉTODOS: Os dados populacionais foram obtidos no banco de dados do Instituto Brasileiro de Geografia e Estatística e os óbitos, por meio do Sistema de Informações sobre Mortalidade da Secretaria de Vigilância em Saúde do Ministério da Saúde, sendo incluídos os códigos I60 a I69 de acordo com a 10ª Classificação Internacional de Doenças. Foi calculada a incidência de óbitos/1.000 habitantes, as taxas de mortalidade bruta e padronizada/100.000 habitantes. A modelagem da tendência das taxas foi feita com modelos de regressão. RESULTADOS: Observou-se um aumento na incidência de óbitos até 2006, seguindo-se um declínio até 2009, quando ocorreu a incidência mínima. Comparando os anos 2000 e 2009, nota-se uma tendência de queda da taxa de mortalidade padronizada em ambos os sexos (masculino = -14,69%; feminino = -17%) e no total (-14,99%), com oscilações no período. Entre 30 e 49 anos em ambos os sexos, houve uma tendência de redução contínua e linear da taxa de mortalidade, enquanto os demais grupos etários apresentaram uma função curvilínea, culminando com uma efetiva diminuição dos valores. CONCLUSÃO: Houve uma tendência de queda na taxa de mortalidade em todas as faixas etárias e em ambos os sexos. A redução da taxa de mortalidade bruta foi mais acentuada no sexo masculino, enquanto a taxa de mortalidade padronizada mostrou uma maior redução no sexo feminino.

Macroecologia de carnívoros do Novo Mundo (Mammalia): envelopes de restrição e análise de padrões filogenéticos

The relationship between body size and geographic range was analyzed for 70 species of terrestrial Carnivora ("fissipeds") of the New World, after the control of phylogenetic patterns in the data using phylogenetic eigenvector regression. The analysis from EcoSim software showed that the variables are related as a triangular envelope. Phylogenetic patterns in data were detected by means of phylogenetic correlograms, and 200 simulations of the phenotypic evolution were also performed over the phylogeny. For body size, the simulations suggested a non-linear relationship for the evolution of this character along the phylogeny. For geographic range size, the correlogram showed no phylogenetic patterns. A phylogenetic eigenvector regression was performed on original data and on data simulated under Ornstein-Uhlenbeck process. Since both characters did not evolve under a simple Brownian motion process, the Type I errors should be around 10%, compatible with other methods to analyze correlated evolution. The significant correlation of the original data (r = 0.38; P < 0.05), as well as the triangular envelope, then indicate ecological and adaptive processes connecting the two variables, such as those proposed in minimum viable population models.