Avaliamos a incidência e o clareamento para o vírus do papiloma humano (HPV) dos tipos oncogênicos e não-oncogênicos em uma coorte aberta com 202 mulheres portadoras do vírus da imunodeficiência humana (HIV), e identificamos alguns fatores de risco e proteção associados a cada desfecho utilizando modelos de fragilidade Gama. No modelo de incidência, foram estudados os tempos até incidência de HPV oncogênicos e não-oncogênicos para cada mulher; no modelo de clareamento, foram estudados os correspondentes tempos até clareamento. Comparamos os erros-padrões estimados pela matriz de informação observada com os erros-padrões bootstrap para ambos os modelos e verificamos que a proposta de Verweij & Houwelingen (1994) para a matriz de variâncias e covariâncias dos parâmetros é a mais apropriada. Para a incidência de HPV oncogênicos, identificamos como fator de risco o uso de drogas em que a taxa de incidência para as pacientes que usam drogas é 1.88 (IC 90%, 1.01; 3.5) vezes aquela correspondente a mulheres que não usam e como fator de proteção a renda em que a taxa de incidência de pacientes com renda igual ou superior a 3 salários mínimos é 0.62 (IC 90%, 0.38; 1.00) vezes a taxa referente àquelas com renda menor que 3 salários mínimos. Para a incidência de HPV não-oncogênicos identificamos como fatores de risco a escolaridade e o total de gestações, em que, para a última, a taxa de incidência para as mulheres que tiveram mais do que uma gestação é 1.76 (IC 90%, 1.09; 2.86) vezes a taxa referente àquelas que tiveram uma ou nenhuma. Para o clareamento de HPV oncogênicos identificamos como fatores que indicam um clareamento mais rápido a renda, a idade e o tratamento antirretroviral (ARV), em que, para a última, supondo mulheres com fragilidades iguais, a taxa de clareamento para as pacientes que eram tratadas com o esquema inibidor de protease (IP) é 1.79 (IC 90%, 1.1; 2.9) vezes aquela correspondente a mulheres que não foram tratadas com nenhum tratamento ARV e como fator que indicam um clareamento mais lento o número de parceiros sexuais no último ano, em que, as pacientes com mais de um parceiro tiveram taxa de clareamento 0.39 (IC 90%, 0.16; 0.98) vezes a taxa de clareamento referente à uma mulher que teve um parceiro ou menos. Para o clareamento de HPV não-oncogênicos tivemos como fator que indica um clareamento mais lento o hábito tabagista em que, supondo fragilidades iguais, pacientes fumantes tem a taxa de clareamento 0.53 (IC 90%, 0.32; 0.87) vezes a taxa referente à uma mulher que não fuma. / We evaluated the incidence and clearance for oncogenic and non-oncogenic human papilloma virus (HPV) in an open cohort of 202 women infected with human immunodeficiency virus (HIV), and we identified some risk factors and protective factors for each outcome using Gamma frailty models. In the incidence model, we studied the incidence of stroke by oncogenic and non-oncogenic HPV for each woman; in the clearance model, the corresponding times to clearance were studied. We compared the standard errors estimated by the observed information matrix with bootstrap standard errors for both models and found that the variance and covariance matrix of the parameters proposed by Verweij & Houwelingen (1994) is more appropriate. For the incidence of oncogenic HPV, identified as a risk factor drug use and the incidence rate for patients who use drugs is 1.88 (90% CI, 1.01; 3.5) times the rate for those who do not use and as a protective factor income where the incidence rate is 0.62 (90% CI, 0.38; 1.00) times the rate for those earning less than 3 minimum wages. For the incidence of non-oncogenic HPV identified as risk factors schooling and total pregnancies, in which, for the latter, the incidence rate for women who had more than one pregnancy is 1.76 (90% CI, 1.09; 2.86) times the rate for those which have one or none. For clearance of oncogenic HPV identified as factors that indicate a faster clearance income, age and antiretroviral therapy (ART), in which, to the last, with women assuming equal frailties, the rate of clearance for patients who were treated with the protease inhibitor (IP) regimen is 1.79 (90% CI, 1.1; 2.9) times the rate for those who were not treated with any antiretroviral regimen and as a factor that indicates slower clearance the number of sexual partners in the last year, and for patients with more than one partner the clearance rate 0.39 (IC 90%, 0.16; 0.98) times the rate referring to a woman who had up to a partner. For the clearance of non-oncogenic HPV had a factor which indicates a slower clearance smoking habit, assuming equal frailties, smokers have the clearance rate 0.53 (90% CI, 0.32; 0.87) times the rate referring to a woman who does not smoke.

