Sound Texture Modelling

In the PhD thesis “Sound Texture Modeling” we deal with statistical modelling or textural sounds like water, wind, rain, etc. For synthesis and classification. Our initial model is based on a wavelet tree signal decomposition and the modeling of the resulting sequence by means of a parametric probabilistic model, that can be situated within the family of models trainable via expectation maximization (hidden Markov tree model ). Our model is able to capture key characteristics of the source textures (water, rain, fire, applause, crowd chatter ), and faithfully reproduces some of the sound classes. In terms of a more general taxonomy of natural events proposed by Graver, we worked on models for natural event classification and segmentation. While the event labels comprise physical interactions between materials that do not have textural propierties in their enterity, those segmentation models can help in identifying textural portions of an audio recording useful for analysis and resynthesis. Following our work on concatenative synthesis of musical instruments, we have developed a pattern-based synthesis system, that allows to sonically explore a database of units by means of their representation in a perceptual feature space. Concatenative syntyhesis with “molecules” built from sparse atomic representations also allows capture low-level correlations in perceptual audio features, while facilitating the manipulation of textural sounds based on their physical and perceptual properties. We have approached the problem of sound texture modelling for synthesis from different directions, namely a low-level signal-theoretic point of view through a wavelet transform, and a more high-level point of view driven by perceptual audio features in the concatenative synthesis setting. The developed framework provides unified approach to the high-quality resynthesis of natural texture sounds. Our research is embedded within the Metaverse 1 European project (2008-2011), where our models are contributting as low level building blocks within a semi-automated soundscape generation system.

Multiple-scale hydrological and geophysical data integration through non-linear Bayesian sequential simulation

Significant progress has been made with regard to the quantitative integration of geophysical and hydrological data at the local scale. However, extending the corresponding approaches to the scale of a field site represents a major, and as-of-yet largely unresolved, challenge. To address this problem, we have developed downscaling procedure based on a non-linear Bayesian sequential simulation approach. The main objective of this algorithm is to estimate the value of the sparsely sampled hydraulic conductivity at non-sampled locations based on its relation to the electrical conductivity logged at collocated wells and surface resistivity measurements, which are available throughout the studied site. The in situ relationship between the hydraulic and electrical conductivities is described through a non-parametric multivariatekernel density function. Then a stochastic integration of low-resolution, large-scale electrical resistivity tomography (ERT) data in combination with high-resolution, local-scale downhole measurements of the hydraulic and electrical conductivities is applied. The overall viability of this downscaling approach is tested and validated by comparing flow and transport simulation through the original and the upscaled hydraulic conductivity fields. Our results indicate that the proposed procedure allows obtaining remarkably faithful estimates of the regional-scale hydraulic conductivity structure and correspondingly reliable predictions of the transport characteristics over relatively long distances.

Fourth order pseudo maximum likelihood methods

We extend PML theory to account for information on the conditional moments up to order four, but without assuming a parametric model, to avoid a risk of misspecification of the conditional distribution. The key statistical tool is the quartic exponential family, which allows us to generalize the PML2 and QGPML1 methods proposed in Gourieroux et al. (1984) to PML4 and QGPML2 methods, respectively. An asymptotic theory is developed. The key numerical tool that we use is the Gauss-Freud integration scheme that solves a computational problem that has previously been raised in several fields. Simulation exercises demonstrate the feasibility and robustness of the methods [Authors]

Changes associated with laboratory rearing in antennal sensilla patterns of Triatoma infestans, Rhodnius prolixus, and Rhodnius pallescens (Hemiptera, Reduviidae, Triatominae)

