An EM-based Semi-Parametric Mixture Model Approach to the Regression Analysis of Competing-Risks Data

We consider a mixture model approach to the regression analysis of competing-risks data. Attention is focused on inference concerning the effects of factors on both the probability of occurrence and the hazard rate conditional on each of the failure types. These two quantities are specified in the mixture model using the logistic model and the proportional hazards model, respectively. We propose a semi-parametric mixture method to estimate the logistic and regression coefficients jointly, whereby the component-baseline hazard functions are completely unspecified. Estimation is based on maximum likelihood on the basis of the full likelihood, implemented via an expectation-conditional maximization (ECM) algorithm. Simulation studies are performed to compare the performance of the proposed semi-parametric method with a fully parametric mixture approach. The results show that when the component-baseline hazard is monotonic increasing, the semi-parametric and fully parametric mixture approaches are comparable for mildly and moderately censored samples. When the component-baseline hazard is not monotonic increasing, the semi-parametric method consistently provides less biased estimates than a fully parametric approach and is comparable in efficiency in the estimation of the parameters for all levels of censoring. The methods are illustrated using a real data set of prostate cancer patients treated with different dosages of the drug diethylstilbestrol. Copyright (C) 2003 John Wiley Sons, Ltd.

Fuzzy linguistic model for evaluating the risk of neonatal death

OBJECTIVE: To introduce a fuzzy linguistic model for evaluating the risk of neonatal death. METHODS: The study is based on the fuzziness of the variables newborn birth weight and gestational age at delivery. The inference used was Mamdani's method. Neonatologists were interviewed to estimate the risk of neonatal death under certain conditions and to allow comparing their opinions and the model values. RESULTS: The results were compared with experts' opinions and the Fuzzy model was able to capture the expert knowledge with a strong correlation (r=0.96). CONCLUSIONS: The linguistic model was able to estimate the risk of neonatal death when compared to experts' performance.

Tsunamigenic earthquakes in the Gulf of Cadiz: fault model and recurrence

The Gulf of Cadiz, as part of the Azores-Gibraltar plate boundary, is recognized as a potential source of big earthquakes and tsunamis that may affect the bordering countries, as occurred on 1 November 1755. Preparing for the future, Portugal is establishing a national tsunami warning system in which the threat caused by any large-magnitude earthquake in the area is estimated from a comprehensive database of scenarios. In this paper we summarize the knowledge about the active tectonics in the Gulf of Cadiz and integrate the available seismological information in order to propose the generation model of destructive tsunamis to be applied in tsunami warnings. The fault model derived is then used to estimate the recurrence of large earthquakes using the fault slip rates obtained by Cunha et al. (2012) from thin-sheet neotectonic modelling. Finally we evaluate the consistency of seismicity rates derived from historical and instrumental catalogues with the convergence rates between Eurasia and Nubia given by plate kinematic models.

Memória de longo prazo nos retornos acionistas dos índices de referência da euronext, implicações para a hipótese de mercados eficientes e contributo fractal para aperfeiçoamento do capital asset pricing model