Parametrizações de modelos e estruturas de correlações dos parâmetros no âmbito agropecuário são importantes por caracterizarem o comportamento de um sistema em resposta a variações de múltiplos cenários (clima, genótipos, dietas nutricionais, dentre outros fatores) que existem em escalas globais. O objetivo foi contribuir com inferências estatísticas na produção de gases CO2 [um potente Gás de Efeito Estufa (GEE)] nas fermentações in vitro de feno de alfafa, comparando métodos frequentistas com novas metodologias surgidas na literatura científica como a combinação dos métodos de Rejeição por Atraso e o Metropólis Adaptativo (RAMA), até então não testados para predições de gases de fermentação in vitro. Além disso, modelos de séries temporais foram usados para previsão da produção de CO2 nas fermentações de gases in vitro de feno de alfafa. Dentro do contexto de crescimento de gado de corte, foi realizada pela primeira vez uma abordagem para predições individuais dos animais para taxa de ganho de peso e a necessidade de energia para mantença baseada na dinâmica de crescimento e composição química corporal do Modelo de Crescimento de Davis (MCD), com comparação de análise de covariância multivariada entre diferentes cenários (gêneros, sistemas e genótipo cruzados), em um experimento a campo no Brasil. Adicionalmente calibrações dos parâmetros baseadas na amostra de cada cenário, pelos ajustes do MCD e usando análise frequentista, bootstrap não-paramétrico e simulações Monte Carlo foram realizadas com os dados nacionais (raça cruzada) e comparada com as estimativas originais do modelo obtido com raças Britânicas (Bos taurus). Os principais critérios adotados para avaliar os ajustes dos modelos foram o Erro Quadrático Médio de Predição (EQMP), o Critério de Informação Akaike (AIC) e o Critério de Informação Bayesiano (BIC). Os resultados não só contribuirão para o avanço da literatura existente, mas também auxiliarão a indústria de carne bovina e produtores rurais a encontrar especificações do mercado de carne, tanto a nível nacional e internacional. Concluiu-se que i) na produção de gases: o modelo ARIMA (1, 1, 2) ajustou a produção acumulativa de CO2, atingindo o valor máximo de 1,1066 (mL) no tempo de 47,5 h e a equação é indicada para estimar a produção de gases; ii) no crescimento de gados de corte usando as estimativas individuais do MCD, os vetores de efeitos de energia de mantença e o acréscimo de proteína possuem efeitos pronunciados quanto as interações entre sistemas e gêneros; iii) no crescimento de gados de corte usando as estimativas da amostra total com MCD, os genótipos cruzados tiveram maior gasto de energia de mantença e foram mais rápidos de maturação em comparação tanto com os animais de genótipos Britânicos (Bos taurus) e touros Nelores. A técnica de bootstrap não-paramétrica estimou com sucesso as distribuições dos parâmetros (que tiveram distribuição probabilidade normal para maioria dos cenários). Correlação negativa entre os parâmetros de acréscimo de DNA e energia de mantença foram encontrados para animais machos não castrados do sistema extensivo, indicando que foram mais eficientes no uso da energia. A generalização de tal relação ainda demanda estudos mais abrangentes e aprofundados. / Model parameter fitting and parameter correlation structures are important for characterize a system\'s behaviour in response to multiple scenarios variations (climate, genotypes, nutritional diet and other factors). The aim was to contribute to statistical inferences in the production of CO2 [a potent greenhouse gas (GHG)] in vitro fermentation of alfalfa hay, comparing frequentist methods with new methodologies that emerged in the scientific literature, such as the combination of a delay Rejection and the Adaptive Metropolis methods (RAMA), not yet tested for in vitro fermentation gases predictions. In addition, time series models were used to predict CO2 production in the in vitro fermentation of alfalfa hay. For the first time, individual predictions of animal weight gain rate and energy of maintenance based on the growth dynamics and body composition Davis Growth Model (DGM) was carried out besides multivariate covariance analysis of different scenarios (genres, systems and crossed genotype). Additionally, parameter estimation based on sample of each scenario, using frequentist analysis, nonparametric bootstrap and Monte Carlo simulations were performed with national data (cross breed) and compared to the original estimates of the model obtained with British breeds (Bos taurus). The main criteria used to evaluate the model accuracy were the Mean Square Error of Prediction (MSEP), the Akaike Information Criterion (AIC) and the Bayesian Information Criterion (BIC). The results not only contribute to the scientific progress in modeling animal growth and composition, but also assist the beef industry and farmers to adjust the production process to the meat market specifications, both nationally and internationally. For in vitro gas production, we concluded that the ARIMA (1, 1, 2) model presented the highest accuracy in predicting cumulative CO2 production and the respective equation is recommended for estimating CO2 production. In the growth beef cattle using the individual estimates DGM, average vectors from maintenance of energy and protein accretion showed pronounced effects as the interactions between systems and genres. Also the total sample with DGM, cross-breed genotypes had higher maintenance energy expenditure and were faster-maturing compared with British genotypes animals(Bos taurus) and Nellore bulls estimates. Bootstrap nonparametric with downhill simplex optimization method successfully estimated the distributions of the parameters (that had normal probability distribution for most scenarios). Uncastrated male animals of the extensive system showed negative correlation between the protein deposition rate and requirement for energy maintenance, indicating that animals with faster lean tissue deposition were also more efficient in energy usage. We warn that the generalization of this finding demands studies with larger populations.

