Exploring the Origin of Power Law Distribution in Single-Molecule Conformation Dynamics: Energy Landscape Perspectives

We explored the origin of power law distribution observed in single-molecule conformational dynamics experiments. By establishing a kinetic master equation approach to study statistically the microscopic state dynamics, we show that the underlying landscape with exponentially distributed density of states leads to power law distribution of kinetics. The exponential density of states emerges when the system becomes glassy and landscape becomes rough with significant trapping.

Probabilistic Solution of Inverse Problems

In this thesis we study the general problem of reconstructing a function, defined on a finite lattice from a set of incomplete, noisy and/or ambiguous observations. The goal of this work is to demonstrate the generality and practical value of a probabilistic (in particular, Bayesian) approach to this problem, particularly in the context of Computer Vision. In this approach, the prior knowledge about the solution is expressed in the form of a Gibbsian probability distribution on the space of all possible functions, so that the reconstruction task is formulated as an estimation problem. Our main contributions are the following: (1) We introduce the use of specific error criteria for the design of the optimal Bayesian estimators for several classes of problems, and propose a general (Monte Carlo) procedure for approximating them. This new approach leads to a substantial improvement over the existing schemes, both regarding the quality of the results (particularly for low signal to noise ratios) and the computational efficiency. (2) We apply the Bayesian appraoch to the solution of several problems, some of which are formulated and solved in these terms for the first time. Specifically, these applications are: teh reconstruction of piecewise constant surfaces from sparse and noisy observationsl; the reconstruction of depth from stereoscopic pairs of images and the formation of perceptual clusters. (3) For each one of these applications, we develop fast, deterministic algorithms that approximate the optimal estimators, and illustrate their performance on both synthetic and real data. (4) We propose a new method, based on the analysis of the residual process, for estimating the parameters of the probabilistic models directly from the noisy observations. This scheme leads to an algorithm, which has no free parameters, for the restoration of piecewise uniform images. (5) We analyze the implementation of the algorithms that we develop in non-conventional hardware, such as massively parallel digital machines, and analog and hybrid networks.

Number of degrees of freedom in inverse diffraction

The problem of inverse diffraction from plane to plane is considered in the case where a finite aperture exists in the boundary plane. Singular values and singular functions for the problem are introduced, and the number of degrees of freedom is defined in terms of the distribution of the singular values. Numerical computations are presented for the one-dimensional problem, and it is shown that the effect of evanescent waves disappears at a distance of approximately one wavelength from the boundary plane, even when the dimension of the slit is comparable with the wavelength of the diffracted field. © 1983 Taylor & Francis Group, LLC.

Stability Problems in Inverse Diffraction

Inverse diffraction consists in determining the field distribution on a boundary surface from the knowledge of the distribution on a surface situated within the domain where the wave propagates. This problem is a good example for illustrating the use of least-squares methods (also called regularization methods) for solving linear ill-posed inverse problem. We focus on obtaining error bounds For regularized solutions and show that the stability of the restored field far from the boundary surface is quite satisfactory: the error is proportional to ∊(ðŗ‚ ≃ 1) ,ðŗœ being the error in the data (Hölder continuity). However, the error in the restored field on the boundary surface is only proportional to an inverse power of │In∊│ (logarithmic continuity). Such a poor continuity implies some limitations on the resolution which is achievable in practice. In this case, the resolution limit is seen to be about half of the wavelength. Copyright © 1981 by The Institute of Electrical and Electronics Engineers, Inc.

A new parameterization of Johnson's SB distribution with application to fitting forest tree diameter data

The SB distributional model of Johnson's 1949 paper was introduced by a transformation to normality, that is, z ~ N(0, 1), consisting of a linear scaling to the range (0, 1), a logit transformation, and an affine transformation, z = γ + δu. The model, in its original parameterization, has often been used in forest diameter distribution modelling. In this paper, we define the SB distribution in terms of the inverse transformation from normality, including an initial linear scaling transformation, u = γ′ + δ′z (δ′ = 1/δ and γ′ = �γ/δ). The SB model in terms of the new parameterization is derived, and maximum likelihood estimation schema are presented for both model parameterizations. The statistical properties of the two alternative parameterizations are compared empirically on 20 data sets of diameter distributions of Changbai larch (Larix olgensis Henry). The new parameterization is shown to be statistically better than Johnson's original parameterization for the data sets considered here.

