Open Access. Powered by Scholars. Published by Universities.®

Numerical Analysis and Computation Commons™

Open Access. Powered by Scholars. Published by Universities.®

Discipline
Institution
Keyword
Publication Year
Publication
Publication Type
File Type

Articles 1291 - 1320 of 1368

Full-Text Articles in Numerical Analysis and Computation

Estimating The Accuracy Of Polymerase Chain Reaction-Based Tests Using Endpoint Dilution, Jim Hughes, Patricia Totten Mar 2003

Estimating The Accuracy Of Polymerase Chain Reaction-Based Tests Using Endpoint Dilution, Jim Hughes, Patricia Totten

UW Biostatistics Working Paper Series

PCR-based tests for various microorganisms or target DNA sequences are generally acknowledged to be highly "sensitive" yet the concept of sensitivity is ill-defined in the literature on these tests. We propose that sensitivity should be expressed as a function of the number of target DNA molecules in the sample (or specificity when the target number is 0). However, estimating this "sensitivity curve" is problematic since it is difficult to construct samples with a fixed number of targets. Nonetheless, using serially diluted replicate aliquots of a known concentration of the target DNA sequence, we show that it is possible to disentangle …

Go to article

Simple Parallel Statistical Computing In R, Anthony Rossini, Luke Tierney, Na Li Mar 2003

Simple Parallel Statistical Computing In R, Anthony Rossini, Luke Tierney, Na Li

UW Biostatistics Working Paper Series

Theoretically, many modern statistical procedures are trivial to parallelize. However, practical deployment of a parallelized implementation which is robust and reliably runs on different computational cluster configurations and environments is far from trivial. We present a framework for the R statistical computing language that provides a simple yet powerful programming interface to a computational cluster. This interface allows the development of R functions that distribute independent computations across the nodes of the computational cluster. The resulting framework allows statisticians to obtain significant speed-ups for some computations at little additional development cost. The particular implementation can be deployed in heterogeneous computing …

Go to article

Literate Statistical Practice, Anthony Rossini, Friedrich Leisch Mar 2003

Literate Statistical Practice, Anthony Rossini, Friedrich Leisch

UW Biostatistics Working Paper Series

