Direct Computation of Stochastic Flow in Reservoirs with Uncertain Parameters

A direct method is presented for determining the uncertainty in reservoir pressure, flow, and net present value (NPV) using the time-dependent, one phase, two- or three-dimensional equations of flow through a porous medium. The uncertainty in the solution is modelled as a probability distribution function and is computed from given statistical data for input parameters such as permeability. The method generates an expansion for the mean of the pressure about a deterministic solution to the system equations using a perturbation to the mean of the input parameters. Hierarchical equations that define approximations to the mean solution at each point and to the field covariance of the pressure are developed and solved numerically. The procedure is then used to find the statistics of the flow and the risked value of the field, defined by the NPV, for a given development scenario. This method involves only one (albeit complicated) solution of the equations and contrasts with the more usual Monte-Carlo approach where many such solutions are required. The procedure is applied easily to other physical systems modelled by linear or nonlinear partial differential equations with uncertain data.

Probability density estimation with tunable kernels using orthogonal forward regression

A generalized or tunable-kernel model is proposed for probability density function estimation based on an orthogonal forward regression procedure. Each stage of the density estimation process determines a tunable kernel, namely, its center vector and diagonal covariance matrix, by minimizing a leave-one-out test criterion. The kernel mixing weights of the constructed sparse density estimate are finally updated using the multiplicative nonnegative quadratic programming algorithm to ensure the nonnegative and unity constraints, and this weight-updating process additionally has the desired ability to further reduce the model size. The proposed tunable-kernel model has advantages, in terms of model generalization capability and model sparsity, over the standard fixed-kernel model that restricts kernel centers to the training data points and employs a single common kernel variance for every kernel. On the other hand, it does not optimize all the model parameters together and thus avoids the problems of high-dimensional ill-conditioned nonlinear optimization associated with the conventional finite mixture model. Several examples are included to demonstrate the ability of the proposed novel tunable-kernel model to effectively construct a very compact density estimate accurately.

Estimators in capture–recapture studies with two sources

This paper investigates the applications of capture–recapture methods to human populations. Capture–recapture methods are commonly used in estimating the size of wildlife populations but can also be used in epidemiology and social sciences, for estimating prevalence of a particular disease or the size of the homeless population in a certain area. Here we focus on estimating the prevalence of infectious diseases. Several estimators of population size are considered: the Lincoln–Petersen estimator and its modified version, the Chapman estimator, Chao’s lower bound estimator, the Zelterman’s estimator, McKendrick’s moment estimator and the maximum likelihood estimator. In order to evaluate these estimators, they are applied to real, three-source, capture-recapture data. By conditioning on each of the sources of three source data, we have been able to compare the estimators with the true value that they are estimating. The Chapman and Chao estimators were compared in terms of their relative bias. A variance formula derived through conditioning is suggested for Chao’s estimator, and normal 95% confidence intervals are calculated for this and the Chapman estimator. We then compare the coverage of the respective confidence intervals. Furthermore, a simulation study is included to compare Chao’s and Chapman’s estimator. Results indicate that Chao’s estimator is less biased than Chapman’s estimator unless both sources are independent. Chao’s estimator has also the smaller mean squared error. Finally, the implications and limitations of the above methods are discussed, with suggestions for further development.

Estimators in capture-recapture studies with two sources

This paper investigates the applications of capture-recapture methods to human populations. Capture-recapture methods are commonly used in estimating the size of wildlife populations but can also be used in epidemiology and social sciences, for estimating prevalence of a particular disease or the size of the homeless population in a certain area. Here we focus on estimating the prevalence of infectious diseases. Several estimators of population size are considered: the Lincoln-Petersen estimator and its modified version, the Chapman estimator, Chao's lower bound estimator, the Zelterman's estimator, McKendrick's moment estimator and the maximum likelihood estimator. In order to evaluate these estimators, they are applied to real, three-source, capture-recapture data. By conditioning on each of the sources of three source data, we have been able to compare the estimators with the true value that they are estimating. The Chapman and Chao estimators were compared in terms of their relative bias. A variance formula derived through conditioning is suggested for Chao's estimator, and normal 95% confidence intervals are calculated for this and the Chapman estimator. We then compare the coverage of the respective confidence intervals. Furthermore, a simulation study is included to compare Chao's and Chapman's estimator. Results indicate that Chao's estimator is less biased than Chapman's estimator unless both sources are independent. Chao's estimator has also the smaller mean squared error. Finally, the implications and limitations of the above methods are discussed, with suggestions for further development.

