Métodos robustos em geoestatística

O objectivo principal da presente tese consiste no desenvolvimento de estimadores robustos do variograma com boas propriedades de eficiência. O variograma é um instrumento fundamental em Geoestatística, pois modela a estrutura de dependência do processo em estudo e influencia decisivamente a predição de novas observações. Os métodos tradicionais de estimação do variograma não são robustos, ou seja, são sensíveis a pequenos desvios das hipóteses do modelo. Essa questão é importante, pois as propriedades que motivam a aplicação de tais métodos, podem não ser válidas nas vizinhanças do modelo assumido. O presente trabalho começa por conter uma revisão dos principais conceitos em Geoestatística e da estimação tradicional do variograma. De seguida, resumem-se algumas noções fundamentais sobre robustez estatística. No seguimento, apresenta-se um novo método de estimação do variograma que se designou por estimador de múltiplos variogramas. O método consiste em quatro etapas, nas quais prevalecem, alternadamente, os critérios de robustez ou de eficiência. A partir da amostra inicial, são calculadas, de forma robusta, algumas estimativas pontuais do variograma; com base nessas estimativas pontuais, são estimados os parâmetros do modelo pelo método dos mínimos quadrados; as duas fases anteriores são repetidas, criando um conjunto de múltiplas estimativas da função variograma; por fim, a estimativa final do variograma é definida pela mediana das estimativas obtidas anteriormente. Assim, é possível obter um estimador que tem boas propriedades de robustez e boa eficiência em processos Gaussianos. A investigação desenvolvida revelou que, quando se usam estimativas discretas na primeira fase da estimação do variograma, existem situações onde a identificabilidade dos parâmetros não está assegurada. Para os modelos de variograma mais comuns, foi possível estabelecer condições, pouco restritivas, que garantem a unicidade de solução na estimação do variograma. A estimação do variograma supõe sempre a estacionaridade da média do processo. Como é importante que existam procedimentos objectivos para avaliar tal condição, neste trabalho sugere-se um teste para validar essa hipótese. A estatística do teste é um estimador-MM, cuja distribuição é desconhecida nas condições de dependência assumidas. Tendo em vista a sua aproximação, apresenta-se uma versão do método bootstrap adequada ao estudo de observações dependentes de processos espaciais. Finalmente, o estimador de múltiplos variogramas é avaliado em termos da sua aplicação prática. O trabalho contém um estudo de simulação que confirma as propriedades estabelecidas. Em todos os casos analisados, o estimador de múltiplos variogramas produziu melhores resultados do que as alternativas usuais, tanto para a distribuição assumida, como para distribuições contaminadas.

Estimação robusta em modelos lineares de equações simultâneas

Os Modelos de Equações Simultâneas (SEM) são modelos estatísticos com muita tradição em estudos de Econometria, uma vez que permitem representar e estudar uma vasta gama de processos económicos. Os estimadores mais usados em SEM resultam da aplicação do Método dos Mínimos Quadrados ou do Método da Máxima Verosimilhança, os quais não são robustos. Em Maronna e Yohai (1997), os autores propõem formas de “robustificar” esses estimadores. Um outro método de estimação com interesse nestes modelos é o Método dos Momentos Generalizado (GMM), o qual também conduz a estimadores não robustos. Estimadores que sofrem de falta de robustez são muito inconvenientes uma vez que podem conduzir a resultados enganadores quando são violadas as hipóteses subjacentes ao modelo assumido. Os estimadores robustos são de grande valor, em particular quando os modelos em estudo são complexos, como é o caso dos SEM. O principal objectivo desta investigação foi o de procurar tais estimadores tendo-se construído um estimador robusto a que se deu o nome de GMMOGK. Trata-se de uma versão robusta do estimador GMM. Para avaliar o desempenho do novo estimador foi feito um adequado estudo de simulação e foi também feita a aplicação do estimador a um conjunto de dados reais. O estimador robusto tem um bom desempenho nos modelos heterocedásticos considerados e, nessas condições, comporta-se melhor do que os estimadores não robustos usados no estudo. Contudo, quando a análise é feita em cada equação separadamente, a especificidade de cada equação individual e a estrutura de dependência do sistema são dois aspectos que influenciam o desempenho do estimador, tal como acontece com os estimadores usuais. Para enquadrar a investigação, o texto inclui uma revisão de aspectos essenciais dos SEM, o seu papel em Econometria, os principais métodos de estimação, com particular ênfase no GMM, e uma curta introdução à estimação robusta.