Não existe uma definição única de processo de memória de longo prazo. Esse processo é geralmente definido como uma série que possui um correlograma decaindo lentamente ou um espectro infinito de frequência zero. Também se refere que uma série com tal propriedade é caracterizada pela dependência a longo prazo e por não periódicos ciclos longos, ou que essa característica descreve a estrutura de correlação de uma série de longos desfasamentos ou que é convencionalmente expressa em termos do declínio da lei-potência da função auto-covariância. O interesse crescente da investigação internacional no aprofundamento do tema é justificado pela procura de um melhor entendimento da natureza dinâmica das séries temporais dos preços dos ativos financeiros. Em primeiro lugar, a falta de consistência entre os resultados reclama novos estudos e a utilização de várias metodologias complementares. Em segundo lugar, a confirmação de processos de memória longa tem implicações relevantes ao nível da (1) modelação teórica e econométrica (i.e., dos modelos martingale de preços e das regras técnicas de negociação), (2) dos testes estatísticos aos modelos de equilíbrio e avaliação, (3) das decisões ótimas de consumo / poupança e de portefólio e (4) da medição de eficiência e racionalidade. Em terceiro lugar, ainda permanecem questões científicas empíricas sobre a identificação do modelo geral teórico de mercado mais adequado para modelar a difusão das séries. Em quarto lugar, aos reguladores e gestores de risco importa saber se existem mercados persistentes e, por isso, ineficientes, que, portanto, possam produzir retornos anormais. O objetivo do trabalho de investigação da dissertação é duplo. Por um lado, pretende proporcionar conhecimento adicional para o debate da memória de longo prazo, debruçando-se sobre o comportamento das séries diárias de retornos dos principais índices acionistas da EURONEXT. Por outro lado, pretende contribuir para o aperfeiçoamento do capital asset pricing model CAPM, considerando uma medida de risco alternativa capaz de ultrapassar os constrangimentos da hipótese de mercado eficiente EMH na presença de séries financeiras com processos sem incrementos independentes e identicamente distribuídos (i.i.d.). O estudo empírico indica a possibilidade de utilização alternativa das obrigações do tesouro (OT’s) com maturidade de longo prazo no cálculo dos retornos do mercado, dado que o seu comportamento nos mercados de dívida soberana reflete a confiança dos investidores nas condições financeiras dos Estados e mede a forma como avaliam as respetiva economias com base no desempenho da generalidade dos seus ativos. Embora o modelo de difusão de preços definido pelo movimento Browniano geométrico gBm alegue proporcionar um bom ajustamento das séries temporais financeiras, os seus pressupostos de normalidade, estacionariedade e independência das inovações residuais são adulterados pelos dados empíricos analisados. Por isso, na procura de evidências sobre a propriedade de memória longa nos mercados recorre-se à rescaled-range analysis R/S e à detrended fluctuation analysis DFA, sob abordagem do movimento Browniano fracionário fBm, para estimar o expoente Hurst H em relação às séries de dados completas e para calcular o expoente Hurst “local” H t em janelas móveis. Complementarmente, são realizados testes estatísticos de hipóteses através do rescaled-range tests R/S , do modified rescaled-range test M - R/S e do fractional differencing test GPH. Em termos de uma conclusão única a partir de todos os métodos sobre a natureza da dependência para o mercado acionista em geral, os resultados empíricos são inconclusivos. Isso quer dizer que o grau de memória de longo prazo e, assim, qualquer classificação, depende de cada mercado particular. No entanto, os resultados gerais maioritariamente positivos suportam a presença de memória longa, sob a forma de persistência, nos retornos acionistas da Bélgica, Holanda e Portugal. Isto sugere que estes mercados estão mais sujeitos a maior previsibilidade (“efeito José”), mas também a tendências que podem ser inesperadamente interrompidas por descontinuidades (“efeito Noé”), e, por isso, tendem a ser mais arriscados para negociar. Apesar da evidência de dinâmica fractal ter suporte estatístico fraco, em sintonia com a maior parte dos estudos internacionais, refuta a hipótese de passeio aleatório com incrementos i.i.d., que é a base da EMH na sua forma fraca. Atendendo a isso, propõem-se contributos para aperfeiçoamento do CAPM, através da proposta de uma nova fractal capital market line FCML e de uma nova fractal security market line FSML. A nova proposta sugere que o elemento de risco (para o mercado e para um ativo) seja dado pelo expoente H de Hurst para desfasamentos de longo prazo dos retornos acionistas. O expoente H mede o grau de memória de longo prazo nos índices acionistas, quer quando as séries de retornos seguem um processo i.i.d. não correlacionado, descrito pelo gBm(em que H = 0,5 , confirmando- se a EMH e adequando-se o CAPM), quer quando seguem um processo com dependência estatística, descrito pelo fBm(em que H é diferente de 0,5, rejeitando-se a EMH e desadequando-se o CAPM). A vantagem da FCML e da FSML é que a medida de memória de longo prazo, definida por H, é a referência adequada para traduzir o risco em modelos que possam ser aplicados a séries de dados que sigam processos i.i.d. e processos com dependência não linear. Então, estas formulações contemplam a EMH como um caso particular possível.

Brownian motion with drift threshold model

In this thesis we implement estimating procedures in order to estimate threshold parameters for the continuous time threshold models driven by stochastic di®erential equations. The ¯rst procedure is based on the EM (expectation-maximization) algorithm applied to the threshold model built from the Brownian motion with drift process. The second procedure mimics one of the fundamental ideas in the estimation of the thresholds in time series context, that is, conditional least squares estimation. We implement this procedure not only for the threshold model built from the Brownian motion with drift process but also for more generic models as the ones built from the geometric Brownian motion or the Ornstein-Uhlenbeck process. Both procedures are implemented for simu- lated data and the least squares estimation procedure is also implemented for real data of daily prices from a set of international funds. The ¯rst fund is the PF-European Sus- tainable Equities-R fund from the Pictet Funds company and the second is the Parvest Europe Dynamic Growth fund from the BNP Paribas company. The data for both funds are daily prices from the year 2004. The last fund to be considered is the Converging Europe Bond fund from the Schroder company and the data are daily prices from the year 2005.

analytical model and measurements of the target erosion depth profile of balanced and unbalanced planar magnetron cathodes

