Métodos de adequação e diagnóstico em modelos de sobrevivência dinâmicos

A análise de dados de sobrevivência tem sido tradicionalmente baseada no modelo de regressão de Cox (COX, 1972). No entanto, a suposição de taxas de falha proporcionais assumida para esse modelo pode não ser atendida em diversas situações práticas. Essa restrição do modelo de Cox tem gerado interesse em abordagens alternativas, dentre elas os modelos dinâmicos que permitem efeito das covariáveis variando no tempo. Neste trabalho, foram revisados os principais modelos de sobrevivência dinâmicos com estrutura aditiva e multiplicativa nos contextos não paramétrico e semiparamétrico. Métodos gráficos baseados em resíduos foram apresentados com a finalidade de avaliar a qualidade de ajuste desses modelos. Uma versão tempo-dependente da área sob a curva ROC, denotada por AUC(t), foi proposta com a finalidade de avaliar e comparar a qualidade de predição entre modelos de sobrevivência com estruturas aditiva e multiplicativa. O desempenho da AUC(t) foi avaliado por meio de um estudo de simulação. Dados de três estudos descritos na literatura foram também analisados para ilustrar ou complementar os cenários que foram considerados no estudo de simulação. De modo geral, os resultados obtidos indicaram que os métodos gráficos apresentados para avaliar a adequação dos modelos em conjunto com a AUC(t) se constituem em um conjunto de ferramentas estatísticas úteis para o próposito de avaliar modelos de sobrevivência dinâmicos nos contextos não paramétrico e semiparamétrico. Além disso, a aplicação desse conjunto de ferramentas em alguns conjuntos de dados evidenciou que se, por um lado, os modelos dinâmicos são atrativos por permitirem covariáveis tempo-dependentes, por outro lado podem não ser apropriados para todos os conjuntos de dados, tendo em vista que estimação pode apresentar restrições para alguns deles.

Modelos não lineares resultantes da soma de regressões lineares ponderadas por funções distribuição acumulada

Os controladores eletrônicos de pulverização visam minimizar a variação das taxas de insumos aplicadas no campo. Eles fazem parte de um sistema de controle, e permitem a compensação da variação de velocidade de deslocamento do pulverizador durante a operação. Há vários tipos de controladores eletrônicos de pulverização disponíveis no mercado e uma das formas de selecionar qual o mais eficiente nas mesmas condições, ou seja, em um mesmo sistema de controle, é quantificar o tempo de resposta do sistema para cada controlador específico. O objetivo desse trabalho foi estimar os tempos de resposta para mudanças de velocidade de um sistema eletrônico de pulverização via modelos de regressão não lineares, estes, resultantes da soma de regressões lineares ponderadas por funções distribuição acumulada. Os dados foram obtidos no Laboratório de Tecnologia de Aplicação, localizado no Departamento de Engenharia de Biossistemas da Escola Superior de Agricultura \"Luiz de Queiroz\", Universidade de São Paulo, no município de Piracicaba, São Paulo, Brasil. Os modelos utilizados foram o logístico e de Gompertz, que resultam de uma soma ponderada de duas regressões lineares constantes com peso dado pela função distribuição acumulada logística e Gumbell, respectivamente. Reparametrizações foram propostas para inclusão do tempo de resposta do sistema de controle nos modelos, com o objetivo de melhorar a interpretação e inferência estatística dos mesmos. Foi proposto também um modelo de regressão não linear difásico que resulta da soma ponderada de regressões lineares constantes com peso dado pela função distribuição acumulada Cauchy seno hiperbólico exponencial. Um estudo de simulação foi feito, utilizando a metodologia de Monte Carlo, para avaliar as estimativas de máxima verossimilhança dos parâmetros do modelo.

Efeitos de dois modelos experimentais de dietas hipercalóricas sobre parâmetros comportamentais e biológicos em ratos da linhagem Wistar