Plankton Of The Fladen Ground During Flex-76 .3. Vertical-Distribution Population-Dynamics And Production Of Calanus-Finmarchicus (Crustacea Copepoda)

Samples taken in the northern North Sea with the Continuous Plankton Recorder (CPR), the Undulating Oceanographic Recorder (UOR), the Longhurst Hardy Plankton Recorder (LHPR) and by our colleagues from other participating Institutes during the Fladen Ground Experiment (FLEX 76) were used to describe the vertical distribution and population dynamics of Calanus finmarchicus (Gunnerus) and to provide estimates of the production and carbon budget of the population from 19 March to 3 June, 1976. Total production of the 19 March to 3 June, 1976. Total production of the nauplii and copepodite stages (including adults), during the exponential growth phase in May, was estimated to be in the range of 0.49 to 0.91 g C m-2 d-1 or 29.0 to 55 g dry wt m-2 (14.5 to 27.8 g C m-2) for the three successive 10 d periods in May. Two gross growth efficiencies (K 1) (20 and 34%), together with the lower value of C. finmarchicus production, were used to calculate the gross ingestion levels of algae as 2.45 and 1.44 g C m-2 d-1 (73.5 and 43.2 g C m-2 over the May period). These ingestion levels, together with the algae ingested by other zooplankton species, are greater than the estimated total phytoplankton production of 45.9 g C m-2 over the FLEX period. A number of factors are discussed which could explain the discrepancies between the production estimates. One suggestion is that the vertical distribution of the development stages of this herbivorous copepod and their diel and ontogenetic migration patterns enable it to efficiently exploit its food source. Data from the FLEX experiment indicated that the depletion of nutrients limited the size of the spring bloom, but that it was the grazing pressure exerted by C. finmarchicus which was responsible for the control and depletion of the phytoplankton in the spring of 1976 in the northern North Sea.

Controls on the contemporary distribution of lake thecamoebians (testate amoebae) in the Greater Toronto Area and their potential as water quality indicators.

Thecamoebians were examined from 71 surface sediment samples collected from 21 lakes and ponds in the Greater Toronto Area to (1) elucidate the controls on faunal distribution in modern lake environments; and (2) to consider the utility of thecamoebians in quantitative studies of water quality change. This area was chosen because it includes a high density of kettle and other lakes which are threatened by urban development and where water quality has deteriorated locally as a result of contaminant inputs, particularly nutrients. Fifty-eight samples yielded statistically significant thecamoebian populations. The most diverse faunas (highest Shannon Diversity Index values) were recorded in lakes beyond the limits of urban development, although the faunas of all lakes showed signs of sub-optimal conditions. The assemblages were divided into five clusters using Q-mode cluster analysis, supported by Detrended Correspondence Analysis. Canonical Correspondence Analysis (CCA) was used to examine species-environment relationships and to explain the observed clusterings. Twenty-four measured environmental variables were considered, including water property attributes (e.g., pH, conductivity, dissolved oxygen), substrate characteristics, sediment-based phosphorus (Olsen P) and 11 environmentally available metals. The thecamoebian assemblages showed a strong association with phosphorus, reflecting the eutrophic status of many of the lakes, and locally to elevated conductivity measurements, which appear to reflect road salt inputs associated with winter de-icing operations. Substrate characteristics, total organic carbon and metal contaminants (particularly Cu and Mg) also influenced the faunas of some samples. A series of partial CCAs show that of the measured variables, sedimentary phosphorus has the largest influence on assemblage distribution, explaining 6.98% (P < 0.002) of the total variance. A transfer function was developed for sedimentary phosphorus (Olsen P) using 58 samples from 15 of the studied lakes. The best performing model was based on weighted averaging with inverse deshrinking (WA Inv, r jack 2= 0.33, RMSEP = 102.65 ppm). This model was applied to a small modern thecamoebian dataset from a eutrophic lake in northern Ontario to predict phosphorus and performed satisfactorily. This preliminary study confirms that thecamoebians have considerable potential as quantitative water quality indicators in urbanising regions, particularly in areas influenced by nutrient inputs and road salts.

Distribution of Uranium and Thorium in Dolomitic Gravel Fill and Shale Saprolite

