摘要 本文主要考慮群集樣本(clustered samples)的存活分析，而每一群集中又分為兩種組別(groups)。假定同群集同組別內的個體共享相同但不可觀測的隨機脆弱性(frailty)，因此面臨的是雙變量脆弱性變數的多變量存活資料。首先，驗證雙變量脆弱性對雙變量對數存活時間及雙變量存活時間之相關係數所造成的影響。接著，假定雙變量脆弱性服從雙變量對數常態分配，條件存活時間模式為韋伯迴歸模式，我們利用EM法則，推導出雙變量脆弱性之多變量存活模式中母數的估計方法。 關鍵詞：雙變量脆弱性，Weibull迴歸模式，對數常態分配，EM法則 / Abstract Consider survival analysis for clustered samples, where each cluster contains two groups. Assume that individuals within the same cluster and the same group share a common but unobservable random frailty. Hence, the focus of this work is on bivariate frailty model in analysis of multivariate survival data. First, we derive expressions for the correlation between the two survival times to show how the bivariate frailty affects these correlation coefficients. Then, the bivariate log-normal distribution is used to model the bivariate frailty. We modified EM algorithm to estimate the parameters for the Weibull regression model with bivariate log-normal frailty. Key words：bivariate frailty, Weibull regression model, log-normal distribution, EM algorithm.

This paper deals with the computation of OC- and ASN-function of sequential probability ratio tests in the multi-parameter case. In generalization of the method of conjugated parameter pairs Wald-like approximations are presented for the OC- and ASN-function. These characteristics can be used describing robustness properties of a sequential test in case of nuisance parameters. As examples tests are considered for the mean and the variance of a normal distribution.

This thesis is focused on mathematical modelling of exchange rate CZK/USD in 1991 - 2014. Time series was divided into 5 parts. First Box-Jenkins methodology models were examined, especially ARIMA model. Unfortunately, the model could not be used because none of the time series showed correlation. The time series is considered as a white noise. The data appear to be completely random and unpredictable. The time series have not constant variance neither normal distribution and therefore GARCH volatility model was used as the second model. It is better not to divide time series when using model of volatility. Volatility model contributes to more accurate prediction than the standard deviation. Results were calculated in RStudio software and MS Excel.

Made available in DSpace on 2016-06-02T20:04:53Z (GMT). No. of bitstreams: 1 6424.pdf: 5426262 bytes, checksum: 82bb9386f85845b0d3db787265ea8236 (MD5) Previous issue date: 2014-11-21 / The Generalized Normal (GN) and Generalized lognormal (logGN) distributions are flexible for accommodating features present in the data that are not captured by traditional distribution, such as the normal and the lognormal ones, respectively. These distributions are considered to be tools for the reduction of outliers and for the obtention of robust estimates. However, computational problems have always been the major obstacle to obtain the effective use of these distributions. This paper proposes the Bayesian reference analysis methodology to estimate the GN and logGN. The reference prior for a possible order of the model parameters is obtained. It is shown that the reference prior leads to a proper posterior distribution for all the proposed model. The development of Monte Carlo Markov Chain (MCMC) is considered for inference purposes. To detect possible influential observations in the models considered, the Bayesian method of influence analysis on a case based on the Kullback-Leibler divergence is used. In addition, a scale mixture of uniform representation of the GN and logGN distributions are exploited, as an alternative method in order, to allow the development of efficient Gibbs sampling algorithms. Simulation studies were performed to analyze the frequentist properties of the estimation procedures. Real data applications demonstrate the use of the proposed models. / As distribuições normal generalizada (NG) e lognormal generalizada (logNG) são flexíveis por acomodarem características presentes nos dados que não são capturadas por distribuições tradicionais, como a normal e a lognormal, respectivamente. Essas distribuições são consideradas ferramentas para reduzir as observações aberrantes e obter estimativas robustas. Entretanto o maior obstáculo para a utilização eficiente dessas distribuições tem sido os problemas computacionais. Este trabalho propõe a metodologia da análise de referência Bayesiana para estimar os parâmetros dos modelos NG e logNG. A função a priori de referência para uma possível ordem dos parâmetros do modelo é obtida. Mostra-se que a função a priori de referência conduz a uma distribuição a posteriori própria, em todos os modelos propostos. Para fins de inferência, é considerado o desenvolvimento de métodos Monte Carlo em Cadeias de Markov (MCMC). Para detectar possíveis observações influentes nos modelos considerados, é utilizado o método Bayesiano de análise de influência caso a caso, baseado na divergência de Kullback-Leibler. Além disso, uma representação de mistura de escala uniforme para as distribuições NG e logNG é utilizada, como um método alternativo, para permitir o desenvolvimento de algoritmos de amostrador de Gibbs. Estudos de simulação foram desenvolvidos para analisar as propriedades frequentistas dos processos de estimação. Aplicações a conjuntos de dados reais mostraram a aplicabilidade dos modelos propostos.

