Statistical modeling of bowing control applied to violin sound synthesis

Excitation-continuous music instrument control patterns are often not explicitly represented in current sound synthesis techniques when applied to automatic performance. Both physical model-based and sample-based synthesis paradigmswould benefit from a flexible and accurate instrument control model, enabling the improvement of naturalness and realism. Wepresent a framework for modeling bowing control parameters inviolin performance. Nearly non-intrusive sensing techniques allow for accurate acquisition of relevant timbre-related bowing control parameter signals.We model the temporal contour of bow velocity, bow pressing force, and bow-bridge distance as sequences of short Bézier cubic curve segments. Considering different articulations, dynamics, and performance contexts, a number of note classes are defined. Contours of bowing parameters in a performance database are analyzed at note-level by following a predefined grammar that dictates characteristics of curve segment sequences for each of the classes in consideration. As a result, contour analysis of bowing parameters of each note yields an optimal representation vector that is sufficient for reconstructing original contours with significant fidelity. From the resulting representation vectors, we construct a statistical model based on Gaussian mixtures suitable for both the analysis and synthesis of bowing parameter contours. By using the estimated models, synthetic contours can be generated through a bow planning algorithm able to reproduce possible constraints caused by the finite length of the bow. Rendered contours are successfully used in two preliminary synthesis frameworks: digital waveguide-based bowed stringphysical modeling and sample-based spectral-domain synthesis.

Spatial statistical analysis and selection of genotypes in plant breeding

The objective of this study was to evaluate the efficiency of spatial statistical analysis in the selection of genotypes in a plant breeding program and, particularly, to demonstrate the benefits of the approach when experimental observations are not spatially independent. The basic material of this study was a yield trial of soybean lines, with five check varieties (of fixed effect) and 110 test lines (of random effects), in an augmented block design. The spatial analysis used a random field linear model (RFML), with a covariance function estimated from the residuals of the analysis considering independent errors. Results showed a residual autocorrelation of significant magnitude and extension (range), which allowed a better discrimination among genotypes (increase of the power of statistical tests, reduction in the standard errors of estimates and predictors, and a greater amplitude of predictor values) when the spatial analysis was applied. Furthermore, the spatial analysis led to a different ranking of the genetic materials, in comparison with the non-spatial analysis, and a selection less influenced by local variation effects was obtained.

Multiway Array Decomposition Analysis of EEGs in Alzheimer’s Disease

Methods for the extraction of features from physiological datasets are growing needs as clinical investigations of Alzheimer’s disease (AD) in large and heterogeneous population increase. General tools allowing diagnostic regardless of recording sites, such as different hospitals, are essential and if combined to inexpensive non-invasive methods could critically improve mass screening of subjects with AD. In this study, we applied three state of the art multiway array decomposition (MAD) methods to extract features from electroencephalograms (EEGs) of AD patients obtained from multiple sites. In comparison to MAD, spectral-spatial average filter (SSFs) of control and AD subjects were used as well as a common blind source separation method, algorithm for multiple unknown signal extraction (AMUSE). We trained a feed-forward multilayer perceptron (MLP) to validate and optimize AD classification from two independent databases. Using a third EEG dataset, we demonstrated that features extracted from MAD outperformed features obtained from SSFs AMUSE in terms of root mean squared error (RMSE) and reaching up to 100% of accuracy in test condition. We propose that MAD maybe a useful tool to extract features for AD diagnosis offering great generalization across multi-site databases and opening doors to the discovery of new characterization of the disease.

Concrete Pavement Mixture Design and Analysis (MDA): An Innovative Approach To Proportioning Concrete Mixtures, TPF-5(205), 2015

Mixture proportioning is routinely a matter of using a recipe based on a previously produced concrete, rather than adjusting the proportions based on the needs of the mixture and the locally available materials. As budgets grow tighter and increasing attention is being paid to sustainability metrics, greater attention is beginning to be focused on making mixtures that are more efficient in their usage of materials yet do not compromise engineering performance. Therefore, a performance-based mixture proportioning method is needed to provide the desired concrete properties for a given project specification. The proposed method should be user friendly, easy to apply in practice, and flexible in terms of allowing a wide range of material selection. The objective of this study is to further develop an innovative performance-based mixture proportioning method by analyzing the relationships between the selected mix characteristics and their corresponding effects on tested properties. The proposed method will provide step-by-step instructions to guide the selection of required aggregate and paste systems based on the performance requirements. Although the provided guidance in this report is primarily for concrete pavements, the same approach can be applied to other concrete applications as well.

