Software reliability estimation method based on Markov usage models using Importance Sampling

Importance sampling is a change-of-measure technique for speeding up the simulation of rare events in stochastic systems. In this paper the authors establish a number of properties characterizing optimal importance sampling measures for Markovian systems. The authors show how these properties may be used to compute the optimal measure and give specific results for a tandem queueing system. The authors' approach has no apparent computational advantage over other direct methods, but it suggests a number of heuristic approximations which may lead to computationally attractive methods.

We present several numerical approaches to investigate rare events in stochastic systems, with a specific focus on application to biological models. We first present several aspects concerning variance reduction of Monte-Carlo methods, with a focus on importance sampling. We show how these methods can be applied to basic biological models. We show that these techniques can be useful in dealing with multiscale continuous-time Markov chains that are important in the context of biochemical reaction networks. We treat with more detail the problem of first passage time for a linear diffusion process, arising from the integrate-and-fire neuron model, to show the kind of mathematical problems that may arise when trying to design an optimal importance sampling scheme. This leads to the observation that large deviation theory can be helpful to solve these questions. We also review a numerical method designed to estimate large deviations quantities such as the quasipotential and the optimal path and apply this method to estimate numerically the quasipotential of the Morris–Lecar neuron model perturbed by small noise.

Statistical testing based on a Markov chain usage model, as a rigorous testing method, has been around for more than two decades. Through the comprehensive application of statistical science to the testing of software, it provides audit trails of evidence to support correctness arguments for a software-intensive system as well as a decision that the system is of requisite quality for its intended use. This paper reports a real-world case study in which we applied standard statistical testing practices to the phone flag copy testing problem at the site of our industrial collaborator, and presents our solution from problem formalization, usage modeling and model analysis to test case generation and analysis. Our results helped evaluate the coverage of a heuristically generated test suite, and shed light on what other test cases to craft/generate to improve test coverage and adequacy, as well as reliability estimates both at the arc (usage event) level and at the path (system) level.

Rare events are an important limitation of discrete-event simulation in some areas of applications. The existing rare-event simulation (RES) techniques, such as importance sampling (IS), splitting, and RESTART, have been used for RES so far. However, these techniques not only have some limitations, but also are not general or easy to use. In this paper, the partition of the region (PoR) method is used as a variance reduction technique in a new RES technique for usage with a wide range of modeling and simulation languages. Three variants of the proposed technique will be introduced: simple PoR, recursive PoR, and advanced PoR. Simple and recursive PoR methods are powerful techniques that have some important features such as bounded relative error; however, they are restricted to some specific models. Advanced PoR is a more general technique that integrates the PoR method and the IS technique. To evaluate the efficiency of the proposed technique and its variants, the comparative simulation results of some illustrative examples using stochastic activity networks will be presented.

We review two types of adaptive Monte Carlo methods for rare event simulations. These methods are based on importance sampling. The first approach selects importance sampling distributions by minimizing the variance of importance sampling estimator. The second approach selects importance sampling distributions by minimizing the cross entropy to the optimal importance sampling distribution. We also review the basic concepts of importance sampling in the rare event simulation context. To make the basic concepts concrete, we introduce these ideas via the study of rare events of M/M/1 queues.

When a population is too large for study, as is the case for all possible uses of a software system, a statistically correct sample must be drawn as a basis for inferences about the population. In statistical testing of software based on a Markov chain usage model, the rich body of analytical results available for Markov chains provides numerous insights that can be used in test planning. Further, the connection between Markov chains and operations research techniques permits a Markov usage model to be expressed as a system of constraints, with mathematical programming used to generate the optimal model for a particular objective function. Since a software usage model is based on the specification, all analyses may be performed early in the development cycle and used as a quantitative basis for management decisions. These techniques have been reduced to engineering practice and used in large projects by IBM, Ericsson, all branches of the US military, and others. In this paper, statistical experiments, Markov models, and optimization techniques are shown to provide a sound theoretical and practical basis far quantifying the reliability of software.