The objectives of this study were to examine (1) the distribution of U and Th in dolomitic gravel fill and shale saprolite, and (2) the removal of uranium from acidic groundwater by dolomitic gravel through precipitation with amorphous basaluminite at the U.S. DOE Oak Ridge Integrated Field Research Challenge (ORIFRC) field site west of the Oak Ridge Y-12 National Security Complex in East Tennessee. Media reactivity and sustainability are a technical concern with the deployment of any subsurface reactive media. Because the gravel was placed in the subsurface and exposed to contaminated groundwater for over 20 years, it provided a unique opportunity to study the solid and water phase geochemical conditions within the media after this length of exposure. This study illustrates that dolomite gravel can remove U from acidic contaminated groundwater with high levels of Al3+, Ca2+, NO3−, and SO42− over the long term. As the groundwater flows through high pH carbonate gravel, U containing amorphous basaluminite precipitates as the pH increases. This is due to an increase in groundwater pH from 3.2 to ∼6.5 as it comes in contact with the gravel. Therefore, carbonate gravel could be considered as a possible treatment medium for removal and sequestration of U and other pH sensitive metals from acidic contaminated groundwater. Thorium concentrations are also high in the carbonate gravel. Thorium generally shows an inverse relationship with U from the surface down into the deeper saprolite. Barite precipitated in the shallow saprolite directly below the dolomitic gravel from barium present in the acidic contaminated groundwater.

Abrasion, transport and distribution of sediment in selected streams of Southern Ontario and Western New York /

The rate of decrease in mean sediment size and weight per square metre along a 54 km reach of the Credit River was found to depend on variations in the channel geometry. The distribution of a specific sediment size consist of: (1) a transport zone; (2) an accumulation zone; and (3) a depletion zone. These zones shift downstream in response to downcurrent decreases in stream competence. Along a .285 km man-made pond, within the Credit River study area, the sediment is also characterized by downstream shifting accumulation zones for each finer clast size. The discharge required to initiate movement of 8 cm and 6 cm blocks in Cazenovia Creek is closely approximated by Baker and Ritter's equation. Incipient motion of blocks in Twenty Mile Creek is best predicted by Yalin's relation which is more efficient in deeper flows. The transport distance of blocks in both streams depends on channel roughness and geometry. Natural abrasion and distribution of clasts may depend on the size of the surrounding sediment and variations in flow competence. The cumulative percent weight loss with distance of laboratory abraded dolostone is defined by a power function. The decrease in weight of dolostone follows a negative exponential. In the abrasion mill, chipping causes the high initial weight loss of dolostone; crushing and grinding produce most of the subsequent weight loss. Clast size was found to have little effect on the abrasion of dolostone within the diameter range considered. Increasing the speed of the mill increased the initial amount of weight loss but decreased the rate of abrasion. The abrasion mill was found to produce more weight loss than stream action. The maximum percent weight loss determined from laboratory and field abrasion data is approximately 40 percent of the weight loss observed along the Credit River. Selective sorting of sediment explains the remaining percentage, not accounted for by abrasion.

Distribution-Free Bounds for Serial Correlation Coefficients in Heteroskedastic Symmetric Time Series

We consider the problem of testing whether the observations X1, ..., Xn of a time series are independent with unspecified (possibly nonidentical) distributions symmetric about a common known median. Various bounds on the distributions of serial correlation coefficients are proposed: exponential bounds, Eaton-type bounds, Chebyshev bounds and Berry-Esséen-Zolotarev bounds. The bounds are exact in finite samples, distribution-free and easy to compute. The performance of the bounds is evaluated and compared with traditional serial dependence tests in a simulation experiment. The procedures proposed are applied to U.S. data on interest rates (commercial paper rate).

Strategies for Inverse Profiling of Two Dimensional Dielectric Scatterers Using Electromagnetic Illumination

