Markov random fields and spatial smoothing in improving of forest inventory estimates

Markov random fields (MRF) are popular in image processing applications to describe spatial dependencies between image units. Here, we take a look at the theory and the models of MRFs with an application to improve forest inventory estimates. Typically, autocorrelation between study units is a nuisance in statistical inference, but we take an advantage of the dependencies to smooth noisy measurements by borrowing information from the neighbouring units. We build a stochastic spatial model, which we estimate with a Markov chain Monte Carlo simulation method. The smooth values are validated against another data set increasing our confidence that the estimates are more accurate than the originals.

Study of the random vibration of nonlinear systems by the Gaussian closure technique

A technique is developed to study random vibration of nonlinear systems. The method is based on the assumption that the joint probability density function of the response variables and input variables is Gaussian. It is shown that this method is more general than the statistical linearization technique in that it can handle non-Gaussian excitations and amplitude-limited responses. As an example a bilinear hysteretic system under white noise excitation is analyzed. The prediction of various response statistics by this technique is in good agreement with other available results.

The non-linear response of a gyrostabilized platform: Effects of correction torque time delay and of random inputs

First, the non-linear response of a gyrostabilized platform to a small constant input torque is analyzed in respect to the effect of the time delay (inherent or deliberately introduced) in the correction torque supplied by the servomotor, which itself may be non-linear to a certain extent. The equation of motion of the platform system is a third order nonlinear non-homogeneous differential equation. An approximate analytical method of solution of this equation is utilized. The value of the delay at which the platform response becomes unstable has been calculated by using this approximate analytical method. The procedure is illustrated by means of a numerical example. Second, the non-linear response of the platform to a random input has been obtained. The effects of several types of non-linearity on reducing the level of the mean square response have been investigated, by applying the technique of equivalent linearization and solving the resulting integral equations by using laguerre or Gaussian integration techniques. The mean square responses to white noise and band limited white noise, for various values of the non-linear parameter and for different types of non-linearity function, have been obtained. For positive values of the non-linear parameter the levels of the non-linear mean square responses to both white noise and band-limited white noise are low as compared to the linear mean square response. For negative values of the non-linear parameter the level of the non-linear mean square response at first increases slowly with increasing values of the non-linear parameter and then suddenly jumps to a high level, at a certain value of the non-linearity parameter.

Algorithms for Colouring Random k-colourable Graphs

The k-colouring problem is to colour a given k-colourable graph with k colours. This problem is known to be NP-hard even for fixed k greater than or equal to 3. The best known polynomial time approximation algorithms require n(delta) (for a positive constant delta depending on k) colours to colour an arbitrary k-colourable n-vertex graph. The situation is entirely different if we look at the average performance of an algorithm rather than its worst-case performance. It is well known that a k-colourable graph drawn from certain classes of distributions can be ii-coloured almost surely in polynomial time. In this paper, we present further results in this direction. We consider k-colourable graphs drawn from the random model in which each allowed edge is chosen independently with probability p(n) after initially partitioning the vertex set into ii colour classes. We present polynomial time algorithms of two different types. The first type of algorithm always runs in polynomial time and succeeds almost surely. Algorithms of this type have been proposed before, but our algorithms have provably exponentially small failure probabilities. The second type of algorithm always succeeds and has polynomial running time on average. Such algorithms are more useful and more difficult to obtain than the first type of algorithms. Our algorithms work as long as p(n) greater than or equal to n(-1+is an element of) where is an element of is a constant greater than 1/4.

Generation of a Manganese Specific Restriction Endonuclease with Nicking Activity

A typical feature of type II restriction endonucleases (REases) is their obligate sequence specificity and requirement for Mg2+ during catalysis. R.KpnI is an exception. Unlike most other type II REases, the active site of this enzyme can accommodate Mg2+, Mn2+, Ca2+, or Zn2+ and cleave DNA. The enzyme belongs to the HNH superfamily of nucleases and is characterized by the presence of a beta beta alpha-Me finger motif. Residues D148, H149, and Q175 together form the HNH active site and are essential for Mg2+ binding and catalysis. The unique ability of the enzyme to cleave DNA in the presence of different metal ions is exploited to generate mutants that are specific to one particular metal ion. We describe the generation of a Mn2+-dependent sequence specific endonuclease, defective in DNA cleavage with Mg2+ and other divalent metal ions. In the engineered mutant, only Mn2+ is selectively bound at the active site, imparting Mn2+-mediated cleavage. The mutant is impaired in concerted double-stranded DNA cleavage, leading to accumulation of nicked intermediates. The nicking activity of the mutant enzyme is further enhanced by altered reaction conditions. The active site fluidity of R Eases allowing flexible accommodation of catalytic cofactors thus forms a basis for engineering selective metal ion-dependent REase additionally possessing nicking activity.