O objetivo deste estudo foi investigar os efeitos de dois modelos experimentais de dietas hipercalóricas em comportamentos de ansiedade, processos de aprendizagem e memória e alterações metabólicas. Os animais foram divididos em seis grupos experimentais, de acordo com a condição nutricional. 1) Controle (C); 2) Dieta de Cafeteria (DC); 3) Dieta Hiperlipídica (DH); 4) Controle AIN-93 (C/AIN-93); 5) Dieta de Cafeteria AIN-93 (DC/AIN-93), e 6) Dieta Hiperlipídica AIN-93 (DH/AIN-93). Posteriormente, os grupos foram subdivididos em dois grupos independentes, conforme a tarefa à qual foram submetidos. Pesagens foram realizadas semanalmente até os 98 dias de vida; foram verificados os pesos do fígado, do coração e o peso de tecido adiposo retroperitoneal e epididimal e foram realizadas dosagens de glicose, triglicérides, TGO e TGP no soro e gordura total, colesterol total e triglicérides no fígado. Os testes utilizados: Labirinto em T Elevado (LTE), Caixa Claro/Escuro e Labirinto Aquático de Morris (LAM). Os resultados de peso corporal, os dados comportamentais do LAM, do LTE e os dados de peso dos tecidos extraídos no dia do sacrifício e as análises bioquímicas foram submetidos a uma Análise de Variância (ANOVA). Quando apropriado, foi utilizado o teste de comparações múltiplas de Newman-Keuls (p< 0,05). Os dados comportamentais do teste claro/escuro foram submetidos ao teste t-Student (p< 0,05). Animais tratados com dieta hiperlipídica apresentaram maiores medidas de peso e ganho de peso comparados aos animais controle e dieta de cafeteria, tratados com pellet e com dieta AIN-93. Animais DH1, DC1, DH1 AIN-93, DH2 AIN-93 e DH2 apresentaram maior peso no dia do sacrifício. Animais DH1, DH1 AIN-93, DH2 e DH2 AIN-93 apresentaram maior acúmulo dos tecidos adiposos retroperitoneal e epididimal. Animais DH1 AIN-93 e DC2 AIN-93 apresentaram maiores níveis de glicose. Animais C2, DH2 e DC2 apresentaram maiores níveis de triglicérides. Animais DH1 e C1 apresentaram menores valores de TGO. Animais C2 e C2 AIN-93 apresentaram maiores níveis de TGO. Animais C1, DH1, C2 e DH2 apresentaram maiores níveis de TGP. Animais DH1 AIN-93, DH1, DH2 e DH2 AIN-93 apresentaram maiores valores de gordura total no fígado. Animais DH1 AIN-93 e DH2 apresentaram maiores níveis de colesterol no fígado. Animais DH1, DC1, DH2 e DH2 AIN-93 apresentaram maiores níveis de triglicérides no fígado. Com relação ao consumo alimentar, animais DH apresentaram maior consumo calórico e maior consumo lipídico quando comparados aos animais C e DC, com ração em pellet ou dieta AIN-93. Com relação ao LTE, não foram verificadas diferenças nas esquivas e na fuga. Animais DC1, DH1 e DH1 AIN-93 apresentaram menores níveis de ansiedade verificados a partir dos dados do teste da caixa claro-escuro. Animais DC2 AIN-93 apresentaram pior desempenho em tarefa de memória. Os dados obtidos a partir deste estudo demonstraram que as dietas utilizadas foram capazes de acarretar ganho de peso, acúmulo de tecido adiposo, alterações metabólicas, diminuição da ansiedade nos animais e pior desempenho em uma tarefa de memória em um dos grupos nutricionais.

Imputação de dados em experimentos multiambientais: novos algoritmos utilizando a decomposição por valores singulares

As análises biplot que utilizam os modelos de efeitos principais aditivos com inter- ação multiplicativa (AMMI) requerem matrizes de dados completas, mas, frequentemente os ensaios multiambientais apresentam dados faltantes. Nesta tese são propostas novas metodologias de imputação simples e múltipla que podem ser usadas para analisar da- dos desbalanceados em experimentos com interação genótipo por ambiente (G×E). A primeira, é uma nova extensão do método de validação cruzada por autovetor (Bro et al, 2008). A segunda, corresponde a um novo algoritmo não-paramétrico obtido por meio de modificações no método de imputação simples desenvolvido por Yan (2013). Também é incluído um estudo que considera sistemas de imputação recentemente relatados na literatura e os compara com o procedimento clássico recomendado para imputação em ensaios (G×E), ou seja, a combinação do algoritmo de Esperança-Maximização com os modelos AMMI ou EM-AMMI. Por último, são fornecidas generalizações da imputação simples descrita por Arciniegas-Alarcón et al. (2010) que mistura regressão com aproximação de posto inferior de uma matriz. Todas as metodologias têm como base a decomposição por valores singulares (DVS), portanto, são livres de pressuposições distribucionais ou estruturais. Para determinar o desempenho dos novos esquemas de imputação foram realizadas simulações baseadas em conjuntos de dados reais de diferentes espécies, com valores re- tirados aleatoriamente em diferentes porcentagens e a qualidade das imputações avaliada com distintas estatísticas. Concluiu-se que a DVS constitui uma ferramenta útil e flexível na construção de técnicas eficientes que contornem o problema de perda de informação em matrizes experimentais.

Expressão de metaloproteinases de matriz (MMPS) e de seus inibidores (TIMPS e RECK) em modelo de progressão tumoral de Câncer de mama e sua correlação com dados clínicos-patológicos