Literate Statistical Practice (LSP, Rossini, 2001) describes an approach for creating self-documenting statistical results. It applies literate programming (Knuth, 1992) and related techniques in a natural fashion to the practice of statistics. In particular, documentation, specification, and descriptions of results are written concurrently with writing and evaluation of statistical programs. We discuss how and where LSP can be integrated into practice and illustrate this with an example derived from an actual statistical consulting project. The approach is simplified through the use of a comprehensive, open source toolset incorporating Noweb, Emacs Speaks Statistics (ESS), Sweave (Ramsey, 1994; Rossini, et al, 2002; …

Go to article

Checking Assumptions In Latent Class Regression Models Via A Markov Chain Monte Carlo Estimation Approach: An Application To Depression And Socio-Economic Status, Elizabeth Garrett, Richard Miech, Pamela Owens, William W. Eaton, Scott L. Zeger Jan 2003

Checking Assumptions In Latent Class Regression Models Via A Markov Chain Monte Carlo Estimation Approach: An Application To Depression And Socio-Economic Status, Elizabeth Garrett, Richard Miech, Pamela Owens, William W. Eaton, Scott L. Zeger

Johns Hopkins University, Dept. of Biostatistics Working Papers

Latent class regression models are useful tools for assessing associations between covariates and latent variables. However, evaluation of key model assumptions cannot be performed using methods from standard regression models due to the unobserved nature of latent outcome variables. This paper presents graphical diagnostic tools to evaluate whether or not latent class regression models adhere to standard assumptions of the model: conditional independence and non-differential measurement. An integral part of these methods is the use of a Markov Chain Monte Carlo estimation procedure. Unlike standard maximum likelihood implementations for latent class regression model estimation, the MCMC approach allows us to …

Go to article

Impulse Control Of Stochastic Navier-Stokes Equations, J. L. Menaldi, S. S. Sritharan Jan 2003

Impulse Control Of Stochastic Navier-Stokes Equations, J. L. Menaldi, S. S. Sritharan

Mathematics Faculty Research Publications

In this paper we study stopping time and impulse control problems for stochastic Navier-Stokes equation. Exploiting a local monotonicity property of the nonlinearity, we establish existence and uniqueness of strong solutions in two dimensions which gives a Markov-Feller process. The variational inequality associated with the stopping time problem and the quasi-variational inequality associated with the impulse control problem are resolved in a weak sense, using semigroup approach with a convergence uniform over path.

Go to article

Simulation Of Dynamic Electrochemical Processes, John Cassidy Jan 2003

Simulation Of Dynamic Electrochemical Processes, John Cassidy

Articles

This work is designed to introduce electrochemists in a tutorial manner to the basics of modeling of electrochemical systems based primarily on diffusion equations. There is an introduction to analytical and numerical methods with examples taken from typical electrochemical experiments. The Laplace transform is used to derive the Cottrell equation and chronopotentiometry. The response of an electrode to a Gaussian concentration profile is detailed. Laplace’s equation is solved for a simple cell to determine the potential distribution. Discrete methods are employed to calculate the current time behavior following a potential step using the explicit finite difference method. Cyclic voltammetry is …

Go to article

An Empirical Study Of Marginal Structural Models For Time-Independent Treatment, Tanya A. Henneman, Mark J. Van Der Laan Oct 2002

An Empirical Study Of Marginal Structural Models For Time-Independent Treatment, Tanya A. Henneman, Mark J. Van Der Laan

U.C. Berkeley Division of Biostatistics Working Paper Series

In non-randomized treatment studies a significant problem for statisticians is determining how best to adjust for confounders. Marginal structural models (MSMs) and inverse probability of treatment weighted (IPTW) estimators are useful in analyzing the causal effect of treatment in observational studies. Given an IPTW estimator a doubly robust augmented IPTW (AIPTW) estimator orthogonalizes it resulting in a more e±cient estimator than the IPTW estimator. One purpose of this paper is to make a practical comparison between the IPTW estimator and the doubly robust AIPTW estimator via a series of Monte- Carlo simulations. We also consider the selection of the optimal …

Go to article

A Review Of Selected Works On Crack Indentification, Kurt M. Bryan Oct 2002

A Review Of Selected Works On Crack Indentification, Kurt M. Bryan

Mathematical Sciences Technical Reports (MSTR)

We give a short survey of some of the results obtained within the last 10 years or so concerning crack identification using impedance imaging techniques. We touch upon uniqueness results, continuous dependence results, and computational algorithms.

Go to article

Characterizing A Defect In A One-Dimensional Bar, Cynthia Gangi, Sameer Shah Sep 2002

Characterizing A Defect In A One-Dimensional Bar, Cynthia Gangi, Sameer Shah

Mathematical Sciences Technical Reports (MSTR)

We examine the inverse problem of locating and describing an internal point defect in a one­ dimensional rod W by controlling the heat inputs and measuring the subsequent temperatures at the boundary of W. We use a variation of the forward heat equation to model heat flow through W, then propose algorithms for locating an internal defect and quantifying the effect the defect has on the heat flow. We implement these algorithms, analyze the stability of the procedures, and provide several computational examples.

Go to article

Reconstruction Of Cracks With Unknown Transmission Condition From Boundary Data, F Ronald Ogborne Iii, Melissa E. Vellela Sep 2002

Reconstruction Of Cracks With Unknown Transmission Condition From Boundary Data, F Ronald Ogborne Iii, Melissa E. Vellela

Mathematical Sciences Technical Reports (MSTR)

We examine the problem of Identifying both the location and constitutive law governing electrical current flow across a one-dimensional linear crack in a two dimensional region when the crack only partially blocks the flow of current. We develop a a constructive numerical procedure for solving the inverse problem and provide computational examples.

Go to article

Bivariate Current Status Data, Mark J. Van Der Laan, Nicholas P. Jewell Sep 2002

Bivariate Current Status Data, Mark J. Van Der Laan, Nicholas P. Jewell

U.C. Berkeley Division of Biostatistics Working Paper Series

In many applications, it is often of interest to estimate a bivariate distribution of two survival random variables. Complete observation of such random variables is often incomplete. If one only observes whether or not each of the individual survival times exceeds a common observed monitoring time C, then the data structure is referred to as bivariate current status data (Wang and Ding, 2000). For such data, we show that the identifiable part of the joint distribution is represented by three univariate cumulative distribution functions, namely the two marginal cumulative distribution functions, and the bivariate cumulative distribution function evaluated on the …

Go to article

A Posteriori Error Estimates Based On Polynomial Preserving Recovery, Zhimin Zhang, Ahmed Naga Sep 2002

A Posteriori Error Estimates Based On Polynomial Preserving Recovery, Zhimin Zhang, Ahmed Naga

Mathematics Research Reports

Superconvergence of order O(h1+rho), for some rho is greater than 0, is established for gradients recovered using Polynomial Preserving Recovery technique when the mesh is mildly structured. Consequently this technique can be used in building a posteriori error estimator that is asymptotically exact.

Go to article

Gradient Recovery And A Posteriori Estimate For Bilinear Element On Irregular Quadrilateral Meshes, Zhimin Zhang Aug 2002

Gradient Recovery And A Posteriori Estimate For Bilinear Element On Irregular Quadrilateral Meshes, Zhimin Zhang

Mathematics Research Reports

A polynomial preserving gradient recovery method is proposed and analyzed for bilinear element under general quadrilateral meshes. It has been proven that the recovered gradient converges at a rate O(h1+rho) for rho = min(alpha, 1) when the mesh is distorted O(h1+alpha) (alpha > 0) from a regular one. Consequently, the a posteriori error estimator based on the recovered gradient is asymptotically exact.

Go to article

Fast Reconstruction Of Cracks Using Boundary Measurements, Nicholas A. Trainor, Rachel M. Krieger Jul 2002

Fast Reconstruction Of Cracks Using Boundary Measurements, Nicholas A. Trainor, Rachel M. Krieger

Mathematical Sciences Technical Reports (MSTR)

This paper develops a fast algorithm for locating one or more perfectly insulating, pair-wise disjoint, linear cracks in a homogeneous two-dimensional electrical conductor, using boundary measurements.

Go to article

Analysis Of Recovery Type A Posteriori Error Estimators For Mildly Structured Grids, Jinchao Xu, Zhimin Zhang Jul 2002

Analysis Of Recovery Type A Posteriori Error Estimators For Mildly Structured Grids, Jinchao Xu, Zhimin Zhang

Mathematics Research Reports

Some recovery type error estimators for linear finite element method are analyzed under O(h1+alpha) (alpha greater than 0) regular grids. Superconvergence is established for recovered gradients by three different methods when solving general non-self-adjoint second-order elliptic equations. As a consequence, a posteriori error estimators based on those recovery methods are asymptotically exact.

Go to article

A Meshless Gradient Recovery Method Part I: Superconvergence Property, Zhiming Zhang, Ahmed Naga Feb 2002

A Meshless Gradient Recovery Method Part I: Superconvergence Property, Zhiming Zhang, Ahmed Naga

Mathematics Research Reports

A new gradient recovery method is introduced and analyzed. It is proved that the method is superconvergent for translation invariant finite element spaces of any order. The method maintains the simplicity, efficiency, and superconvergence properties of the Zienkiewicz-Zhu patch recovery method. In addition, under uniform triangular meshes, the method is superconvergent for the Chevron pattern, and ultraconvergence at element edge centers for the regular pattern.

Go to article

A New Partitioning Around Medoids Algorithm, Mark J. Van Der Laan, Katherine S. Pollard, Jennifer Bryan Feb 2002

A New Partitioning Around Medoids Algorithm, Mark J. Van Der Laan, Katherine S. Pollard, Jennifer Bryan

U.C. Berkeley Division of Biostatistics Working Paper Series

Kaufman & Rousseeuw (1990) proposed a clustering algorithm Partitioning Around Medoids (PAM) which maps a distance matrix into a specified number of clusters. A particularly nice property is that PAM allows clustering with respect to any specified distance metric. In addition, the medoids are robust representations of the cluster centers, which is particularly important in the common context that many elements do not belong well to any cluster. Based on our experience in clustering gene expression data, we have noticed that PAM does have problems recognizing relatively small clusters in situations where good partitions around medoids clearly exist. In this …

Go to article

Dimensionality-Reducing Expansion, Boundary Type Quadrature Formulas, And The Boundary Element Method, Tian-Xiao He Jan 2002

Dimensionality-Reducing Expansion, Boundary Type Quadrature Formulas, And The Boundary Element Method, Tian-Xiao He

Scholarship

This paper discusses the connection between boundary quadrature formulas constructed by using solutions of partial differential equations and boundary element schemes.

Go to article

Discrete Maximum Principle For Nonsmooth Optimal Control Problems With Delays, Boris S. Mordukhovich, Ilya Shvartsman Dec 2001

Discrete Maximum Principle For Nonsmooth Optimal Control Problems With Delays, Boris S. Mordukhovich, Ilya Shvartsman

Mathematics Research Reports

We consider optimal control problems for discrete-time systems with delays. The main goal is to derive necessary optimality conditions of the discrete maximum principle type in the case of nonsmooth minimizing functions. We obtain two independent forms of the discrete maximum principle with transversality conditions described in terms of subdifferentials and superdifferentials, respectively. The superdifferential form is new even for non-delayed systems and may be essentially stronger than a more conventional subdifferential form in some situations.

Go to article

Maximum Likelihood Estimation Of Ordered Multinomial Parameters, Nicholas P. Jewell, John D. Kalbfleisch Oct 2001

Maximum Likelihood Estimation Of Ordered Multinomial Parameters, Nicholas P. Jewell, John D. Kalbfleisch

U.C. Berkeley Division of Biostatistics Working Paper Series

The pool-adjacent violator-algorithm (Ayer, et al., 1955) has long been known to give the maximum likelihood estimator of a series of ordered binomial parameters, based on an independent observation from each distribution (see Barlow et al., 1972). This result has immediate application to estimation of a survival distribution based on current survival status at a set of monitoring times. This paper considers an extended problem of maximum likelihood estimation of a series of ‘ordered’ multinomial parameters. By making use of variants of the pool adjacent violator algorithm, we obtain a simple algorithm to compute the maximum likelihood estimator and demonstrate …

Go to article

Ultraconvergence Of Zz Patch Recovery At Mesh Symmetry Points, Zhimin Zhang, Runchang Lin Sep 2001

Ultraconvergence Of Zz Patch Recovery At Mesh Symmetry Points, Zhimin Zhang, Runchang Lin

Mathematics Research Reports

Ultraconvergence property of the Zienkiewicz-Zhu gradient patch recovery technique based on local discrete least squares fitting is established for a large class of even-order finite elements. The result is valid at all rectangular mesh symmetry points. Different smoothing strategies are discussed. Superconvergence recovery for the Q8 element is proved and ultraconvergence numerical examples are demonstrated.

Go to article

Emacs Speaks Statistics: A Universal Interface For Statistical Analysis, Anthony Rossini, Martin Maechler, Kurt Hornik, Richard M. Heiberger, Rodney Sparapani Jul 2001

Emacs Speaks Statistics: A Universal Interface For Statistical Analysis, Anthony Rossini, Martin Maechler, Kurt Hornik, Richard M. Heiberger, Rodney Sparapani

UW Biostatistics Working Paper Series

Emacs Speaks Statistics (ESS) is a user interface for developing statistical applications and performing data analysis using any of several common statistical programming languages. ESS falls in the programming tools category of Integrated Development Environments (IDEs), which are approaches for developing and visualizing computer programs. We discuss how it works, the advantages of using it, and extensions for increasing statistical programming efficiency.

Go to article

Alternative Principal Components Regression Procedures For Dendrohydrologic Reconstructions, Hugo G. Hidalgo, Thomas C. Piechota, John A. Dracup Nov 2000

Alternative Principal Components Regression Procedures For Dendrohydrologic Reconstructions, Hugo G. Hidalgo, Thomas C. Piechota, John A. Dracup

Civil and Environmental Engineering and Construction Faculty Research

Streamflow reconstruction using tree ring information (dendrohydrology) has traditionally used principal components analysis (PCA) and stepwise regression to form a transfer function. However, PCA has several procedural choices that may result in very different reconstructions. This study assesses the different procedures in PCA-based regression and suggests alternative procedures for selection of variables and principal components. Cross-validation statistics are presented as an alternative for independently testing and identifying the optimal model. The objective is to use these statistics as a measure of the model's performance to find a conceptually acceptable model with a low prediction error and the fewest number of …

Go to article

Stochastic Hybrid Control, A. Bensoussan, J. L. Menaldi Sep 2000

Stochastic Hybrid Control, A. Bensoussan, J. L. Menaldi

Mathematics Faculty Research Publications

The objective of this paper is to study the stochastic version of a previous paper of the authors, in which hybrid control for deterministic systems was considered. The modelling is quite similar to the deterministic case. We have a system whose state is composed of a continuous part and a discrete part. They are affected by a continuous type control and an impulse control. The dynamics is moreover perturbed by noise, also a continuous and a discrete noise process. The Markovian character of the state process is preserved. We develop the model and show how the dynamic programming approach leads …

Go to article

Collected Papers Vol. Iii, Florentin Smarandache Jan 2000

Collected Papers Vol. Iii, Florentin Smarandache

Branch Mathematics and Statistics Faculty and Staff Publications

No abstract provided.

Go to article

Mathematical Modelling Of Road Bumps, E. S. Massawe Dec 1999

Mathematical Modelling Of Road Bumps, E. S. Massawe

Tanzania Journal of Engineering and Technology (TJET)

In this paper, a general mathematical model that may be used for construction of a road bump on an inclined road is developed. It is shown that if the parameters for the inclination are set to zero, the resulting model reflects the special case of a flat road. Existing experiences are briefly highlighted and a case study was done on different locations along the University Road to demonstrate the developed model.

Go to article

Application Of Vibrational Techniques In Determination Of Dynamic Properties Of Agricultural Products-State Of The Arton Of Vibrational Techniques In Determination Of Dynamic Properties Of Agricultural Products-State Of The Art, Silas Kajuna Dec 1999

Application Of Vibrational Techniques In Determination Of Dynamic Properties Of Agricultural Products-State Of The Arton Of Vibrational Techniques In Determination Of Dynamic Properties Of Agricultural Products-State Of The Art, Silas Kajuna

Tanzania Journal of Engineering and Technology (TJET)

Vibration is one of the techniques employed in the determination of dynamic properties of fruits and vegetables. It entails generation of a mechanical or acoustic vibrational signal which is propagated through the flesh of the agricultural material. A transducer is either attached or held close to the specimen to monitor the propagation of the signal through the specimen. The manner in which the signal is transmitted through the material is analyzed, and the dynamic properties of the specimen which relate to its firmness or its internal being are derived. The technique has been around for the past 30 years or …

Go to article

Development And Application Of Adaptive Dynamic Grids For Wake Vortex/Ground Interaction, Ilker Emin Türkgeldi Jul 1999

Development And Application Of Adaptive Dynamic Grids For Wake Vortex/Ground Interaction, Ilker Emin Türkgeldi

Mechanical & Aerospace Engineering Theses & Dissertations

An adaptive and dynamic grid approach is developed for the numerical computations of wake vortex flow problems with ground interaction. The unnecessary grid points at the far field of the original, large, single zone grid are removed for computational efficiency. Next, the grid is modified into an adaptive one, where the fine core region keeps track of the vortices and can be split into two regions when the vortex separation distance exceeds a specified limit length. In the final step a multi zonal feature is added to this grid by means of two embedded grids placed at the vortex centers …

Go to article

Invariant Measure For Diffusions With Jumps, Jose-Luis Menaldi, Maurice Robin Feb 1999

Invariant Measure For Diffusions With Jumps, Jose-Luis Menaldi, Maurice Robin

Mathematics Faculty Research Publications

Our purpose is to study an ergodic linear equation associated to diffusion processes with jumps in the whole space. This integro-differential equation plays a fundamental role in ergodic control problems of second order Markov processes. The key result is to prove the existence and uniqueness of an invariant density function for a jump diffusion, whose lower order coefficients are only Borel measurable. Based on this invariant probability, existence and uniqueness (up to an additive constant) of solutions to the ergodic linear equation are established.

Go to article

Stability And Reconstruction For An Inverse Problem For The Heat Equations, Kurt M. Bryan, Lester Caudill May 1998

Stability And Reconstruction For An Inverse Problem For The Heat Equations, Kurt M. Bryan, Lester Caudill

Mathematical Sciences Technical Reports (MSTR)

We examine the inverse problem of determining the shape of some unknown portion of the boundary of a region W from measurements of the Cauchy data for solutions to the heat equation on W. By suitably linearizing the inverse problem we obtain uniqueness and continuous dependence results. We propose an algorithm for recovering estimates of the unknown portion of the surface and use the insight gained from a detailed analysis of the inverse problem to regularize the inversion. Several computational examples are presented.

Go to article

First Prev Next Last