Scaling Multiple Point Statistics for Non-Stationary Geostatistical Modeling
|
|
- Kory Ford
- 5 years ago
- Views:
Transcription
1 Scaling Multile Point Statistics or Non-Stationary Geostatistical Modeling Julián M. Ortiz, Steven Lyster and Clayton V. Deutsch Centre or Comutational Geostatistics Deartment o Civil & Environmental Engineering University o Alberta Multile-oint statistics are used in geostatistical simulation to imrove orecasting o resonses that are highly deendent on the reroduction o comlex eatures o the henomenon that cannot be catured by conventional two-oint simulation methods. Inerence o multile-oint statistics is oten based on a training image that deicts the eatures that rovide the character to the geological event being modelled. One limitation o this aroach is that the univariate distribution o categories (acies or rock tyes) in the training image may not match the statistics o the inal model. The question o scaling multile-oint statistics arises, the idea being to take the statistics rom the training image and scale them in a reeatable manner and honouring the univariate roortions o categories. An iterative scaling aroach based on the exression or scaling multile-oint statistics in a urely random case is roosed. A multi-gaussian alternative is discussed, which shows some roblems due to the symmetry o the multivariate Gaussian distribution. The imlementation is illustrated through an examle where it is shown that the roosed method lies in between two extreme cases or a Boolean simulation, namely, the change in size o the obects and the change in their number o occurrences. A second examle is resented to illustrate the otential use o this scaling rocedure or non-stationary multile-oint geostatistical simulation. Introduction Multile-oint statistics can be used in an attemt to imrove geostatistical modelling o variables distributed in sace. The relationshis between several oints at a time are estimated rom training data and imosed during the simulation rocess, in order to achieve a numerical model that correctly reresents these intricate relationshis that conventional simulation cannot cature. Figure shows a reerence image where the relationshis between the dierent geological units or acies cannot be correctly catured by conventional simulation methods. Exhaustive images such as the one deicted here can be used as an analogue to the henomenon that is being modelled; multile-oint statistics can be extracted by scanning the image and comuting the requency with which acies arranged in seciic atterns occur. The training inormation requires abundant data located over a regular grid o oints, in order to have enough relicates o each articular multile-oint event. This is oten solved by utilizing a training image (Guardiano and Srivastava, 992; Deutsch, 992) or by using available seudoregularly saced roduction data, as is the case in mining alications (Ortiz, 2003; Ortiz and Deutsch, 2004; Ortiz and Emery, 2005). Multile-oint geostatistical simulation can be erormed using any o the available methods. The single normal equation simulation roosed by Strebelle and Journel (2000) estimates the -
2 conditional distribution at every location given a multile-oint coniguration by calculating the requency with which the indicator at the location being simulated is one given that the multileoint event occurs in the training image. Alternatively, simulated annealing (Deutsch, 992) can be used to match the multile-oint requencies extracted rom a training image into a simulated numerical model. Again, the multile-oint statistics are read rom the training image as requencies o articular events occurring. Other methods such as neural networks also rely on the use o a training image to extract and reroduce the multile-oint statistics (Caers and Journel, 998). There is an imlicit assumtion o stationarity during the exorting o multile-oint statistics rom the training image to the simulated model. A stationarity assumtion stricter than the usual second order stationarity is required, since the use o conigurations o several oints also locks all lower order statistics. For instance, i our-oint conigurations are matched during the multile-oint simulation rocess, then all three-oint and two-oint statistics whose coniguration is included in the our-oint coniguration originally used are imlicitly matched ( Figure 2). Most imortantly, the histogram (one-oint statistic) is also locked when higher order statistics are honoured. This aer addresses the issue o scaling multile-oint statistics, which may be required in two distinct cases. First, consistency roblems may arise when the training image does not have exactly the same one-oint and two-oint distribution as the underlying henomenon. This consistency roblem occurs requently, since training images usually are not available with the exact same histogram o the variable that is being modelled (roortions o acies). Direct use o the multile-oint statistics extracted rom the training image will distort the histogram o the simulated realizations generating a bias in the roortions. A basic requirement or the simulated model to be acceted as a lausible reresentation o the true henomenon is that it honours a given histogram. Scaling o multile-oint statistics is then required to allow the use o the training image in a consistent manner with the data and reserving its character, hence generating accetable numerical models. Another relevant use o a scaling rocedure or multile-oint statistics is to imose nonstationary eatures, but reserving the character rovided by the training image. One could devise the use o a single training image to model a ield with locally varying roortions o the acies. We resent a methodology or scaling multile-oint statistics to generate consistent results rom simulation methods that account or this inormation and to consider non-stationary eatures during the modeling rocess. The methodology is general and could be adated to be used with continuous variables, as long as the data are coded as indicators by deining classes through a set o thresholds. In what ollows, the roblem is resented in the case o a categorical variable, which is where simulation accounting or multile-oint statistics has seen a aster develoment. Problem setting Consider a training image that deicts the satial arrangement o K categories or acies. The global roortions with which these categories are resent in the training image are denoted: k, k =,..., K. The general aearance o the training image is deemed aroriate to model a given geological setting, hence the modeler decides the training image is to be used or inerence o the multileoint statistics that a given simulation algorithm will imose to a set o realizations. The multile-oint statistics to be considered are deined by a satial arrangement o nodes and by the combination o acies values in these nodes. I we consider the case, where an N-oints -2
3 statistic is considered, then K N ossible combinations are available. Each o these combinations occurs in the training image with a given requency. In act, as soon as N or K become relatively large, there is big chance that many o the K N combinations do not occur in the training image. Each one o the ossible combinations is identiied with an index that comletely deines the acies values within the N oints, however the ordering to identiy the oints in the attern must be deined rior to the calculation o the index o each multile-oint event (Figure 3). The index o each multile-oint coniguration is calculated as: = + N n= ( i n ) K where i n is the code o the n th node o the attern that identiies its acies. Here, it is assumed that acies are denoted by consecutive integers starting with code. N The requency o each multile-oint event in the training image is denoted:, =,..., K. Knowing the indexing and the requencies with which each multile-oint coniguration occurs, ermits calculation o the acies roortions k, k =,..., K.. Denoting by,..., N k, ; k =,..., K; = K, the roortion o acies k in the multile-oint arrangement identiied with the index, the roortions o the acies in the training image can be retrieved rom the multile-oint statistics as ollows. Consider the multile-oint index o interest is. We can calculate the value o the n th node, by taking: = i N int N K where int reresents the integer art o the division. Subsequently, the indexes o the n remaining nodes o the attern can be calculated recursively using the residual o the division (ractional art): n rac n K i = n int n K n =,..., N Knowing the acies values o the N nodes o each multile-oint index and their requencies, the univariate roortions can be easily calculated. Now, the simulated model must honor a set o statistics inerred rom a set o samles. For instance, the available data may show that the global declustered roortions o the acies in the domain are: k, k =,..., K, with k not necessarily equal to k or some k =,..., K. The goal o this aer is to roose a methodology to calculate the corrected requencies o multile-oint events that will honor the roortions, reserving the character o the training image, that is, keeing the eatures that make it distinct. These corrected requencies are denoted: *, =,..., K N. Scaling aroaches or multile-oint statistics To understand the concet o scaling multile-oint statistics, the ollowing examle illustrates ossible outcomes rom an increase in the roortion o a acies in a binary case and where the -3
4 scaled models reserve the character o the original training image. A ield o 000 by 000 ixels is oulated with 0 by 0 ixels squares, where the centre o the squares are randomly located in the ield, that is they are generated through a Poisson rocess ( Figure 4). There are enough squares to cover 20% o the domain, leaving the remaining 80% as background. Scaling this training image so that the roortion o squares goes u to 40% may generate two equally valid outcomes: we can have more 0 by 0 ixels squares or have a smaller number o larger squares, say 20 by 20 ixels squares. These situations are o course extremes and since we do not know the exact multile-oint statistics o the variable and we borrow this inormation rom a training image that does not exactly match the roortions, a scaling rocedure that lays somewhere in between the two extremes resented above can be used as a modelling decision to handle the roblem. Two aroaches could be devised to scale multile-oint statistics when the univariate distribution changes: Scaling based on how a variable randomly distributed changes its multile-oint statistics when the univariate roortions o the acies change. Scaling based on how a categorical variable obtained by multile truncations o a correlated multigaussian continuous variable changes as the thresholds that deine the acies are modiied. Dimensionality is always a roblem when dealing with multile-oint statistics since the combinatorial becomes raidly very large as the number o oints in the multile-oint coniguration N or the number o acies K increase. For examle, 0 acies can be arranged on a 9 oints attern in billion ossible combinations. O course, this examle is extreme, but attern size oten increases in an exonential ashion: 4, 9, 6, 25 oints in 2D, or 8, 27, 64, 25 oints in 3D. This roblem may be artially solved by roceeding airwise, that is, considering the most relevant acies and all other acies together. Then, reezing the most imortant acies already simulated, one can simulate another relevant acies against all remaining acies on a reduced domain, and so on, until all acies have been individually taken into account. Scaling MPS based on the random case The irst aroach roosed is based on how multile-oint statistics change when the acies codes o the oints in the attern are randomly distributed. In this case, the requency with which a multile-oint event occurs can be calculated as the roduct o the robabilities o occurrence (roortions) o each acies values in its nodes. Since the acies values are considered uncorrelated, it is straightorward to scale the requency o occurrence o multile-oint statistics. The requency o a multile-oint coniguration can be calculated by: = K k = N k, k Since k, is the roortion o acies k in class, N k, is the number o occurrences o acies k in class. For examle (Figure 5), considering a case where two acies are available and the robability o acies revailing at a given location is 0.25, then the robability o having a ouroints coniguration where acies revails in two nodes and acies 2 revails in the remaining nodes would be calculated as: = =
5 Scaling the random case is quite simle; multilying the original requency by a series roduct o the desired requency over the old requency results in a simle equation: * = K k = k It is easy to see that the scaled multile-oint requencies honour the univariate roortions: * = K k = k k N k, ( K k k= = K N k, ) ( k ) k= k N N k, k, * = = Following u on the revious examle, we can see that a change in the global roortions o the acies will change the robability o occurrence o the multile-oint event deicted in Figure 5. Since the uncorrelated case is simle, we can easily calculate this robability or the new global roortions. Consider the case, the robability o acies revailing at a location goes u to 0.4, leaving a robability o 0.6 or acies 2. Considering the same our-oints coniguration, we can calculate its new robability o occurrence: new = = It can be checked that the exression rovided or scaling multile-oint statistics in the random case rovides the same result calculated above: * 3 = 3 = N,3 N , = This exression is valid only in the case o a multile-oint statistics rom an uncorrelated variable (ure nugget eect). As soon as the variable is satially correlated, multile-oint roortions cannot be directly calculated. An iterative aroach is roosed next to calculate the multile-oint requencies that honour the roortions, rom the initial multile-oint requencies inerred rom the training image. The ormula above overcorrects the multile-oint requencies. Convergence can be achieved by iterating using the ollowing modiied exression: * * = K k = k k k, N Multile-oint requencies or all indexes =,..., K must be udated. The new global roortions k, k =,..., K are recalculated or each iteration. This ormula does not require a large number o iterations, and usually, the desired multile-oint requencies can be obtained with a nearly erect match o the roortions, with less than 50 iterations, which takes only a ew seconds. -5
6 Scaling MPS based on the multigaussian case An alternative method consists on considering how a categorical variable deined by multile truncations o a multigaussian ield changes as the thresholds or truncations change. The idea can be understood quite simly, by considering a two-oint attern ( Figure 6). A generalization to N oints can be achieved, associating a multivariate Gaussian distribution o dimension N. In this simle examle, two oints searated by a vector h constitute the multile-oint attern. Initially, we can consider a binary variable, hence it can be deined by a single truncation o the Gaussian variable with a threshold y C. The two categories (acies) have robabilities: = G( y ) 2 C = C = G( y An indicator variogram can be calculated and associated with a variogram or the underlying multigaussian variable (Kyriakidis and others, 999). Given the satial arrangement o the two oints, Their satial correlation can be calculated rom the variogram model and assuming multivariate Gaussianity, all conditional robabilities can be retrieved. The robabilities o all our ossible multile-oint events can be calculated, by using conditional robabilities that only deend on two-oint statistics, that is, the variogram model. For examle, assuming acies is associated with the continuous Gaussian variable when the value is below the threshold y C : Prob ( yc i = ; i2 = 2) = g( y, y2 ) dy2 dy y where g(y,y 2 ) is the bivariate Gaussian robability density unction. This simle case can be easily extended to multile oints (N > 2), however, calculating the oint robabilities becomes cumbersome as the dimensionality increases. A simle numerical aroximation can be used to calculate the multile-oint event requencies. In the multile-oint case, all correlations between the nodes in the attern must be taken into account in order to aroximate the integration o the multivariate Gaussian distribution using the aroriate thresholds to deine the acies. A Monte-Carlo aroach to simulate multivariate Gaussian vectors using the matrix decomosition method is roosed (Davis, 987). The general rocedure is:. Given the training image and roortions k, k =,..., K and k, k =,..., K, calculate the thresholds that give those robabilities when truncating a standard Gaussian distribution: y y C, N = G M = G N i= i C C ) C, N = G M C, ( ) C, ( ) 2. Comute the indicator variograms and iner the variogram o the Gaussian variable that better its the calculated indicator variograms. 3. Given the searation vectors between the N oints in the multile-oint attern, build the N by N matrix o covariances between oints in the attern ( y y = G N i= i -6
7 4. Figure 7): C( h) = C( h 2) C M C( h N) C( h C( h C( h 5. Using the matrix decomosition simulation algorithm generate a large number o simulated multivariate vectors that honor the covariance between the oints. The algorithm roceeds as ollows: a. Decomose the covariance matrix as a roduct o a lower (L) and an uer (U) triangular matrix such that L = U T (Cholesky decomosition). b. Draw a vector o indeendent standard Gaussian deviates w. c. Generate a vector o correlated Gaussian values by multilying the lower matrix by the vector o random normal deviates: y = L w. 6. Identiy the multile-oint index o each simulated vector by alying the truncation rule deined by the thresholds, that is transorm the vectors o correlated Gaussian values to vectors o acies values, by coding the values as acies using the thresholds. Then, calculate the multile-oint index corresonding to each coniguration o acies. 7. For each multile-oint coniguration that existed in the training image, calculate the roortion o simulated vectors that match that multile-oint index when alying the thresholds calculated or the training image roortions and when the roortions are used. 8. Correct the multile-oint statistics o the training image, using a ratio o those roortions. 9. Standardize to ensure the sum o roortions rom the scaled multile-oint requencies add to one. This rocedure accounts or the correlation between the oints. It is an aroximation due to the act that the Gaussian variable is continuous, hence all multile-oint combination o indexes in the attern has a robability o occurrence, which is not the case in the training image, since this is a discrete reresentation. A correction is required to match exactly the roortions, but this cannot always be achieved, since using the multivariate Gaussian case imoses constraints with resect to the occurrence o symmetric events: or instance the resence o acies one on the right o acies two imlies that acies two must exist also on the right o acies one. For this reason, this aroach is generally not recommended. Examles First, to illustrate the use o these scaling aroaches, the small examle resented in Figure 4 is ollowed u. The multile-oint histogram or a two by two oints attern is calculated. The requency o occurrence o all 6 multile-oint events is comuted and a histogram o the requency with which each indexed event occurs is lotted as a summary. The 6 ossible events can be easily interreted in terms o the geometry o the obects (squares), as shown in Figure 8. The reerence image (ma on the let hand side o M 2 22 N2 ) ) ) L L O L C( h C( h M C( h N 2N NN ) ) ) -7
8 Figure 4) is used to calculate the multile-oint requencies in a 2 by 2 ixels attern and these are scaled to reach a global roortion o squares o These scaled statistics are then comared with the two extreme cases resented on the right hand side o Figure 4. Figure 9 shows the three multile-oint histograms suerimosed. It can be seen that the scaling using the random aroach rovides statistics that lie in between the two extreme cases. We can exect that these scaled u statistics reresent a case where there are more squares in the domain, but these are also slightly larger than the ones in the training image. The relevant eatures o the reerence image are then catured in the statistics. A second examle has been designed to illustrate the use o the scaling rocedure roosed. A binary two dimensional training image has been created with black horizontal stries two units thick in a white background. The training image is shown in Figure 0. The original roortions are black = 0.25 and white = Multile-oint statistics have been inerred rom the training image and unconditional realizations have been comuted or dierent roortions o black (stries) and white (background) acies. A simulated annealing rogram to match the multile-oint statistics was reared. Dierent attern sizes imly control over dierent scales o the henomenon. For examle, Figure shows two realizations considering a 4 by 4 nodes attern. Notice that with this attern size the general aearance o the resulting models is close to the training image. I a dierent attern is considered, the long range eatures may not be roerly catured. Figure 2 shows realizations made with dierent attern sizes to illustrate this eect. The scaling rocedure was used to scale u the roortion o the black acies. The initial roortions (25% black, 75% white) was changed and the multile-oint statistics were scaled. These were then used to simulated realizations with the intent o generating realizations that ket the essential eatures o the training image, but honouring the acies roortions imosed to the new models. Figure 3 shows three realizations constructed using the multile-oint statistics or a 3 by 3 nodes attern, scaled to match the ollowing acies roortions: 30% black and 70% white. These roortions were then changed to 50% or the black acies and 50% or the white background, and realizations with statistics rom 2 by 2, 3 by 3, and 4 by 4 nodes atterns considered. These results are shown in Figure 4. It can be seen rom the realizations dislayed above that, deending on the attern size, two situations can arise: the stries may become wider i the attern is small enough not to cature its entire width, or the stries may become more abundant i the attern is big enough to cature the relationshis between the width o the black acies obects and the background, as is the case when 3 by 3 or 4 by 4 atterns are used. Conclusions Scaling multile-oint statistics can imrove the current use o multile-oint geostatistical simulation techniques, by allowing locally varying roortions o acies to be reroduced still honouring the relationshis catured by multile-oint statistics. Another imortant use o a multile-oint scaling technique is the use o a reresentative training image that does not exactly match the simulation roortions. This article resents an aroach to scale multile-oint statistics by iteratively changing them in a manner similar to the change that occurs in multile-oint statistics o a urely random ield where the roortions o the acies are modiied. This aroach can be shown to rovide -8
9 reasonable values o the statistics that lie in between two extreme cases that are easy to show in the case o obects: an increase in the number o obects or an increase in the size o the obects. The scaling rocedure has the otential o being incororated in a non-stationary simulation algorithm that honors the eatures rovided by a training image, through the local udating o multile-oint statistics. Acknowledgements The authors would like to acknowledge the industry sonsors o the Centre or Comutational Geostatistics at the University o Alberta and the Chair in Ore Reserve Evaluation at the University o Chile sonsored by Codelco Chile or suorting this research. Reerences Caers, J. and Journel, A. G., Stochastic reservoir simulation using neural networks trained on outcro data, in 998 SPE Annual Technical Conerence and Exhibition, New Orleans, LA. Society o Petroleum Engineers. SPE aer # 49026, 998, Davis, M. W., Production o conditional simulations via the LU triangular decomosition o the covariance matrix. Mathematical Geology, vol. 9, no. 2, 987, Deutsch, C. V., Annealing Techniques Alied to Reservoir Modeling and the Integration o Geological and Engineering (Well Test) Data. PhD thesis, Stanord University, Stanord, CA, 992. Guardiano, F. and Srivastava, M., Multivariate geostatistics: Beyond bivariate moments, in A. Soares, editor, Geostatistics Tróia 92. Kluwer, Dordrecht, vol., 993, Kyriakidis, P. C., Deutsch C. V., and Grant M. L., Calculation o the normal scores variogram used or truncated Gaussian lithoacies simulation: theory and FORTRAN code. Comuters & Geosciences 25 (2): 6-69, March 999. Ortiz, J. M., Characterization o High Order Correlation or Enhanced Indicator Simulation. PhD thesis, University o Alberta, Edmonton, AB, Canada, Ortiz, J. M. and Deutsch, C. V., Indicator simulation accounting or multile-oint statistics. Mathematical Geology, vol. 36, no. 6, 2004, Ortiz, J. M. and Emery, X., Integrating multile-oint statistics into sequential simulation algorithms, in O. Leuangthong and C. V. Deutsch, editors, Geostatistics Ban Kluwer Academic, Dordrecht, The Netherlands, Vol. 2, 2005, Strebelle, S. and Journel, A. G., Sequential simulation drawing structures rom training images, in 6th International Geostatistics Congress, Cae Town, South Arica, Geostatistical Association o Southern Arica,
10 Figure : an exhaustive image showing intrincated relationshis o our acies in multile-oint atterns that cannot be easily catured by conventional simulation techniques. Figure 2: A our-oints attern and all lower order conigurations that are imlicitly matched by honouring the our-oint statistics: three-oints statistics in our conigurations; two-oints statistics or the corresonding lag distances searating the centres o the nodes (these reresent indicator variogram values); and the one-oint statistics that corresonds to the histogram. Nodes are reresented by their surrounding squares or illustration. -0
11 Figure 3: A our-oint coniguration, the order or considering the nodes and codes o the acies. Three examles o calculation o the multile-oints indexes are illustrated. Figure 4: Two ossible uscaled models rom the training image on the let. The to right ma shows more squares o the same dimension as in the reerence image; the bottom right ma shows larger squares. Both mas have a roortion o ixels belonging to squares o
12 Figure 5: examle o our-oints coniguration with two acies to calculate the robability o occurrence o an event in the random case. Figure 6: An illustration o the relationshi between global roortions, deined by the thresholds that truncate the multigaussian distribution, and the multile-oint requencies. A two-oints attern o nodes searated by a lag vector h allows calculation o the correlation coeicient rom the normal scores variogram γ(h). This deines the bivariate Gaussian distribution that is truncated so that the robability o being below y C is and comletely deines the multile-oint robabilities. These robabilities change in a tractable manner when the global roortions are changed by modiying the threshold. -2
13 Figure 7: Four-oints coniguration and the corresonding lag searation vectors used to build the covariance matrix to simulate multivariate Gaussian vectors. Figure 8: Geometric interretation o the dierent multile-oint conigurations or the examle o drawing squares randomly in the domain. -3
14 Frequency Multile-Point Index scaled 0 x 0 20 x 20 Figure 9: Multile-oint histogram or the scaled statistics and the two extreme cases: increased number o 0 by 0 squares and increased size o the squares to 20 by 20. A logarithmic scale has been used or the requency axis to better show that he scaled statistics lie in between the requencies or the extreme cases. Figure 0: Reerence training image. Black stries are two nodes wide and searated by 6 nodes in the north direction. -4
15 Figure : Two realizations that honor the original roortions, made considering a 4 by 4 nodes attern. Figure 2: Two realizations that honor the original roortions, made considering statistics extracted rom a 2 by 2 (let) and a 3 by 3 nodes attern (right). -5
16 Figure 3: Three realizations considering the statistics rom a 3 by 3 nodes attern and scaled roortion o 30% or the black acies and 70% or the white background. Figure 4: Three realizations with multile-oint statistics scaled to 50% black acies and 50% white background, with increasing attern size: 2 by 2 (let), 3 by 3 (middle), and 4 by 4 (right). -6
Combining Logistic Regression with Kriging for Mapping the Risk of Occurrence of Unexploded Ordnance (UXO)
Combining Logistic Regression with Kriging for Maing the Risk of Occurrence of Unexloded Ordnance (UXO) H. Saito (), P. Goovaerts (), S. A. McKenna (2) Environmental and Water Resources Engineering, Deartment
More informationDeriving Indicator Direct and Cross Variograms from a Normal Scores Variogram Model (bigaus-full) David F. Machuca Mory and Clayton V.
Deriving ndicator Direct and Cross Variograms from a Normal Scores Variogram Model (bigaus-full) David F. Machuca Mory and Clayton V. Deutsch Centre for Comutational Geostatistics Deartment of Civil &
More informationA MultiGaussian Approach to Assess Block Grade Uncertainty
A MultiGaussian Approach to Assess Block Grade Uncertainty Julián M. Ortiz 1, Oy Leuangthong 2, and Clayton V. Deutsch 2 1 Department of Mining Engineering, University of Chile 2 Department of Civil &
More informationANALYSIS OF ULTRA LOW CYCLE FATIGUE PROBLEMS WITH THE BARCELONA PLASTIC DAMAGE MODEL
XII International Conerence on Comutational Plasticity. Fundamentals and Alications COMPLAS XII E. Oñate, D.R.J. Owen, D. Peric and B. Suárez (Eds) ANALYSIS OF ULTRA LOW CYCLE FATIGUE PROBLEMS WITH THE
More informationA BSS-BASED APPROACH FOR LOCALIZATION OF SIMULTANEOUS SPEAKERS IN REVERBERANT CONDITIONS
A BSS-BASED AROACH FOR LOCALIZATION OF SIMULTANEOUS SEAKERS IN REVERBERANT CONDITIONS Hamid Reza Abutalebi,2, Hedieh Heli, Danil Korchagin 2, and Hervé Bourlard 2 Seech rocessing Research Lab (SRL), Elec.
More informationState Estimation with ARMarkov Models
Deartment of Mechanical and Aerosace Engineering Technical Reort No. 3046, October 1998. Princeton University, Princeton, NJ. State Estimation with ARMarkov Models Ryoung K. Lim 1 Columbia University,
More informationHotelling s Two- Sample T 2
Chater 600 Hotelling s Two- Samle T Introduction This module calculates ower for the Hotelling s two-grou, T-squared (T) test statistic. Hotelling s T is an extension of the univariate two-samle t-test
More informationSolved Problems. (a) (b) (c) Figure P4.1 Simple Classification Problems First we draw a line between each set of dark and light data points.
Solved Problems Solved Problems P Solve the three simle classification roblems shown in Figure P by drawing a decision boundary Find weight and bias values that result in single-neuron ercetrons with the
More information8 STOCHASTIC PROCESSES
8 STOCHASTIC PROCESSES The word stochastic is derived from the Greek στoχαστικoς, meaning to aim at a target. Stochastic rocesses involve state which changes in a random way. A Markov rocess is a articular
More informationEstimation of the large covariance matrix with two-step monotone missing data
Estimation of the large covariance matrix with two-ste monotone missing data Masashi Hyodo, Nobumichi Shutoh 2, Takashi Seo, and Tatjana Pavlenko 3 Deartment of Mathematical Information Science, Tokyo
More informationUsing a Computational Intelligence Hybrid Approach to Recognize the Faults of Variance Shifts for a Manufacturing Process
Journal of Industrial and Intelligent Information Vol. 4, No. 2, March 26 Using a Comutational Intelligence Hybrid Aroach to Recognize the Faults of Variance hifts for a Manufacturing Process Yuehjen E.
More informationCheng, N. S., and Law, A. W. K. (2003). Exponential formula for computing effective viscosity. Powder Technology. 129(1-3),
THIS PAPER SHOULD BE CITED AS Cheng, N. S., and Law, A. W. K. (2003). Exonential ormula or comuting eective viscosity. Powder Technology. 129(1-3), 156 160. EXPONENTIAL FORMULA FOR COMPUTING EFFECTIVE
More informationCHAPTER-II Control Charts for Fraction Nonconforming using m-of-m Runs Rules
CHAPTER-II Control Charts for Fraction Nonconforming using m-of-m Runs Rules. Introduction: The is widely used in industry to monitor the number of fraction nonconforming units. A nonconforming unit is
More informationarxiv: v1 [physics.data-an] 26 Oct 2012
Constraints on Yield Parameters in Extended Maximum Likelihood Fits Till Moritz Karbach a, Maximilian Schlu b a TU Dortmund, Germany, moritz.karbach@cern.ch b TU Dortmund, Germany, maximilian.schlu@cern.ch
More informationLower Confidence Bound for Process-Yield Index S pk with Autocorrelated Process Data
Quality Technology & Quantitative Management Vol. 1, No.,. 51-65, 15 QTQM IAQM 15 Lower onfidence Bound for Process-Yield Index with Autocorrelated Process Data Fu-Kwun Wang * and Yeneneh Tamirat Deartment
More information4. Score normalization technical details We now discuss the technical details of the score normalization method.
SMT SCORING SYSTEM This document describes the scoring system for the Stanford Math Tournament We begin by giving an overview of the changes to scoring and a non-technical descrition of the scoring rules
More informationAn Improved Generalized Estimation Procedure of Current Population Mean in Two-Occasion Successive Sampling
Journal of Modern Alied Statistical Methods Volume 15 Issue Article 14 11-1-016 An Imroved Generalized Estimation Procedure of Current Poulation Mean in Two-Occasion Successive Samling G. N. Singh Indian
More informationLogistics Optimization Using Hybrid Metaheuristic Approach under Very Realistic Conditions
17 th Euroean Symosium on Comuter Aided Process Engineering ESCAPE17 V. Plesu and P.S. Agachi (Editors) 2007 Elsevier B.V. All rights reserved. 1 Logistics Otimization Using Hybrid Metaheuristic Aroach
More informationGeneral Linear Model Introduction, Classes of Linear models and Estimation
Stat 740 General Linear Model Introduction, Classes of Linear models and Estimation An aim of scientific enquiry: To describe or to discover relationshis among events (variables) in the controlled (laboratory)
More informationUncorrelated Multilinear Principal Component Analysis for Unsupervised Multilinear Subspace Learning
TNN-2009-P-1186.R2 1 Uncorrelated Multilinear Princial Comonent Analysis for Unsuervised Multilinear Subsace Learning Haiing Lu, K. N. Plataniotis and A. N. Venetsanooulos The Edward S. Rogers Sr. Deartment
More informationPaper C Exact Volume Balance Versus Exact Mass Balance in Compositional Reservoir Simulation
Paer C Exact Volume Balance Versus Exact Mass Balance in Comositional Reservoir Simulation Submitted to Comutational Geosciences, December 2005. Exact Volume Balance Versus Exact Mass Balance in Comositional
More informationBayesian Spatially Varying Coefficient Models in the Presence of Collinearity
Bayesian Satially Varying Coefficient Models in the Presence of Collinearity David C. Wheeler 1, Catherine A. Calder 1 he Ohio State University 1 Abstract he belief that relationshis between exlanatory
More informationRadial Basis Function Networks: Algorithms
Radial Basis Function Networks: Algorithms Introduction to Neural Networks : Lecture 13 John A. Bullinaria, 2004 1. The RBF Maing 2. The RBF Network Architecture 3. Comutational Power of RBF Networks 4.
More informationMODELING THE RELIABILITY OF C4ISR SYSTEMS HARDWARE/SOFTWARE COMPONENTS USING AN IMPROVED MARKOV MODEL
Technical Sciences and Alied Mathematics MODELING THE RELIABILITY OF CISR SYSTEMS HARDWARE/SOFTWARE COMPONENTS USING AN IMPROVED MARKOV MODEL Cezar VASILESCU Regional Deartment of Defense Resources Management
More informationTowards understanding the Lorenz curve using the Uniform distribution. Chris J. Stephens. Newcastle City Council, Newcastle upon Tyne, UK
Towards understanding the Lorenz curve using the Uniform distribution Chris J. Stehens Newcastle City Council, Newcastle uon Tyne, UK (For the Gini-Lorenz Conference, University of Siena, Italy, May 2005)
More informationA Recursive Block Incomplete Factorization. Preconditioner for Adaptive Filtering Problem
Alied Mathematical Sciences, Vol. 7, 03, no. 63, 3-3 HIKARI Ltd, www.m-hiari.com A Recursive Bloc Incomlete Factorization Preconditioner for Adative Filtering Problem Shazia Javed School of Mathematical
More informationMany spatial attributes are classified into mutually exclusive
Published online May 16, 2007 Transiograms for Characterizing Satial Variability of Soil Classes Weidong Li* De. of Geograhy Kent State Univ. Kent, OH 44242 The characterization of comlex autocorrelations
More informationShadow Computing: An Energy-Aware Fault Tolerant Computing Model
Shadow Comuting: An Energy-Aware Fault Tolerant Comuting Model Bryan Mills, Taieb Znati, Rami Melhem Deartment of Comuter Science University of Pittsburgh (bmills, znati, melhem)@cs.itt.edu Index Terms
More informationUse of Transformations and the Repeated Statement in PROC GLM in SAS Ed Stanek
Use of Transformations and the Reeated Statement in PROC GLM in SAS Ed Stanek Introduction We describe how the Reeated Statement in PROC GLM in SAS transforms the data to rovide tests of hyotheses of interest.
More informationModeling and Estimation of Full-Chip Leakage Current Considering Within-Die Correlation
6.3 Modeling and Estimation of Full-Chi Leaage Current Considering Within-Die Correlation Khaled R. eloue, Navid Azizi, Farid N. Najm Deartment of ECE, University of Toronto,Toronto, Ontario, Canada {haled,nazizi,najm}@eecg.utoronto.ca
More informationJohn Weatherwax. Analysis of Parallel Depth First Search Algorithms
Sulementary Discussions and Solutions to Selected Problems in: Introduction to Parallel Comuting by Viin Kumar, Ananth Grama, Anshul Guta, & George Karyis John Weatherwax Chater 8 Analysis of Parallel
More informationq-ary Symmetric Channel for Large q
List-Message Passing Achieves Caacity on the q-ary Symmetric Channel for Large q Fan Zhang and Henry D Pfister Deartment of Electrical and Comuter Engineering, Texas A&M University {fanzhang,hfister}@tamuedu
More informationEstimation of Separable Representations in Psychophysical Experiments
Estimation of Searable Reresentations in Psychohysical Exeriments Michele Bernasconi (mbernasconi@eco.uninsubria.it) Christine Choirat (cchoirat@eco.uninsubria.it) Raffaello Seri (rseri@eco.uninsubria.it)
More informationAI*IA 2003 Fusion of Multiple Pattern Classifiers PART III
AI*IA 23 Fusion of Multile Pattern Classifiers PART III AI*IA 23 Tutorial on Fusion of Multile Pattern Classifiers by F. Roli 49 Methods for fusing multile classifiers Methods for fusing multile classifiers
More informationChapter 7 Sampling and Sampling Distributions. Introduction. Selecting a Sample. Introduction. Sampling from a Finite Population
Chater 7 and s Selecting a Samle Point Estimation Introduction to s of Proerties of Point Estimators Other Methods Introduction An element is the entity on which data are collected. A oulation is a collection
More informationApproximating min-max k-clustering
Aroximating min-max k-clustering Asaf Levin July 24, 2007 Abstract We consider the roblems of set artitioning into k clusters with minimum total cost and minimum of the maximum cost of a cluster. The cost
More informationA Comparison between Biased and Unbiased Estimators in Ordinary Least Squares Regression
Journal of Modern Alied Statistical Methods Volume Issue Article 7 --03 A Comarison between Biased and Unbiased Estimators in Ordinary Least Squares Regression Ghadban Khalaf King Khalid University, Saudi
More informationMorten Frydenberg Section for Biostatistics Version :Friday, 05 September 2014
Morten Frydenberg Section for Biostatistics Version :Friday, 05 Setember 204 All models are aroximations! The best model does not exist! Comlicated models needs a lot of data. lower your ambitions or get
More informationA Qualitative Event-based Approach to Multiple Fault Diagnosis in Continuous Systems using Structural Model Decomposition
A Qualitative Event-based Aroach to Multile Fault Diagnosis in Continuous Systems using Structural Model Decomosition Matthew J. Daigle a,,, Anibal Bregon b,, Xenofon Koutsoukos c, Gautam Biswas c, Belarmino
More informationTime Frequency Aggregation Performance Optimization of Power Quality Disturbances Based on Generalized S Transform
Time Frequency Aggregation Perormance Otimization o Power Quality Disturbances Based on Generalized S Transorm Mengda Li Shanghai Dianji University, Shanghai 01306, China limd @ sdju.edu.cn Abstract In
More informationResearch of power plant parameter based on the Principal Component Analysis method
Research of ower lant arameter based on the Princial Comonent Analysis method Yang Yang *a, Di Zhang b a b School of Engineering, Bohai University, Liaoning Jinzhou, 3; Liaoning Datang international Jinzhou
More informationSystem Reliability Estimation and Confidence Regions from Subsystem and Full System Tests
009 American Control Conference Hyatt Regency Riverfront, St. Louis, MO, USA June 0-, 009 FrB4. System Reliability Estimation and Confidence Regions from Subsystem and Full System Tests James C. Sall Abstract
More informationInformation collection on a graph
Information collection on a grah Ilya O. Ryzhov Warren Powell February 10, 2010 Abstract We derive a knowledge gradient olicy for an otimal learning roblem on a grah, in which we use sequential measurements
More informationPositive decomposition of transfer functions with multiple poles
Positive decomosition of transfer functions with multile oles Béla Nagy 1, Máté Matolcsi 2, and Márta Szilvási 1 Deartment of Analysis, Technical University of Budaest (BME), H-1111, Budaest, Egry J. u.
More informationMATHEMATICAL MODELLING OF THE WIRELESS COMMUNICATION NETWORK
Comuter Modelling and ew Technologies, 5, Vol.9, o., 3-39 Transort and Telecommunication Institute, Lomonosov, LV-9, Riga, Latvia MATHEMATICAL MODELLIG OF THE WIRELESS COMMUICATIO ETWORK M. KOPEETSK Deartment
More informationSAS for Bayesian Mediation Analysis
Paer 1569-2014 SAS for Bayesian Mediation Analysis Miočević Milica, Arizona State University; David P. MacKinnon, Arizona State University ABSTRACT Recent statistical mediation analysis research focuses
More informationA SIMPLE PLASTICITY MODEL FOR PREDICTING TRANSVERSE COMPOSITE RESPONSE AND FAILURE
THE 19 TH INTERNATIONAL CONFERENCE ON COMPOSITE MATERIALS A SIMPLE PLASTICITY MODEL FOR PREDICTING TRANSVERSE COMPOSITE RESPONSE AND FAILURE K.W. Gan*, M.R. Wisnom, S.R. Hallett, G. Allegri Advanced Comosites
More informationMULTIVARIATE STATISTICAL PROCESS OF HOTELLING S T CONTROL CHARTS PROCEDURES WITH INDUSTRIAL APPLICATION
Journal of Statistics: Advances in heory and Alications Volume 8, Number, 07, Pages -44 Available at htt://scientificadvances.co.in DOI: htt://dx.doi.org/0.864/jsata_700868 MULIVARIAE SAISICAL PROCESS
More informationAlgorithms for Air Traffic Flow Management under Stochastic Environments
Algorithms for Air Traffic Flow Management under Stochastic Environments Arnab Nilim and Laurent El Ghaoui Abstract A major ortion of the delay in the Air Traffic Management Systems (ATMS) in US arises
More informationDETERMINISTIC STOCHASTIC SUBSPACE IDENTIFICATION FOR BRIDGES
DEERMINISIC SOCHASIC SBSPACE IDENIFICAION FOR BRIDGES H. hai, V. DeBrunner, L. S. DeBrunner, J. P. Havlicek 2, K. Mish 3, K. Ford 2, A. Medda Deartment o Electrical and Comuter Engineering Florida State
More informationBrownian Motion and Random Prime Factorization
Brownian Motion and Random Prime Factorization Kendrick Tang June 4, 202 Contents Introduction 2 2 Brownian Motion 2 2. Develoing Brownian Motion.................... 2 2.. Measure Saces and Borel Sigma-Algebras.........
More informationApproximate Dynamic Programming for Dynamic Capacity Allocation with Multiple Priority Levels
Aroximate Dynamic Programming for Dynamic Caacity Allocation with Multile Priority Levels Alexander Erdelyi School of Oerations Research and Information Engineering, Cornell University, Ithaca, NY 14853,
More informationDETC2003/DAC AN EFFICIENT ALGORITHM FOR CONSTRUCTING OPTIMAL DESIGN OF COMPUTER EXPERIMENTS
Proceedings of DETC 03 ASME 003 Design Engineering Technical Conferences and Comuters and Information in Engineering Conference Chicago, Illinois USA, Setember -6, 003 DETC003/DAC-48760 AN EFFICIENT ALGORITHM
More informationSupplementary Materials for Robust Estimation of the False Discovery Rate
Sulementary Materials for Robust Estimation of the False Discovery Rate Stan Pounds and Cheng Cheng This sulemental contains roofs regarding theoretical roerties of the roosed method (Section S1), rovides
More informationSpin as Dynamic Variable or Why Parity is Broken
Sin as Dynamic Variable or Why Parity is Broken G. N. Golub golubgn@meta.ua There suggested a modification of the Dirac electron theory, eliminating its mathematical incomleteness. The modified Dirac electron,
More informationBayesian inference & Markov chain Monte Carlo. Note 1: Many slides for this lecture were kindly provided by Paul Lewis and Mark Holder
Bayesian inference & Markov chain Monte Carlo Note 1: Many slides for this lecture were kindly rovided by Paul Lewis and Mark Holder Note 2: Paul Lewis has written nice software for demonstrating Markov
More informationLinear diophantine equations for discrete tomography
Journal of X-Ray Science and Technology 10 001 59 66 59 IOS Press Linear diohantine euations for discrete tomograhy Yangbo Ye a,gewang b and Jiehua Zhu a a Deartment of Mathematics, The University of Iowa,
More informationDistributed Rule-Based Inference in the Presence of Redundant Information
istribution Statement : roved for ublic release; distribution is unlimited. istributed Rule-ased Inference in the Presence of Redundant Information June 8, 004 William J. Farrell III Lockheed Martin dvanced
More informationPassive Identification is Non Stationary Objects With Closed Loop Control
IOP Conerence Series: Materials Science and Engineering PAPER OPEN ACCESS Passive Identiication is Non Stationary Obects With Closed Loo Control To cite this article: Valeriy F Dyadik et al 2016 IOP Con.
More informationAnalysis of M/M/n/K Queue with Multiple Priorities
Analysis of M/M/n/K Queue with Multile Priorities Coyright, Sanjay K. Bose For a P-riority system, class P of highest riority Indeendent, Poisson arrival rocesses for each class with i as average arrival
More informationCMSC 425: Lecture 4 Geometry and Geometric Programming
CMSC 425: Lecture 4 Geometry and Geometric Programming Geometry for Game Programming and Grahics: For the next few lectures, we will discuss some of the basic elements of geometry. There are many areas
More informationStatics and dynamics: some elementary concepts
1 Statics and dynamics: some elementary concets Dynamics is the study of the movement through time of variables such as heartbeat, temerature, secies oulation, voltage, roduction, emloyment, rices and
More informationOptimal array pattern synthesis with desired magnitude response
Otimal array attern synthesis with desired magnitude resonse A.M. Pasqual a, J.R. Arruda a and P. erzog b a Universidade Estadual de Caminas, Rua Mendeleiev, 00, Cidade Universitária Zeferino Vaz, 13083-970
More informationPlotting the Wilson distribution
, Survey of English Usage, University College London Setember 018 1 1. Introduction We have discussed the Wilson score interval at length elsewhere (Wallis 013a, b). Given an observed Binomial roortion
More informationLower bound solutions for bearing capacity of jointed rock
Comuters and Geotechnics 31 (2004) 23 36 www.elsevier.com/locate/comgeo Lower bound solutions for bearing caacity of jointed rock D.J. Sutcliffe a, H.S. Yu b, *, S.W. Sloan c a Deartment of Civil, Surveying
More informationPrinciples of Computed Tomography (CT)
Page 298 Princiles of Comuted Tomograhy (CT) The theoretical foundation of CT dates back to Johann Radon, a mathematician from Vienna who derived a method in 1907 for rojecting a 2-D object along arallel
More informationFig. 4. Example of Predicted Workers. Fig. 3. A Framework for Tackling the MQA Problem.
217 IEEE 33rd International Conference on Data Engineering Prediction-Based Task Assignment in Satial Crowdsourcing Peng Cheng #, Xiang Lian, Lei Chen #, Cyrus Shahabi # Hong Kong University of Science
More informationEfficient Approximations for Call Admission Control Performance Evaluations in Multi-Service Networks
Efficient Aroximations for Call Admission Control Performance Evaluations in Multi-Service Networks Emre A. Yavuz, and Victor C. M. Leung Deartment of Electrical and Comuter Engineering The University
More informationGenetic Algorithms, Selection Schemes, and the Varying Eects of Noise. IlliGAL Report No November Department of General Engineering
Genetic Algorithms, Selection Schemes, and the Varying Eects of Noise Brad L. Miller Det. of Comuter Science University of Illinois at Urbana-Chamaign David E. Goldberg Det. of General Engineering University
More informationInformation collection on a graph
Information collection on a grah Ilya O. Ryzhov Warren Powell October 25, 2009 Abstract We derive a knowledge gradient olicy for an otimal learning roblem on a grah, in which we use sequential measurements
More informationAggregate Prediction With. the Aggregation Bias
100 Aggregate Prediction With Disaggregate Models: Behavior of the Aggregation Bias Uzi Landau, Transortation Research nstitute, Technion-srael nstitute of Technology, Haifa Disaggregate travel demand
More informationAsymptotically Optimal Simulation Allocation under Dependent Sampling
Asymtotically Otimal Simulation Allocation under Deendent Samling Xiaoing Xiong The Robert H. Smith School of Business, University of Maryland, College Park, MD 20742-1815, USA, xiaoingx@yahoo.com Sandee
More informationLending Direction to Neural Networks. Richard S. Zemel North Torrey Pines Rd. La Jolla, CA Christopher K. I.
Lending Direction to Neural Networks Richard S Zemel CNL, The Salk Institute 10010 North Torrey Pines Rd La Jolla, CA 92037 Christoher K I Williams Deartment of Comuter Science University of Toronto Toronto,
More informationHidden Predictors: A Factor Analysis Primer
Hidden Predictors: A Factor Analysis Primer Ryan C Sanchez Western Washington University Factor Analysis is a owerful statistical method in the modern research sychologist s toolbag When used roerly, factor
More informationPairwise active appearance model and its application to echocardiography tracking
Pairwise active aearance model and its alication to echocardiograhy tracking S. Kevin Zhou 1, J. Shao 2, B. Georgescu 1, and D. Comaniciu 1 1 Integrated Data Systems, Siemens Cororate Research, Inc., Princeton,
More informationKEY ISSUES IN THE ANALYSIS OF PILES IN LIQUEFYING SOILS
4 th International Conference on Earthquake Geotechnical Engineering June 2-28, 27 KEY ISSUES IN THE ANALYSIS OF PILES IN LIQUEFYING SOILS Misko CUBRINOVSKI 1, Hayden BOWEN 1 ABSTRACT Two methods for analysis
More informationPER-PATCH METRIC LEARNING FOR ROBUST IMAGE MATCHING. Sezer Karaoglu, Ivo Everts, Jan C. van Gemert, and Theo Gevers
PER-PATCH METRIC LEARNING FOR ROBUST IMAGE MATCHING Sezer Karaoglu, Ivo Everts, Jan C. van Gemert, and Theo Gevers Intelligent Systems Lab, Amsterdam, University of Amsterdam, 1098 XH Amsterdam, The Netherlands
More informationAn Analysis of Reliable Classifiers through ROC Isometrics
An Analysis of Reliable Classifiers through ROC Isometrics Stijn Vanderlooy s.vanderlooy@cs.unimaas.nl Ida G. Srinkhuizen-Kuyer kuyer@cs.unimaas.nl Evgueni N. Smirnov smirnov@cs.unimaas.nl MICC-IKAT, Universiteit
More informationA MIXED CONTROL CHART ADAPTED TO THE TRUNCATED LIFE TEST BASED ON THE WEIBULL DISTRIBUTION
O P E R A T I O N S R E S E A R C H A N D D E C I S I O N S No. 27 DOI:.5277/ord73 Nasrullah KHAN Muhammad ASLAM 2 Kyung-Jun KIM 3 Chi-Hyuck JUN 4 A MIXED CONTROL CHART ADAPTED TO THE TRUNCATED LIFE TEST
More informationUnobservable Selection and Coefficient Stability: Theory and Evidence
Unobservable Selection and Coefficient Stability: Theory and Evidence Emily Oster Brown University and NBER August 9, 016 Abstract A common aroach to evaluating robustness to omitted variable bias is to
More informationImproved Capacity Bounds for the Binary Energy Harvesting Channel
Imroved Caacity Bounds for the Binary Energy Harvesting Channel Kaya Tutuncuoglu 1, Omur Ozel 2, Aylin Yener 1, and Sennur Ulukus 2 1 Deartment of Electrical Engineering, The Pennsylvania State University,
More informationA generalization of Amdahl's law and relative conditions of parallelism
A generalization of Amdahl's law and relative conditions of arallelism Author: Gianluca Argentini, New Technologies and Models, Riello Grou, Legnago (VR), Italy. E-mail: gianluca.argentini@riellogrou.com
More informationCHAPTER 5 STATISTICAL INFERENCE. 1.0 Hypothesis Testing. 2.0 Decision Errors. 3.0 How a Hypothesis is Tested. 4.0 Test for Goodness of Fit
Chater 5 Statistical Inference 69 CHAPTER 5 STATISTICAL INFERENCE.0 Hyothesis Testing.0 Decision Errors 3.0 How a Hyothesis is Tested 4.0 Test for Goodness of Fit 5.0 Inferences about Two Means It ain't
More informationAn Improved Calibration Method for a Chopped Pyrgeometer
96 JOURNAL OF ATMOSPHERIC AND OCEANIC TECHNOLOGY VOLUME 17 An Imroved Calibration Method for a Choed Pyrgeometer FRIEDRICH FERGG OtoLab, Ingenieurbüro, Munich, Germany PETER WENDLING Deutsches Forschungszentrum
More informationEcological Resemblance. Ecological Resemblance. Modes of Analysis. - Outline - Welcome to Paradise
Ecological Resemblance - Outline - Ecological Resemblance Mode of analysis Analytical saces Association Coefficients Q-mode similarity coefficients Symmetrical binary coefficients Asymmetrical binary coefficients
More informationCharacteristics of Beam-Based Flexure Modules
Shorya Awtar e-mail: shorya@mit.edu Alexander H. Slocum e-mail: slocum@mit.edu Precision Engineering Research Grou, Massachusetts Institute of Technology, Cambridge, MA 039 Edi Sevincer Omega Advanced
More informationAdaptive estimation with change detection for streaming data
Adative estimation with change detection for streaming data A thesis resented for the degree of Doctor of Philosohy of the University of London and the Diloma of Imerial College by Dean Adam Bodenham Deartment
More informationCharacterizing the Behavior of a Probabilistic CMOS Switch Through Analytical Models and Its Verification Through Simulations
Characterizing the Behavior of a Probabilistic CMOS Switch Through Analytical Models and Its Verification Through Simulations PINAR KORKMAZ, BILGE E. S. AKGUL and KRISHNA V. PALEM Georgia Institute of
More informationute measures of uncertainty called standard errors for these b j estimates and the resulting forecasts if certain conditions are satis- ed. Note the e
Regression with Time Series Errors David A. Dickey, North Carolina State University Abstract: The basic assumtions of regression are reviewed. Grahical and statistical methods for checking the assumtions
More informationLOGISTIC REGRESSION. VINAYANAND KANDALA M.Sc. (Agricultural Statistics), Roll No I.A.S.R.I, Library Avenue, New Delhi
LOGISTIC REGRESSION VINAANAND KANDALA M.Sc. (Agricultural Statistics), Roll No. 444 I.A.S.R.I, Library Avenue, New Delhi- Chairerson: Dr. Ranjana Agarwal Abstract: Logistic regression is widely used when
More informationSCHUR S LEMMA AND BEST CONSTANTS IN WEIGHTED NORM INEQUALITIES. Gord Sinnamon The University of Western Ontario. December 27, 2003
SCHUR S LEMMA AND BEST CONSTANTS IN WEIGHTED NORM INEQUALITIES Gord Sinnamon The University of Western Ontario December 27, 23 Abstract. Strong forms of Schur s Lemma and its converse are roved for mas
More informationCorrespondence Between Fractal-Wavelet. Transforms and Iterated Function Systems. With Grey Level Maps. F. Mendivil and E.R.
1 Corresondence Between Fractal-Wavelet Transforms and Iterated Function Systems With Grey Level Mas F. Mendivil and E.R. Vrscay Deartment of Alied Mathematics Faculty of Mathematics University of Waterloo
More informationObserver/Kalman Filter Time Varying System Identification
Observer/Kalman Filter Time Varying System Identification Manoranjan Majji Texas A&M University, College Station, Texas, USA Jer-Nan Juang 2 National Cheng Kung University, Tainan, Taiwan and John L. Junins
More informationEstimating function analysis for a class of Tweedie regression models
Title Estimating function analysis for a class of Tweedie regression models Author Wagner Hugo Bonat Deartamento de Estatística - DEST, Laboratório de Estatística e Geoinformação - LEG, Universidade Federal
More informationChapter 7 Rational and Irrational Numbers
Chater 7 Rational and Irrational Numbers In this chater we first review the real line model for numbers, as discussed in Chater 2 of seventh grade, by recalling how the integers and then the rational numbers
More informationOn the Relationship Between Packet Size and Router Performance for Heavy-Tailed Traffic 1
On the Relationshi Between Packet Size and Router Performance for Heavy-Tailed Traffic 1 Imad Antonios antoniosi1@southernct.edu CS Deartment MO117 Southern Connecticut State University 501 Crescent St.
More informationOptimal Design of Truss Structures Using a Neutrosophic Number Optimization Model under an Indeterminate Environment
Neutrosohic Sets and Systems Vol 14 016 93 University of New Mexico Otimal Design of Truss Structures Using a Neutrosohic Number Otimization Model under an Indeterminate Environment Wenzhong Jiang & Jun
More informationSampling. Inferential statistics draws probabilistic conclusions about populations on the basis of sample statistics
Samling Inferential statistics draws robabilistic conclusions about oulations on the basis of samle statistics Probability models assume that every observation in the oulation is equally likely to be observed
More informationPartial Identification in Triangular Systems of Equations with Binary Dependent Variables
Partial Identification in Triangular Systems of Equations with Binary Deendent Variables Azeem M. Shaikh Deartment of Economics University of Chicago amshaikh@uchicago.edu Edward J. Vytlacil Deartment
More informationProbability Estimates for Multi-class Classification by Pairwise Coupling
Probability Estimates for Multi-class Classification by Pairwise Couling Ting-Fan Wu Chih-Jen Lin Deartment of Comuter Science National Taiwan University Taiei 06, Taiwan Ruby C. Weng Deartment of Statistics
More information