We examined changes in the array of antennal sensilla of three species of Triatominae (Triatoma infestans, Rhodnius prolixus, and R. pallescens) following their establishment for different periods in laboratory culture. In each case, the laboratory colonies were compared with conspecific samples taken directly from the field, by quantitative analysis of the sensilla arrays on the three distal segments of the antenna in terms of the densities of three types of chemoreceptors (basiconics and thick and thin walled trichoids) and one type of mechanoreceptor (bristles). Sensilla densities were compared by ANOVA or non-parametric tests, and by multivariate discriminant analysis. Strains of the same species reared in different laboratories showed significant differences in their sensilla arrays, especially when compared to field-collected material from the same geographic origin. A Bolivian strain of T. infestans reared in the laboratory for 15 years and fed at monthly intervals, showed greatest differences from its conspecific wild forms, especially in terms of reductions in the number of chemoreceptors. By contrast, an Argentine strain of T. infestans reared for 25 years in the laboratory and fed weekly, showed a relative increase in the density of mechanoreceptors. A Colombian strain of R. prolixus reared for 20 years and fed weekly or fortnightly, showed only modest differences in the sensilla array when compared to its wild populations from the same area. However, a Colombian strain of R. pallescens reared for 12 years and fed fortnightly, did show highly significant reductions in one form of chemoreceptor compared to its conspecific wild populations. For all populations, multivariate analysis clearly discriminated between laboratory and field collected specimens, suggesting that artificial rearing can lead to modifications in the sensory array. This not only supports the idea of morphological plasticity in these species, but also suggests caution in the use of long-established laboratory material for experimental studies designed to extrapolate the natural behaviour and physiology of these species.

An empirical hierarchical Bayesian unification of occupational exposure assessment methods.

In occupational exposure assessment of airborne contaminants, exposure levels can either be estimated through repeated measurements of the pollutant concentration in air, expert judgment or through exposure models that use information on the conditions of exposure as input. In this report, we propose an empirical hierarchical Bayesian model to unify these approaches. Prior to any measurement, the hygienist conducts an assessment to generate prior distributions of exposure determinants. Monte-Carlo samples from these distributions feed two level-2 models: a physical, two-compartment model, and a non-parametric, neural network model trained with existing exposure data. The outputs of these two models are weighted according to the expert's assessment of their relevance to yield predictive distributions of the long-term geometric mean and geometric standard deviation of the worker's exposure profile (level-1 model). Bayesian inferences are then drawn iteratively from subsequent measurements of worker exposure. Any traditional decision strategy based on a comparison with occupational exposure limits (e.g. mean exposure, exceedance strategies) can then be applied. Data on 82 workers exposed to 18 contaminants in 14 companies were used to validate the model with cross-validation techniques. A user-friendly program running the model is available upon request.

HIV behavioural surveillance in Switzerland: trends in drug consumption, injection and syringe distribution, 1993-2011

Background: As part of the second generation surveillance system for HIV/Aids in Switzerland, repeated cross-sectional surveys were conducted in 1993, 1994, 1996, 2000, 2006 and 2011 among attenders of all low threshold facilities (LTFs) with needle exchange programmes and/or supervised drug consumption rooms for injection or inhalation. The number of syringes distributed to the injectors has also been measured annually since 2000. Distribution in other settings, such as pharmacies, is also monitored nationally. Methods: Periodic surveys of LTFs have been conducted using an interviewer/self-administered questionnaire structured along four themes: socio-demographic characteristics, drug consumption, risk/preventive behaviour and health. Analysis is restricted to attenders who had injected drugs during their lifetime (IDU´s). Pearson's chi-square test and trend analysis were conducted on annual aggregated data. Trend significance was assessed using Stata's non parametric test nptrend. Results: Median age of IDU´s increased from 26 years in 1993 to 40 in 2011; most are men (78%). Total yearly number of syringes distributed by LTFs has decreased by 44% in 10 years. Use of cocaine has increased (Table 1). Injection, regular use of heroin and borrowing of syringes/needles have decreased, while sharing of other material remains stable. There are fewer new injectors; more IDU´s report substitution treatment. Most attenders had ever been tested for HIV (90% in 1993, 94% in 2011). Reported prevalence of HIV remained stable around 10%; that of HCV decreased from 62% in 2000 to 42% in 2011. Conclusions: Overall, findings indicate a decrease in injection as a means of drug consumption in that population. This interpretation is supported by data from other sources, such as a national decrease in distribution from other delivery points. Switzerland's behavioural surveillance system is sustainable and allows the HIV epidemic to be monitored among this hard-to-reach population, providing information for planning and evaluation.

Early white matter alterations in men predisposed to FXTAS revealed by MT imaging.