Orientador: Luiz Koodi Hotta / Dissertação (mestrado) - Universidade Estadual de Campinas, Instituto de Matemática, Estatística e Computação Científica / Made available in DSpace on 2018-08-20T22:42:22Z (GMT). No. of bitstreams: 1 TruciosMaza_CarlosCesar_M.pdf: 13820849 bytes, checksum: 0cc000af0d7cb7cb6ee6c05ef9f3afbd (MD5) Previous issue date: 2012 / Resumo: Mercados financeiros têm mostrado um grande interesse em intervalos de previsão como uma medida de incerteza. Além das previsões do nível, a previsão da volatilidade é importante em várias aplicações em finanças. O modelo GARCH tem sido bastante utilizado na modelagem da volatilidade. A partir deste modelo, outros modelos foram propostos para incorporar outros fatos estilizados, como o efeito de alavancagem. Neste sentido, temos os modelos EGARCH e GJR-GARCH. Os métodos tradicionais de construção de intervalos de previsão para séries temporais geralmente assumem que os parâmetros do modelo são conhecidos e os erros normais. Quando estas suposições não são verdadeiras, o que costuma acontecer na prática, o intervalo de previsão obtido tenderá a ter uma cobertura abaixo da nominal. Nesta dissertação propomos uma adaptação do algoritmo PRR (Pascual, Romo e Ruiz) desenvolvido para obter intervalos de previsão em modelos GARCH, para obter intervalos de previsão em modelos EGARCH e GJR-GARCH. As adaptações feitas são analisadas através de experimentos Monte Carlo e verifica-se que tiveram bom desempenho apresentando valores da cobertura estimada próximos da cobertura nominal. As adaptações propostas assim como o algoritmo PRR são aplicadas para obter intervalos de previsão dos retornos e das volatilidades para a série de retornos da Ibovespa e para a série NYSE COMPOSITE(DJ) da bolsa de valores de Nova Iorque, obtendo em ambos os casos resultados satisfatórios / Abstract: Financial Markets have shown a big interest in forecast intervals (prediction intervals) as a uncertain measure. Besides the level prediction, the prediction of the volatility is very important in many financial applications. The GARCH model, has been very used in volatility modeling. From this model, other have been proposed to incorporate other stylized facts, such as the leverage effect. In this sense, we have the EGARCH and GJR-GARCH models. Traditional methods for constructing predictions intervals for time series generally assume that the model parameters are known and the erros are normal. When these assumptions are not true, that it is very often in practice, the obtained prediction interval, will tend to have a cover under the nominal. In this theses we propose an adaptation of the PRR (Pascual, Romo and Ruiz) algorithm developed to obtain prediction intervals in GARCH models, to obtain prediction intervals in EGARCH and GJR-GARCH models. These adaptations are analized through Monte-Carlo experiments and It was verified that they have a good performance showing estimated cover values close to the nominal cover. The proposed adaptations, such as the PRR algorithm are applied to obtain prediction intervals from the returns and volatilities for the Ibovespa return series and for the New York Stock Exchange NYSE COMPOSITE(DJ) series, obtaining satisfactory results in both cases / Mestrado / Estatistica / Mestre em Estatística