Made available in DSpace on 2016-06-02T20:06:02Z (GMT). No. of bitstreams: 1 2390.pdf: 3256865 bytes, checksum: 7ad1acbefc5f29dddbaad3f14dbcef7c (MD5) Previous issue date: 2009-04-17 / Financiadora de Estudos e Projetos / Microarrays technologies are used to measure the expression levels of a large amount of genes or fragments of genes simultaneously in diferent situations. This technology is useful to determine genes that are responsible for genetic diseases. A common statistical methodology used to determine whether a gene g has evidences to diferent expression levels is the t-test which requires the assumption of normality for the data (Saraiva, 2006; Baldi & Long, 2001). However this assumption sometimes does not agree with the nature of the analyzed data. In this work we use the skew-normal distribution described formally by Azzalini (1985), which has the normal distribution as a particular case, in order to relax the assumption of normality. Considering a frequentist approach we made a simulation study to detect diferences between the gene expression levels in situations of control and treatment through the t-test. Another simulation was made to examine the power of the t-test when we assume an asymmetrical model for the data. Also we used the likelihood ratio test to verify the adequability of an asymmetrical model for the data. / Os microarrays são ferramentas utilizadas para medir os níveis de expressão de uma grande quantidade de genes ou fragmentos de genes simultaneamente em situações variadas. Com esta ferramenta é possível determinar possíveis genes causadores de doenças de origem genética. Uma abordagem estatística comumente utilizada para determinar se um gene g apresenta evidências para níveis de expressão diferentes consiste no teste t, que exige a suposição de normalidade aos dados (Saraiva, 2006; Baldi & Long, 2001). No entanto, esta suposição pode não condizer com a natureza dos dados analisados. Neste trabalho, será utilizada a distribuição normal assimétrica descrita formalmente por Azzalini (1985), que tem a distribuição normal como caso particular, com o intuito de flexibilizar a suposição de normalidade. Considerando a abordagem clássica, é realizado um estudo de simulação para detectar diferenças entre os níveis de expressão gênica em situações de controle e tratamento através do teste t, também é considerado um estudo de simulação para analisar o poder do teste t quando é assumido um modelo assimétrico para o conjunto de dados. Também é realizado o teste da razão de verossimilhança, para verificar se o ajuste de um modelo assimétrico aos dados é adequado.