Electromagnetic tomography has been applied to problems in nondestructive evolution, ground-penetrating radar, synthetic aperture radar, target identification, electrical well logging, medical imaging etc. The problem of electromagnetic tomography involves the estimation of cross sectional distribution dielectric permittivity, conductivity etc based on measurement of the scattered fields. The inverse scattering problem of electromagnetic imaging is highly non linear and ill posed, and is liable to get trapped in local minima. The iterative solution techniques employed for computing the inverse scattering problem of electromagnetic imaging are highly computation intensive. Thus the solution to electromagnetic imaging problem is beset with convergence and computational issues. The attempt of this thesis is to develop methods suitable for improving the convergence and reduce the total computations for tomographic imaging of two dimensional dielectric cylinders illuminated by TM polarized waves, where the scattering problem is defmed using scalar equations. A multi resolution frequency hopping approach was proposed as opposed to the conventional frequency hopping approach employed to image large inhomogeneous scatterers. The strategy was tested on both synthetic and experimental data and gave results that were better localized and also accelerated the iterative procedure employed for the imaging. A Degree of Symmetry formulation was introduced to locate the scatterer in the investigation domain when the scatterer cross section was circular. The investigation domain could thus be reduced which reduced the degrees of freedom of the inverse scattering process. Thus the entire measured scattered data was available for the optimization of fewer numbers of pixels. This resulted in better and more robust reconstructions of the scatterer cross sectional profile. The Degree of Symmetry formulation could also be applied to the practical problem of limited angle tomography, as in the case of a buried pipeline, where the ill posedness is much larger. The formulation was also tested using experimental data generated from an experimental setup that was designed. The experimental results confirmed the practical applicability of the formulation.

Bayesian Inference in Exponential and Pareto populations in the presence of Outliers

This thesis Entitled Bayesian inference in Exponential and pareto populations in the presence of outliers. The main theme of the present thesis is focussed on various estimation problems using the Bayesian appraoch, falling under the general category of accommodation procedures for analysing Pareto data containing outlier. In Chapter II. the problem of estimation of parameters in the classical Pareto distribution specified by the density function. In Chapter IV. we discuss the estimation of (1.19) when the sample contain a known number of outliers under three different data generating mechanisms, viz. the exchangeable model. Chapter V the prediction of a future observation based on a random sample that contains one contaminant. Chapter VI is devoted to the study of estimation problems concerning the exponential parameters under a k-outlier model.

Optimization of inverse reflector design

Aquesta tesi presenta un nou mètode pel disseny invers de reflectors. Ens hem centrat en tres temes principals: l’ús de fonts de llum reals i complexes, la definició d’un algoritme ràpid pel càlcul de la il•luminació del reflector, i la definició d’un algoritme d’optimització per trobar més eficientment el reflector desitjat. Les fonts de llum estan representades per models near-field, que es comprimeixen amb un error molt petit, fins i tot per fonts de llum amb milions de raigs i objectes a il•luminar molt propers. Llavors proposem un mètode ràpid per obtenir la distribució de la il•luminació d’un reflector i la seva comparació amb la il•luminació desitjada, i que treballa completament en la GPU. Finalment, proposem un nou mètode d’optimització global que permet trobar la solució en menys passos que molts altres mètodes d’optimització clàssics, i alhora evitant mínims locals.

Bayesian analysis of an inverse Gaussian correlated frailty model

In survival analysis frailty is often used to model heterogeneity between individuals or correlation within clusters. Typically frailty is taken to be a continuous random effect, yielding a continuous mixture distribution for survival times. A Bayesian analysis of a correlated frailty model is discussed in the context of inverse Gaussian frailty. An MCMC approach is adopted and the deviance information criterion is used to compare models. As an illustration of the approach a bivariate data set of corneal graft survival times is analysed. (C) 2006 Elsevier B.V. All rights reserved.

Assembling spatially explicit landscape models of pollen and spore dispersal by wind for risk assessment

Models of windblown pollen or spore movement are required to predict gene flow from genetically modified (GM) crops and the spread of fungal diseases. We suggest a simple form for a function describing the distance moved by a pollen grain or fungal spore, for use in generic models of dispersal. The function has power-law behaviour over sub-continental distances. We show that air-borne dispersal of rapeseed pollen in two experiments was inconsistent with an exponential model, but was fitted by power-law models, implying a large contribution from distant fields to the catches observed. After allowance for this 'background' by applying Fourier transforms to deconvolve the mixture of distant and local sources, the data were best fit by power-laws with exponents between 1.5 and 2. We also demonstrate that for a simple model of area sources, the median dispersal distance is a function of field radius and that measurement from the source edge can be misleading. Using an inverse-square dispersal distribution deduced from the experimental data and the distribution of rapeseed fields deduced by remote sensing, we successfully predict observed rapeseed pollen density in the city centres of Derby and Leicester (UK).