[pt] O bootstrap de B. Efron, que não poderia ser imaginado sem os computadores de hoje, pode resolver vários problemas livre da suposição de Gaussianidade para os dados. Este trabalho tem o objetivo de apresentar essa técnica computacionalmente intensiva no contexto de Séries temporais - Metodologia Box and Jenkins. Como se sabe essa Metodologia possui alguns resultados assintóticos. Então, na fase da identificação da estrutura do modelo, pode apresentar problemas em regiões do espaço paramétrico aqui determinadas,. O bootstrap é proposto como opção e um estudo de simulação, comparativo, é apresentado. Constrói- se a distribuição bootstrap da autocorrelação e autocorrelação parcial, amostrais, e ainda a distribuição bootstrap do estimador de MQNL dos coeficientes de modelos ARMA (p, q). consequentemente, fica disponí­vel medida não- paramétrica da precisão da estimativa. O estudo de simulação que aborda o estimador de MQNL dos coeficientes enfoca, basicamente, a região de fronteira da estacionariedade e inversibilidade. / [en] The bootstrap of B. Efron, what should not be imagined without fast andcheaper computation, can solve several problems free from assumption that the data conform to a bell-shaped curve. This work has the aim to present this computer-intensive technics in the context of Time Series - Box and Jenkins´s Methodology. As we know this methodology own some asymptotic results. Then in the identification stage of the structure of the model it may present some troubles on regions of the parametric space, as we show later on the bootstrap is proposed as an aption and a comparative simulation study is pointed out. We build up the bootstrap distribution of the sample autocorrelation and sample partial autocorrelation, and yet a bootstrap distribution to the non-linear LS estimator of the coefficients to the ARMA (p,q) model. As a consequence we get the non- parametric measure of the accuracy of the estimates. The study of simulation wich takes into account the non-linear LS estimato to the coefficients, actually focalize the borden of the stationarity and invertibility region.

Une classe générale de statistiques de rangs basées sur la fonction caractéristique est introduite afin de tester l'hypothèse composite d'appartenance à une famille de copules multidimensionnelles. Ces statistiques d'adéquation sont définies comme des distances fonctionnelles de type L_2 pondérées entre une version non paramétrique et une version semi-paramétrique de la fonction caractéristique que l'on peut associer à une copule. Il est démontré que ces statistiques de test se comportent asymptotiquement comme des V-statistiques dégénérées d'ordre quatre et que leurs lois limites s'expriment en termes de sommes pondérées de variables khi-deux indépendantes. La convergence des tests sous des alternatives générales est établie, de même que la validité du bootstrap paramétrique pour le calcul de valeurs critiques. Le comportement des nouveaux tests sous des tailles d'échantillons faibles et modérées est étudié à l'aide de simulations et est comparé à celui d'un test concurrent fondé sur la copule empirique. La méthodologie est finalement illustrée sur un jeu de données à plusieurs dimensions.