Introduction: The Fragile X - associated Tremor Ataxia Syndrome (FXTAS) is a recently described, and under-diagnosed, late onset (≈ 60y) neurodegenerative disorder affecting male carriers of a premutation in the Fragile X Mental Retardation 1 (FMR1) gene. The premutation is an CGG (Cytosine-Guanine-Guanine) expansion (55 to 200 CGG repeats) in the proximal region of the FMR1 gene. Patients with FXTAS primarily present with cerebellar ataxia and intention tremor. Neuroradiological features of FXTAS include prominent white matter disease in the periventricular, subcortical, middle cerebellar peduncles and deep white matter of the cerebellum on T2-weighted or FLAIR MR imaging (Jacquemmont 2007, Loesch 2007, Brunberg 2002, Cohen 2006). We hypothesize that a significant white matter alteration is present in younger individuals many years prior to clinical symptoms and/or the presence of visible lesions on conventional MR sequences and might be detectable by magnetization transfer (MT) imaging. Methods: Eleven asymptomatic premutation carriers (mean age = 55 years) and seven intra-familial controls participated to the study. A standardized neurological examination was performed on all participants and a neuropsychological evaluation was carried out before MR scanning performed on a 3T Siemens Trio. The protocol included a sagittal T1-weighted 3D gradient-echo sequence (MPRAGE, 160 slices, 1 mm^3 isotropic voxels) and a gradient-echo MTI (FA 30, TE 15, matrix size 256*256, pixel size 1*1 mm, 36 slices (thickness 2mm), MT pulse duration 7.68 ms, FA 500, frequency offset 1.5 kHz). MTI was performed by acquiring consecutively two set of images; first with and then without the MT saturation pulse. MT images were coregistered to the T1 acquisition. The MTR for every intracranial voxel was calculated as follows: MTR = (M0 - MS)/M0*100%, creating a MTR map for each subject. As first analysis, the whole white matter (WM) was used to mask the MTR image in order to create an histogram of the MTR distribution in the whole tissue class over the two groups examined. Then, for each subject, we performed a segmentation and parcellation of the brain by means of Freesurfer software, starting from the high resolution T1-weighted anatomical acquisition. Cortical parcellations was used to assign a label to the underlying white matter by the construction of a Voronoi diagram in the WM voxels of the MR volume based on distance to the nearest cortical parcellation label. This procedure allowed us to subdivide the cerebral WM in 78 ROIs according to the cortical parcellation (see example in Fig 1). The cerebellum, by the same procedure, was subdivided in 5 ROIs (2 per each hemisphere and one corresponding to the brainstem). For each subject, we calculated the mean value of MTR within each ROI and averaged over controls and patients. Significant differences between the two groups were tested using a two sample T-test (p<0.01). Results: Neurological examination showed that no patient met the clinical criteria of Fragile X Tremor and Ataxia Syndrome yet. Nonetheless, premutation carriers showed some subtle neurological signs of the disorder. In fact, premutation carriers showed a significant increase of tremor (CRST, T-test p=0.007) and increase of ataxia (ICARS, p=0.004) when compared to controls. The neuropsychological evaluation was normal in both groups. To obtain general characterizations of myelination for each subject and premutation carriers, we first computed the distribution of MTR values across the total white matter volume and averaged for each group. We tested the equality of the two distributions with the non parametric Kolmogorov-Smirnov test and we rejected the null-hypothesis at a p=0.03 (fig. 2). As expected, when comparing the asymptomatic permutation carriers with control subjects, the peak value and peak position of the MTR values within the whole WM were decreased and the width of the distribution curve was increased (p<0.01). These three changes point to an alteration of the global myelin status of the premutation carriers. Subsequently, to analyze the regional myelination and white matter integrity of the same group, we performed a ROI analysis of MTR data. The ROI-based analysis showed a decrease of mean MTR value in premutation carriers compared to controls in bilateral orbito-frontal and inferior frontal WM, entorhinal and cingulum regions and cerebellum (Fig 3). The detection of these differences in these regions failed with other conventional MR techniques. Conclusions: These preliminary data confirm that in premutation carriers, there are indeed alterations in "normal appearing white matter" (NAWM) and these alterations are visible with the MT technique. These results indicate that MT imaging may be a relevant approach to detect both global and local alterations within NAWM in "asymptomatic" carriers of premutations in the Fragile X Mental Retardation 1 (FMR1) gene. The sensitivity of MT in the detection of these alterations might point towards a specific physiopathological mechanism linked to an underlying myelin disorder. ROI-based analyses show that the frontal, parahippocampal and cerebellar regions are already significantly affected before the onset of symptoms. A larger sample will allow us to determine the minimum CGG expansion and age associated with these subclinical white matter alterations.