The erosion depth profile of planar targets in balanced and unbalanced magnetron cathodes with cylindrical symmetry is measured along the target radius. The magnetic fields have rotational symmetry. The horizontal and vertical components of the magnetic field B are measured at points above the cathode target with z = 2 x 10(-3) m. The experimental data reveal that the target erosion depth profile is a function of the angle. made by B with a horizontal line defined by z = 2 x 10(-3) m. To explain this dependence a simplified model of the discharge is developed. In the scope of the model, the pathway lengths of the secondary electrons in the pre-sheath region are calculated by analytical integration of the Lorentz differential equations. Weighting these lengths by using the distribution law of the mean free path of the secondary electrons, we estimate the densities of the ionizing events over the cathode and the relative flux of the sputtered atoms. The expression so deduced correlates for the first time the erosion depth profile of the target with the angle theta. The model shows reasonably good fittings to the experimental target erosion depth profiles confirming that ionization occurs mainly in the pre-sheath zone.

Credit default swap (CDS) prediction model & trading strategy

This project focuses on the study of different explanatory models for the behavior of CDS security, such as Fixed-Effect Model, GLS Random-Effect Model, Pooled OLS and Quantile Regression Model. After determining the best fitness model, trading strategies with long and short positions in CDS have been developed. Due to some specifications of CDS, I conclude that the quantile regression is the most efficient model to estimate the data. The P&L and Sharpe Ratio of the strategy are analyzed using a backtesting analogy, where I conclude that, mainly for non-financial companies, the model allows traders to take advantage of and profit from arbitrages.

Low coronary perfusion pressure is associated with endocardial fibrosis in a rat model of volume overload cardiac hypertrophy

Left ventricular hypertrophy following volume overload is regarded as an example of cardiac remodeling without increased fibrosis accumulation. However, infarction is associated with increased fibrosis within the noninfarcted, hypertrophied myocardium, particularly in the subendocardial regions. It is conceivable to suppose that, as also occurs postinfarction, low coronary driving pressure may also interfere with accumulation of myocardial fibrosis following aortocaval fistula. PURPOSE: To investigate the role of acute hemodynamic changes in subsequent deposition of cardiac fibrosis in response to aortocaval fistula. METHOD: Aortocaval fistula were created in 4 groups of Wistar rats that were followed over 4 and 8 weeks: aortocaval fistula 4 and aortocaval fistula 8 (10 rats each) and their respective controls (sham-operated controls - Sh), Sh4 and Sh8 (8 rats each). Hemodynamic measurements were performed 1 week after surgery. Hypertrophy and fibrosis were quantified by myocyte diameter and collagen volume fraction at the end of follow up. RESULT: Compared with Sh4 and Sh8, pulse pressure, left ventricular end-diastolic pressure, and +dP/dt were higher in aortocaval fistula 4 and aortocaval fistula 8, but -dP/dt was similar. Coronary driving pressure (mm Hg), used as an estimate of perfusion pressure, was lower in aortocaval fistula 8 (52.6 ± 4.1) than in Sh8 (100.8 ± 1.3), but comparable between aortocaval fistula 4 (50.0 ± 8.9) and Sh4 (84.8 ± 2.3). Myocyte diameter was greater in aortocaval fistula 8, whereas interstitial and subendocardial fibrosis were greater in aortocaval fistula 4 and aortocaval fistula 8. Coronary driving pressure correlated inversely and independently with subendocardial fibrosis (r² = .86, P <.001), whereas left ventricular systolic pressure (r² = 0.73, P = .004) and end-diastolic pressure (r² = 0.55, P = 012) correlated positively and independently with interstitial fibrosis. CONCLUSION: Coronary driving pressure falls and ventricular pressures increase early after aortocaval fistula and are associated with subsequent myocardial fibrosis deposition.

A bio-inspired computational model for motion detection

Tese de Doutoramento (Programa Doutoral em Engenharia Biomédica)

A nutrient balance model for typical Irish farming systems

The twin objectives of the work described were to construct nutrient balance models (NBM) for a range of Irish animal production systems and to evaluate their potential as a means of estimating the nutrient composition of farm wastes. The NBM has three components. The first is the intake of nutrients in the animal's diet. The second is retention or the nutrients the animal retains for the production of milk, meat or eggs. The third is the balance or the difference between the nutrient intake and retention. Data on the intake levels and their nutrient value for dairy cows, beef cattle, pigs and poultry systems were assembled. Literature searches and interviews with National experts were the primary sources of information. NBMs were then constructed for each production system. Summary tables of the nutrient values for the common diet constituents used in Irish animal production systems, the nutrient composition of the animal products and the NBMs (nutrient intake, retention and excretion) for a range of production systems were assembled. These represent the first comprehensive data set of this type for Irish animal production systems. There was generally good agreement between the derived NBMs values and those published in the literature. The NBMs were validated on a number of farms. Data on animal numbers, fertiliser use, concentrates inputs and production output were recorded on seven farms. Using the data a nutrient input/output balance was constructed for each farm. This was compared with the NBM estimate of the farm nutrient balance. The results showed good agreement between the measured balance and the NBM estimate particularly for the pig and poultry farms. However, the validation emphasised the inherent risks associated with NBMs. The average values used for feed intake and production parameters in the NEMs may result in the under or over estimate of actual nutrient balances on individual farms where these variables are substantially different. On the grassland farms there was a poor correlation between the input/output estimate and the NBM. This possibly results from the omission of the soil's contribution to the nutrient balance. However, the results indicate that the NBMs developed are a potentially useful tool for estimating nutrient balances. They also will serve to highlight the significant fraction of the nutrient inputs into farming systems that are retained on the farm. The potential of the NBM as a means of estimating the nutrient composition of farm wastes was evaluated on two farms. Feed intake and composition, animal production, slurry production was monitored during the indoor winter feeding period. Slurry samples were taken for analysis. The appropriates NBMs were used to estimate the nutrient balance for each farm. The nutrient content of the slurry produced was calculated. There was a good agreement between the NBM estimate and the measured values. This preliminary evaluation suggests that the NBM has a potential to provide the farmer with a simple means of estimating the nutrient value of his slurry.

