Applications of the bootstrap to survival analysis

En este trabajo se revisan algunas de las aplicaciones clásicas del bootstrap al análisis de la supervivencia. Se consideran en primer lugar el estimador bootstrap de la varianza y el estimador de la mediana corregido para el sesgo del estimador de Kaplan-Meier de la función de supervivencia. A continuación se consideran algunos aspectos mas recientes, tales como métodos para construir bandas de confianza para el estimador de la funcidn de supervivencia y contrastes aproximados para la comparación de funciones de supervivencia. En ambas situaciones el bootstrap resulta de gran utilidad para la aproximación de 10s valores críticos necesarios.

Aplicació web responsiva amb CodeIgniter i Bootstrap

Els objectius que es pretenen aconseguir mitjançant la realització d'aquest treball són els següents: 1. Posar en pràctica coneixements adquirits al llarg de la carrera sobre enginyeria del programari i tècniques de desenvolupament. 2. Aprofundir en tecnologies responsives per al desenvolupament de solucions de gestió, prenent com a referència una plantilla basada en Bootstrap. 3. Demostrar la validesa del paradigma MVC com a fonament arquitectural d'una aplicació de gestió. 4. Demostrar la solvència dels "frameworks" CodeIgniter i Bootstrap per a la implantació d'aquest tipus de solucions distribuïdes.

Creating and using a non-dedicated HPC cluster with ParallelKnoppix

This note describes ParallelKnoppix, a bootable CD that allows econometricians with average knowledge of computers to create and begin using a high performance computing cluster for parallel computing in very little time. The computers used may be heterogeneous machines, and clusters of up to 200 nodes are supported. When the cluster is shut down, all machines are in their original state, so their temporary use in the cluster does not interfere with their normal uses. An example shows how a Monte Carlo study of a bootstrap test procedure may be done in parallel. Using a cluster of 20 nodes, the example runs approximately 20 times faster than it does on a single computer.

Price level convergence, purchasing power parity and multiple structural breaks: an application to US cities

This article provides a fresh methodological and empirical approach for assessing price level convergence and its relation to purchasing power parity (PPP) using annual price data for seventeen US cities. We suggest a new procedure that can handle a wide range of PPP concepts in the presence of multiple structural breaks using all possible pairs of real exchange rates. To deal with cross-sectional dependence, we use both cross-sectional demeaned data and a parametric bootstrap approach. In general, we find more evidence for stationarity when the parity restriction is not imposed, while imposing parity restriction provides leads toward the rejection of the panel stationar- ity. Our results can be embedded on the view of the Balassa-Samuelson approach, but where the slope of the time trend is allowed to change in the long-run. The median half-life point estimate are found to be lower than the consensus view regardless of the parity restriction.

Inference of distributional parameters from compositional samples containing nondetects

Low concentrations of elements in geochemical analyses have the peculiarity of beingcompositional data and, for a given level of significance, are likely to be beyond thecapabilities of laboratories to distinguish between minute concentrations and completeabsence, thus preventing laboratories from reporting extremely low concentrations of theanalyte. Instead, what is reported is the detection limit, which is the minimumconcentration that conclusively differentiates between presence and absence of theelement. A spatially distributed exhaustive sample is employed in this study to generateunbiased sub-samples, which are further censored to observe the effect that differentdetection limits and sample sizes have on the inference of population distributionsstarting from geochemical analyses having specimens below detection limit (nondetects).The isometric logratio transformation is used to convert the compositional data in thesimplex to samples in real space, thus allowing the practitioner to properly borrow fromthe large source of statistical techniques valid only in real space. The bootstrap method isused to numerically investigate the reliability of inferring several distributionalparameters employing different forms of imputation for the censored data. The casestudy illustrates that, in general, best results are obtained when imputations are madeusing the distribution best fitting the readings above detection limit and exposes theproblems of other more widely used practices. When the sample is spatially correlated, itis necessary to combine the bootstrap with stochastic simulation

KAM theory for conformally symplectic systems