O câncer de mama é o tipo de câncer mais comumente detectado em mulheres de todo o mundo. Na maioria das pacientes, a causa de morte se deve, principalmente, à doença metastática que pode se desenvolver a partir do tumor primário. O processo metastático envolve uma complexa cascata de eventos, incluindo a quebra organizada dos componentes da matriz extracelular por metaloproteinases de matriz (MMPs). A atividade das MMPs é precisamente regulada por inibidores específicos, os inibidores teciduais das MMPs (TIMPs). Dado seu papel na progressão tumoral, níveis elevados de MMPs têm sido associados com prognóstico desfavorável para pacientes com câncer. Por outro lado, sendo os TIMPs proteínas multifuncionais, níveis elevados de TlMP-1 e de TIMP-2 correlacionam com agressividade do tumor e prognóstico ruim em diferentes tipos de câncer, incluindo o câncer de mama. O gene supressor de metástase RECK codifica uma glicoproteína de membrana capaz de inibir a invasão e a metástase tumoral através da regulação negativa da atividade de MMPs envolvidas em carcinogênese: MMP-2, MMP-9 e MMP-14 (MT1-MMP). A fim de analisar o papel das MMPs e de seus inibidores (TIMPs e RECK) na progressão tumoral do câncer de mama, o perfil de expressão destes genes foi detectado, através de ensaios de Real-Time PCR, em um painel de cinco linhagens celulares de carcinoma de mama humano com diferentes potenciais invasivos e metastáticos e em 72 amostras teciduais de tumores primários de mama e 30 amostras teciduais de borda normal adjacente ao tumor. O perfil de expressão protéica de RECK foi avaliado em 236 amostras de tumores primários de mama através de ensaios de Tissue Microarray. Além disso, a atividade proteolítica das MMPs foi detectada em ensaios de Zimografia. Os resultados obtidos indicam que a progressão do câncer de mama humano está relacionada com um aumento dos níveis de expressão das MMPs e de seus inibidores específicos. O aumento dos níveis de expressão dos TIMPs parece estar relacionado ao seu papel como proteína multifuncional que pode estar funcionando de maneira a promover, mais do que suprimir, a progressão tumoral. Níveis elevados da expressão protéica de RECK estão associados com pior prognóstico. No entanto, para pacientes em estádios clínicos avançados, altos níveis de expressão de RECK podem estar correlacionados com melhor prognóstico, dependendo do balanço MMP/inibidor. Os níveis de expressão das MMPs apresentaram correlação positiva em relação aos níveis de expressão de seus inibidores específicos, sugerindo a existência de fatores e vias de sinalização comuns envolvidas na regulação coordenada destes genes. Além disso, a síntese do inibidor pode estar relacionada a uma resposta celular ao aumento da expressão e atividade de proteases. O balanço transcricional enzima/inibidor favorece a enzima nas amostras tumorais e, de modo contrário, o inibidor específico nas amostras de borda normal, sugerindo o balanço como o principal fator na determinação da degradação da MEC em processos invasivos e metastáticos. Os resultados obtidos podem contribuir para um melhor entendimento da complexidade dos mecanismos envolvidos na metástase do câncer de mama.

Simulação e modelagem computacional de dados de espalhamento à baixos ângulos - enfoque em estruturas de alta simetria

Esta tese apresenta uma abordagem para a criação rápida de modelos em diferentes geometrias (complexas ou de alta simetria) com objetivo de calcular a correspondente intensidade espalhada, podendo esta ser utilizada na descrição de experimentos de es- palhamento à baixos ângulos. A modelagem pode ser realizada com mais de 100 geome- trias catalogadas em um Banco de Dados, além da possibilidade de construir estruturas a partir de posições aleatórias distribuídas na superfície de uma esfera. Em todos os casos os modelos são gerados por meio do método de elementos finitos compondo uma única geometria, ou ainda, compondo diferentes geometrias, combinadas entre si a partir de um número baixo de parâmetros. Para realizar essa tarefa foi desenvolvido um programa em Fortran, chamado de Polygen, que permite modelar geometrias convexas em diferentes formas, como sólidos, cascas, ou ainda com esferas ou estruturas do tipo DNA nas arestas, além de usar esses modelos para simular a curva de intensidade espalhada para sistemas orientados e aleatoriamente orientados. A curva de intensidade de espalhamento é calculada por meio da equação de Debye e os parâmetros que compõe cada um dos modelos, podem ser otimizados pelo ajuste contra dados experimentais, por meio de métodos de minimização baseados em simulated annealing, Levenberg-Marquardt e algorítmicos genéticos. A minimização permite ajustar os parâmetros do modelo (ou composição de modelos) como tamanho, densidade eletrônica, raio das subunidades, entre outros, contribuindo para fornecer uma nova ferramenta para modelagem e análise de dados de espalhamento. Em outra etapa desta tese, é apresentado o design de modelos atomísticos e a sua respectiva simulação por Dinâmica Molecular. A geometria de dois sistemas auto-organizado de DNA na forma de octaedro truncado, um com linkers de 7 Adeninas e outro com linkers de ATATATA, foram escolhidas para realizar a modelagem atomística e a simulação por Dinâmica Molecular. Para este sistema são apresentados os resultados de Root Mean Square Deviations (RMSD), Root Mean Square Fluctuations (RMSF), raio de giro, torção das hélices duplas de DNA além da avaliação das ligações de Hidrogênio, todos obtidos por meio da análise de uma trajetória de 50 ns.