Low latitude Paleocene to early Eocene radiolarian systematics, biochronology and paleoenvironmental analysis

Highly diverse radiolarian faunas of latest Maastrichtian to early Eocene age have been recovered from the low latitude realm in order to contribute to the clarification of radiolarian taxonomy, construct a zonation based on a discrete sequence of co-existence intervals of species ranging from the late Paleocene to early Eocene and to describe a rich low latitude latest Cretaceous to late Paleocene fauna. 225 samples of late Paleocene to early Eocene age have been collected from ODP Leg 171 B-Hole 1051 A (Blake Nose), DSDP Leg 43-Site 384 (Northwest Atlantic) and DSDP Leg 10-Sites 86, 94, 95, 96. Sequences consist of mainly pelagic oozes and chalks, with some clay and ash layers. A new imaging technique is devised to perform (in particular on topotypic material) both transmitted light microscopy and SEM imaging on individual radiolarian specimens. SEM precedes transmitted light imaging. Radiolarians are adhered to a cover slip (using nail varnish) which is secured to a stub using conductive levers. Specimens are then photographed in low vacuum (40-50Pa; 0.5mbar), which enables charge neutralization by ionized molecules of the chamber atmosphere. Thus gold coating is avoided and subsequently this allows transmitted light imaging to follow. The conductive levers are unscrewed and the cover slip is simply overturned and mounted with Canada balsam. In an attempt towards a post-Haeckelian classification, the initial spicule (Entactinaria), micro- or macrosphere (Spumellaria) and initial spicule and cephalis (Nassellaria) have been studied by slicing Entactinaria and Spumellaria, and by tilting Nassellaria in the SEM chamber. A new genus of the family Coccodiscidae is erected and Spongatractus HAECKEL is re-located to the subfamily Axopruinae. The biochronology has been carried out using the Unitary Association Method (Guex 1977, 1991). A database recording the occurrences of 112 species has been used to establish a succession of 22 Unitary Associations. Each association is correlated to chronostratigraphy via calcareous microfossils that were previously studied by other authors. The 22 UAs have been united into seven Unitary Associations Zones (UAZones) (JP10- JE4). The established zones permit to distinguish supplementary subdivisions within the existing zonation. The low-latitude Paleocene radiolarian zonation established by Sanfilippo and Nigrini (1998a) is incomplete due to the lack of radiolarian-bearing early Paleocene sediments. In order to contribute to the study of sparsely known low latitude early Paleocene faunas, 80 samples were taken from the highly siliceous Guayaquil Formation (Ecuador). The sequence consists of black cherts, shales, siliceous limestones and volcanic ash layers. The carbonate content increases up section. Age control is supplied by sporadic occurrences of silicified planktonic foraminifera casts. One Cretaceous zone and seven Paleocene zones have been identified. The existing zonation for the South Pacific can be applied to the early-early late Paleocene sequence, although certain marker species have significantly shorter ranges (notably Buryella foremanae and B. granulata). Despite missing marker species in the late Paleocene, faunal distribution correlates reasonably to the Low-Latitude zonation. An assemblage highly abundant in Lithomelissa, Lophophaena and Cycladophora in the upper RP6 zone (correlated by the presence of Pterocodon poculum, Circodiscus circularis, Pterocodon? sp. aff. P. tenellus and Stylotrochus nitidus) shows a close affinity to contemporaneous faunas reported from Site 1121, Campbell Plateau. Coupled with a high diatom abundance (notably Aulacodiscus spp. and Arachnoidiscus spp.), these faunas are interpreted as reflecting a period of enhanced biosiliceous productivity during the late Paleocene. The youngest sample is void of radiolarians, diatoms and sponge spicules yet contains many pyritized infaunal benthic foraminifera which are akin to the midway-type fauna. The presence of this fauna suggests deposition in a neritic environment. This is in contrast to the inferred bathyal slope depositional environment of the older Paleocene sediments and suggests a shoaling of the depositional environment which may be related to a coeval major accretionary event. RESUME DE LA THESE Des faunes de radiolaires de basses latitudes très diversifiées d'âge Maastrichtien terminal à Eocène inférieur, ont été étudiées afin de contribuer à la clarification de leur taxonomie, de construire une biozonation basée sur une séquence discrète d'intervalles de coexistence des espèces d'age Paléocène supérieur à Eocène inférieur et de décrire une riche faune de basse latitude allant du Crétacé terminal au Paléocène supérieur. L'étude de cette faune contribue particulièrement à la connaissance des insaisissables radiolaires de basses latitudes du Paléocène inférieur. 225 échantillons d'âge Paléocène supérieur à Eocène inférieur provenant des ODP Leg 171B-Site 1051A (Blake Nose), Leg DSDP 43-Site 384 (Atlantique Nord -Ouest) et des DSDP Leg 10 -Sites 86, 94, 95, 96, ont été étudiés. Ces séquences sont constituées principalement de « ooze » et de « chalks »pélagiques ainsi que de quelques niveaux de cendres et d'argiles. Une nouvelle technique d'imagerie a été conçue afin de pouvoir prendre conjointement des images en lumière transmise et au Microscope Electronique à Balayage (MEB) de spécimens individuels. Ceci à été particulièrement appliqué à l'étude des topotypes. L'imagerie MEB précède l'imagerie en lumière transmise. Les radiolaires sont collés sur une lame pour micropaléontologie (au moyen de vernis à ongles) qui est ensuite fixée à un porte-objet à l'aide de bras métalliques conducteurs. Les spécimens sont ensuite photographiés en vide partiel (40-50Pa; 0.5mbar), ce qui permet la neutralisation des charges électrostatiques dues à la présence de molécules ionisées dans l'atmosphère de la chambre d'observation. Ainsi la métallisation de l'échantillon avec de l'or n'est plus nécessaire et ceci permet l'observation ultérieure en lumière transmise. Les bras conducteurs sont ensuite dévissés et la lame est simplement retournée et immergée dans du baume du Canada. Dans une approche de classification post Haeckelienne, le spicule initial (Entactinaires), la micro- ou macro -sphère (Spumellaires) et le spicule initial et cephalis (Nassellaires) ont été étudiés. Ceci a nécessité le sectionnement d'Entactinaires et de Spumellaires, et de pivoter les Nassellaires dans la chambre d'observation du MEB. Un nouveau genre de la Famille des Coccodiscidae a été érigé et Spongatractus HAECKEL à été réassigné à la sous-famille des Axopruninae. L'analyse biostratigraphique à été effectuée à l'aide de la méthode des Associations Unitaires {Guex 1977, 1991). Une base de données enregistrant les présences de 112 espèces à été utilisée poux établir une succession de 22 Associations Unitaires. Chaque association est corrélée à la chronostratigraphie au moyen de microfossiles calcaires précédemment étudiés par d'autres auteurs. Les 22 UAs ont été combinées en sept Zones d'Associations Unitaires (UAZones) (JP10- JE4). Ces Zones permettent d'insérer des subdivisions supplémentaires dans la zonation actuelle. La zonation de basses latitudes du Paléocène établie par Sanfilippo et Nigrini (1998a) est incomplète due au manque de sédiments du Paléocène inférieur contenant des radiolaires. Afin de contribuer à l'étude des faunes peu connues des basses latitudes du Paléocène inférieur, 80 échantillons ont été prélevés d'une section siliceuse de la Formation de Guayaquil (Equateur). La séquence est composée de cherts noirs, de shales, de calcaires siliceux et de couches de cendres volcaniques. La fraction carbonatée augmente vers le haut de la section. Des contraintes chronologiques sont fournies par la présence sporadique de moules de foraminifères planctoniques. Une zone d'intervalles du Crétacé et sept du Paléocène ont été mises en évidence. Bien que certaines espèces marqueur ont des distributions remarquablement plus courtes (notamment Buryella foremanae et B. granulata), la zonation existante pour le Pacifique Sud est applicable à la séquence d'age Paléocène inférieure à Paléocène supérieur basal étudiée. Malgré l'absence d'espèces marqueur du Paléocène supérieur, la succession faunistique se corrèle raisonnablement avec la zonation pour les basses latitudes. Un assemblage contenant d'abondants représentant du genre Lithomelissa, Lophophaena et Cycladophora dans la zone RP6 (correlée par la présence de Pterocodon poculum, Circodiscus circularis, Pterocodon? sp. aff. P. tenellus et Stylotrochus nitidus) montre une grande similitude avec certaines faunes issues des hauts latitudes et d'age semblable décrites par Hollis (2002, Site 1121, Campbell Plateau). Ceci, en plus d'une abondance importante en diatomés (notamment Aulacodiscus spp. et Arachnoidiscus spp.) nous mènent à interpréter cette faune comme témoin d'un épisode de productivité biosiliceuse accrue dans le Paléocène supérieur. L'échantillon le plus jeune, dépourvu de radiolaires, de diatomés et de spicules d'éponge contient de nombreux foraminifères benthiques infaunaux pyritisés. Les espèces identifiées sont caractéristiques d'une faune de type midway. La présence de ces foraminifères suggère un environnement de type néritique. Ceci est en contraste avec l'environnement de pente bathyale caractérisent les sédiments sous-jacent. Cette séquence de diminution de la tranche d'eau peut être associée à un événement d'accrétion majeure. RESUME DE LA THESE (POUR LE GRAND PUBLIC) Les radiolaires constituent le groupe de plancton marin le plus divers et le plus largement répandu de l'enregistrement fossile. Un taux d'évolution rapide et une variation géographique considérable des populations font des radiolaires un outil de recherche sans égal pour la biostratigraphie et la paléocéanographie. Néanmoins, avant de pouvoir les utiliser comme outils de travail, il est essentiel d'établir une solide base taxonomique. L'étude des Radiolaires peut impliquer plusieurs techniques d'extraction, d'observation et d'imagerie qui sont dépendantes du degré d'altération diagénétique des spécimens. Le squelette initial, qu'il s'agisse d'un spicule initial (Entactinaria), d'une micro- ou macro -sphère (Spumellaria) ou d'un spicule initial et d'un cephalis (Nassellaria), est l'élément le plus constant au cours de l'évolution et devrait représenter le fondement de la systématique. Des échantillons provenant de carottes de basses latitudes du Deep Sea Drilling Project et de l' Ocean Drilling ont été étudiés. De nouvelles techniques d'imagerie et de sectionnement ont été développées sur des topotypes de radiolaires préservés en opale, dans le but d'étudier les caractéristiques de leur squelette initial qui n'étaient pas visibles dans leur illustration originale. Ceci aide entre autre à comparer des spécimens recristallisés en quartz, provenant de terrains accrétés, avec les holotypes en opale de la littérature. La distribution des espèces étudiés a fourni des données biostratigraphiques qui ont été compilées à l'aide de la méthode des Associations Unitaires (Guez 1977, 1991). Il s'agit d'un modèle mathématique déterministe conçu pour exploiter la totalité de l'assemblage plutôt que de se confiner à l'utilisation de taxons marqueurs individuels. Une séquence de 22 Associations Unitaires a été établie pour la période allant du Paléocène supérieur à l'Éocène inférieur. Chaque Association Unitaire a été corrélée à l'échelle de temps absolue à l'aide de microfossiles calcaires. Les 22 UAs ont été combinées en sept Zones d'Associations Unitaires (JP10- JE4). Ces Zones permettent d'insérer des subdivisions supplémentaires dans la zonation actuelle. Les radiolaires du Paléocène inférieur à moyen des basses latitudes sont rares. Les meilleures sections connues se trouvent dans les hautes latitudes (Nouvelle Zélande). Quelques assemblages épars ont été mentionnés par le passé en Californie, en Équateur et en Russie. Une séquence siliceuse de 190 mètres dans la Formation de Guayaquil (Équateur), s'étendant du Maastrichtien supérieur au Paléocène supérieur, a fourni des faunes relativement bien préservées. L'étude de ces faunes a permis de mettre en évidence la première séquence complète de radiolaires de basses latitudes dans le Paléocène inférieure. Huit zones allant du Crétacé terminal au Paléocène supérieur ont pu être appliqués et la présence de foraminifères planctoniques a fournie plusieurs points d'attache chronologiques. Dans le Paléocène supérieur, un riche assemblage contenant d'abondants diatomés et radiolaires ayant des similitudes faunistiques marquantes avec des assemblages de hautes latitudes de Nouvelle Zélande, témoigne d'un épisode de productivité biosiliceuse accrue pendant cette période. Étant donné que la pointe du continent sud-américain et l'Antarctique étaient plus proches au cours du Paléocène, ce phénomène peut être expliqué par le transport, le long de la côte ouest de l'Amérique du Sud, d'eaux riches en nutriments en provenance de l'Océan Antarctique. Suite à cet épisode, l'enregistrement en radiolaires est interrompu. Ceci peut être associé à des événements tectoniques régionaux qui ont eu pour effet de diminuer la tranche d'eau relative, rendant l'environnement plus favorable aux foraminifères benthiques qui sont abondamment présents dans l'échantillon le plus jeune de la séquence.