Can a Lucas model with habit generate realistic conditional volatility in exchange rate returns?

In this paper, we attempt to give a theoretical underpinning to the well established empirical stylized fact that asset returns in general and the spot FOREX returns in particular display predictable volatility characteristics. Adopting Moore and Roche s habit persistence version of Lucas model we nd that both the innovation in the spot FOREX return and the FOREX return itself follow "ARCH" style processes. Using the impulse response functions (IRFs) we show that the baseline simulated FOREX series has "ARCH" properties in the quarterly frequency that match well the "ARCH" properties of the empirical monthly estimations in that when we scale the x-axis to synchronize the monthly and quarterly responses we find similar impulse responses to one unit shock in variance. The IRFs for the ARCH processes we estimate "look the same" with an approximately monotonic decreasing fashion. The Lucas two-country monetary model with habit can generate realistic conditional volatility in spot FOREX return.

Technology shocks and aggregate fluctuations in an estimated hybrid RBC model

This paper contributes to the on-going empirical debate regarding the role of the RBC model and in particular of technology shocks in explaining aggregate fluctuations. To this end we estimate the model’s posterior density using Markov-Chain Monte-Carlo (MCMC) methods. Within this framework we extend Ireland’s (2001, 2004) hybrid estimation approach to allow for a vector autoregressive moving average (VARMA) process to describe the movements and co-movements of the model’s errors not explained by the basic RBC model. The results of marginal likelihood ratio tests reveal that the more general model of the errors significantly improves the model’s fit relative to the VAR and AR alternatives. Moreover, despite setting the RBC model a more difficult task under the VARMA specification, our analysis, based on forecast error and spectral decompositions, suggests that the RBC model is still capable of explaining a significant fraction of the observed variation in macroeconomic aggregates in the post-war U.S. economy.

Realised and Optimal Monetary Policy Rules in an Estimated Markov-Switching DSGE Model of the United Kingdom

This paper investigates underlying changes in the UK economy over the past thirtyfive years using a small open economy DSGE model. Using Bayesian analysis, we find UK monetary policy, nominal price rigidity and exogenous shocks, are all subject to regime shifting. A model incorporating these changes is used to estimate the realised monetary policy and derive the optimal monetary policy for the UK. This allows us to assess the effectiveness of the realised policy in terms of stabilising economic fluctuations, and, in turn, provide an indication of whether there is room for monetary authorities to further improve their policies.

A bayesian spatial individual effects probit model of the 2010 UK general election

The Conservative Party emerged from the 2010 United Kingdom General Election as the largest single party, but their support was not geographically uniform. In this paper, we estimate a hierarchical Bayesian spatial probit model that tests for the presence of regional voting effects. This model allows for the estimation of individual region-specic effects on the probability of Conservative Party success, incorporating information on the spatial relationships between the regions of the mainland United Kingdom. After controlling for a range of important covariates, we find that these spatial relationships are significant and that our individual region-specic effects estimates provide additional evidence of North-South variations in Conservative Party support.

Productivity shocks and aggregate fluctuations in an estimated endogenous growth model with human capital

Employing an endogenous growth model with human capital, this paper explores how productivity shocks in the goods and human capital producing sectors contribute to explaining aggregate fluctuations in output, consumption, investment and hours. Given the importance of accounting for both the dynamics and the trends in the data not captured by the theoretical growth model, we introduce a vector error correction model (VECM) of the measurement errors and estimate the model’s posterior density function using Bayesian methods. To contextualize our findings with those in the literature, we also assess whether the endogenous growth model or the standard real business cycle model better explains the observed variation in these aggregates. In addressing these issues we contribute to both the methods of analysis and the ongoing debate regarding the effects of innovations to productivity on macroeconomic activity.