New dioxadiaza- and trioxadiaza-macrocycles containing one dibenzofuran unit with two amino pendant arms: synthesis, protonation and complexation studies

New dioxadiaza- and trioxadiaza-macrocycles containing one rigid dibenzofuran unit (DBF) and N-(2-aminoethyl) pendant arms were synthesized, N,N'-bis(2-aminoethyl)-[17]( DBF) N2O2 (L-1) and N,N'-bis(2-aminoethyl)-[22](DBF)N2O3 (L-2), respectively. The binding properties of both macrocycles to metal ions and structural studies of their metal complexes were carried out. The protonation constants of both compounds and the stability constants of their complexes with Co2+, Ni2+, Cu2+, Zn2+, Cd2+, and Pb2+ were determined at 298.2 K, in aqueous solutions, and at ionic strength 0.10 mol dm(-3) in KNO3. Mononuclear complexes with both ligands were formed, and dinuclear complexes were only found for L-2. The thermodynamic binding affinities of the metal complexes of L-2 are lower than those of L-1 as expected, but the Pb2+ complexes of both macrocycles exhibit close stability constant values. On the other hand, the binding affinities of Cd2+ and Pb2+ for L-1 are very high, when compared to those of Co2+, Ni2+ and Zn2+. These interesting properties were explained by the presence of the rigid DBF moiety in the backbone of the macrocycle and to the special match between the macrocyclic cavity size and the studied larger metal ions. To elucidate the adopted structures of complexes in solution, the nickel(II) and copper( II) complexes with both ligands were further studied by UV-vis-MR spectroscopy in DMSO-H2O 1 : 1 (v/v) solution. The copper(II) complexes were also studied by EPR spectroscopy in the same mixture of solvents. The crystal structure of the copper complex of L-1 was also determined. The copper(II) displays an octahedral geometry, the four nitrogen atoms forming the equatorial plane and two oxygen atoms, one from the DBF unit and the other one from the ether oxygen, in axial positions. One of the ether oxygens of the macrocycle is out of the coordination sphere. Our results led us to suggest that this geometry is also adopted by the Co2+ to Zn2+ complexes, and only the larger Cd2+ and Pb2+ manage to form complexes with the involvement of all the oxygen atoms of the macrocyclic backbone.

Dimetallic complexes of macrocycles with two rigid dibenzofuran units as receptors for detection of anionic substrates