Integrative analysis of the Regulatory Region of the FGFR3 Oncogene

The study of transcriptional regulation often needs the integration of diverse yet independent data. In the present work, sequence conservation, predic-tion of transcription factor binding sites (TFBS) and gene expression analysis have been applied to the detection of putative transcription factor (TF) modules in the regulatory region of the FGFR3 oncogene. Several TFs with conserved binding sites in the FGFR3 regulatory region have shown high positive or negative corre-lation with FGFR3 expression both in urothelial carcinoma and in benign nevi. By means of conserved TF cluster analysis, two different TF modules have been iden-tified in the promoter and first intron of FGFR3 gene. These modules contain acti-vating AP2, E2F, E47 and SP1 binding sites plus motifs for EGR with possible repressor function.

EEG signal analysis via a cleaning procedure based on multivariate empirical mode decomposition

Artifacts are present in most of the electroencephalography (EEG) recordings, making it difficult to interpret or analyze the data. In this paper a cleaning procedure based on a multivariate extension of empirical mode decomposition is used to improve the quality of the data. This is achieved by applying the cleaning method to raw EEG data. Then, a synchrony measure is applied on the raw and the clean data in order to compare the improvement of the classification rate. Two classifiers are used, linear discriminant analysis and neural networks. For both cases, the classification rate is improved about 20%.