The serine repeat antigen (SERA) gene family phylogeny in Plasmodium: the impact of GC content and reconciliation of gene and species trees.

Plasmodium falciparum is the parasite responsible for the most acute form of malaria in humans. Recently, the serine repeat antigen (SERA) in P. falciparum has attracted attention as a potential vaccine and drug target, and it has been shown to be a member of a large gene family. To clarify the relationships among the numerous P. falciparum SERAs and to identify orthologs to SERA5 and SERA6 in Plasmodium species affecting rodents, gene trees were inferred from nucleotide and amino acid sequence data for 33 putative SERA homologs in seven different species. (A distance method for nucleotide sequences that is specifically designed to accommodate differing GC content yielded results that were largely compatible with the amino acid tree. Standard-distance and maximum-likelihood methods for nucleotide sequences, on the other hand, yielded gene trees that differed in important respects.) To infer the pattern of duplication, speciation, and gene loss events in the SERA gene family history, the resulting gene trees were then "reconciled" with two competing Plasmodium species tree topologies that have been identified by previous phylogenetic studies. Parsimony of reconciliation was used as a criterion for selecting a gene tree/species tree pair and provided (1) support for one of the two species trees and for the core topology of the amino acid-derived gene tree, (2) a basis for critiquing fine detail in a poorly resolved region of the gene tree, (3) a set of predicted "missing genes" in some species, (4) clarification of the relationship among the P. falciparum SERA, and (5) some information about SERA5 and SERA6 orthologs in the rodent malaria parasites. Parsimony of reconciliation and a second criterion--implied mutational pattern at two key active sites in the SERA proteins-were also seen to be useful supplements to standard "bootstrap" analysis for inferred topologies.

Genetic diversity of environmental Aspergillus flavus strains in the state of São Paulo, Brazil by random amplified polymorphic DNA

Aspergillus flavus is a very important toxigenic fungus that produces aflatoxins, a group of extremely toxic substances to man and animals. Toxigenic fungi can grow in feed crops, such as maize, peanuts, and soybeans, being thus of high concern for public health. There are toxigenic and non-toxigenic A. flavus variants, but the necessary conditions for expressing the toxigenic potential are not fully understood. Therefore, we have studied total-DNA polymorphism from toxigenic and non toxigenic A. flavus strains isolated from maize crops and soil at two geographic locations, 300 km apart, in the Southeast region of Brazil. Total DNA from each A. flavus isolate was extracted and subjected to polymerase chain reaction amplification with five randomic primers through the RAPD (random amplified polymorphic DNA) technique. Phenetic and cladistic analyses of the data, based on bootstrap analyses, led us to conclude that RAPD was not suitable to discriminate toxigenic from non toxigenic strains. But the present results support the use of RAPD for strain characterization, especially for preliminary evaluation over extensive collections.

Effect of inspiratory threshold loading on ventilatory kinetics during constant-load exercise.

Humoral factors play an important role in the control of exercise hyperpnea. The role of neuromechanical ventilatory factors, however, is still being investigated. We tested the hypothesis that the afferents of the thoracopulmonary system, and consequently of the neuromechanical ventilatory loop, have an influence on the kinetics of oxygen consumption (VO2), carbon dioxide output (VCO2), and ventilation (VE) during moderate intensity exercise. We did this by comparing the ventilatory time constants (tau) of exercise with and without an inspiratory load. Fourteen healthy, trained men (age 22.6 +/- 3.2 yr) performed a continuous incremental cycle exercise test to determine maximal oxygen uptake (VO2max = 55.2 +/- 5.8 ml x min(-1) x kg(-1)). On another day, after unloaded warm-up they performed randomized constant-load tests at 40% of their VO2max for 8 min, one with and the other without an inspiratory threshold load of 15 cmH2O. Ventilatory variables were obtained breath by breath. Phase 2 ventilatory kinetics (VO2, VCO2, and VE) could be described in all cases by a monoexponential function. The bootstrap method revealed small coefficients of variation for the model parameters, indicating an accurate determination for all parameters. Paired Student's t-tests showed that the addition of the inspiratory resistance significantly increased the tau during phase 2 of VO2 (43.1 +/- 8.6 vs. 60.9 +/- 14.1 s; P < 0.001), VCO2 (60.3 +/- 17.6 vs. 84.5 +/- 18.1 s; P < 0.001) and VE (59.4 +/- 16.1 vs. 85.9 +/- 17.1 s; P < 0.001). The average rise in tau was 41.3% for VO2, 40.1% for VCO2, and 44.6% for VE. The tau changes indicated that neuromechanical ventilatory factors play a role in the ventilatory response to moderate exercise.