We present a KAM theory for some dissipative systems (geometrically, these are conformally symplectic systems, i.e. systems that transform a symplectic form into a multiple of itself). For systems with n degrees of freedom depending on n parameters we show that it is possible to find solutions with n-dimensional (Diophantine) frequencies by adjusting the parameters. We do not assume that the system is close to integrable, but we use an a-posteriori format. Our unknowns are a parameterization of the solution and a parameter. We show that if there is a sufficiently approximate solution of the invariance equation, which also satisfies some explicit non–degeneracy conditions, then there is a true solution nearby. We present results both in Sobolev norms and in analytic norms. The a–posteriori format has several consequences: A) smooth dependence on the parameters, including the singular limit of zero dissipation; B) estimates on the measure of parameters covered by quasi–periodic solutions; C) convergence of perturbative expansions in analytic systems; D) bootstrap of regularity (i.e., that all tori which are smooth enough are analytic if the map is analytic); E) a numerically efficient criterion for the break–down of the quasi–periodic solutions. The proof is based on an iterative quadratically convergent method and on suitable estimates on the (analytical and Sobolev) norms of the approximate solution. The iterative step takes advantage of some geometric identities, which give a very useful coordinate system in the neighborhood of invariant (or approximately invariant) tori. This system of coordinates has several other uses: A) it shows that for dissipative conformally symplectic systems the quasi–periodic solutions are attractors, B) it leads to efficient algorithms, which have been implemented elsewhere. Details of the proof are given mainly for maps, but we also explain the slight modifications needed for flows and we devote the appendix to present explicit algorithms for flows.

CIAU : Aplicació web per a gestió d'incidències en entorns virtuals

Aplicació web per entrar i gestionar incidències desenvolupada amb CodeIgniter (PHP) i Bootstrap.

Stepwise multiple testing as formalized data snooping

It is common in econometric applications that several hypothesis tests arecarried out at the same time. The problem then becomes how to decide whichhypotheses to reject, accounting for the multitude of tests. In this paper,we suggest a stepwise multiple testing procedure which asymptoticallycontrols the familywise error rate at a desired level. Compared to relatedsingle-step methods, our procedure is more powerful in the sense that itoften will reject more false hypotheses. In addition, we advocate the useof studentization when it is feasible. Unlike some stepwise methods, ourmethod implicitly captures the joint dependence structure of the teststatistics, which results in increased ability to detect alternativehypotheses. We prove our method asymptotically controls the familywise errorrate under minimal assumptions. We present our methodology in the context ofcomparing several strategies to a common benchmark and deciding whichstrategies actually beat the benchmark. However, our ideas can easily beextended and/or modied to other contexts, such as making inference for theindividual regression coecients in a multiple regression framework. Somesimulation studies show the improvements of our methods over previous proposals. We also provide an application to a set of real data.

Improved nonparametric confidence intervals in time series regressions

Condence intervals in econometric time series regressions suffer fromnotorious coverage problems. This is especially true when the dependencein the data is noticeable and sample sizes are small to moderate, as isoften the case in empirical studies. This paper suggests using thestudentized block bootstrap and discusses practical issues, such as thechoice of the block size. A particular data-dependent method is proposedto automate the method. As a side note, it is pointed out that symmetricconfidence intervals are preferred over equal-tailed ones, since theyexhibit improved coverage accuracy. The improvements in small sampleperformance are supported by a simulation study.

Exact and approximate stepdown methods for multiple hypothesis testing

Consider the problem of testing k hypotheses simultaneously. In this paper,we discuss finite and large sample theory of stepdown methods that providecontrol of the familywise error rate (FWE). In order to improve upon theBonferroni method or Holm's (1979) stepdown method, Westfall and Young(1993) make eective use of resampling to construct stepdown methods thatimplicitly estimate the dependence structure of the test statistics. However,their methods depend on an assumption called subset pivotality. The goalof this paper is to construct general stepdown methods that do not requiresuch an assumption. In order to accomplish this, we take a close look atwhat makes stepdown procedures work, and a key component is a monotonicityrequirement of critical values. By imposing such monotonicity on estimatedcritical values (which is not an assumption on the model but an assumptionon the method), it is demonstrated that the problem of constructing a validmultiple test procedure which controls the FWE can be reduced to the problemof contructing a single test which controls the usual probability of a Type 1error. This reduction allows us to draw upon an enormous resamplingliterature as a general means of test contruction.

Statistics in archaeology: New directions

Connections between Statistics and Archaeology have always appeared veryfruitful. The objective of this paper is to offer an outlook of somestatistical techniques that are being developed in the most recentyears and that can be of interest for archaeologists in the short run.

Validation procedures in radiological diagnostic models. Neural network and logistic regression

The objective of this paper is to compare the performance of twopredictive radiological models, logistic regression (LR) and neural network (NN), with five different resampling methods. One hundred and sixty-seven patients with proven calvarial lesions as the only known disease were enrolled. Clinical and CT data were used for LR and NN models. Both models were developed with cross validation, leave-one-out and three different bootstrap algorithms. The final results of each model were compared with error rate and the area under receiver operating characteristic curves (Az). The neural network obtained statistically higher Az than LR with cross validation. The remaining resampling validation methods did not reveal statistically significant differences between LR and NN rules. The neural network classifier performs better than the one based on logistic regression. This advantage is well detected by three-fold cross-validation, but remains unnoticed when leave-one-out or bootstrap algorithms are used.