RNA interference: a new and powerful tool for functional genomic analysis.

In many species, the introduction of double-stranded RNA induces potent and specific gene silencing, referred to as RNA interference. This phenomenon, which is based on targeted degradation of mRNAs and occurs in almost any eukaryote, from trypanosomes to mice including plants and fungi, has sparked general interest from both applied and fundamental standpoints. RNA interference, which is currently used to investigate gene function in a variety of systems, is linked to natural resistance to viruses and transposon silencing, as if it were a primitive immune system involved in genome surveillance. Here, we review the mechanism of RNA interference in post-transcriptional gene silencing, its function in nature, its value for functional genomic analysis, and the modifications and improvements that may make it more efficient and inheritable. We also discuss the future directions of this versatile technique in both fundamental and applied science.

High-resolution mass spectrometry applied to the study of metabolome modifications in various chicken tissues after amoxicillin administration

The performance of high resolution accurate mass spectrometry (HRMS) operating in full scan MS mode was investigated for the quantitative determination of amoxicillin (AMX) as well as qualitative analysis of metabolomic profiles in tissues of medicated chickens. The metabolomic approach was exploited to compile analytical information on changes in the metabolome of muscle, kidney and liver from chickens subjected to a pharmacological program with AMX. Data consisting of m/z features taken throughout the entire chromatogram were extracted and filtered to be treated by Principal Component Analysis. As a result, it was found that medicated and non-treated animals were clearly clustered in distinct groups. Besides, the multivariate analysis revealed some relevant mass features contributing to this separation. In this context, recognizing those potential markers of each chicken class was a priority research for both metabolite identification and, obviously, evaluation of food quality and health effects associated to food consumption.