Utilização da estatística univariada e multivariada na interpretação de dados obtidos em sedimentos de linhas de água e solos na área Ferreira do Zêzere, Vila do Rei Santiago de Montalegre

Resumo indisponível

Contributions to the theory of maximum entropy estimation for ill-posed models

As técnicas estatísticas são fundamentais em ciência e a análise de regressão linear é, quiçá, uma das metodologias mais usadas. É bem conhecido da literatura que, sob determinadas condições, a regressão linear é uma ferramenta estatística poderosíssima. Infelizmente, na prática, algumas dessas condições raramente são satisfeitas e os modelos de regressão tornam-se mal-postos, inviabilizando, assim, a aplicação dos tradicionais métodos de estimação. Este trabalho apresenta algumas contribuições para a teoria de máxima entropia na estimação de modelos mal-postos, em particular na estimação de modelos de regressão linear com pequenas amostras, afetados por colinearidade e outliers. A investigação é desenvolvida em três vertentes, nomeadamente na estimação de eficiência técnica com fronteiras de produção condicionadas a estados contingentes, na estimação do parâmetro ridge em regressão ridge e, por último, em novos desenvolvimentos na estimação com máxima entropia. Na estimação de eficiência técnica com fronteiras de produção condicionadas a estados contingentes, o trabalho desenvolvido evidencia um melhor desempenho dos estimadores de máxima entropia em relação ao estimador de máxima verosimilhança. Este bom desempenho é notório em modelos com poucas observações por estado e em modelos com um grande número de estados, os quais são comummente afetados por colinearidade. Espera-se que a utilização de estimadores de máxima entropia contribua para o tão desejado aumento de trabalho empírico com estas fronteiras de produção. Em regressão ridge o maior desafio é a estimação do parâmetro ridge. Embora existam inúmeros procedimentos disponíveis na literatura, a verdade é que não existe nenhum que supere todos os outros. Neste trabalho é proposto um novo estimador do parâmetro ridge, que combina a análise do traço ridge e a estimação com máxima entropia. Os resultados obtidos nos estudos de simulação sugerem que este novo estimador é um dos melhores procedimentos existentes na literatura para a estimação do parâmetro ridge. O estimador de máxima entropia de Leuven é baseado no método dos mínimos quadrados, na entropia de Shannon e em conceitos da eletrodinâmica quântica. Este estimador suplanta a principal crítica apontada ao estimador de máxima entropia generalizada, uma vez que prescinde dos suportes para os parâmetros e erros do modelo de regressão. Neste trabalho são apresentadas novas contribuições para a teoria de máxima entropia na estimação de modelos mal-postos, tendo por base o estimador de máxima entropia de Leuven, a teoria da informação e a regressão robusta. Os estimadores desenvolvidos revelam um bom desempenho em modelos de regressão linear com pequenas amostras, afetados por colinearidade e outliers. Por último, são apresentados alguns códigos computacionais para estimação com máxima entropia, contribuindo, deste modo, para um aumento dos escassos recursos computacionais atualmente disponíveis.

Assessing the combined effect of land management and wildfire on runoff and soil erosion in north central Portugal