Uppsatsen undersöker vilka mål Växjö Transportkyla har vid dimensionering av deras säkerhetslager. Genom att studera litteratur och genomföra intervjuer på företaget sammanställs målen till följande: Hög servicenivå Låg kapitalbindning Kontroll över storleken på säkerhetslager Differentierade servicenivåer Användarvänlig beräkning av säkerhetslager baserad på teori Den andra frågeställningen för uppsatsen syftar till att bestämma vilken beräkningsmetod som Växjö Transportkyla bör välja för att uppnå de mål som satts upp. Fem olika beräkningsmetoder testas och det undersöks även vilken statistisk fördelning produkternas efterfråga har. Efter att ha testat de fem beräkningsmetoderna på data insamlad från Växjö Transportkyla blir slutsatsen att för de 426 poissonfördelade produkterna bör SERV2 poissonfördelning användas. För de övriga ca 700 artiklarna finns inget resultat då efterfrågefördelningen endast gick att bestämma för ett fåtal av dem. En rekommendation till företaget blir att testa fler fördelningar och på så sätt komma fram till bästa resultat på övriga produkter. Det poängteras även att det är bra att mäta servicenivåer och differentiera servicenivåmålet för olika produkter. Uppsatsen ger sitt teoretiska bidrag i form av en generalisering av vilka produktegenskaper som påverkar valet av beräkningsmetod vid säkerhetslageroptimering och hur. De fyra övriga beräkningsmetoderna jämförs här med SERV2 poissonfördelning för att se om det vid vissa produktegenskaper passar med en annan beräkningsmetod. Detta leder fram till slutsatsen att SERV2 normalfördelning var den av de fyra beräkningsmetoderna som hade mest korrekt resultat och det var vid medellånga ledtider och låg efterfrågan. / The thesis examines which goals Växjö Transportkyla has with calculating their safety stock. By analysing literature and responses from interviews at the company, the following goals have been put together: High service level Low tied-up costs Control over the size of the safety stock Differentiated service levels User-friendly calculation of safety stock based on theory The purpose of the second research question is to decide which safety stock calculation should be used by Växjö Transportkyla to reach the goals. Five different safety stock calculations are tested and the statistical distribution of the products’ demand is examined. After testing the safety stock calculations on data from Växjö Transportkyla, the conclusion is that for the 426 poisson distributed products, SERV2 poisson distribution should be used. For the rest of the products, which are around 700, there is no result since the correct distribution only was found for a few of them. A recommendation for the company is to test more distributions to be able to calculate the safety stock for the rest of their products. It is also pointed out that the company should measure their service levels and differentiate the service level for the different products. The thesis provides a theoretical contribution by examining which product features affect the choice of safety stock calculation and how. The four remaining safety stock calculations are compared to SERV2 normal distribution to examine if it is plausible for certain product features to use another safety stock calculation. This leads to the conclusion that SERV2 normal distribution was the one out of the four safety stock calculations that had the most correct result in relation to SERV2 poisson distribution and that was for the product features of middle-long lead times and low demand.

Este trabalho tem como objetivo principal apresentar uma abordagem clássica e Bayesiana dos modelos lineares com observações censuradas, que é uma nova área de pesquisa com grandes possibilidades de aplicações. Aqui, substituimos o uso convencional da distribuição normal para os erros por uma família de distribuições mais flexíveis, o que nos permite lidar de forma mais adequada com observações censuradas na presença de outliers. Esta família é obtida através de um mecanismo de fácil construção e possui como casos especiais as distribuições t de Student, Pearson tipo VII, slash, normal contaminada e, obviamente, a normal. Para o caso de respostas correlacionadas e censuradas propomos um modelo de regressão linear robusto baseado na distribuição t de Student, desenvolvendo um algoritmo tipo EM que depende dos dois primeiros momentos da distribuição t de Student truncada. / This work aims to present a classical and Bayesian approach to linear models with censored observations, which is a new area of research with great potential for applications. Here, we replace the conventional use of the normal distribution for the errors of a more flexible family of distributions, which deal in more appropriately with censored observations in the presence of outliers. This family is obtained through a mechanism easy to construct and has as special cases the distributions Student t, Pearson type VII, slash, contaminated normal, and obviously normal. For the case of correlated and censored responses we propose a model of robust linear regression based on Student\'s t distribution and we developed an EM type algorithm based on the first two moments of the truncated Student\'s t distribution.