Tsoi Yuen Resistance and Post-80s Generation: A new wave of domestic social movement

Hong Kong was once a British colony and has been under the sovereignty of People’s Republic of China (PRC) since 1997. However, some of the unjust practices and colonial legacies are infiltrated into the development ideology as well as the social structures. The construction of intercity express railway project announced in 2008 causing the demolishment of Tsoi Yuen Tsuen, a “non-indigenous” agricultural village in Hong Kong, was one of the current examples. Tsoi Yuen village was established under the former colonial sovereignty sixty years ago. Approximately 450 populations were affected that they had to relocate their homeland involuntarily. However, these villagers were very attached to their homelands and were unwilling to move, and meanwhile they found that they were absent in the government’s consultation and decision-making process. Soon they began their resistance and demanded for “No Move! No Demolish!”. Their movement was strongly supported by a group of “Post-80s generation” and turned into the most important social movement of the city in recent years. In fact, demolition of Tsoi Yuen Village for city development is not an isolated case in the city. Meanwhile the situation is getting worse in Mainland China. I chose the case study of Tsoi Yuen Resistance from 2008 to 2011 for revelation of the complicated colonial history and postcolonial era of Hong Kong. I focused on discussing the Tsoi Yuen Resistance and the Post-80s movement, and how they have exposed the tension between top-down urban planning and development and public movements fighting for a more democratic process in choosing their way of living. Through the study of a village movement which as well as the rationale behind the Post-80s’ support, I hoped to illustrate how this movement has awaken a different sense of living for the new generations in the midst of the high-sounding urban development. It is an opportunity to examine Hong Kong’s colonial epoch in a different perspective: through studying the Tsoi Yuen Village, let them (subalterns) speak for themselves. Furthermore, the significance of this resistance, taking place eleven years after the handover to the PRC, is an important fact that I shall not miss in later discussion. Last but not least, during the resistance, advanced technology and social networks such as Facebook, Twitter, iPhone were used by Post 80s generation to spread the latest information in order to attract public’s concern and participation. Therefore, apart from studying Tsoi Yuen Resistance as a local social movement, I also regard it as a part of the global movement in perusing ecological lifestyle and civil society. How Post 80s’ generation manipulates the global idea in a local context will also be examined.

Probability Distribution of the Eigenvalues of the Random String Equation

The probability distribution of the eigenvalues of a second-order stochastic boundary value problem is considered. The solution is characterized in terms of the zeros of an associated initial value problem. It is further shown that the probability distribution is related to the solution of a first-order nonlinear stochastic differential equation. Solutions of this equation based on the theory of Markov processes and also on the closure approximation are presented. A string with stochastic mass distribution is considered as an example for numerical work. The theoretical probability distribution functions are compared with digital simulation results. The comparison is found to be reasonably good.

Search on a hypercubic lattice using a quantum random walk. I. d > 2

Random walks describe diffusion processes, where movement at every time step is restricted to only the neighboring locations. We construct a quantum random walk algorithm, based on discretization of the Dirac evolution operator inspired by staggered lattice fermions. We use it to investigate the spatial search problem, that is, to find a marked vertex on a d-dimensional hypercubic lattice. The restriction on movement hardly matters for d > 2, and scaling behavior close to Grover's optimal algorithm (which has no restriction on movement) can be achieved. Using numerical simulations, we optimize the proportionality constants of the scaling behavior, and demonstrate the approach to that for Grover's algorithm (equivalent to the mean-field theory or the d -> infinity limit). In particular, the scaling behavior for d = 3 is only about 25% higher than the optimal d -> infinity value.

Search on a hypercubic lattice using a quantum random walk. II. d = 2

We investigate the spatial search problem on the two-dimensional square lattice, using the Dirac evolution operator discretized according to the staggered lattice fermion formalism. d = 2 is the critical dimension for the spatial search problem, where infrared divergence of the evolution operator leads to logarithmic factors in the scaling behavior. As a result, the construction used in our accompanying article A. Patel and M. A. Rahaman, Phys. Rev. A 82, 032330 (2010)] provides an O(root N ln N) algorithm, which is not optimal. The scaling behavior can be improved to O(root N ln N) by cleverly controlling the massless Dirac evolution operator by an ancilla qubit, as proposed by Tulsi Phys. Rev. A 78, 012310 (2008)]. We reinterpret the ancilla control as introduction of an effective mass at the marked vertex, and optimize the proportionality constants of the scaling behavior of the algorithm by numerically tuning the parameters.