Cette these examine deux aspects des theories conformes des champs (TCC) en d dimensions.Sa premiere parti est dediee aux blocs conformes, des fonctions speciales qui contribuent au developpement en ondes partielles des fonctions a quatre points dans les TCC. On montre que ces blocs admettent un developpement en coordonnees polaires dont les coecients se calculent par une recurrence. Les blocs conformes sont naturellement denis sur le plan complexe : on considere alors leur restriction a l'axe r eel, an de montrer qu'ils obeissent une equation dierentielle sur ce domaine, ce qui mene a un algorithme ecace pour calculer les blocs conformes et leurs derivees pour tout d. Quelques applications au programme de bootstrap sont developpees. La seconde partie de cette these examine les perturbations d'une TCC par des operateurs pertinents. On etudie de tels ots du groupe de renormalisation en utilisant la Methode de Troncature Conforme (MTC) de Yurov et Zamolodchikov, une methode numerique qui permet de faire des calculs non-perturbatifs en theorie quantique des champs. Deux theories derentes sont considerees : le boson libre avec un terme de masse, et la theorie 4. Pour le dernier cas, les resultats de la MTC mettent en evidence la brisure de symetrie Z2. Finalement, on developpe une methode pour reduire les erreurs de troncature en ajoutant des contre-termes a l'action \nue" de la MTC, suivant des travaux anterieurs en d = 2 dimensions. / This thesis investigates two aspects of Conformal Field Theories (CFTs) in d dimensions. Its rst part is devoted to conformal blocks, special functions that arise in the partial wave expansion of CFT four-point functions. We prove that these conformal blocks admit an expansion in terms of polar coordinates and show that the expansion coecients are determined by recursion relations. Conformal blocks are naturally dened on the complex plane: we study their restriction to the real line, and show that they obey a fourth-order dierential equation there. This ODE can be used to eciently compute conformal blocks and their derivatives in general d. Several applications to the conformal bootstrap program are mentioned. The second half of this thesis investigates RG ows that are dened by perturbing a CFT by a number of relevant operators. We study such ows using the Truncated Conformal Space Approach (TCSA) of Yurov and Zamolodchikov, a numerical method that allows for controlled computations in strongly coupled QFTs. Two dierent RG ows are considered: the free scalar feld deformed by a mass term, and 4 theory. The former is used as a benchmark, in order to compare numerical TCSA results to exact predictions. TCSA results for 4 theory display spontaneous Z2 symmetry breaking at strong coupling: we study the spectrum of this theory both in the Z2-broken and preserved phase, and we compare the critical exponents governing the phase transition to known values. In a separate chapter, we show how truncation errors can be reduced by adding suitable counterterms to the bare TCSA action, following earlier work in d = 2 dimensions.

Finding a well-predicting model is one of the main goals of regression analysis. However, to evaluate a model's prediction abilities, it is a normal practice to use criteria which either do not serve this purpose, or criteria of insufficient reliability. As an alternative, there are relatively new methods which use repeated simulations for estimating an appropriate loss function -- prediction error. Cross-validation and bootstrap belong to this category. This thesis describes how to utilize these methods in order to select a regression model that best predicts new values of the response variable.

La symétrie conforme impose de très fortes contraintes sur les théories quantiques des champs. En deux dimensions, l’algèbre des symétries conformes est infinie, et les théories conformes bidimensionnelles peuvent être complètement résolubles, dans le sens où toutes leurs fonctions de corrélation peuvent être calculées. Ces théories ont un grand domaine d'application, de la théorie des cordes jusqu'aux systèmes critiques en physique statistique, et elles ont été largement étudiées pendant les dernières décennies.Dans cette thèse nous étudions les théories conformes bidimensionnelles dont l’algèbre de symétrie est celle de Virasoro, en suivant l'approche connue sous le nom de bootstrap conforme. Sous l'hypothèse de l'existence de champs dégénérés, nous généralisons le bootstrap conforme analytique aux théories avec des spectres non-diagonaux. Nous écrivons les équations qui déterminent les constantes de structure, et nous trouvons des solutions explicites en termes de fonctions spéciales. Nous validons ces résultats en faisant des calculs numériques des fonctions de corrélation à quatre points dans des modèles minimaux diagonaux et non-diagonaux, et en vérifiant que la symétrie de croisement est respectée.En outre, nous construisons une proposition pour une famille de théories conformes non-diagonales et non-rationnelles pour toute charge centrale telle que Re(c) < 13. Cette proposition est motivée par les limites des spectres des modèles minimaux de la série D. Nous réalisons des calculs numériques des fonctions à quatre points dans ces théories, et nous trouvons qu'elles obéissent à la symétrie de croisement. Ces théories peuvent être interprétées comme des extensions non-diagonales de la théorie de Liouville. / Conformal symmetry imposes very strong constraints on quantum field theories. In two dimensions, the conformal symmetry algebra is infinite-dimensional, and two-dimensional conformal field theories can be completely solvable, in the sense that all their correlation functions may be computed. These theories have an ample range of applications, from string theory to critical phenomena in statistical physics, and they have been widely studied during the last decades.In this thesis we study two-dimensional conformal field theories with Virasoro algebra symmetry, following the conformal bootstrap approach. Under the assumption that degenerate fields exist, we provide an extension of the analytic conformal bootstrap method to theories with non-diagonal spectrums. We write the equations that determine structure constants, and find explicit solutions in terms of special functions. We validate this results by numerically computing four-point functions in diagonal and non-diagonal minimal models, and verifying that crossing symmetry is satisfied.In addition, we build a proposal for a family of non-diagonal, non-rational conformal field theories for any central charges such that Re(c) < 13. This proposal is motivated by taking limits of the spectrum of D-series minimal models. We perform numerical computations of four-point functions in these theories, and find that they satisfy crossing symmetry. These theories may be understood as non-diagonal extensions of Liouville theory.