Rare event simulation in the context of queueing networks has been an active area of research for more than two decades. A commonly used technique to increase the efficiency of Monte Carlo simulation is importance sampling. However, there are few rigorous results on the design of efficient or asymptotically optimal importance sampling schemes for queueing networks. Using a recently developed game/subsolution approach, we construct simple and efficient state-dependent importance sampling schemes for simulating buffer overflows in stable open Jackson networks. The sampling distributions do not depend on the particular event of interest, and hence overflow probabilities for different events can be estimated simultaneously. A by-product of the analysis is the identification of the minimizing trajectory for the calculus of variation problem that is associated with the sample-path large deviation rate function.

The paper concerns automatic generation of usage models from real-time software UML models. Firstly, we define the reasonably constrained real-time software UML artifacts, which include use case diagrams, timed sequence diagrams and the execution probability of each sequence diagram in its associated use case. Secondly, the paper presents a method that derives the software usage model from the constrained UML artifacts. The method elicits the messages associated with the objects under testing and their occurrence probabilities to generate the usage model of each use case. Timing constraints in sequence diagrams are considered during usage model generation. Then the usage models of use cases are integrated into the software usage model by utilizing the execution sequence relations between use cases. The usage models can be used to generate real-time software statistical test cases and facilitate real-time software statistical testing.

A new type of exponential asymptotical expansions in mixed large deviation and quasi-ergodic theorems and asymptotical expansions for quasistationary distributions are presented for nonlinearly perturbed Markov chains. Applications to analysis of rare events in stochastic systems are discussed.

Markovian jumping systems have arisen naturally in the mathematical modeling of phenomena spanning disciplines in the social sciences, natural sciences, and engineering. This kind of stochastic dynamical systems can be employed to model the dynamics when parameters are subject to random abrupt changes due to sudden environment changes, subsystem switching, system noises, executor faults, and so forth. Much attention has been given to modeling, optimization, and real applications of such stochastic dynamical systems in the literature in recent years. As the advanced control and optimization will provide a basis for the design and application of such stochastic systems, these advanced techniques would result in substantial and sustainable benefits. The accepted papers in this special issue include stochastic stability, stabilization, stochastic control optimization, system modeling and identification methods, predictive control, signal processing, robust filtering, multiagent systems, networked control systems, time-delayed systems, neural networks, the Takagi-Sugeno fuzzy systems, simulated annealing, and fault detection methods.