Strong and sometimes extreme responses in runoff and soil erosion following wildfires have been reported worldwide. However, in the case of North-Central Portugal, little research had been carried out regarding the hydrologic and erosive impacts of several land management activities in recently burnt areas (such as ground preparation, post-fire logging or post-fire mitigation treatments). This study aims to assess post-fire runoff and soil erosion response on Eucalypt and Maritime pine plantations during the first, second and third years following wildfires. The effect of several pre-fire ground preparation operations (ploughed down-slope, contour ploughed and inclined terraces), post-fire logging activities (on both the eucalypt and pine plantations), as well as the application of hydromulch (a post-fire emergency treatment) on overland flow and soil erosion were compared to burnt but undisturbed and untreated areas. The intensive monitoring of runoff, soil erosion and selected soil properties served to determine the main factors involved in post-fire runoff and soil erosion and their spatial and temporal variation. Soil water repellency deserved special attention, due to its supposed important role for overland flow generation. Repeated rainfall simulation experiments (RSE’s), micro-scale runoff plots and bounded sediment fences were carried out and/or installed immediately after the wildfire on seven burnt slopes. Micro-scale runoff plots results under natural rainfall conditions were also compared to the RSE’s results, which was useful for assessing the representativeness of the data obtained with artificial rainfall. The results showed comparable runoff coefficient (20-60%) but lower sediment losses (125-1000 g m-2) than prior studies in Portugal, but especially outside Portugal. Lower sediment losses were related with the historic intensive land use in the area. In evaluating these losses, however, the shallowness and stoniness of the soils, as well as the high organic matter fraction of the eroded sediments (50%) must not be overlooked. Sediment limited erosion was measured in all the ploughed sites, probably due to the time since ploughing (several years). The disturbance of the soil surface cover due to post-fire logging and wood extraction substantially increased sediment losses at both the pine and eucalypt sites. Hydromulch effectiveness in reducing the runoff (70%) and sediment losses (83%) was attributed to the protective high coverage provided by hydromulch. The hydromulch significantly affected the soil cover and other soil properties and these changes also reduced the soil erosion risk. The rainfall amount was the main factor explaining the variance in runoff. However, a shift from rainfall amount to rainfall intensity was detected when either the surface cover or the infiltration capacity (hydrophilic conditions) increased. Sediment losses were controlled by rainfall intensity and surface cover. The role of soil water repellency on runoff generation was not consistent; the overall repellency levels alone were not enough to assess its hydrological impact. Soil water repellency explained runoff generation in the specific-sites model better than in the overall model. Additionally, soil moisture content was a better predictor for soil water repellency than antecedent rainfall. The natural rainfall results confirmed that RSE’s were able to capture the specific sediment losses and its organic matter content as well as the differences between the ploughed and unploughed sites. Repeated RSE’s also captured the seasonal variations in runoff and sediment losses attributed to soil water repellency. These results have implications for post-fire soil erosion modelling and soil conservation practices in the region, or areas with the same land use, climate and soil characteristics. The measured sediment loss, as well as the increasing frequency of ploughing in recently burnt and unburnt eucalypt stands, suggests ploughing is not an effective as a soil conservation measure. Logging activities with less impact are recommended in order to maintain the forest litter protecting the soil surface. Due to its high effectiveness in reducing runoff and soil erosion, hydromulch is recommended for highly sensitive and vulnerable areas.

Channel estimation and parameters acquisition systems employing cooperative diversity

This work investigates new channel estimation schemes for the forthcoming and future generation of cellular systems for which cooperative techniques are regarded. The studied cooperative systems are designed to re-transmit the received information to the user terminal via the relay nodes, in order to make use of benefits such as high throughput, fairness in access and extra coverage. The cooperative scenarios rely on OFDM-based systems employing classical and pilot-based channel estimators, which were originally designed to pointto-point links. The analytical studies consider two relaying protocols, namely, the Amplifyand-Forward and the Equalise-and-Forward, both for the downlink case. The relaying channels statistics show that such channels entail specific characteristics that comply to a proper filter and equalisation designs. Therefore, adjustments in the estimation process are needed in order to obtain the relay channel estimates, refine these initial estimates via iterative processing and obtain others system parameters that are required in the equalisation. The system performance is evaluated considering standardised specifications and the International Telecommunication Union multipath channel models.

Liderança, virtuosidade e criatividade explicando o desempenho de equipes