The hexaazamacrocycles [28](DBF)2N6 {cyclo[bis(4,6-dimethyldibenzo[b,d]furaniminoethyleneiminoethylene]} and [32](DBF)2N6 {cyclo[bis(4,6-dimethyldibenzo[b,d]furaniminopropyleneiminopropylene]} form stable dinuclear copper(II) complexes suitable to behave as receptors for several anionic substrates. These two receptors were used to study the binding interactions with several substrates, such as imidazole (Him) and some carboxylates [benzoate (bz−), oxalate (ox2−), malonate (mal2−), phthalate (ph2−), isophthalate (iph2−), and terephthalate (tph2−)] by spectrophotometric titrations and EPR spectroscopy in MeOH (or H2O):DMSO (1:1 v/v) solution. The largest association constant was found for ox2− with Cu2[32](DBF)2N64+, whereas for the aromatic dicarboxylate anions the binding constants follow the trend ph2− > iph2− > tph2−, i.e. decrease with the increase of the distance of the two binding sites of the substrate. On the other hand, the large blue shift of 68 nm observed by addition of Him to Cu2[32](DBF)2N64+ points out for the formation of the bridged CuimCu cascade complex, indicating this receptor as a potential sensor for the detection and determination of imidazole in solution. The X-band EPR spectra of the Cu2[28](DBF)2N64+ and Cu2[32](DBF)2N6]4+ complexes and the cascade complexes with the substrates, performed in H2O:DMSO (1:1 v/v) at 5 to 15 K, showed that the CuCu distance is slightly larger than the one found in crystal state and that this distance increases when the substrate is accommodated between the two copper centres. The crystal structure of [Cu2[28](DBF)2N6(ph)2]·CH3OH was determined by X-ray diffraction and revealed the two copper centres bridged by two ph2− anions at a Cu···Cu distance of 5.419(1) Å. Each copper centre is surrounded by three carboxylate oxygen atoms from two phthalate anions and three contiguous nitrogen atoms of the macrocycle in a pseudo octahedral coordination environment.

Emotion processing in infancy: specificity in risk for social anxiety and associations with two year outcomes

The current study examined the specificity of patterns of responding to high and low intensity negative emotional expressions of infants of mothers with social phobia, and their association with child outcomes at two years of age. Infants of mothers with social phobia, generalised anxiety disorder (GAD) or no history of anxiety were shown pairs of angry and fearful emotional expressions at 10 weeks of age. Symptoms of social withdrawal, anxiety and sleep problems were assessed at two years of age. Only infants of mothers with social phobia showed a tendency to look away from high intensity fear faces; however infants of mothers with both social phobia and GAD showed a bias towards high intensity angry faces. Among the offspring of mothers with social phobia, anxiety symptoms at two years of age were associated with a preference for high intensity fear faces in infancy. The reverse pattern was found amongst the offspring of non-anxious mothers. These findings suggest a possible specific response to emotional expressions among the children of mothers with social phobia.

Knowledge of A/A'-dependencies on subject extraction with two types of infinitives in non-native Portuguese adult bilingualism

Within generative L2 acquisition research there is a longstanding debate as to what underlies observable differences in L1/L2 knowledge/ performance. On the one hand, Full Accessibility approaches maintain that target L2 syntactic representations (new functional categories and features) are acquirable (e.g., Schwartz & Sprouse, 1996). Conversely, Partial Accessibility approaches claim that L2 variability and/or optionality, even at advanced levels, obtains as a result of inevitable deficits in L2 narrow syntax and is conditioned upon a maturational failure in adulthood to acquire (some) new functional features (e.g., Beck, 1998; Hawkins & Chan, 1997; Hawkins & Hattori, 2006; Tsimpli & Dimitrakopoulou, 2007). The present study tests the predictions of these two sets of approaches with advanced English learners of L2 Brazilian Portuguese (n = 21) in the domain of inflected infinitives. These advanced L2 learners reliably differentiate syntactically between finite verbs, uninflected and inflected infinitives, which, as argued, only supports Full Accessibility approaches. Moreover, we will discuss how testing the domain of inflected infinitives is especially interesting in light of recent proposals that Brazilian Portuguese colloquial dialects no longer actively instantiate them (Lightfoot, 1991; Pires, 2002, 2006; Pires & Rothman, 2009; Rothman, 2007).

[C7H10N][In3Se5]: a layered selenide with two indium coordination environments

A new layered indium selenide, [C7H10N][In3Se5], was prepared under solvothermal conditions using 3,5-dimethylpyridine as a structure-directing agent. The crystal structure contains anionic layers of stoichiometry [In3Se5]− in which indium atoms with octahedral and tetrahedral coordination coexist. This material represents the first occurrence of octahedrally coordinated indium in a solvothermally-prepared indium selenide.

Recovering the counterfactual wage distribution with selective return migration

This paper recovers the distribution of wages for Mexican-born workers living in the U.S. if no return migration of Mexican-born workers occurred. Because migrants self-select in the decision to return, the overarching problem addressed by this study is the use of an estimator that also accounts for selection on unobservables. I find that Mexican returnees are middle- to high-wage earners at all levels of educational attainment. Taking into account self-selection in return migration, wages would be approximately 7.7% higher at the median and 4.5% higher at the mean. Owing to positive self-selection, the immigrant-native wage gap would, therefore, partially close if there was no return migration.

Sparse density estimator with tunable kernels

A new sparse kernel density estimator with tunable kernels is introduced within a forward constrained regression framework whereby the nonnegative and summing-to-unity constraints of the mixing weights can easily be satisfied. Based on the minimum integrated square error criterion, a recursive algorithm is developed to select significant kernels one at time, and the kernel width of the selected kernel is then tuned using the gradient descent algorithm. Numerical examples are employed to demonstrate that the proposed approach is effective in constructing very sparse kernel density estimators with competitive accuracy to existing kernel density estimators.

Unimodality of wave amplitude in the northern Hemisphere

A novel statistic for local wave amplitude of the 500-hPa geopotential height field is introduced. The statistic uses a Hilbert transform to define a longitudinal wave envelope and dynamical latitude weighting to define the latitudes of interest. Here it is used to detect the existence, or otherwise, of multimodality in its distribution function. The empirical distribution function for the 1960-2000 period is close to a Weibull distribution with shape parameters between 2 and 3. There is substantial interdecadal variability but no apparent local multimodality or bimodality. The zonally averaged wave amplitude, akin to the more usual wave amplitude index, is close to being normally distributed. This is consistent with the central limit theorem, which applies to the construction of the wave amplitude index. For the period 1960-70 it is found that there is apparent bimodality in this index. However, the different amplitudes are realized at different longitudes, so there is no bimodality at any single longitude. As a corollary, it is found that many commonly used statistics to detect multimodality in atmospheric fields potentially satisfy the assumptions underlying the central limit theorem and therefore can only show approximately normal distributions. The author concludes that these techniques may therefore be suboptimal to detect any multimodality.

Rapid radiation in spiny lobsters (Palinurus spp) as revealed by classic and ABC methods using mtDNA and microsatellite data

Background: Molecular tools may help to uncover closely related and still diverging species from a wide variety of taxa and provide insight into the mechanisms, pace and geography of marine speciation. There is a certain controversy on the phylogeography and speciation modes of species-groups with an Eastern Atlantic-Western Indian Ocean distribution, with previous studies suggesting that older events (Miocene) and/or more recent (Pleistocene) oceanographic processes could have influenced the phylogeny of marine taxa. The spiny lobster genus Palinurus allows for testing among speciation hypotheses, since it has a particular distribution with two groups of three species each in the Northeastern Atlantic (P. elephas, P. mauritanicus and P. charlestoni) and Southeastern Atlantic and Southwestern Indian Oceans (P. gilchristi, P. delagoae and P. barbarae). In the present study, we obtain a more complete understanding of the phylogenetic relationships among these species through a combined dataset with both nuclear and mitochondrial markers, by testing alternative hypotheses on both the mutation rate and tree topology under the recently developed approximate Bayesian computation (ABC) methods. Results: Our analyses support a North-to-South speciation pattern in Palinurus with all the South-African species forming a monophyletic clade nested within the Northern Hemisphere species. Coalescent-based ABC methods allowed us to reject the previously proposed hypothesis of a Middle Miocene speciation event related with the closure of the Tethyan Seaway. Instead, divergence times obtained for Palinurus species using the combined mtDNA-microsatellite dataset and standard mutation rates for mtDNA agree with known glaciation-related processes occurring during the last 2 my. Conclusion: The Palinurus speciation pattern is a typical example of a series of rapid speciation events occurring within a group, with very short branches separating different species. Our results support the hypothesis that recent climate change-related oceanographic processes have influenced the phylogeny of marine taxa, with most Palinurus species originating during the last two million years. The present study highlights the value of new coalescent-based statistical methods such as ABC for testing different speciation hypotheses using molecular data.