Point estimates of the Standard Chain Ladder method (CLM) and of the more complex Munich Chain Ladder method (MCL) are compared to real data on 38 different datasets in order to evaluate if MCL produces better predictions on average with a dataset from an arbitrary insurance portfolio. MCL is also examined to determine if the future paid and incurred claims converge as time progresses. A bootstrap model based on MCL (BMCL) is examined in order to evaluate its possibility to estimate the probability density function (PDF) of future claims and observable claim development results (OCDR). The results show that the paid and incurred predictions by MCL converge. The results also show that when considering all datasets MCL produce on average better estimations than CLM with paid data but no improvement can be seen with incurred data. Further the results show that by considering a subset of datasets which fulfil certain criteria, or by only considering accident years after 1999 the percentage of datasets in which MCL produce superior estimations increases. When examining BMCL one finds that it can produce estimated PDFs of ultimate reserves and OCDRs, however the mean of estimate of ultimate reserves does not converge to the MCL estimates nor do the mean of the OCDRs converge to zero. In order to get the right convergence the estimated OCDR PDFs are centered and the mean of the BMCL estimated ultimate reserve is set to the MCL estimate by multiplication. / Punktskattningar gjorda med Standard Chain Ladder (CLM) och den mer komplexa Munich Chain Ladder-metoden (MCL) jämförs med verklig data för 38 olika dataset för att evaluera om MCL ger bättre prediktioner i genomsnitt än CLM för en godtycklig försäkringsportfölj. MCLs prediktioner undersöks också för att se om de betalda och de kända skadekostnaderna konvergerar. En bootstrapmodell baserad på MCL (BMCL) undersöks för att utvärdera om möjligheterna att estimera täthetsfunktionen (probability density function, PDF) av framtida skadekostnader och av ”observable claim development results (OCDR)”. Resultaten visar att MCLs estimerade betalda och kända skadekostnader konvergerar. Resultaten visar även att när man evaluerar alla dataseten så ger MCL i genomsnitt bättre prediktioner än CLM med betald data, men ingen förbättring kan ses med CLM med känd skadekostnadsdata. Vidare visar resultaten även att genom att bara titta på dataset som uppfyller vissa krav, eller genom att bara använda olycksår efter 1999, så ökar andelen dataset där MCL ger bättre prediktioner än CLM.Vid evaluering av BMCL ser man att den kan producera estimerade PDF:er för ultimo-reserver och OCDR:er, men att medelvärdet av ultimo-reserv prediktionerna från BMCL inte konvergerar mot MCL-prediktionerna och att medelvärdet av OCDR:erna inte konvergerar mot noll. För att få rätt konvergens så centreras OCDR PDF:erna och ultimo-reservernas medelvärden sätts till motsvarande MCL-prediktionens värde genom multiplikation.

Skewness and kurtosis are measures used to describe shape characteristics of distributions. In this thesis, we examine the interval estimates about the skewness and kurtosis by using jackknife empirical likelihood (JEL), adjusted JEL, extended JEL, traditional bootstrap, percentile bootstrap, and BCa bootstrap methods. The limiting distribution of the JEL ratio is the standard chi-squared distribution. The simulation study of this thesis makes a comparison of different methods in terms of the coverage probabilities and interval lengths under the standard normal distribution and exponential distribution. The proposed adjusted JEL and extended JEL perform better than the other methods. Finally we illustrate the proposed JEL methods and different bootstrap methods with three real data sets.