Addition of inspiratory resistance increases the amplitude of the slow component of O2 uptake kinetics.

The contribution of respiratory muscle work to the development of the O(2) consumption (Vo(2)) slow component is a point of controversy because it has been shown that the increased ventilation in hypoxia is not associated with a concomitant increase in Vo(2) slow component. The first purpose of this study was thus to test the hypothesis of a direct relationship between respiratory muscle work and Vo(2) slow component by manipulating inspiratory resistance. Because the conditions for a Vo(2) slow component specific to respiratory muscle can be reached during intense exercise, the second purpose was to determine whether respiratory muscles behave like limb muscles during heavy exercise. Ten trained subjects performed two 8-min constant-load heavy cycling exercises with and without a threshold valve in random order. Vo(2) was measured breath by breath by using a fast gas exchange analyzer, and the Vo(2) response was modeled after removal of the cardiodynamic phase by using two monoexponential functions. As anticipated, when total work was slightly increased with loaded inspiratory resistance, slight increases in base Vo(2), the primary phase amplitude, and peak Vo(2) were noted (14.2%, P < 0.01; 3.5%, P > 0.05; and 8.3%, P < 0.01, respectively). The bootstrap method revealed small coefficients of variation for the model parameter, including the slow-component amplitude and delay (15 and 19%, respectively), indicating an accurate determination for this critical parameter. The amplitude of the Vo(2) slow component displayed a 27% increase from 8.1 +/- 3.6 to 10.3 +/- 3.4 ml. min(-1). kg(-1) (P < 0.01) with the addition of inspiratory resistance. Taken together, this increase and the lack of any differences in minute volume and ventilatory parameters between the two experimental conditions suggest the occurrence of a Vo(2) slow component specific to the respiratory muscles in loaded condition.

Echovirus 30 associated with cases of aseptic meningitis in state of Pará, Northern Brazil

Investigation of the aetiology of viral meningitis in Brazil is most often restricted to cases that occur in the Southern and Southeastern Regions; therefore, the purpose of this study is to describe the viral meningitis cases that occurred in state of Pará, Northern Brazil, from January 2005-December 2006. The detection of enterovirus (EV) in cerebrospinal fluid was performed using cell culture techniques, RT-PCR, nested PCR and nucleotide sequencing. The ages of the 91 patients ranged from < one year old to > 60 years old (median age 15.90 years). Fever (87.1%), headache (77.0%), vomiting (61.5%) and stiffness (61.5%) were the most frequent symptoms. Of 91 samples analyzed, 18 (19.8%) were positive for EV. Twelve were detected only by RT- PCR followed by nested PCR, whereas six were found by both cell culture and RT-PCR. From the last group, five were sequenced and classified as echovirus 30 (Echo 30). Phylogenetic analyses revealed that Echo 30 detected in Northern Brazil clustered within a unique group with a bootstrap value of 100% and could constitute a new subgroup (4c) according to the phylogenetic tree described by Oberste et al. (1999). This study described the first molecular characterization of Echo 30 in Brazil and this will certainly contribute to future molecular analyses involving strains detected in other regions of Brazil.

The agreement between ipsative and normative questionnaires using compositional data analysis techniques

The main instrument used in psychological measurement is the self-report questionnaire. One of its majordrawbacks however is its susceptibility to response biases. A known strategy to control these biases hasbeen the use of so-called ipsative items. Ipsative items are items that require the respondent to makebetween-scale comparisons within each item. The selected option determines to which scale the weight ofthe answer is attributed. Consequently in questionnaires only consisting of ipsative items everyrespondent is allotted an equal amount, i.e. the total score, that each can distribute differently over thescales. Therefore this type of response format yields data that can be considered compositional from itsinception.Methodological oriented psychologists have heavily criticized this type of item format, since the resultingdata is also marked by the associated unfavourable statistical properties. Nevertheless, clinicians havekept using these questionnaires to their satisfaction. This investigation therefore aims to evaluate bothpositions and addresses the similarities and differences between the two data collection methods. Theultimate objective is to formulate a guideline when to use which type of item format.The comparison is based on data obtained with both an ipsative and normative version of threepsychological questionnaires, which were administered to 502 first-year students in psychology accordingto a balanced within-subjects design. Previous research only compared the direct ipsative scale scoreswith the derived ipsative scale scores. The use of compositional data analysis techniques also enables oneto compare derived normative score ratios with direct normative score ratios. The addition of the secondcomparison not only offers the advantage of a better-balanced research strategy. In principle it also allowsfor parametric testing in the evaluation