Directing role of functional groups in selective generation of C-H center dot center dot center dot pi interactions: In situ cryo-crystallographic studies on benzyl derivatives

The in situ cryo-crystallization study of benzyl derivatives reveals that the molecular packing in these compounds is either through methylene (sp(3)) C-H center dot center dot center dot pi or aromatic (sp(2)) C-H center dot center dot center dot pi interactions depending on the level of acidity of the benzyl proton. These studies of low melting compounds bring out the subtle features of such weak interactions and point to the directional preferences depending on the nature (electron withdrawing, polarizability) of the neighbouring functional group.

Criticality of the exponential rate of decay for the largest nearest-neighbor link in random geometric graphs

Let n points be placed independently in d-dimensional space according to the density f(x) = A(d)e(-lambda parallel to x parallel to alpha), lambda, alpha > 0, x is an element of R-d, d >= 2. Let d(n) be the longest edge length of the nearest-neighbor graph on these points. We show that (lambda(-1) log n)(1-1/alpha) d(n) - b(n) converges weakly to the Gumbel distribution, where b(n) similar to ((d - 1)/lambda alpha) log log n. We also prove the following strong law for the normalized nearest-neighbor distance (d) over tilde (n) = (lambda(-1) log n)(1-1/alpha) d(n)/log log n: (d - 1)/alpha lambda <= lim inf(n ->infinity) (d) over tilde (n) <= lim sup(n ->infinity) (d) over tilde (n) <= d/alpha lambda almost surely. Thus, the exponential rate of decay alpha = 1 is critical, in the sense that, for alpha > 1, d(n) -> 0, whereas, for alpha <= 1, d(n) -> infinity almost surely as n -> infinity.

Narrowband random excitation of a limit cycle system

The response of the Van der Pol oscillator to stationary narrowband Gaussian excitation is considered. The central frequency of excitation is taken to be in the neighborhood of the system limit cycle frequency. The solution is obtained using a non-Gaussian closure approximation on the probability density function of the response. The validity of the solution is examined with the help of a stochastic stability analysis. Solution based on Stratonovich''s quasistatic averaging technique is also obtained. The comparison of the theoretical solutions with the digital simulations shows that the theoretical estimates are reasonably good.

Rocking response of rectangular rigid blocks under random noise base excitations

The response of a rigid rectangular block resting on a rigid foundation and acted upon simultaneously by a horizontal and a vertical random white-noise excitation is considered. In the equation of motion, the energy dissipation is modeled through a viscous damping term. Under the assumption that the body does not topple, the steady-state joint probability density function of the rotation and the rotational velocity is obtained using the Fokker-Planck equation approach. Closed form solution is obtained for a specific combination of system parameters. A more general but approximate solution to the joint probability density function based on the method of equivalent non-linearization is also presented. Further, the problem of overturning of the block is approached in the framework of the diffusion methods for first passage failure studies. The overturning of the block is deemed incipient when the response trajectories in the phase plane cross the separatrix of the conservative unforced system. Expressions for the moments of first passage time are obtained via a series solution to the governing generalized Pontriagin-Vitt equations. Numerical results illustra- tive of the theoretical solutions are presented and their validity is examined through limited amount of digital simulations.

In-Network Computation in Random Wireless Networks: A PAC Approach to Constant Refresh Rates with Lower Energy Costs

We propose a method to compute a probably approximately correct (PAC) normalized histogram of observations with a refresh rate of Theta(1) time units per histogram sample on a random geometric graph with noise-free links. The delay in computation is Theta(root n) time units. We further extend our approach to a network with noisy links. While the refresh rate remains Theta(1) time units per sample, the delay increases to Theta(root n log n). The number of transmissions in both cases is Theta(n) per histogram sample. The achieved Theta(1) refresh rate for PAC histogram computation is a significant improvement over the refresh rate of Theta(1/log n) for histogram computation in noiseless networks. We achieve this by operating in the supercritical thermodynamic regime where large pathways for communication build up, but the network may have more than one component. The largest component however will have an arbitrarily large fraction of nodes in order to enable approximate computation of the histogram to the desired level of accuracy. Operation in the supercritical thermodynamic regime also reduces energy consumption. A key step in the proof of our achievability result is the construction of a connected component having bounded degree and any desired fraction of nodes. This construction may also prove useful in other communication settings on the random geometric graph.