Clinical practice involves measuring quantities for a variety of purposes, such as aiding diagnosis, predicting future patient outcomes, and serving as endpoints in studies or randomized trials. Measurements are almost always prone to various sorts of errors, which cause the measured value to differ from the true value; accordingly, studies investigating measurement error frequently appear in this and other journals. The importance of measurement error depends upon the context in which the measurements in question are to be used. For example, a certain degree of measurement error may be acceptable if measurements are to be used as an outcome in a comparative study such as a clinical trial, but the same measurement errors may be unacceptably large to make measurements usable in individual patient management, such as screening or risk prediction. In the past 20 years many papers have been published advocating how studies of measurement error should be analyzed, with a paper by Bland and Altman1 being one of the most cited and well known examples. There has been much controversy concerning the choice of parameter to be estimated and reported, and consequently confusion surrounding the meaning and interpretation of results from studies investigating measurement error. In this paper we first distinguish between the general concepts of agreement and reliability to aid researchers in considering which are relevant for their particular application. We then review the statistical methods that can be used to investigate and quantify agreement and reliability, dealing separately with the different types of measurement error study, while emphasizing the largely common techniques that should be used for data analysis. We reiterate that the judgment of whether agreement or reliability are acceptable must be related to the clinical application, and cannot be proven by a statistical test. We highlight the fact that reliability depends on the population in which measurements are made, and not just on the measurement errors of the measurement method. We discuss the advantages of method comparison studies making at least two measurements with each measurement method on each subject. A key advantage is that the cause of a correlation between paired differences and means in the so-called Bland–Altman plot can be determined, in contrast to when only a single measurement is made with each method. Throughout the paper, we try to emphasize that calculated values of agreement and reliability from measurement error studies are estimates of parameters, and as such we should report such estimates with CIs to indicate the uncertainty with which they have been estimated. We restrict our attention to measurements of a continuous quantity; alternative methods are required for categorical data2. One difficulty in the measurement error field is the number of different terms used to describe studies of measurement error. The terms 'agreement', 'reliability', 'reproducibility' and 'repeatability' are used with varying degrees of consistency in the medical literature. We first make clear the distinction between the statistical concepts of agreement and reliability3. Agreement quantifies how close two measurements made on the same subject are, and is measured on the same scale as the measurements themselves. Two measurements of the same subject may differ for a number of reasons, depending on the conditions under which the measurements were made. In a method comparison study there will be differences because of inherent variability in each of the measurement methods, as well as potentially a bias between the measurements from the methods. If the measurements are made by different observers or raters, differences may be due to biases between the observers. Agreement between measurements is a characteristic of the measurement method(s) involved, which does not depend on the population in which measurements are made, unless bias or measurement precision varies with the true value being measured. One popular way of quantifying agreement is to estimate the 95% limits of agreement, as proposed by Bland and Altman1. These limits are defined such that we expect that, in the long run, 95% of future differences between measurements made on the same subject will lie within the limits. If reliability is high, measurement errors are small in comparison to the true differences between subjects, so that subjects can be relatively well distinguished (in terms of the quantity being measured) on the basis of the error-prone measurements. Conversely, if measurement errors tend to be large compared with the true differences between subjects, reliability will be low because differences between measurements of two subjects could be due purely to error rather than to a genuine difference in their true values. The reliability parameter is also known as an intraclass correlation (ICC), as it equals the correlation between any two measurements made on the same subject. Reliability takes values between zero and one, with a value of one corresponding to zero measurement error and a value of zero meaning that all the variability in measurements is due to measurement error. As a dimensionless quantity, it is arguably quite difficult to interpret, and deciding what value constitutes sufficiently high reliability is often made in a subjective fashion. Repeatability of measurements refers to the variation in repeat measurements made on the same subject under identical conditions4. This means that measurements are made by the same instrument or method, the same observer (or rater) if human input is required, and that the measurements are made over a short period of time, over which the underlying value can be considered to be constant. Variability in measurements made on the same subject in a repeatability study can then be ascribed only to errors due to the measurement process itself. Reproducibility refers to the variation in measurements made on a subject under changing conditions4. The changing conditions may be due to different measurement methods or instruments being used, measurements being made by different observers or raters, or measurements being made over a period of time, within which the 'error-free' level of the variable could undergo non-negligible change. The first type of study we consider is a repeatability study, in which we investigate and quantify the repeatability of measurements made by a single instrument or method, and in which the conditions of measurement remain constant. The second type of study we consider is a method comparison study, in which measurements are made using two measurement methods on a sample of subjects. The use of two different methods means that this is a reproducibility study, but the term method comparison is used as it clearly communicates the changing conditions under which measurements have been made. In contrast to a repeatability study, systematic bias may exist between measurements made by two different methods, and their measurement errors may have different SDs. The last type of study we consider is one in which measurements are made by different observers or raters. Again this is a type of reproducibility study. As with method comparison studies, biases may exist between observers, and their measurement SDs may differ. As we discuss below (see Measurements with observers or raters), if interest lies in quantifying the measurement error characteristics of the particular observers in one's study, exactly the same analysis methods should be used as for a method comparison study. In order to investigate the repeatability of measurements, a repeatability study must, for an appropriately selected sample, make at least two measurements per subject under identical conditions. This means that the measurements must be made by the same measurement method, or the same observer or rater. The objective is to then quantify the agreement and reliability of measurements made by that particular method or observer. If the differences between two measurements made on a subject are approximately Normally distributed, in the long run we expect the absolute difference between two measurements on a subject to differ by no more than the repeatability coefficient on 95% of occasions. To estimate the within-subject SD, we can fit a one-way analysis of variance (ANOVA) model to the data containing the repeat measurements made on subjects. ANOVAs can be fitted in all modern statistical packages. ANOVA models partition variability in data into that which can be ascribed to differences between groups, and that remaining to within groups. For a repeatability study, the groups are defined by the subjects under measurement, and this must be specified in the statistical package used. The ANOVA model estimates how much variation in measurements can be attributed to differences in the true, or 'error-free' values of subjects, with the remainder constituting measurement error. Fitting the ANOVA model results in estimates of the between-subject and within-subject SDs (or alternatively the corresponding variances, which are the SDs squared). The estimate of the within-subject SD can be used in the above formula to give an estimate of the repeatability coefficient. We illustrate with a study by Järvelä et al., who investigated the repeatability of measurements of flow volume made by one observer using three-dimensional (3D) power Doppler ultrasonography from 22 ovaries5. The estimated repeatability coefficient was 3.97, meaning that the absolute difference between any two future measurements made by that particular observer on a particular subject/unit are estimated to be no greater than 3.97 on 95% of occasions. It is important to note that the repeatability of another observer may be different, because of differences in the training and ability of observers. Because the repeatability coefficient calculated is an estimate, it is important to calculate a CI for it to indicate how precisely it has been estimated. A CI may be given automatically by statistical software, but in the Appendix we review how a 95% CI for the within-subject SD can be calculated. Such a CI can be used to find a CI for the repeatability coefficient by multiplying the CI limits by . If the CI for the within-subject variance is given by software instead, the limits must first be square rooted to give a CI for the within-subject SD. The ANOVA model assumes that the measurement errors are statistically independent of the true 'error-free' value, and that the SD of the errors is constant throughout the range of 'error-free' values. Sometimes the SD of errors increases with the true value being measured. This should be checked by plotting paired differences between measurements against their mean, the so called Bland–Altman plot. We illustrate this later in the context of a method comparison study (see Plotting the data) and describe how such heterogeneous errors can often be dealt with by making a log transformation (see Non-uniform variability in measurement of the repeatability coefficient on the differences between measurements being approximately Normally This can be checked by a or plot of the paired differences in measurements on each subject. made in the of the CI for the within-subject SD is that the measurement errors are If this is in the may be to The reliability of a measurement method is often of interest when measurements are to be used to between subjects or groups of subjects. For example, if we have a choice of two measurement methods that could be used to an outcome in a clinical or study, using the method with reliability will give greater statistical power to differences between groups, for a given sample In this we describe how measurements from a single method can be used to estimate the reliability in a given we discuss the use of reliability to measurement methods (see Reliability in method comparison and different observers (see in The same ANOVA model as above (see can be used to estimate statistical will automatically give the estimated which is the of we can the estimates of the and within-subject SDs into the formula given above (see Agreement and to estimate the For the flow data of Järvelä et an estimated of was This means that of the variability in measurements of the flow was estimated to be due to genuine differences in flow between with the remaining being due to errors in the measurement process and the observer Because the measurements were all made by one the reliability may be to as As with the repeatability the calculated is an estimate, and a 95% CI should be statistical give a but in the Appendix we review the of a 95% CI for the For the flow data of Järvelä et the estimate of a 95% of to Because we use the same ANOVA model to estimate reliability as we describe for agreement, the same of 95% CIs for the estimate on of the true 'error-free' values. The reliability of a measurement method depends upon the of the population in which the measurements are made. the of reliability given above (see Agreement and we that the of true 'error-free' values in the measured by the between-subject SD, the value of agreement between repeat measurements is a characteristic of the method or instrument the of measurement errors is the range of true reliability depends on the of measurement errors and the true in the population in which measurements are made. The is with a a method for measuring has a within-subject SD of which does not with the underlying value being measured. we a study to estimate the reliability, and we sample from a heterogeneous population in which the between-subject SD in true is 20 This give a reliability or of . we sample our subjects from a more in which the SD of true values is the same as the variability of the within-subject error. In this the reliability is to . If studies report only an estimate of the reliability (ICC), can only make use of the estimate if the population in which the to use such measurements has We that report estimates of and within-subject SD, in to the In this can whether the measurement method will be sufficiently for their application, in which the between subjects may be we use a measurement method or in clinical we must that the measurements it are to by the measurement method used, that measurements made with the method are using the method. If the measurements from the two methods are sufficiently close and the of a patient on their basis be the the method could the method in clinical because the method is or to The as to what constitutes depends on how the measurements are to be used. If the two measurements are made on the same or we can quantify their To investigate and quantify the agreement between measurements made by two methods, we must at a a sample of subjects using the measurement method and the proposed method. The data from such a study of of measurements from each with the containing the measurements from the two measurement methods. we making two measurements on each subject using each method because the one measurement per method is the most common we by analysis. The first to such a is to plot the The plot is of measurements from the method against from the method (or If measurements were from we expect the to lie on the of the of can be used to if there is if more data are above the than this that the method on the measurements on the of fit for this plot and report a statistical for whether the of the plot from the of As by Bland and in a paper in this we expect the to have a than one if the method on the any measurement and so a of the that the true is to one is not a it the same as a plot with the of of the between the measurements from two methods is often more by plotting the difference in a measurements from the two methods against the of their measurements, as first by and this is to as the Bland–Altman and is frequently in of measurement error We illustrate the plot using an from the literature. et volume measurement in with and volume was measured using and also by We have approximately their plot of differences against means and it is as using only data from the in volume measured by three-dimensional and by against their from study by et the SD and the SD, each with 95% CI The plot the difference between calculated using and against the of two values. et not their plot or of results to indicate which measurements were from It is important to be and always the same measurement from the The the of the paired from zero an estimate of the bias between the two methods. which measurements have been from which to which method, on measurements. The indicate the estimated limits of agreement their CI which we describe The plot can be used to the differences between measurements made by the two methods. The variability of the differences between the two methods how well the methods If the variability of the paired differences is the range of measurements, and there is no between the difference and mean, we can quantify the agreement between the methods by the limits of the of the limits of agreement, we illustrate how the for limits of agreement may be and how they can be dealt on we expect the difference in volume as measured by the two methods to lie between and for 95% of future measurements. It is important to that the limits of agreement calculated in this are just and as with any type of estimate it is to quantify and report how precisely the limits are estimated of 95% The estimated limits of agreement indicate how large the between the measurements will be on 95% of occasions. The whether the methods sufficiently well must then be made depending on the context in which the measurements will be used. In contrast to the repeatability which assumes no bias between measurements, the limits of agreement method this The of the paired differences whether on one method to or measurements to the measurements of the second method, which we to as a bias between the methods. In the data in from et al., there was of a bias between the two methods, because the difference between the methods was close to at We can a statistical to whether there is of bias by a of the differences paired of the measurements from each we the that the true of the differences is corresponding to no bias between the methods. For the data in this that there was no of a bias between and made using The limits of agreement method assumes that the SD of the differences is throughout the range of measurements. Sometimes this will not be the In frequently the SD of the differences will with the In this for small values the limits of agreement will be than for values the limits will be of the plot of differences against from the study by et that the variability of the differences may be when the value being measured is The SD of the paired differences below the of the means was it was for above the This can often be by of the measurements of the two If the difference in the of the measurements results in constant we can calculate the limits of agreement and CIs for the limits in the Because the difference of two is to the of the we can all the estimates made on the log scale by Because of the limits of agreement to the of one measurements to the For the data in the of the difference in measurements was with SD The limits of agreement on the log scale are to the of with 95% limits of agreement for the of to The 95% CIs can be calculated on the log scale to the which measurements were made by which method, we cannot whether the is for measurements by to by or a plot of the of the measurements by the two methods to their mean, using a with the estimated 95% limits of agreement and their 95% This that the of a SD for the paired differences is more of measured by three-dimensional (3D) and by against of measurements, on data published by et the SD and the SD, each with 95% CI to log transformation can also be used to try to For example, et estimate limits of agreement using paired differences in measurements as a of the of the two measurements. The key is that a transformation should be used which that are approximately Normally and have approximately SD the range of It is for there to be an between the paired differences and of this is in a study by et al., who the agreement between as measured using and the calculated using in We have the data from their Bland–Altman plot of the difference between the volume and volume against the mean, and it in There to be a between the paired differences and differences between volume measured by three-dimensional and by from study by et correlation the SD and the SD, each with 95% CI We can a statistical to the for a whether the correlation coefficient between the paired differences and means from zero or by of the differences against the The is known as in the context of the SD of measurements from two methods, which for the data in is statistically for an between paired differences and how should we than models for the difference as a of the we should be given to the cause of the A between the paired differences and means when the SD of measurements from one method from the SD of measurements from the other method. There are at least two of an between paired differences and The first is that there is between the difference in measurements from the two methods and true value being the bias between methods over the range of true values. The second is that the within-subject SDs of the two methods differ. This will in the of changing bias if a method has or measurement errors than the method. This is and when a or measurement method is compared with a If the cause of the correlation is within-subject SDs and we were to plot the paired differences against the true value being there be no The between the Bland–Altman plot and a plot against true values because the paired means measurement errors, in contrast to the true values. as by and with only one measurement per method per we cannot which of two is the To we must make an that the within-subject SDs of the two methods are so that the correlation is attributed to changing or that the two within-subject SDs are and that there in no between bias and true value we may have a of For the of we with the from et by that the correlation is due to a difference in the within-subject SDs of the two methods. A correlation between differences and means will when the method with measurement error is from measurements of the method with error. For the data in the correlation is the errors than the measured by which and to the volume measured by as the true that such measurements were to no an that there is no changing we can calculate 95% limits of agreement, as by et al., because the limits of agreement method does not that the two within-subject SDs are In order to from the data whether changing bias is we must make at least two measurements with each method on each subject. This study has often been but is we describe how such studies can be to whether a correlation between paired differences and means is due to changing bias between the methods (see with two measurements from each As reliability may be a parameter with which to two different measurement methods. To estimate each reliability, we must make at least two measurements of each subject with each of the two methods. The repeat measurements from each method can then be as two repeatability studies (see estimates of each reliability, which can be advantage of using reliability to measurement methods is that it can be used to methods when their measurements are given on different or as the reliability is a dimensionless Because reliability depends on the of the true values in the population (see Reliability depends on population it is that reliability are compared only if they have been estimated from the same report a single reliability estimate from method comparison studies in which it that only a single measurement per subject per method is It appear that the estimates are using a one-way ANOVA as above (see The one-way ANOVA model assumes that there are no systematic biases between measurements within subjects, and that the within-subject SDs are for all measurements measurements are from two different methods, of may be bias may exist between the methods and their within-subject SDs may differ. If is the model is and the estimates of within and between subject SD are estimates of different The reliability estimate then has a different interpretation from the We against the estimate from the one-way ANOVA to data from two different measurement methods. As making two measurements with each method an of whether bias between the methods is constant or whether the measurement error SDs differ. the measurements from each method can be separately as two repeatability studies, using the methods estimates of the repeatability coefficient and reliability for each method. The Bland–Altman plot of paired differences against means can be by using the of a two measurements from each method in of the single As with the Bland–Altman plot on a single measurement from each method, an correlation between differences and means could be due to changing bias or a difference in measurement error between the methods. We describe a to whether any bias between measurements made by a method with the as measured by the method. We by the of the two measurements made using the method and we by the measurement made on subject by the measurement method. A plot of against can then be used to investigate whether the bias between the and methods with the as measured by the method. It is that is against rather than against because the a to the common measurement error in and It is for the same that plotting paired differences between measurements from two methods against the measurements from one method is in method comparison between values of and that the bias between the methods with the true value as measured by the method. A statistical of this can be by a for with as the The for the that the in this from zero the of against the of constant If of bias changing with true value is we may be in the at which the bias between the methods We can estimate this by the estimate from the above model by our estimate of the This quantity is an estimate of the in bias method for a in the true value as measured by the method. The choice to against is We could have against and have a different and estimate of the at which the bias This choice is because we will two different and because it means the is statistically The model underlying is a type of and as such can be fitted using software for The is to and for measurements are made by observers or raters. measurements made by two different observers are than are two measurements made by the same observer. as with two methods, measurements from two observers may differ due to bias between the observers observer and their measurement errors may also have different SDs. For example, measurements from an observer who can make more measurements will have a SD than made by a observer. If measurements in the future are to be made by different observers, we to describe and quantify the differences between such measurements in order to whether differences are genuine or may be due to measurement error. As with a method comparison study, the way to study this is for each observer to make at least two measurements of a sample of subjects. The of such a study and the type of statistical analysis that are should be by whether interest lies in a particular of observers, or whether we are in a population of observers or raters. If future measurements are to be made by a of observers, the measurement error study should each of observers making measurements of a sample of subjects, with at least two measurements per observer per subject. We can then use exactly the same methods as for method different observers as different measurement methods. If each observer two or more measurements on each we can whether there is bias or changing bias between observers. We can the bias of future measurements (in by measurements using the corresponding estimated We can also estimate repeatability and reliability for each using the one-way ANOVA model (see It may be of interest to which observers are more and if differences in reliability can be related to observer such as of or If we are to that biases between observers are we can fit a so-called model to such a for a subject and observer The estimates indicate the and for bias between observers. such models can be to the measurement error SD to differ between observers. The observers in a measurement error study can often be considered as a sample of observers from a population of observers who may be used in future studies or clinical In this we are not in the particular observers in the measurement error study, but only in the that they the population of observers. In this it is important that a number of observers is used in the measurement error study. For example, if the measurement error study only two observers, can be the population of observers, because we have a sample of just Because the observers in the measurement error study are considered a sample from the population of observers, we such a study using a model that the observer as a We use a with subject and observer Fitting a model estimates of the between-subject SD, SD and the measurement error SD to be the same for different The SD variability in measurements due to observers, or biases between the observers. such estimates one can calculate the SD of the difference in two measurements, made by the same observer or by two different observers. the SD is the will be than the This is because biases between observers to make measurements from different observers it is more difficult to distinguish between subjects on the basis of measurements made by two different observers than if the subjects been measured by the same observer. In this paper we have distinguished between the concepts of agreement and parameter is to the as they describe different characteristics of the measurement The choice of what to report in a particular study should be by how measurements are to be used in the and also by the fact that may to use a measurement method in a different We have the fact that the reliability of a method depends not only on the of the measurement errors but on the of true values in the population in which measurements are made. As measurement techniques potentially may be used in a variety of clinical and different it is to report estimates of and within-subject SDs. We have which methods we are for the analysis of repeatability studies, method comparison studies, and studies with measurements made by different observers or raters. In we not a single reliability coefficient should be used for method comparison If the reliability of two methods are to be each reliability should be estimated by making at least two measurements on each subject with each measurement method. We have how an between paired differences and means may not be by changing bias between two methods. Such an may also be by a difference in the measurement error but with only one measurement per subject per method it is not to which of is the comparison studies should make at least two measurements per subject per method. This an of the of any between paired differences and means for measurements made by the two methods, and also the repeatability and reliability of each method to be estimated. measurements an observer or measurement error studies must use an number of observers if interest lies in making a population of observers. we results for the of CIs for two in the one-way ANOVA in which two measurements are from each

