The Limits of Multiplexing
|
|
- Loreen Dennis
- 6 years ago
- Views:
Transcription
1 Article type: Focus Article The Limits of Multiplexing Dan Shen,, D.P. Dittmer 2, an J. S. Marron 3 Keywors Multiplexing, Genomics, nanostring TM, DASL TM, Probability Abstract We were motivate by three novel technologies, which exemplify a new esign paraigm in high throughput genomics: nanostring TM, DNA-meiate Annealing, Selection, extension, an Ligation (DASL TM an multiplex real-time quantitative polymerase chain reaction (QPCR. All three are solution hybriization base, an all three employ on - DNA sequence probes in a small volume, each probe specific for a particular sequence in a ifferent human gene. Nanostring TM uses 5-mer, DASL an multiplex QPCR use 2-mer probes. Assuming a nm probe concentration in a µl volume, there are 9 x 9 x 6.23 x 23 or 6.23 x 5 molecules of each probe present in the reaction compare to -, target molecules. Excess probe rives the sensitivity of the reaction. We are intereste in the limits of multiplexing, i.e. the probability that in such a esign a particular probe woul bin to any other, sequence-relate probe rather than the intene, specific target. If this were to happen with appreciable frequency, this woul result in much reuce sensitivity an potential failure of this esign. We establishe upper an lower bouns for the probability that in a multiplex assay at least one probe woul bin to another sequencerelate probe rather than its cognate target. These bouns are reassuring, because for reasonable egrees of multiplexing ( 3 probes the probability for such an event is practically negligible. As the egree of multiplexing increases to 6 probes, our theoretical bounaries gain practical importance an establish a principal upper limit for the use of highly multiplexe solution-base assays vis-à-vis soli-support anchore esigns. Recently solution-base multiplex hybriization base methos have been evelope an use for messenger ribonucleic aci (mrna profiling experiments that were previously the purview of soli-state, anchore methos the so-calle microarrays or chips. By most practical accounts their performance seems * Correspon to: anshen@usf.eu; ittmer@me.unc.eu Interisciplinary Data Sciences Consortium, Department of Mathematics an Statistics, University of South Floria. 2 Department of Microbiology an Immunology an Center for AIDS Research, Comprehensive Cancer Center, University of North Carolina at Chapel Hill. 3 Department of Statistics an Operations Research, Comprehensive Cancer Center, University of North Carolina at Chapel Hill.
2 equal, but practical experiments represent examples with a bias towars reporting positive ata. They are not exhaustive an o not represent a general solution. Because massive multiplexing involves - 6 probes, iniviual experimental valiation is no longer feasible. At the time the reaction volume of polymerase chain reaction (PCR an hybriization assays has been reuce ue to nanotechnology. Conventional PCR instruments now can use the 536-well format e.g. the Roche system with µl. Newer microfluiics-base machines perform up to 2, iniviual reactions on the same chip e.g. Fluiigm system with.85 - nanoliter. Recently, a picoliter evice has been escribe by White et al. (2. As the engineering pushes the technical bounaries of miniaturization, it becomes important to efine the statistical bounaries of experimental esigns. Problem We were concerne that the multiplex esign introuces the potential for cross-hybriization among probe molecules, resulting in a loss of sensitivity or etection failure. This possibility is not present in soli-state, anchore esigns, since in this case there are no free probes available in solution only the target molecules. Traitionally crosshybriization refers to a scenario where the probe bins to a secon, relate but not the intene target. In solution-base multiplex esigns there exists in aition the possibility that the probe cross-hybriizes to another probe rather than any target at all. Since every probe has to be in excess over any potential target in orer to rive the hybriization reaction to completion, cross-hybriization to an unrelate probe woul be favore an woul prevent etection of the cognate target. Here we etermine theoretical bouns for this cross-hybriization problem. Figure lays out the problem. A microarray can be seen as a set S of probes in a special orientation. Each probe of length (soli arrow is physically attache to the support surface. Two probes never touch each other. These are then hybriize to a mixture of target mrnas. The probes are in molar excess compare to the target. Each probe is an oligonucleotie of length, i.e. a -mer. The -mer is mae up of the four bases A, T, C, G. We assume that A an T have the same frequency, as o C an G. Furthermore. we parameterize the CG-ratio (i.e. the relative frequency of A or T as p, where p. Each -mer bins to the target mrna with perfect complementarity an we assume percent efficiency. We also assume that a given -mer oes not bin non-specifically to non-target mrnas. These assumptions are suppose to hol in maximum efficient instruments or assays. In a soli-state anchore microarray, this is the only interaction that can take place (Figure panel A. There is no bining of -mer probes to each other. The situation is ifferent in a solution-hybriization base multiplex esign, e.g. a bea array (Figure panel B. In a bea array probes are couple to iniviual beas an in principle any two beas with complementary probes (blue an re arrows in Figure can hybriize to each other. Other esign rely on free oligonucloties/ probes with no 2
3 bea array Probe: -mer microarray Set of S probes Probe: -mer with m mismatches target mrna A. B. Figure : Conceptual illustrations of bining possibility for (A soli-state anchore micro array, (B multiplex solution array. Soli arrows inicate probes, the otte arrow the correct target. The blue arrow refers to a probe or oligonucleotie of length, which is esire, perfectly complementary to the target mrna (ashe arrow. The re arrow refers to an oligonucleotie of length, which is similar to the target mrna an thus can bin the blue probe except for m mismatches. The blue-re interactions are possible ue to sequence complementarity. beas attache. In aition to each -mer bining is cognate target mrna, each -mer can also bin to any other -mer in the probe set. If the -mer bins to another -mer rather than the target mrna, the assay fails. It is important to keep in min that for solution hybriization base esigns the concentration of probe is orers of magnitue large than the target mrna that is being etecte. Otherwise, the assay woul not be quantitative. Since hybriization efficiency is a function of probe concentration, unwante probe -mer to probe -mer hybriizations poses a novel problem for solution-base multiplex approaches. In a complete set that contains all possible -mers, there exists for each -mer one an only one perfectly complementary -mer. For example, for a -mer of sequence ACT G the perfect complement woul be T GAC. In a complete set S, all -mer probes woul hybriize to their complementary -mer probe an none woul bin to the target mrna. This is avoie in traitional microarray esigns, since all -mers are spatially separate by anchoring them to a soli support matrix, i.e. a slie microarray or chip. The size of the complete set, i.e. the number of all possible -mers is 4. In praxi, one woul never multiplex the complete set, but only a subset of all possible -mers. This subset S is much smaller than the complete set S. In praxi, not only perfectly completely complementary -mers woul hybriize, but also those with m mismatches. The number of mismatches being etermine by the stringency of the hybriization reaction. The more mismatches are tolerate by the reaction conitions the more the size of S approaches the size of the complete set S. Suppose = 4, then the size of S = 64. Assume we only have 3
4 Probability Probability Probability Panel A (= 4, s=, p=.5 Panel C (= 4, s=, p= Lower boun Upper boun Panel E (= 4, s=, p= m, number of mismatches log ( Probability log ( Probability log ( Probability Panel B (= 4, s=, p=.5 5 Lower boun Upper boun Panel D (= 4,s=,p= Panel F (= 4, s=, p= m, number of mismatches Figure 2: Upper (soli an lower (ashe bouns on the probability that there exist no m mismatche -mers in a subset uner the ifferent CG-ratio: Show a absolute probability scale in the left panels. To epict closeness to at higher resolution, the right panels are plotte on the log(-probability scale. Rows compare affect of the CG-ratio p. Show upper an lower bouns are equally very close. Dashe re an pink lines give specific interesting comparison. one probe of length = 4, then size of S =, e.g. ACT G. S /S = /64. Now we allow one mismatch at the en to yiel : ACT g, ACT t, ACT c, ACT a. The size of S = 4 an S /S = 4/64. If we allow 4 mismatches S = S. This paper gives useful theoretical bouns on how many -mers can be multiplexe an how these epen on the length, the number of mismatches an the CG-ratio p. There exists a vibrant literature on the probability of an iniviual -mers an the complete set S of all possible sequence permutation, D yachkov et al. (25; Bishop et al. (27; Dyachkov an Voronina (29. Fewer stuies have investigate this problem in the context of subsets S of S an how subset size influences the probability of annealing. Results We efine no m mismatche -mers in a subset S if the number of mismatches between any two -mers in a subset S oesn t equal to m. The probability of no m mismatche -mers in a subset S of all possible -mers is a function of s (the size of the subset S, of (the length of the -mer, of m (the number of mismatches an of p (the CG-ratio. Panel A of the Figure 2 show our lower an upper bouns on the probability of no m mismatche -mers in the subset S for s =, = 4 an p =.5. Lower an upper bouns on this probability are close to when m is small, then ecrease to for m between 6 an 24 an increase again for m near = 4. 4
5 Panels C an E show the lower an upper bouns plot with the same s = an = 4 but a ifferent CG-ratio p =.28 an p =.2 respectively. The boun curves in Panels C an E have a similar tren as in Panel A. For smaller values of the CGratio, the sharp ecrease from to happens for smaller m. Also the increase happens for larger m, an oes not occur for p =.2 in Panel C. Panels of Figure 2 also show that the lower an upper bouns on the probability of no m mismatche -mers in the subset S is almost equal to when m is small. For example, the vertical scale in the panel A oes not effectively istinguish the bouns from, for m < 5. A better visualization of this practically important range is achieve by applying the log function to minus the bouns as shown in the secon column panels of Figure 2. These transforme plots clearly show the orer of magnitue of the ifference between the probability an. For example, the ashe re an pink lines in the plot show that to have a probability within. = 3 of, we nee m 3, 8, or 6, for the CG-ratios p =.5,.28, or.2 respectively, an to be within 6, m, 5, or 3 respectively. As mentione above, the combinatorial problem (Graham, 995; Rioran, 22 of fining a general close form for the probability of no m mismatche -mers in the subset S is very challenging. This motivate us to instea fin close forms for the lower an upper bouns on the probability of no m mismatche -mers in the subset S. Without loss of generality, we assume that the CG-ratio p.5. First, we nee to introuce some notation. Given, m an s, we efine Ni l an N i u, i =,, s, as If i 2, Ni l = ip an Ni u = i( p. If 2 l k= < i 2 l k=, where l, then N l i = 2 l k= p k ( p k + p l ( p l (i 2 l k= ( p l p l (i 2 l k=. an N u i = 2 l k= ( p k p k + In aition, we nee to efine Mi l an M i u, i =,, s, as If i ( m < 2, then Mi l = i ( m ( + p m p m an Mi u p m. = i ( m (2 p m ( If 2 l ( k= < i m 2 l k=, where l m, then M l i = 2 l k= ( + p m ( p k p m k + [i ( m 2 l k= ]( + p m ( p l p m l an Mi u = 2 l k= (2 p m ( p m k p k + [i ( m 2 l k= ](2 p m ( p m l p l. If 2 l k= 2 m k= ( < i m 2 l k=, where m < l, then M l i = (+p m ( p k p m k +2 l k= m (2 p k +m (+ 5
6 p k ( p m + [i ( m 2 l k= ](2 p l +m ( + p l ( p m an Mi l = 2 m k= (2 p m p k ( p m k + 2 l k= m ( + p k +m (2 p k p m +[i ( m 2 l k= ](+p l +m (2 p l p m. Note that Ni l N i u an Mi l M i u, i =,, s an equality hols when the CG-ratio p =.5. A convenient notation is: The number of possible results for the first -mer is M = (2p + 2 2p = 2. Given i -mers with no m mismatches, the number of possible -mers that are m-mismatche with the given -mers is greater than or equal to Mi l an less than or equal to Mi u. Given two ifferent -mers that are not m mismatche, the number of possible -mers H that are m-mismatche with one of them is { M l H + ( m 2 m p, m < M l + ( + p p, m =. Let the lower boun be H l = M l + ( m 2 m p for m < an H l = M l + ( + p p for m =. Lower an upper bouns on the probability of no m mismatche -mers in the subset S are: If M M u s N u s <, then If M M u s N u s, then P [no m mismatches] =. ( s i= (M Mi u Ni u (M Ni l P [no m mismatches] ( M M l s ( M max { } H l, Ni u i=2 s i= (M N u i. (2 The erivation of equations ( an (2 is shown in the supplementary material. 6
7 Panel A ( = 4, p=.5 Panel B ( =4, p=.5 Lower boun.5 Upper boun Panel C ( = 4, p=.28 Panel D ( =4, p=.28 Lower boun.5 Upper boun Lower boun.5 5 Panel E ( = 4, p=.2 Upper boun.5 5 Panel F ( =4, p=.2 s m s m Figure 3: Bouns on the probability of no m mismatche = 4-mers in subsets with ifferent sizes: Extension of Figure 2 to a range of ifferent values of s (the size of the subset. This continues to show large probability for small m (more so for small s, for both lower (Panel A, C an E an upper (Panel B, D an F bouns. The boun remain close, inicating goo approximation quality, over a range of ifferent CG-ratios, p =.5 (A an B,.28 (C an D an.2 (E an F. Discussion Given the CG-ratio p, equations ( an (2 give bouns on the probability of no m mismatche -mers in the subset S as functions of s an m. As an example we moele the bouns for s =, the size of the subset S, in Figure 2. This showe that for probes (-mers with 5% (p =.5 CG- content even if we allow as many as mismatches in a probe of length = 4 the probability that any two probes in this set anneal to each other is in,,, i.e. very unlikely. The situation becomes less favorable as the CG- content ratio becomes more skewe. At 2% (p =.2 CG-content, such as experience in certain microorganisms (mycoplasma has 24% CG-content allowing for 8 mismatches yiels a chance of in in, that any two probes woul anneal to each other. Because our moel is symmetric aroun 5% CG-content the same reasoning applies to positively skewe ration such as foun in streptomyces species, which average 72% CG-content. Hence, multiplex assays with > probes are limite to organisms with balance CG- content. Next, we explore how the bouns change when both subset size s an mismatch number m change. We use Figure 3 to illustrate the lower an upper bouns as functions of both the size of the subset s an the number of allowe mismatches m, for the CGratios p =.5,.28,.2. For practical applications, we want to maximize the size of the subset s, which increases the egree of multiplexing, an we want to minimize the number of allowe mismatches m, which increases specificity. As seen from Figure 3, a practical limit of the egree of multiplexing again epens on the CG-content 7
8 of the target organism. Up to a set size of s = probes, though it is extremely unlikely that any two probes in a multiplex assay woul bin to each other. This assumes large probes of length = 4 or longer as use in the Nanostring TM assay. The chance of unwante cross-hybriization increases as the probe length ecreases. At a probe length of = 2, such as use in multiplex PCR an applie to the worst case scenario of a microorganism with heavily skewe CG-content, allowing as little as 2 mismatches per probe may result in cross-hybriization between probes in the probe set. Luckily homo sapiens has a balance CG-content, which allows the use of highly multiplexe assays for clinical applications. Current soli-state microarrays can achieve a size of s =.8 6 ifferent probes per chip an m =, since they can etect single nucleotie polymorphisms (SNPs. Base on our calculation, we can answer the question: Can a solution-base, multiplex esign reach or excee this performance? For = 4, s =.8 6 an m =, we have the probability (of no m mismatches within 9, 3 an 2 2 of for the CG-ratios p =.5,.28, or.2 respectively. Hence, solution-base SNP arrays base on probe sizes of 4 or longer have comparable performance to soli-state microarrays only for balance (CG-ratio p =.5 probes. If the CG-ratio rops, as is known for many microbial genomes, solution-base SNP arrays unerperform ue to cross-hybriization among probes. Mismatch an mismatch probability have a concrete biophysical meaning, see Cantor an Schimmel (98. Every match lowers the free energy G of the probe-target uplex an every mismatch m increases G. Every probe-target uplex has a characteristic melting temperature T, which is a function of G. We show here that it is extremely unlikely that in a set S of size s < we woul encounter any pair of probes of length = 4 with m < 3, 8, or 6 (corresponing to the CG-ratio p =.5,.28 or.2 mismatches between them. In sum, the current multiplex assays ( e.g. nanostring TM, DASL TM are expecte to work an have a large margin of error built in before they encounter the theoretical bounaries, which we erive here. As we move into higher an higher moes of multiplexing, it is important to know the principal bounaries of each esign. As it is no longer possible to experimentally test all possible failure scenarios or experimentally valiate the performance for each an every probe our theoretical unerstaning nees to improve to near certainty. Otherwise the true potential of highly multiplexe methos cannot be realize. Acknowlegements. This work was partially supporte by the Startup Fun of University of South Floria, an public health service grants CA94 an AI78 to DPD. References Bishop, M. A., D yachkov, A. G., Macula, A. J., Renz, T. E. an Rykov, V. V. (27 Free energy gap an statistical thermoynamic fielity of na coes. Journal of Computational Biology 4(8,
9 Cantor, C. R. an Schimmel, P. R. (98 Biophysical chemistry: Part III: the behavior of biological macromolecules. Macmillan. D yachkov, A. G., Vilenkin, P. A., Ismagilov, I. K., Sarbaev, R. S., Macula, A., Torney, D. an White, S. (25 On na coes. Problems of Information Transmission 4(4, Dyachkov, A. G. an Voronina, A. N. (29 Dna coes for aitive stem similarity. Problems of Information Transmission 45(2, Graham, R. L. (995 Hanbook of combinatorics, Volume. Elsevier. Rioran, J. (22 Introuction to combinatorial analysis. Courier Corporation. White, A. K., VanInsberghe, M., Petriv, I., Hamii, M., Sikorski, D., Marra, M. A., Piret, J., Aparicio, S. an Hansen, C. L. (2 High-throughput microfluiic single-cell rt-qpcr. Proceeings of the National Acaemy of Sciences 8(34,
'HVLJQ &RQVLGHUDWLRQ LQ 0DWHULDO 6HOHFWLRQ 'HVLJQ 6HQVLWLYLW\,1752'8&7,21
Large amping in a structural material may be either esirable or unesirable, epening on the engineering application at han. For example, amping is a esirable property to the esigner concerne with limiting
More informationTable of Common Derivatives By David Abraham
Prouct an Quotient Rules: Table of Common Derivatives By Davi Abraham [ f ( g( ] = [ f ( ] g( + f ( [ g( ] f ( = g( [ f ( ] g( g( f ( [ g( ] Trigonometric Functions: sin( = cos( cos( = sin( tan( = sec
More informationEE 370L Controls Laboratory. Laboratory Exercise #7 Root Locus. Department of Electrical and Computer Engineering University of Nevada, at Las Vegas
EE 370L Controls Laboratory Laboratory Exercise #7 Root Locus Department of Electrical an Computer Engineering University of Nevaa, at Las Vegas 1. Learning Objectives To emonstrate the concept of error
More informationSYNCHRONOUS SEQUENTIAL CIRCUITS
CHAPTER SYNCHRONOUS SEUENTIAL CIRCUITS Registers an counters, two very common synchronous sequential circuits, are introuce in this chapter. Register is a igital circuit for storing information. Contents
More informationRFSS: Lecture 4 Alpha Decay
RFSS: Lecture 4 Alpha Decay Reaings Nuclear an Raiochemistry: Chapter 3 Moern Nuclear Chemistry: Chapter 7 Energetics of Alpha Decay Geiger Nuttall base theory Theory of Alpha Decay Hinrance Factors Different
More informationExperimental Robustness Study of a Second-Order Sliding Mode Controller
Experimental Robustness Stuy of a Secon-Orer Sliing Moe Controller Anré Blom, Bram e Jager Einhoven University of Technology Department of Mechanical Engineering P.O. Box 513, 5600 MB Einhoven, The Netherlans
More informationELEC3114 Control Systems 1
ELEC34 Control Systems Linear Systems - Moelling - Some Issues Session 2, 2007 Introuction Linear systems may be represente in a number of ifferent ways. Figure shows the relationship between various representations.
More informationDesign and performance of in vitro transcription rate regulatory circuits
Conference on Decision an Control (CDC), 2008 http://www.cs.caltech.eu/~murray/papers/2008j_fm08-cc.html Design an performance of in vitro transcription rate regulatory circuits Elisa Franco an Richar
More informationMath 342 Partial Differential Equations «Viktor Grigoryan
Math 342 Partial Differential Equations «Viktor Grigoryan 6 Wave equation: solution In this lecture we will solve the wave equation on the entire real line x R. This correspons to a string of infinite
More informationCONTROL CHARTS FOR VARIABLES
UNIT CONTOL CHATS FO VAIABLES Structure.1 Introuction Objectives. Control Chart Technique.3 Control Charts for Variables.4 Control Chart for Mean(-Chart).5 ange Chart (-Chart).6 Stanar Deviation Chart
More informationExperiment 2, Physics 2BL
Experiment 2, Physics 2BL Deuction of Mass Distributions. Last Upate: 2009-05-03 Preparation Before this experiment, we recommen you review or familiarize yourself with the following: Chapters 4-6 in Taylor
More informationLATTICE-BASED D-OPTIMUM DESIGN FOR FOURIER REGRESSION
The Annals of Statistics 1997, Vol. 25, No. 6, 2313 2327 LATTICE-BASED D-OPTIMUM DESIGN FOR FOURIER REGRESSION By Eva Riccomagno, 1 Rainer Schwabe 2 an Henry P. Wynn 1 University of Warwick, Technische
More informationEvaluation of Column Breakpoint and Trajectory for a Plain Liquid Jet Injected into a Crossflow
ILASS Americas, 1 st Annual Conference on Liqui Atomization an Spray Systems, Orlano, Floria, May 008 Evaluation of Column Breakpoint an Trajectory for a Plain Liqui Jet Injecte into a Crossflow S.M. Thawley,
More informationAnalytic Scaling Formulas for Crossed Laser Acceleration in Vacuum
October 6, 4 ARDB Note Analytic Scaling Formulas for Crosse Laser Acceleration in Vacuum Robert J. Noble Stanfor Linear Accelerator Center, Stanfor University 575 San Hill Roa, Menlo Park, California 945
More informationThe Exact Form and General Integrating Factors
7 The Exact Form an General Integrating Factors In the previous chapters, we ve seen how separable an linear ifferential equations can be solve using methos for converting them to forms that can be easily
More informationRole of parameters in the stochastic dynamics of a stick-slip oscillator
Proceeing Series of the Brazilian Society of Applie an Computational Mathematics, v. 6, n. 1, 218. Trabalho apresentao no XXXVII CNMAC, S.J. os Campos - SP, 217. Proceeing Series of the Brazilian Society
More informationII. First variation of functionals
II. First variation of functionals The erivative of a function being zero is a necessary conition for the etremum of that function in orinary calculus. Let us now tackle the question of the equivalent
More informationPolynomial Inclusion Functions
Polynomial Inclusion Functions E. e Weert, E. van Kampen, Q. P. Chu, an J. A. Muler Delft University of Technology, Faculty of Aerospace Engineering, Control an Simulation Division E.eWeert@TUDelft.nl
More informationConstruction of the Electronic Radial Wave Functions and Probability Distributions of Hydrogen-like Systems
Construction of the Electronic Raial Wave Functions an Probability Distributions of Hyrogen-like Systems Thomas S. Kuntzleman, Department of Chemistry Spring Arbor University, Spring Arbor MI 498 tkuntzle@arbor.eu
More informationFinal Exam Study Guide and Practice Problems Solutions
Final Exam Stuy Guie an Practice Problems Solutions Note: These problems are just some of the types of problems that might appear on the exam. However, to fully prepare for the exam, in aition to making
More informationComputing Exact Confidence Coefficients of Simultaneous Confidence Intervals for Multinomial Proportions and their Functions
Working Paper 2013:5 Department of Statistics Computing Exact Confience Coefficients of Simultaneous Confience Intervals for Multinomial Proportions an their Functions Shaobo Jin Working Paper 2013:5
More informationQubit channels that achieve capacity with two states
Qubit channels that achieve capacity with two states Dominic W. Berry Department of Physics, The University of Queenslan, Brisbane, Queenslan 4072, Australia Receive 22 December 2004; publishe 22 March
More informationImplicit Differentiation
Implicit Differentiation Thus far, the functions we have been concerne with have been efine explicitly. A function is efine explicitly if the output is given irectly in terms of the input. For instance,
More information6. Friction and viscosity in gasses
IR2 6. Friction an viscosity in gasses 6.1 Introuction Similar to fluis, also for laminar flowing gases Newtons s friction law hols true (see experiment IR1). Using Newton s law the viscosity of air uner
More informationinflow outflow Part I. Regular tasks for MAE598/494 Task 1
MAE 494/598, Fall 2016 Project #1 (Regular tasks = 20 points) Har copy of report is ue at the start of class on the ue ate. The rules on collaboration will be release separately. Please always follow the
More information7.1 Support Vector Machine
67577 Intro. to Machine Learning Fall semester, 006/7 Lecture 7: Support Vector Machines an Kernel Functions II Lecturer: Amnon Shashua Scribe: Amnon Shashua 7. Support Vector Machine We return now to
More informationWhy Bernstein Polynomials Are Better: Fuzzy-Inspired Justification
Why Bernstein Polynomials Are Better: Fuzzy-Inspire Justification Jaime Nava 1, Olga Kosheleva 2, an Vlaik Kreinovich 3 1,3 Department of Computer Science 2 Department of Teacher Eucation University of
More informationLinear First-Order Equations
5 Linear First-Orer Equations Linear first-orer ifferential equations make up another important class of ifferential equations that commonly arise in applications an are relatively easy to solve (in theory)
More informationShort Intro to Coordinate Transformation
Short Intro to Coorinate Transformation 1 A Vector A vector can basically be seen as an arrow in space pointing in a specific irection with a specific length. The following problem arises: How o we represent
More informationLecture 2 Lagrangian formulation of classical mechanics Mechanics
Lecture Lagrangian formulation of classical mechanics 70.00 Mechanics Principle of stationary action MATH-GA To specify a motion uniquely in classical mechanics, it suffices to give, at some time t 0,
More informationTechnion - Computer Science Department - M.Sc. Thesis MSC Constrained Codes for Two-Dimensional Channels.
Technion - Computer Science Department - M.Sc. Thesis MSC-2006- - 2006 Constraine Coes for Two-Dimensional Channels Keren Censor Technion - Computer Science Department - M.Sc. Thesis MSC-2006- - 2006 Technion
More informationTHE VAN KAMPEN EXPANSION FOR LINKED DUFFING LINEAR OSCILLATORS EXCITED BY COLORED NOISE
Journal of Soun an Vibration (1996) 191(3), 397 414 THE VAN KAMPEN EXPANSION FOR LINKED DUFFING LINEAR OSCILLATORS EXCITED BY COLORED NOISE E. M. WEINSTEIN Galaxy Scientific Corporation, 2500 English Creek
More informationThermal conductivity of graded composites: Numerical simulations and an effective medium approximation
JOURNAL OF MATERIALS SCIENCE 34 (999)5497 5503 Thermal conuctivity of grae composites: Numerical simulations an an effective meium approximation P. M. HUI Department of Physics, The Chinese University
More informationIERCU. Institute of Economic Research, Chuo University 50th Anniversary Special Issues. Discussion Paper No.210
IERCU Institute of Economic Research, Chuo University 50th Anniversary Special Issues Discussion Paper No.210 Discrete an Continuous Dynamics in Nonlinear Monopolies Akio Matsumoto Chuo University Ferenc
More informationOn the enumeration of partitions with summands in arithmetic progression
AUSTRALASIAN JOURNAL OF COMBINATORICS Volume 8 (003), Pages 149 159 On the enumeration of partitions with summans in arithmetic progression M. A. Nyblom C. Evans Department of Mathematics an Statistics
More information. Using a multinomial model gives us the following equation for P d. , with respect to same length term sequences.
S 63 Lecture 8 2/2/26 Lecturer Lillian Lee Scribes Peter Babinski, Davi Lin Basic Language Moeling Approach I. Special ase of LM-base Approach a. Recap of Formulas an Terms b. Fixing θ? c. About that Multinomial
More informationTest of the validity of the gradient expansion for slowly-varying extended systems and coefficients of the sixth order for the KE in 1d
Test of the valiity of the graient expansion for slowly-varying extene systems an coefficients of the sixth orer for the KE in 1 Joseph B. Dizon Department of Chemistry, University of California, Irvine,
More information19 Eigenvalues, Eigenvectors, Ordinary Differential Equations, and Control
19 Eigenvalues, Eigenvectors, Orinary Differential Equations, an Control This section introuces eigenvalues an eigenvectors of a matrix, an iscusses the role of the eigenvalues in etermining the behavior
More informationNonlinear Adaptive Ship Course Tracking Control Based on Backstepping and Nussbaum Gain
Nonlinear Aaptive Ship Course Tracking Control Base on Backstepping an Nussbaum Gain Jialu Du, Chen Guo Abstract A nonlinear aaptive controller combining aaptive Backstepping algorithm with Nussbaum gain
More informationLevel Construction of Decision Trees in a Partition-based Framework for Classification
Level Construction of Decision Trees in a Partition-base Framework for Classification Y.Y. Yao, Y. Zhao an J.T. Yao Department of Computer Science, University of Regina Regina, Saskatchewan, Canaa S4S
More informationA Path Planning Method Using Cubic Spiral with Curvature Constraint
A Path Planning Metho Using Cubic Spiral with Curvature Constraint Tzu-Chen Liang an Jing-Sin Liu Institute of Information Science 0, Acaemia Sinica, Nankang, Taipei 5, Taiwan, R.O.C., Email: hartree@iis.sinica.eu.tw
More informationTime-of-Arrival Estimation in Non-Line-Of-Sight Environments
2 Conference on Information Sciences an Systems, The Johns Hopkins University, March 2, 2 Time-of-Arrival Estimation in Non-Line-Of-Sight Environments Sinan Gezici, Hisashi Kobayashi an H. Vincent Poor
More informationAn Approach for Design of Multi-element USBL Systems
An Approach for Design of Multi-element USBL Systems MIKHAIL ARKHIPOV Department of Postgrauate Stuies Technological University of the Mixteca Carretera a Acatlima Km. 2.5 Huajuapan e Leon Oaxaca 69000
More informationNew Simple Controller Tuning Rules for Integrating and Stable or Unstable First Order plus Dead-Time Processes
Proceeings of the 3th WSEAS nternational Conference on SYSTEMS New Simple Controller Tuning Rules for ntegrating an Stable or Unstable First Orer plus Dea-Time Processes.G.ARVANTS Department of Natural
More informationChapter 6: Energy-Momentum Tensors
49 Chapter 6: Energy-Momentum Tensors This chapter outlines the general theory of energy an momentum conservation in terms of energy-momentum tensors, then applies these ieas to the case of Bohm's moel.
More informationNonlinear Dielectric Response of Periodic Composite Materials
onlinear Dielectric Response of Perioic Composite aterials A.G. KOLPAKOV 3, Bl.95, 9 th ovember str., ovosibirsk, 639 Russia the corresponing author e-mail: agk@neic.nsk.su, algk@ngs.ru A. K.TAGATSEV Ceramics
More informationA SIMPLE ENGINEERING MODEL FOR SPRINKLER SPRAY INTERACTION WITH FIRE PRODUCTS
International Journal on Engineering Performance-Base Fire Coes, Volume 4, Number 3, p.95-3, A SIMPLE ENGINEERING MOEL FOR SPRINKLER SPRAY INTERACTION WITH FIRE PROCTS V. Novozhilov School of Mechanical
More informationSurvey Sampling. 1 Design-based Inference. Kosuke Imai Department of Politics, Princeton University. February 19, 2013
Survey Sampling Kosuke Imai Department of Politics, Princeton University February 19, 2013 Survey sampling is one of the most commonly use ata collection methos for social scientists. We begin by escribing
More informationLower Bounds for the Smoothed Number of Pareto optimal Solutions
Lower Bouns for the Smoothe Number of Pareto optimal Solutions Tobias Brunsch an Heiko Röglin Department of Computer Science, University of Bonn, Germany brunsch@cs.uni-bonn.e, heiko@roeglin.org Abstract.
More informationIPA Derivatives for Make-to-Stock Production-Inventory Systems With Backorders Under the (R,r) Policy
IPA Derivatives for Make-to-Stock Prouction-Inventory Systems With Backorers Uner the (Rr) Policy Yihong Fan a Benamin Melame b Yao Zhao c Yorai Wari Abstract This paper aresses Infinitesimal Perturbation
More information4. Important theorems in quantum mechanics
TFY4215 Kjemisk fysikk og kvantemekanikk - Tillegg 4 1 TILLEGG 4 4. Important theorems in quantum mechanics Before attacking three-imensional potentials in the next chapter, we shall in chapter 4 of this
More informationExperiment I Electric Force
Experiment I Electric Force Twenty-five hunre years ago, the Greek philosopher Thales foun that amber, the harene sap from a tree, attracte light objects when rubbe. Only twenty-four hunre years later,
More informationSchrödinger s equation.
Physics 342 Lecture 5 Schröinger s Equation Lecture 5 Physics 342 Quantum Mechanics I Wenesay, February 3r, 2010 Toay we iscuss Schröinger s equation an show that it supports the basic interpretation of
More informationQ( t) = T C T =! " t 3,t 2,t,1# Q( t) T = C T T T. Announcements. Bezier Curves and Splines. Review: Third Order Curves. Review: Cubic Examples
Bezier Curves an Splines December 1, 2015 Announcements PA4 ue one week from toay Submit your most fun test cases, too! Infinitely thin planes with parallel sies μ oesn t matter Term Paper ue one week
More informationAn Analytical Expression of the Probability of Error for Relaying with Decode-and-forward
An Analytical Expression of the Probability of Error for Relaying with Decoe-an-forwar Alexanre Graell i Amat an Ingmar Lan Department of Electronics, Institut TELECOM-TELECOM Bretagne, Brest, France Email:
More information3.7 Implicit Differentiation -- A Brief Introduction -- Student Notes
Fin these erivatives of these functions: y.7 Implicit Differentiation -- A Brief Introuction -- Stuent Notes tan y sin tan = sin y e = e = Write the inverses of these functions: y tan y sin How woul we
More informationOptimization of Geometries by Energy Minimization
Optimization of Geometries by Energy Minimization by Tracy P. Hamilton Department of Chemistry University of Alabama at Birmingham Birmingham, AL 3594-140 hamilton@uab.eu Copyright Tracy P. Hamilton, 1997.
More informationState-Space Model for a Multi-Machine System
State-Space Moel for a Multi-Machine System These notes parallel section.4 in the text. We are ealing with classically moele machines (IEEE Type.), constant impeance loas, an a network reuce to its internal
More informationUnified kinetic model of dopant segregation during vapor-phase growth
PHYSICAL REVIEW B 72, 195419 2005 Unifie kinetic moel of opant segregation uring vapor-phase growth Craig B. Arnol 1, * an Michael J. Aziz 2 1 Department of Mechanical an Aerospace Engineering an Princeton
More informationA New Backstepping Sliding Mode Guidance Law Considering Control Loop Dynamics
pp. 9-6 A New Backstepping liing Moe Guiance Law Consiering Control Loop Dynamics V. Behnamgol *, A. Vali an A. Mohammai 3, an 3. Department of Control Engineering, Malek Ashtar University of Technology
More informationθ x = f ( x,t) could be written as
9. Higher orer PDEs as systems of first-orer PDEs. Hyperbolic systems. For PDEs, as for ODEs, we may reuce the orer by efining new epenent variables. For example, in the case of the wave equation, (1)
More informationSlide10 Haykin Chapter 14: Neurodynamics (3rd Ed. Chapter 13)
Slie10 Haykin Chapter 14: Neuroynamics (3r E. Chapter 13) CPSC 636-600 Instructor: Yoonsuck Choe Spring 2012 Neural Networks with Temporal Behavior Inclusion of feeback gives temporal characteristics to
More informationThermal runaway during blocking
Thermal runaway uring blocking CES_stable CES ICES_stable ICES k 6.5 ma 13 6. 12 5.5 11 5. 1 4.5 9 4. 8 3.5 7 3. 6 2.5 5 2. 4 1.5 3 1. 2.5 1. 6 12 18 24 3 36 s Thermal runaway uring blocking Application
More informationPhysics 505 Electricity and Magnetism Fall 2003 Prof. G. Raithel. Problem Set 3. 2 (x x ) 2 + (y y ) 2 + (z + z ) 2
Physics 505 Electricity an Magnetism Fall 003 Prof. G. Raithel Problem Set 3 Problem.7 5 Points a): Green s function: Using cartesian coorinates x = (x, y, z), it is G(x, x ) = 1 (x x ) + (y y ) + (z z
More informationThe Principle of Least Action
Chapter 7. The Principle of Least Action 7.1 Force Methos vs. Energy Methos We have so far stuie two istinct ways of analyzing physics problems: force methos, basically consisting of the application of
More informationarxiv:cond-mat/ v1 12 Oct 2005
Creation an estruction of a spin gap in weakly couple quarter-fille laers B. Eegger,, H.G. Evertz, an R.M. Noack Institut für Theoretische Physik, Technische Universität Graz, A-8 Graz, Austria Institut
More informationON THE OPTIMALITY SYSTEM FOR A 1 D EULER FLOW PROBLEM
ON THE OPTIMALITY SYSTEM FOR A D EULER FLOW PROBLEM Eugene M. Cliff Matthias Heinkenschloss y Ajit R. Shenoy z Interisciplinary Center for Applie Mathematics Virginia Tech Blacksburg, Virginia 46 Abstract
More information1 Introuction In the past few years there has been renewe interest in the nerson impurity moel. This moel was originally propose by nerson [2], for a
Theory of the nerson impurity moel: The Schrieer{Wol transformation re{examine Stefan K. Kehrein 1 an nreas Mielke 2 Institut fur Theoretische Physik, uprecht{karls{universitat, D{69120 Heielberg, F..
More information23 Implicit differentiation
23 Implicit ifferentiation 23.1 Statement The equation y = x 2 + 3x + 1 expresses a relationship between the quantities x an y. If a value of x is given, then a corresponing value of y is etermine. For
More informationCalculus of Variations
16.323 Lecture 5 Calculus of Variations Calculus of Variations Most books cover this material well, but Kirk Chapter 4 oes a particularly nice job. x(t) x* x*+ αδx (1) x*- αδx (1) αδx (1) αδx (1) t f t
More informationA Modification of the Jarque-Bera Test. for Normality
Int. J. Contemp. Math. Sciences, Vol. 8, 01, no. 17, 84-85 HIKARI Lt, www.m-hikari.com http://x.oi.org/10.1988/ijcms.01.9106 A Moification of the Jarque-Bera Test for Normality Moawa El-Fallah Ab El-Salam
More information2Algebraic ONLINE PAGE PROOFS. foundations
Algebraic founations. Kick off with CAS. Algebraic skills.3 Pascal s triangle an binomial expansions.4 The binomial theorem.5 Sets of real numbers.6 Surs.7 Review . Kick off with CAS Playing lotto Using
More informationHomework 7 Due 18 November at 6:00 pm
Homework 7 Due 18 November at 6:00 pm 1. Maxwell s Equations Quasi-statics o a An air core, N turn, cylinrical solenoi of length an raius a, carries a current I Io cos t. a. Using Ampere s Law, etermine
More informationVectors in two dimensions
Vectors in two imensions Until now, we have been working in one imension only The main reason for this is to become familiar with the main physical ieas like Newton s secon law, without the aitional complication
More informationChapter-2. Steady Stokes flow around deformed sphere. class of oblate axi-symmetric bodies
hapter- Steay Stoes flow aroun eforme sphere. class of oblate axi-symmetric boies. General In physical an biological sciences, an in engineering, there is a wie range of problems of interest lie seimentation
More informationA Note on Exact Solutions to Linear Differential Equations by the Matrix Exponential
Avances in Applie Mathematics an Mechanics Av. Appl. Math. Mech. Vol. 1 No. 4 pp. 573-580 DOI: 10.4208/aamm.09-m0946 August 2009 A Note on Exact Solutions to Linear Differential Equations by the Matrix
More informationA Quantitative Analysis of Coupling for a WPT System Including Dielectric/Magnetic Materials
Progress In Electromagnetics Research Letters, Vol. 72, 127 134, 2018 A Quantitative Analysis of Coupling for a WPT System Incluing Dielectric/Magnetic Materials Yangjun Zhang *, Tatsuya Yoshiawa, an Taahiro
More informationMath Notes on differentials, the Chain Rule, gradients, directional derivative, and normal vectors
Math 18.02 Notes on ifferentials, the Chain Rule, graients, irectional erivative, an normal vectors Tangent plane an linear approximation We efine the partial erivatives of f( xy, ) as follows: f f( x+
More informationOn the Aloha throughput-fairness tradeoff
On the Aloha throughput-fairness traeoff 1 Nan Xie, Member, IEEE, an Steven Weber, Senior Member, IEEE Abstract arxiv:1605.01557v1 [cs.it] 5 May 2016 A well-known inner boun of the stability region of
More informationSystems & Control Letters
Systems & ontrol Letters ( ) ontents lists available at ScienceDirect Systems & ontrol Letters journal homepage: www.elsevier.com/locate/sysconle A converse to the eterministic separation principle Jochen
More informationd dx But have you ever seen a derivation of these results? We ll prove the first result below. cos h 1
Lecture 5 Some ifferentiation rules Trigonometric functions (Relevant section from Stewart, Seventh Eition: Section 3.3) You all know that sin = cos cos = sin. () But have you ever seen a erivation of
More informationOptimum design of tuned mass damper systems for seismic structures
Earthquake Resistant Engineering Structures VII 175 Optimum esign of tune mass amper systems for seismic structures I. Abulsalam, M. Al-Janabi & M. G. Al-Taweel Department of Civil Engineering, Faculty
More information24th European Photovoltaic Solar Energy Conference, September 2009, Hamburg, Germany
4th European hotovoltaic Solar Energy Conference, 1-5 September 9, Hamburg, Germany LOCK-IN THERMOGRAHY ON CRYSTALLINE SILICON ON GLASS (CSG) THIN FILM MODULES: INFLUENCE OF ELTIER CONTRIBUTIONS H. Straube,
More informationThe effect of dissipation on solutions of the complex KdV equation
Mathematics an Computers in Simulation 69 (25) 589 599 The effect of issipation on solutions of the complex KV equation Jiahong Wu a,, Juan-Ming Yuan a,b a Department of Mathematics, Oklahoma State University,
More informationInfluence of Radiation on Product Yields in a Film Boiling Reactor
R&D NOTES Influence of Raiation on Prouct Yiels in a Film Boiling Reactor C. Thomas Aveisian, Wing Tsang, Terence Daviovits, an Jonah R. Allaben Sibley School of Mechanical an Aerospace Engineering, Cornell
More informationTime-Optimal Motion Control of Piezoelectric Actuator: STM Application
Time-Optimal Motion Control of Piezoelectric Actuator: STM Application Yongai Xu, Peter H. Mecl Abstract This paper exaes the problem of time-optimal motion control in the context of Scanning Tunneling
More informationParameter estimation: A new approach to weighting a priori information
Parameter estimation: A new approach to weighting a priori information J.L. Mea Department of Mathematics, Boise State University, Boise, ID 83725-555 E-mail: jmea@boisestate.eu Abstract. We propose a
More informationTo understand how scrubbers work, we must first define some terms.
SRUBBERS FOR PARTIE OETION Backgroun To unerstan how scrubbers work, we must first efine some terms. Single roplet efficiency, η, is similar to single fiber efficiency. It is the fraction of particles
More informationInternational Conference KNOWLEDGE-BASED ORGANIZATION Vol. XXIII No
International Conference KNOWLEDGE-BAED ORGANIZATION Vol. XXIII No 3 2017 METHOD FOR DETERMINATION OF THE PARAMETER OF A MACHINE GUN UPENION MOUNTED ON AN ARMOURED VEHICLE vilen PIRDONOV, vilen TEFANOV
More informationHow to Minimize Maximum Regret in Repeated Decision-Making
How to Minimize Maximum Regret in Repeate Decision-Making Karl H. Schlag July 3 2003 Economics Department, European University Institute, Via ella Piazzuola 43, 033 Florence, Italy, Tel: 0039-0-4689, email:
More information28.1 Parametric Yield Estimation Considering Leakage Variability
8.1 Parametric Yiel Estimation Consiering Leakage Variability Rajeev R. Rao, Aniruh Devgan*, Davi Blaauw, Dennis Sylvester University of Michigan, Ann Arbor, MI, *IBM Corporation, Austin, TX {rrrao, blaauw,
More informationA new identification method of the supply hole discharge coefficient of gas bearings
Tribology an Design 95 A new ientification metho of the supply hole ischarge coefficient of gas bearings G. Belforte, F. Colombo, T. Raparelli, A. Trivella & V. Viktorov Department of Mechanics, Politecnico
More informationImproved Rate-Based Pull and Push Strategies in Large Distributed Networks
Improve Rate-Base Pull an Push Strategies in Large Distribute Networks Wouter Minnebo an Benny Van Hout Department of Mathematics an Computer Science University of Antwerp - imins Mielheimlaan, B-00 Antwerp,
More informationLagrangian and Hamiltonian Mechanics
Lagrangian an Hamiltonian Mechanics.G. Simpson, Ph.. epartment of Physical Sciences an Engineering Prince George s Community College ecember 5, 007 Introuction In this course we have been stuying classical
More informationHarmonic Modelling of Thyristor Bridges using a Simplified Time Domain Method
1 Harmonic Moelling of Thyristor Briges using a Simplifie Time Domain Metho P. W. Lehn, Senior Member IEEE, an G. Ebner Abstract The paper presents time omain methos for harmonic analysis of a 6-pulse
More informationSharp Thresholds. Zachary Hamaker. March 15, 2010
Sharp Threshols Zachary Hamaker March 15, 2010 Abstract The Kolmogorov Zero-One law states that for tail events on infinite-imensional probability spaces, the probability must be either zero or one. Behavior
More informationEfficient Macro-Micro Scale Coupled Modeling of Batteries
A00 Journal of The Electrochemical Society, 15 10 A00-A008 005 0013-651/005/1510/A00/7/$7.00 The Electrochemical Society, Inc. Efficient Macro-Micro Scale Couple Moeling of Batteries Venkat. Subramanian,*,z
More informationOn combinatorial approaches to compressed sensing
On combinatorial approaches to compresse sensing Abolreza Abolhosseini Moghaam an Hayer Raha Department of Electrical an Computer Engineering, Michigan State University, East Lansing, MI, U.S. Emails:{abolhos,raha}@msu.eu
More informationarxiv: v1 [math.co] 29 May 2009
arxiv:0905.4913v1 [math.co] 29 May 2009 simple Havel-Hakimi type algorithm to realize graphical egree sequences of irecte graphs Péter L. Erős an István Miklós. Rényi Institute of Mathematics, Hungarian
More informationLight-Soaking Effects on the Open-Circuit Voltage of a-si:h Solar Cells
Syracuse University SURFACE Physics College of Arts an Sciences 2005 Light-Soaking Effects on the Open-Circuit Voltage of a-si:h Solar Cells Jianjun Liang Syracuse University Eric A. Schiff Syracuse University
More information