As organizações estão se tornando cada vez mais competitivas e buscando, constantemente, formas de aumentar o seu desempenho através da compreensão dos fatores que o antecedem. A nossa investigação incorpora dois estudos, ambos incidentes sobre as lojas de uma grande cadeia de varejo brasileira. O primeiro estudo, envolvendo 68 lojas, investiga como a liderança autêntica explica o desempenho das lojas através do efeito mediador da virtuosidade e da potência das equipes (uma espécie de sentido de eficácia coletiva). O segundo estudo, abrangendo 46 lojas, testa como a potência das equipes explica o respetivo desempenho através do efeito mediador da criatividade. Os principais resultados são: (a) a relação entre a liderança autêntica e o desempenho das lojas é totalmente mediada pela virtuosidade e pela potência das equipes; (b) a criatividade medeia a relação entre a potência das equipes e o desempenho das lojas; (c) a criatividade não explica diretamente o desempenho, mas parece influenciá-lo através do efeito mediador da potência. Os cuidados metodológicos destinados a diminuir os riscos de variância do método comum aumentam a validade dos estudos. Ao focar variáveis relacionadas com a virtuosidade, a liderança autêntica, e as forças das equipes, e ao incluir dados de desempenho objetivo, a nossa investigação é um contributo relevante para o movimento dos “Estudos Organizacionais Positivos”.

Modelos autoregressivos por limiares em séries de contagem

A modelação e análise de séries temporais de valores inteiros têm sido alvo de grande investigação e desenvolvimento nos últimos anos, com aplicações várias em diversas áreas da ciência. Nesta tese a atenção centrar-se-á no estudo na classe de modelos basedos no operador thinning binomial. Tendo como base o operador thinning binomial, esta tese focou-se na construção e estudo de modelos SETINAR(2; p(1); p(2)) e PSETINAR(2; 1; 1)T , modelos autorregressivos de valores inteiros com limiares autoinduzidos e dois regimes, admitindo que as inovações formam uma sucessão de variáveis independentes com distribuição de Poisson. Relativamente ao primeiro modelo analisado, o modelo SETINAR(2; p(1); p(2)), além do estudo das suas propriedades probabilísticas e de métodos, clássicos e bayesianos, para estimar os parâmetros, analisou-se a questão da seleção das ordens, no caso de elas serem desconhecidas. Com este objetivo consideraram-se algoritmos de Monte Carlo via cadeias de Markov, em particular o algoritmo Reversible Jump, abordando-se também o problema da seleção de modelos, usando metodologias clássica e bayesiana. Complementou-se a análise através de um estudo de simulação e uma aplicação a dois conjuntos de dados reais. O modelo PSETINAR(2; 1; 1)T proposto, é também um modelo autorregressivo com limiares autoinduzidos e dois regimes, de ordem unitária em cada um deles, mas apresentando uma estrutura periódica. Estudaram-se as suas propriedades probabilísticas, analisaram-se os problemas de inferência e predição de futuras observações e realizaram-se estudos de simulação.

Um modelo de medição e gestão de desempenho para as universidades públicas

Nas últimas décadas as instituições de ensino superior têm sido alvo de uma pressão crescente para aumentar a sua eficiência e a sua eficácia. Fatores como a diversificação da oferta formativa, a massificação, o aumento da internacionalização, entre outros, conduziram a uma maior monitorização das instituições e, por sua vez, geraram o desenvolvimento de novos sistemas de informação. A constante solicitação de informação, quer pelo Estado, quer pelo mercado constitui a base para a definição do objetivo desta investigação: construir um modelo integrado de medição e gestão de desempenho para as universidades públicas e testá-lo no universo português. Para a construção conceptual deste modelo foi realizada uma revisão da literatura baseada em diferentes contextos: organizacional, europeu e nacional. Este modelo foi depois aplicado às universidades públicas portuguesas recorrendo a dados disponibilizados por entidades públicas, tendo em conta a área de educação e formação (CNAEF). Uma análise descritiva aos dados constituiu uma contribuição desta investigação para a prática, no sentido de que permitiu fornecer um conjunto de recomendações às universidades e às entidades oficiais responsáveis pela recolha de dados a nível nacional, relativamente aos sistemas de informação e processos de recolha de dados. O modelo proposto constitui a contribuição teórica desta investigação, no sentido de que integra, no mesmo modelo, as diferentes dimensões de desempenho, a visão dos diferentes stakeholders, quer no contexto do ensino, quer no da investigação, quer no da própria gestão da universidade. A vertente analítica deste modelo, representada pelas diferentes relações entre os grupos de indicadores, foi testada recorrendo à técnica de análise de equações estruturais baseada na variância (nomeadamente Partial Least Squares) em quatro áreas CNAEF. Os resultados demonstraram que o comportamento em termos de medição e gestão de desempenho difere consoante a CNAEF, identificando um maior ajustamento às áreas hard e evidenciando que as áreas soft necessitam de indicadores mais ajustados às suas especificidades.