Network characteristics of film tourism websites in Catalonia: an application of network analysis in tourism

This master thesis presents a research on the analysis of film tourism stakeholders in Catalonia applying the network analysis approach. The research aims to provide an analysis of the relations between local tourism stakeholders with local film offices through their websites. Therefore, the development of the present work involved the review of literature on the themes of film tourism and network analysis. Then the main stakeholders of film and tourism of Catalonia were identified and their websites analyzed. The measures indicators for network analysis such as centrality, closeness and betweenness degree have been applied on the analysis of the websites to determine the extent of the relations of film and tourism stakeholders in Catalonia. Results and conclusions are presented on the referred sections

Fourier transform convolution integrals applied to generalized Born molecular volume.

Generalized Born methods are currently among the solvation models most commonly used for biological applications. We reformulate the generalized Born molecular volume method initially described by (Lee et al, 2003, J Phys Chem, 116, 10606; Lee et al, 2003, J Comp Chem, 24, 1348) using fast Fourier transform convolution integrals. Changes in the initial method are discussed and analyzed. Finally, the method is extensively checked with snapshots from common molecular modeling applications: binding free energy computations and docking. Biologically relevant test systems are chosen, including 855-36091 atoms. It is clearly demonstrated that, precision-wise, the proposed method performs as good as the original, and could better benefit from hardware accelerated boards.