Desenvolupament d'una aplicació per a la creació de pressupostos on-line

En aquests moments, el mercat no es troba en una bona situació, per aquest motiu les empreses han de buscar noves maneres de créixer, expandir-se i noves formes d’interactuar amb els clients. La idea original d’aquest projecte sorgeix de la necessitat de disposar d’una manera diferent de promocionar-se i oferir nous serveis a través d’internet mitjançant una pàgina web. Degut a la situació actual, el preu és un aspecte molt important i influent a l’hora de realitzar una obra. Per aquest motiu es va pensar que seria molt interessant que el client pogués demanar de forma fàcil i ràpida un pressupost, i a l’instant tingués un preu orientatiu del que li pot costar la obra. D’aquesta manera l’interessat s’estalvia i agilitza molts passos previs abans de començar una obra. Després d’analitzar quina era la millor manera de portar a terme el projecte informàtic, s’ha determinat que la pàgina web es desenvoluparà utilitzant els llenguatges HTML i PHP combinant-lo amb el framework CodeIgniter. El disseny de la web es realitzarà mitjançant fulles d’estil CSS conjuntament amb el framework BootStrap. Per realitzar l’aplicació web que realitza els pressupostos s’utilitzarà AJAX i jQuery perquè d’aquesta manera el procés sigui dinàmic. L’entorn de desenvolupament escollit és el NetBeans i per provar el projecte s’utilitza el XAMPP. Un usuari només necessitarà un navegador i connexió a internet per fer servir totes les funcions de la web. Podrà realitzar pressupostos, concertar visites, contactar i per suposat veure tota la part informativa de la pàgina.

Trends in indices of daily precipitation extremes in Catalonia (NE Spain), 1951-2003

The aim of this paper is to quantitatively characterize the climatology of daily precipitation indices in Catalonia (northeastern Iberian Peninsula) from 1951 to 2003. This work has been performed analyzing a subset of the ETCCDI (Expert Team on Climate Change Detection and Indices) precipitation indices calculated from a new interpolated dataset of daily precipitation, namely SPAIN02, regular at 0.2° horizontal resolution (around 20 km) and from two high-quality stations: the Ebro and Fabra observatories. Using a jack-knife technique, we have found that the sampling error of the SPAIN02 regional averaged is relatively low. The trend analysis has been implemented using a Circular Block Bootstrap procedure applicable to non-normal distributions and autocorrelated series. A running trend analysis has been applied to analyze the trend persistence. No general trends at a regional scale are observed, considering the annual or the seasonal regional averaged series of all the indices for all the time windows considered. Only the consecutive dry days index (CDD) at annual scale shows a locally coherent spatial trend pattern; around 30% of the Catalonia area has experienced an increase of around 2¿3 days decade¿1. The Ebro and Fabra observatories show a similar CDD trend, mainly due to the summer contribution. Besides this, a significant decrease in total precipitation (around ¿10 mm decade¿1) and in the index "highest precipitation amount in five-day period" (RX5DAY, around ¿5 mm decade¿1), have been found in summer for the Ebro observatory.

Anàlisi multivariada craniomètrica de les poblacions catalano-balears del passat

S'ha portat a terme una analisi estadística de les dades craniometriques de les principals poblacions catalano-balears del passat, des de J'Edat del Bronze fins a l'Edat Mitjana, i posteriorment s'han integrat les poblacions en un context peninsular més ampli. S'han analitzalles relacions entre poblacions per mitja d'un analisi d'agrupament (cluster), emprant els algoritmes UPGMA i del "veí més proper" i la distancia euclídea al quadrat. Així mateix, s'ha realitzat un re-mostreig (bootstrap) per comprovar la robustes de les agrupacions obtingudes i un test de Mantel per comprovar el grau de concordan~a entre les poblacions femenines i masculines. Els resultats mostren que, a nivell peninsular, els Bascos són la població més diferenciada de la resta des d'un punt de vista morfologic, la qual cosa pot estar relacionada amb una major antiguitat d'aquesta població, potser present des de temps pre-neolítics. La segona font de variació peninsular la constitueixen dues poblacions amb possibles influencies foranies d'origen africa: jueus i musulmans. La resta de poblacions, incloent-hi les de Mallorca, són notablement homogenies a pesar de la seva heterogene'itat cultural, temporal i geogràtfica. Finalment, es conclou que no hi ha cap analisi estadística que objectivament recolzi I'existencia deIs tipus racials tradicionalmcnt descrits a la Península Iberica, i que responen a una visió tipologica i estatica de la diversitat humana.