Sustainable viticulture in Bairrada appellation: vineyard and harvest year effects on grapes oenological potential

The sustainable viticulture of a region passes, among other aspects, for maximizing the varieties potential minimizing subsequent interventions during winemaking, which should contribute to the production of quality wines maintaining their typicity and rationalizing costs. The detailed knowledge of each Appellation specificities, namely vineyard parcel (soil type and topographical peculiarities) and harvest climatic conditions is crucial for sustainability in this sector. Thus, in line with this current trend, the starting point for the development of this PhD thesis was to evaluate the oenological potential of different varieties cultivated throughout Bairrada Appellation (Portugal). During maturation several changes in grape varieties occur, namely berries become sweeter, less acidic, and they develop flavour, aroma and colour properties. The development of these characteristics is essential to define grapes oenological potential, i.e. to estimate the possibility of their usage to produce specific wines. A three years sampling plan was designed to evaluate the effect of harvest year and parcel characteristics on V. vinifera cv. Arinto, Bical, Sauvignon Blanc, Baga, Castelão, Touriga Nacional, and Sousão grapes composition. For each variety, 3 parcels with different characteristics were selected. Several physicochemical parameters were evaluated, during maturation: berry weight, pH, titratable acidity, sugar and phenolic contents, antiradical activity, and volatile composition (free fraction). Special attention was devoted to grapes at technologic maturity, since, besides these parameters, glycosidically-linked fraction was also considered. By using the results obtained at technologic maturity, a comprehensive approach was applied to identify the significance of harvest and parcel characteristics effects on each variety composition. Considering all the parameters under study, it may be highlighted some significant differences. According to the obtained results determined during maturation, it was possible to conclude that Arinto, Bical and Sauvignon Blanc grapes from parcels with clay-sandy and clay-calcareous soils have higher phenolic content and antiradical activity. Otherwise, Sauvignon Blanc presented similar volatile composition for grapes cultivated in the 3 parcels, while Arinto and Bical exhibited higher volatile content in grapes from claysandy and clay-calcareous soils. For Baga, Castelão and Touriga Nacional red varieties, grapes with higher phenolic content, antiradical activity, and volatile content were obtained from clayey and clay-calcareous soils. Furthermore, for Touriga Nacional, parcels altitude seems also to modulate grapes composition. Beyond parcel effect, harvest year conditions also influence grapes composition: 2011 harvest was related with lower phenolic and volatile contents, as well as lower antiradical activity.For grapes collected at technologic maturity, analysis of variance-simultaneous component analysis (ASCA) was applied combining all the parameters under study, in order to assess the influence of harvest and parcel characteristics on each variety oenological potential. The results obtained using this comprehensive approach is closely related with those observed during maturation and revealed that harvest was the main factor that influenced grapes composition (53% to 68% of the total data set variance) followed by parcel characteristics, explaining ca. 15-19% of the total data set variance. The oenological potential of each variety may be different from one parcel to another, i.e., clay-sandy and clay-calcareous related-environments seem to favour Arinto and Bical white grapes composition, but for the red varieties, grapes composition was favoured by clayey and clay-calcareous soils. Besides, also higher altitude seems to favour Touriga Nacional grapes composition. Sauvignon Blanc seems to be a variety well adapted to the different parcel characteristics. In order to go forward in the valuation of these varieties, the aroma properties of 6 monovarietal wines were studied based on an aroma network-approach, linking molecular data related to volatile composition and aroma data about the key odor active molecules. This approach allowed to identify different wine aroma properties and to infer about the consumer’s sensory perception. It was found that aroma properties differ from one wine variety to another: while Arinto and Sauvignon Blanc wine exhibited higher tree fruity, sweety and flowery aromas, related essentially with ester compounds and C13 norisoprenoids, the opposite was obtained for Bical wine, corroborating the aroma sensory perceptions of the trained panel. Sauvignon Blanc also exhibited higher toasted aromas (related with thiols, mainly with 2-methyl-3-furanthiol). Touriga Nacional red wine exhibited higher tree, tropical, and berry fruits notes (sensory described as sweet fruits), toasted and flowery aromas, while these are similar for the other red wines under study. Besides Portuguese Bairrada wines, this aroma network approach is a tool that can be used to explain the aroma properties of wines worldwide. The grape and wine data generated under the present PhD thesis, in the context of Bairrada Appellation, shows the unique character of each variety, and may be used by growers and wine producers as a support for decisionmaking based on objective criteria, increasing the sustainability in this sector. For instance, it is possible to take advantage of the natural resources and produce products with different characteristics obtained from the same variety, minimizing costs during the winemaking process.