Neural Correlates of the Severity of Cocaine, Heroin, Alcohol, MDMA and Cannabis Use in Polysubstance Abusers: A Resting-PET Brain Metabolism Study

INTRODUCTION Functional imaging studies of addiction following protracted abstinence have not been systematically conducted to look at the associations between severity of use of different drugs and brain dysfunction. Findings from such studies may be relevant to implement specific interventions for treatment. The aim of this study was to examine the association between resting-state regional brain metabolism (measured with 18F-fluorodeoxyglucose Positron Emission Tomography (FDG-PET) and the severity of use of cocaine, heroin, alcohol, MDMA and cannabis in a sample of polysubstance users with prolonged abstinence from all drugs used. METHODS Our sample consisted of 49 polysubstance users enrolled in residential treatment. We conducted correlation analyses between estimates of use of cocaine, heroin, alcohol, MDMA and cannabis and brain metabolism (BM) (using Statistical Parametric Mapping voxel-based (VB) whole-brain analyses). In all correlation analyses conducted for each of the drugs we controlled for the co-abuse of the other drugs used. RESULTS The analysis showed significant negative correlations between severity of heroin, alcohol, MDMA and cannabis use and BM in the dorsolateral prefrontal cortex (DLPFC) and temporal cortex. Alcohol use was further associated with lower metabolism in frontal premotor cortex and putamen, and stimulants use with parietal cortex. CONCLUSIONS Duration of use of different drugs negatively correlated with overlapping regions in the DLPFC, whereas severity of cocaine, heroin and alcohol use selectively impact parietal, temporal, and frontal-premotor/basal ganglia regions respectively. The knowledge of these associations could be useful in the clinical practice since different brain alterations have been associated with different patterns of execution that may affect the rehabilitation of these patients.

DAS-28-based EULAR response and HAQ improvement in rheumatoid arthritis patients switching between TNF antagonists

INTRODUCTION No definitive data are available regarding the value of switching to an alternative TNF antagonist in rheumatoid arthritis patients who fail to respond to the first one. The aim of this study was to evaluate treatment response in a clinical setting based on HAQ improvement and EULAR response criteria in RA patients who were switched to a second or a third TNF antagonist due to failure with the first one. METHODS This was an observational, prospective study of a cohort of 417 RA patients treated with TNF antagonists in three university hospitals in Spain between January 1999 and December 2005. A database was created at the participating centres, with well-defined operational instructions. The main outcome variables were analyzed using parametric or non-parametric tests depending on the level of measurement and distribution of each variable. RESULTS Mean (+/- SD) DAS-28 on starting the first, second and third TNF antagonist was 5.9 (+/- 2.0), 5.1 (+/- 1.5) and 6.1 (+/- 1.1). At the end of follow-up, it decreased to 3.3 (+/- 1.6; Delta = -2.6; p > 0.0001), 4.2 (+/- 1.5; Delta = -1.1; p = 0.0001) and 5.4 (+/- 1.7; Delta = -0.7; p = 0.06). For the first TNF antagonist, DAS-28-based EULAR response level was good in 42% and moderate in 33% of patients. The second TNF antagonist yielded a good response in 20% and no response in 53% of patients, while the third one yielded a good response in 28% and no response in 72%. Mean baseline HAQ on starting the first, second and third TNF antagonist was 1.61, 1.52 and 1.87, respectively. At the end of follow-up, it decreased to 1.12 (Delta = -0.49; p < 0.0001), 1.31 (Delta = -0.21, p = 0.004) and 1.75 (Delta = -0.12; p = 0.1), respectively. Sixty four percent of patients had a clinically important improvement in HAQ (defined as > or = -0.22) with the first TNF antagonist and 46% with the second. CONCLUSION A clinically significant effect size was seen in less than half of RA patients cycling to a second TNF antagonist.