Using Digital Video Analysis to Monitor Driver Behavior at Intersections

Commercially available instruments for road-side data collection take highly limited measurements, require extensive manual input, or are too expensive for widespread use. However, inexpensive computer vision techniques for digital video analysis can be applied to automate the monitoring of driver, vehicle, and pedestrian behaviors. These techniques can measure safety-related variables that cannot be easily measured using existing sensors. The use of these techniques will lead to an improved understanding of the decisions made by drivers at intersections. These automated techniques allow the collection of large amounts of safety-related data in a relatively short amount of time. There is a need to develop an easily deployable system to utilize these new techniques. This project implemented and tested a digital video analysis system for use at intersections. A prototype video recording system was developed for field deployment. A computer interface was implemented and served to simplify and automate the data analysis and the data review process. Driver behavior was measured at urban and rural non-signalized intersections. Recorded digital video was analyzed and used to test the system.

Rural Expressway Intersection Synthesis of Practice and Crash Analysis, 2004

This report presents a national synthesis of rural expressway, two-way stop -controlled (TWSC) intersection safety strategies and intersection designs and an analysis of Iowa expressway TWSC intersection crash characteristics. A rural expressway is a multi-lane highway with a divided median and with mostly at -grade intersections, although some intersections may be grade separated. The synthesis of intersection strategies is conducted in two parts. The first is a literature review and the second part is a national survey of strategies currently being applied by state transportation agencies. The characterization of crash patterns at TWSC expressway intersections is examined through the analysis of 5 years of crash data at 644 intersections.

Finite mixture analysis of beauty-contest data using generalised beta distributions

This paper introduces a mixture model based on the beta distribution, without preestablishedmeans and variances, to analyze a large set of Beauty-Contest data obtainedfrom diverse groups of experiments (Bosch-Domenech et al. 2002). This model gives a bettert of the experimental data, and more precision to the hypothesis that a large proportionof individuals follow a common pattern of reasoning, described as iterated best reply (degenerate),than mixture models based on the normal distribution. The analysis shows thatthe means of the distributions across the groups of experiments are pretty stable, while theproportions of choices at dierent levels of reasoning vary across groups.

Analysis of cannabinoids in oral fluid by liquid chromatography-tandem mass spectrometry

A sensitive method was developed for quantifying a wide range of cannabinoids in oral fluid (OF) by liquid chromatography-tandem mass spectrometry (LC-MS/MS). These cannabinoids include a dagger(9)-tetrahydrocannabinol (THC), 11-hydroxy-a dagger(9)-tetrahydrocannabinol (11-OH-THC), 11-nor-9-carboxy-a dagger(9)-tetrahydrocannabinol (THCCOOH), cannabinol (CBN), cannabidiol (CBD), a dagger(9)-tetrahydrocannabinolic acid A (THC-A), 11-nor-9-carboxy-a dagger(9)-tetrahydrocannabinol glucuronide (THCCOOH-gluc), and a dagger(9)-tetrahydrocannabinol glucuronide (THC-gluc). Samples were collected using a Quantisal (TM) device. The advantages of performing a liquid-liquid extraction (LLE) of KCl-saturated OF using heptane/ethyl acetate versus a solid-phase extraction (SPE) using HLB copolymer columns were determined. Chromatographic separation was achieved in 11.5 min on a Kinetex (TM) column packed with 2.6-mu m core-shell particles. Both positive (THC, 11-OH-THC, CBN, and CBD) and negative (THCCOOH, THC-gluc, THCCOOH-gluc, and THC-A) electrospray ionization modes were employed with multiple reaction monitoring using a high-end AB Sciex API 5000 (TM) triple quadrupole LC-MS/MS system. Unlike SPE, LLE failed to extract THC-gluc and THCCOOH-gluc. However, the LLE method was more sensitive for the detection of THCCOOH than the SPE method, wherein the limit of detection (LOD) and limit of quantification (LOQ) decreased from 100 to 50 pg/ml and from 500 to 80 pg/ml, respectively. The two extraction methods were successfully applied to OF samples collected from volunteers before and after they smoked a homemade cannabis joint. High levels of THC were measured soon after smoking, in addition to significant amounts of THC-A. Other cannabinoids were found in low concentrations. Glucuronide conjugate levels were lower than the method's LOD for most samples. Incubation studies suggest that glucuronides could be enzymatically degraded by glucuronidase prior to OF collection