Fiabilidade e validade do instrumento Physiosensing na avaliação do equilíbrio em pessoas com deficiência intelectual

Enquadramento: O ‘Physiosensing’ (‘PhyS’) é um dispositivo médico destinado ao treino do controlo postural nas posições de sentado e de pé, bem como no levantar e sentar, possibilitando também a avaliação do desempenho a este nível. Este trabalho teve como objetivo avaliar a fiabilidade e a validade da plataforma ‘Physiosensing’ na avaliação do equilíbrio em pessoas com deficiência intelectual (PCDI). Métodos: Para o grupo experimental (GE) foram recrutados 47 indivíduos com deficiência intelectual e para o grupo de controlo (GC) 39 indivíduos sem deficiência intelectual, provenientes da região do BaixoMondego. A avaliação da fiabilidade incluiu as análises da concordância entre observadores, reprodutibilidade temporal e consistência interna. A análise fatorial exploratória analisou os pressupostos de subdomínios propostos pelos autores. A validade discriminante foi analisada através da comparação de resultados entre o GE e o GC, e a validade concorrente pela análise dos valores de associação entre os resultados do ‘PhyS’ com a Escala de Equilíbrio de Berg (EEB). Resultados: O subconjunto de exercícios relacionados com a transferência de peso lateralmente (TPL) apresentou os resultados mais elevados a nível da concordância entre observadores (0,40 ≤ CCI > 0,75) e na reprodutibilidade intemporal (CCI ≥ 0,75). O instrumento apresenta uma consistência interna fraca (α = 0,63) quando considerados todos os exercícios, tendo-se obtido o melhor resultado para o subconjunto de exercícios TPL (α = 0,81). A análise fatorial exploratória devolveu quatro fatores, explicando 76,4% da variância, agrupando no primeiro fator o subconjunto de exercícios TPL. Foram encontradas diferenças estatisticamente significativas entre os resultados dos participantes com e sem deficiência intelectual, em dez dos onze exercícios que compõem a configuração da plataforma. Seis exercícios, que incluem os exercícios TPL, apresentam valores de associação estatisticamente significativos com a EEB. Conclusão: Alguns exercícios da plataforma ‘PhyS’ não se mostram adequados para medir o equilíbrio em PCDI, não podendo ser incluídos numa medida global. No entanto, os exercícios TPL poderão constituir um indicador global do equilíbrio em pessoas com PCDI. Recomenda-se a definição de procedimentos de medição de forma a melhorar os índices de fiabilidade, o aprofundamento da configuração de exercícios para a avaliação do equilíbrio e o estudo do potencial da plataforma em programas de intervenção para o treino das funções do equilíbrio.