Avaliações e testes dos métodos MASW e ReMi por meio do tratamento de dados empíricos e sintéticos em um programa de processamento e inversão desenvolvido em MATLAB e sua implicação em um problema geológico na Bacia de Taubaté

Os métodos de ondas superficiais com ênfase nas ondas Rayleigh foram utilizados como o núcleo desse trabalho de Doutorado. Inicialmente, as ondas Rayleigh foram modeladas permitindo o estudo de sensibilidade de suas curvas de dispersão sob diferentes configurações de parâmetros físicos representando diversos modelos de camadas, em que pôde ser observado parâmetros com maior e menor sensibilidade e também alguns efeitos provocados por baixas razões de Poisson. Além disso, na fase de inversão dos dados a modelagem das ondas Rayleigh foi utilizada para a construção da função objeto, que agregada ao método de mínimos quadrados, a partir do método de Levenberg-Marquardt, permitiu a implementação de um algoritmo de busca local responsável pela inversão de dados das ondas superficiais. Por se tratar de um procedimento de busca local, o algoritmo de inversão foi complementado por uma etapa de pré-inversão com a geração de um modelo inicial para que o procedimento de inversão fosse mais rápido e eficiente. Visando uma eficiência ainda maior do procedimento de inversão, principalmente em modelos de camadas com inversão de velocidades, foi implementado um algoritmo de pós-inversão baseado em um procedimento de tentativa e erro minimizando os valores relativos da raiz quadrada do erro quadrático médio (REQMr) da inversão dos dados. Mais de 50 modelos de camadas foram utilizados para testar a modelagem, a pré-inversão, inversão e pós-inversão dos dados permitindo o ajuste preciso de parâmetros matemáticos e físicos presentes nos diversos scripts implementados em Matlab. Antes de inverter os dados adquiridos em campo, os mesmos precisaram ser tratados na etapa de processamento de dados, cujo objetivo principal é a extração da curva de dispersão originada devido às ondas superficiais. Para isso, foram implementadas, também em Matlab, três metodologias de processamento com abordagens matemáticas distintas. Essas metodologias foram testadas e avaliadas com dados sintéticos e reais em que foi possível constatar as virtudes e deficiências de cada metodologia estudada, bem como as limitações provocadas pela discretização dos dados de campo. Por último, as etapas de processamento, pré-inversão, inversão e pós-inversão dos dados foram unificadas para formar um programa de tratamento de dados de ondas superficiais (Rayleigh). Ele foi utilizado em dados reais originados pelo estudo de um problema geológico na Bacia de Taubaté em que foi possível mapear os contatos geológicos ao longo dos pontos de aquisição sísmica e compará-los a um modelo inicial existente baseado em observações geomorfológicas da área de estudos, mapa geológico da região e informações geológicas globais e locais dos movimentos tectônicos na região. As informações geofísicas associadas às geológicas permitiram a geração de um perfil analítico da região de estudos com duas interpretações geológicas confirmando a suspeita de neotectônica na região em que os contatos geológicos entre os depósitos Terciários e Quaternários foram identificados e se encaixaram no modelo inicial de hemi-graben com mergulho para Sudeste.

Modelos de rede de transporte público e individual para estudos de acessibilidade em São Paulo.

O presente estudo tem como objetivo comparar a eficiência dos sistemas de transporte público e individual através da modelagem de redes em SIG. Foram revisados os conceitos de acessibilidade e indicadores de acessibilidade, modelos de rede de transportes, abordadas as etapas para criação dos modelos de rede, bases de dados utilizadas, configurações e atributos no sistema SIG. Os comparativos foram realizados através de dois níveis de detalhamento (simples e avançado) para as redes de transporte e elaborados índices de acessibilidade para os bairros do município de São Paulo através de parâmetros extraídos das redes modeladas. O estudo de caso consiste na medição de um índice de acessibilidade para empregos de baixa renda. Como resultado, os bairros centrais do município apresentam maior acessibilidade, porém, para o transporte público, alguns bairros fora da zona central também apresentaram alta acessibilidade devido à oferta de transporte público (metrô).