nÃo hà / Este trabalho analisa durante o perÃodo de 01/2008 a 12/2011 o risco de mercado de seis Ãndices setoriais da Bolsa de Valores de SÃo Paulo (BM&FBovespa): o Ãndice imobiliÃrio (IMOB), o Ãndice de energia elÃtrica (IEE), o Ãndice de consumo (ICON), o Ãndice do setor industrial (INDX), o Ãndice financeiro (IFNC) e o Ãndice setorial de telecomunicaÃÃes (ITEL). AtravÃs da mÃtrica Value-at-Risk (VaR) estimam-se quatro modelos. Dois desses modelos sÃo ditos incondicionais no que se refere à variÃncia: o VaR Gaussiano Incondicional, admitindo que os retornos seguem uma distribuiÃÃo normal, e o VaR Best Fitting Incondicional, construÃdo a partir da distribuiÃÃo de probabilidades que melhor se ajusta Ãs sÃries de retornos. Os outros dois modelos sÃo chamados de condicionais, assumindo que a volatilidade varia ao longo do tempo. Os modelos autoregressivos do tipo GARCH sÃo utilizados para estimar a variÃncia condicional de cada Ãndice, possibilitando a estimaÃÃo do VaR Gaussiano Incondicional e do VaR Best Fitting Incondicional. Em seguida, realizam-se backtestings dos modelos de VaR, revelando a superioridade dos modelos condicionais. Por fim, atravÃs de grÃficos de Balzer, observou-se a performance dos Ãndices por meio de confrontos entre eles. Foi constatado que, para o perÃodo analisado, o IEE vence todos os embates feitos com os demais Ãndices, apresentando a melhor relaÃÃo risco x retorno. O setor imobiliÃrio, representado pelo IMOB, perde todos os confrontos. / This work analyzes during the period between 2008/01 and 2011/12 the market risk of six sectorial indexes from the SÃo PauloÂs Stock Market (BM&FBovespa): the real state index (IMOB), the eletric power index (IEE), the consumption index (ICON), the industrial sector index (INDX), the financial index (IFNC) and the telecommunications sector index (ITEL). Throughout the Value-at-Risk metric (VaR), four models are estimated. Two of those models are called unconditional, due to its variance: the Unconditional Gaussian VaR, that admits that the returns follow a normal distribution, and the Unconditional Best Fitting VaR, built from the distribution of probabilities that better fits to the returns series. The other two models are called conditionals, assuming that the volatility changes along the time. The GARCH autoregressive models are used to estimate the conditional variance of each index, allowing an estimation of the Unconditional Gaussian VaR and the Unconditional Best Fitting VaR. Afterwards, the VaR models backtestings are realized, revealing the conditional models superiority. Finally, throughout the BalzerÂs graphics, the indexes performances were observed over the confrontations between them. It was found that, for the analyzed period, the IEE wins every confrontation against the all other indexes, showing the best relation risk x return. The real state index sector, represented by the IMOB, lost all the confronts.