In this paper, we present a variance minimization (VM) procedure for rare event simulation in tandem queueing networks. We prove that the VM method can produce a zero variance. The VM method is suitable to compute optimal importance sampling (IS) parameters for small scale tandem networks. For large scale tandem networks we propose a sub-optimal IS (SOIS) method, which projects the optimal biased transition probabilities of the corresponding small scale system into those of a large scale system. In other words, we establish an efficient IS method for a large scale system by zooming into a small scale system and then projecting our findings into the large scale system. The numerical results show that our SOIS method can produce accurate results with very short CPU time, while many other methods often require much longer.

When simulating small probabilities, say of order 10 -6 or less, by importance sampling, an established principle is to choose the importance sampling distribution as close to the conditional distribution given the rare event as possible. Implementing this often leads into large deviations calculations and exponential change of measure. We survey some of the standard examples where this approach works and supplement existing counterexamples with new ones. Difficulties often arise as consequence of reflecting barriers and we present an algorithm which at least in simple cases is able to deal with this problem. Also the case of heavy-tailed distributions is considered

Introduction Static problems Markov chains Algorithms References

Safety-critical software often includes some rarely executed critical operations that are often inadequately tested in statistical testing based reliability estimation. However, it is necessary to assess the reliability of safety-critical software to a high degree of confidence before it is deployed in the field. This paper discusses how to reduce safety-critical software statistical testing cost based on importance sampling technique. When both the critical operations and the entire software are adequately tested, the method can still compute the unbiased software reliability from the test results with much less test cases. Thus, the statistical testing cost of safety-critical software can be reduced effectively The simulated annealing algorithm for calculating optimum transition probabilities of the Markov chain usage model for reducing software statistical testing cost is also presented.