Madeira laminada colada (MLC) é um material de construção muito empregado em estruturas. Este produto é composto de lâminas classificadas de madeira, coladas horizontalmente, para formar peças estruturais de madeira grandes dimensões. A união das lâminas é realizada através de distintos tipos de emendas longitudinais, sendo as emendas dentadas as mais utilizadas. Considerando que os componentes estruturais de MLC, em geral, são solicitadas a carregamentos cíclicos, este trabalho tem por finalidade a proposta de um modelo estatístico para a determinação da vida à fadiga em menos dentadas de madeira. O trabalho foi realizado no Laboratório de Madeiras e de Estruturas de madeira (LaMEM), com o estudo teórico e experimental da fadiga em corpos-de-prova tracionados de ligações com emendas dentadas, baseado em um planejamento estatístico de experimentos e na NBR 7190/97. Os resultados contidos mostram que o modelo Polinomial Ortogonal Múltiplo da distribuição de Birnbaum-Saunders é de grande precisão e o mais adequado ao estudo da fadiga. O uso deste modelo pode ser de grande benefício, já que com ele se poderá estimar e caracterizar com maior confiabilidade e precisão a um menor custo a vida à fadiga em emendas dentadas de madeira, considerando as variáveis independentes de tensão e freqüência. / In the production of structural elements of Glued Laminated Timber (GLULAM), the horizontal union of lumbers are made with finger joints. Considering that the structural components of GLULAM need a great number of finger joints, and some of these structures are subject to cyclic loading, the objective of this work is to present a statistical model to estimate the fatigue life in lumber tension finger joints. The theoretical and experimental work was made in the Laboratory of Wood and Timber Structures (LaMEM), based on a statistical experiment design using tension tests and in agreement with the NBR 7190/97 code. The estimation procedure was based on the multiple orthogonal polynomial Birnbaum-Saunders model and the results show that the parameter estimates of the multiple orthogonal polynomial Birnbaum-Saunders model are obtained with at good accuracy and the most appropriate for the wood fatigue study. The use of the model can be of great benefit, since with it can be fit and characterize with greater reliability and precision at a smaller cost the fatigue life in finger joints of wood, considering the independent variables stress and frequency.

Submitted by Divisão de Documentação/BC Biblioteca Central (ddbc@ufam.edu.br) on 2018-05-29T14:38:33Z No. of bitstreams: 2 license_rdf: 0 bytes, checksum: d41d8cd98f00b204e9800998ecf8427e (MD5) VersaoFinal.pdf: 2882901 bytes, checksum: a35c6d27fe0f9aa61dfe3a96244b3140 (MD5) / Approved for entry into archive by Divisão de Documentação/BC Biblioteca Central (ddbc@ufam.edu.br) on 2018-05-29T14:38:46Z (GMT) No. of bitstreams: 2 license_rdf: 0 bytes, checksum: d41d8cd98f00b204e9800998ecf8427e (MD5) VersaoFinal.pdf: 2882901 bytes, checksum: a35c6d27fe0f9aa61dfe3a96244b3140 (MD5) / Made available in DSpace on 2018-05-29T14:38:46Z (GMT). No. of bitstreams: 2 license_rdf: 0 bytes, checksum: d41d8cd98f00b204e9800998ecf8427e (MD5) VersaoFinal.pdf: 2882901 bytes, checksum: a35c6d27fe0f9aa61dfe3a96244b3140 (MD5) Previous issue date: 2018-03-12 / CAPES - Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / The traditional estimation of mixture regression models is based on the assumption of normality of component errors and thus is sensitive to outliers, heavy-tailed and/or asymmetric errors. Another drawback is that, in general, the analysis is restricted to directly observed predictors. We present a proposal to deal with these issues simultaneously in the context of mixture regression by extending the classic normal model by assuming that, for each mixture component, the random errors and the covariates jointly follow a scale mixture of skew-normal distributions. It is also assumed that the covariates are observed with error. An MCMC-type algorithm to perform Bayesian inference is developed and, in order to show the efficacy of the proposed methods, simulated and real data sets are analyzed. / A estimação tradicional em mistura de modelos de regressão é baseada na suposição de normalidade para os erros aleatórios, sendo assim, sensível a outliers, caudas pesadas e erros assimétricos. Outra desvantagem é que, em geral, a análise é restrita a preditores que são observados diretamente. Apresentamos uma proposta para lidar com estas questões simultaneamente no contexto de mistura de regressões estendendo o modelo normal clássico. Assumimos que, conjuntamente e em cada componente da mistura, os erros aleatórios e as covariáveis seguem uma mistura de escala da distribuição normal assimétrica. Além disso, é feita a suposição de que as covariáveis são observadas com erro aditivo. Um algorítmo do tipo MCMC foi desenvolvido para realizar inferência Bayesiana. A eficácia do modelo proposto é verificada via análises de dados simulados e reais.