Published: 14 October 2013

Size: px
Start display at page:

Download "Published: 14 October 2013"

Transcription

1 Electronic Journal of Alied Statistical Analysis EJASA, Electron. J. A. Stat. Anal. htt://siba-ese.unisalento.it/index.h/ejasa/index e-issn: DOI: /i275948v6n213 Estimation of Parameters of Exonentiated Pareto Distribution for Progressive Tye-II Censored Data with Binomial Random Removals Scheme By Singh et al. Published: 14 October 213 This work is coyrighted by Università del Salento, and is licensed under a Creative Commons Attribuzione - Non commerciale - Non oere derivate 3. Italia License. For more information see: htt://creativecommons.org/licenses/by-nc-nd/3./it/

2 Electronic Journal of Alied Statistical Analysis Vol. 6, Issue 2, 213, DOI: /i275948v6n213 Estimation of Parameters of Exonentiated Pareto Distribution for Progressive Tye-II Censored Data with Binomial Random Removals Scheme Sanjay Kumar Singh a, Umesh Singh a, Manoj Kumar a, and G.P Singh b a Deartment of Statistics and DST-CIMS, Bananas Hindu University,Varanasi-2215 b Deartment of Community Medicine, I.M.S and DST-CIMS, Bananas Hindu University, Varanasi-2215 Published: 14 October 213 In this aer, we roose maximum likelihood estimators and Bayes estimators of arameters of exonentiated Pareto distribution under general Entroy loss function and squared error loss function for Progressive tye-ii censored data with binomial removals. The maximum likelihood estimators and corresonding Bayes estimators are comared in terms of their risks based on simulated samles from exonentiated Pareto distribution keywords: Maximum likelihood estimators, exonentiated Pareto distribution under general Entroy loss function, Progressive tye-ii censored data with binomial removals. 1 Introduction In life testing exeriments, situations do arise when units are lost or removed from the exeriments while they are still alive; i.e, we get censored data from the exeriment. The loss of units may occur due to time constraints giving tye-i censored data. In such censoring scheme, exeriment is terminated at secified time. Sometimes, the exeriment is terminated after a refixed number of observations due to cost constraints Corresonding authors: manustats@gmail.com c Università del Salento ISSN: htt://siba-ese.unisalento.it/index.h/ejasa/index

3 Electronic Journal of Alied Statistical Analysis 131 and we get tye-ii censored data. Besides these two controlled causes, units may dro out of the exeriment randomly due to some uncontrolled causes. For examle, consider that a doctor erform an exeriment with n cancer atients but after the death of first atient, some atient leave the exeriment and go for treatment to other doctor/ hosital. Similarly, after the second death a few more leave and so on. Finally the doctor stos taking observation as soon as the redetermined number of deaths (say m are recorded. It may be assumed here that each stage the articiating atients may indeendently decide to leave the exeriment and the robability ( of leaving the exeriment is same for all the atients. Thus the number of atients who leave the exeriment at a secified stage will follow binomial distribution with robability of success (. The exeriment is similar to a life test exeriment which starts with n units. At the first failure X 1, r 1 (random units are removed from the remaining (n 1 surviving units. At second failure X 2, r 2 unit from remaining n 2 r 1 units are removed, and so on, till m th failure is observed; i.e., at m th failure all the remaining r m = n m r 1 r 2 r m 1 units are removed. It may be re-emhasized that, here, m is re-fixed constant and r, is are random. Such a censoring mechanism is termed as rogressive tye-ii censoring with random removal scheme. As stated above, if we assume that robability of removal of a unit, at every stage, is for each unit then r i can be considered to follow a Binomial distribution; i.e, r i B(n m i 1 l= r l, for i = 1, 2, 3, m 1 with r =. For further details, reader are referred to Balakrishnan (27. In last few years, the estimation of arameters of different life time distribution based on rogressive censored samles have been studied by several authors such as Childs and Balakrishnan (2, Balakrishnan and Kannan (21, Mousa and Jaheen (22, Chan and Balakrishnan (22, Sarhan and Abuamooh (28, Ashour and Afify (27 and Ashour and Afify (28. The rogressive tye-ii censoring with binomial removal has been considered by Tse, Yang and Yuen (2 for Weibull distribution, Wu and Chang (22 for Exonential distribution. Under the rogressive tye-ii censoring with random removals, Wu and Chang (23 and Yuen and Tse (1996 develoed the estimation roblem for the Pareto distribution and Weibull distribution resectively, when the number of units removed at each failure time has a discrete uniform distribution, the exected time of this censoring lan is discussed and comared numerically. Let us assume that the life time of the units follow the exonentiated Pareto distribution (EPD with cumulative distribution function F (x, α, = [ 1 (1 + x α] ; x >, α >, > (1 This distribution was introduced by Guta, Guta and Guta (1998. The robability density function (df of X takes the following form with two shae arameters α and : f(x, α, = α [ 1 (1 + x α] 1 (1 + x (α+1 ; x >, α >, > (2 The corresonding survival function is S(x = 1 F (x, α, = 1 [ 1 (1 + x α] ; x >, α >, > (3

4 132 Singh et al. It may be noted that for = 1, (2 reduces to f 1 (x, α, = α(1 + x (α+1 ; x >, α > (4 which is standard Pareto distribution of second kind Kotz and Balakrishnan (1994. The estimators of the arameters of exonentiated Pareto distribution have been obtained by Shawky, Abu-Zinadah and Hanna (29 under different estimation rocedures for comlete samle case. The estimation of arameters has also been attemted by Afify (21 under tye-i and tye-ii censoring scheme. However, no attemt has been made to develo estimators for the arameters of EPD under rogressive tye-ii censoring with binomial removals. Therefore, we roose to develo such an estimation rocedure. Rest of the aer is organized as follows: Section 2 rovides the likelihood function. In section 3, maximum likelihood estimator (MLE and Bayes estimators have been obtained. An algorithm for simulating the rogressive tye-ii samle with binomial removal is resented in section 4. For illustration urose estimates have been obtained for a simulated data in section 5. The comarison of MLE s and corresonding Bayes estimators are given in section 6. Comarisons are based on simulation studies of risk (average loss over samle sace of the estimators. Finally conclusions are rovided in the last section. 2 Likelihood Function Let (X 1, R 1, (X 2, R 2, (X 3, R 3,, (X m, R m, denote a random rogressive tye II censored samle from (2, where X 1 < X 2 < X 3, < X m. For given number of removals, say R 1 = r 1, R 2 = r 2, R 3 = r 3,, R m = r m, the conditional likelihood function can be written as (see Cohen (1963: m L(α, ; x R = r = c f (x i [S (x i ] r i (5 where c = n(n r 1 1(n r 1 r 2 2(n r 1 r 2 r 3 3 (n r 1 r 2 r 3,, r m m+1, and r i (n m r 1 r 2 r 3 r i 1, for i = 1, 2, 3..., m 1 and r =. Substituting (2 and (3 into (5, we get L(α, ; x R = r = c m α [ 1 (1 + x i α] 1 { 1 [ 1 (1 + x i α] } r i (1 + xi (α+1 (6 We assume that units removed from the test at i th failure, i = 1, 2,, m 1, are indeendent of each other. However, robability of removal of unit remains same (say, for all units at all failures; i.e., R i (the number of units removed at the i th failure i = 1, 2, 3,, m 1 follows binomial distribution with arameters n m i 1 l= r l and. Therefore, ( n m P (R 1 = r 1 = r 1 (1 n m r 1, (7 r 1

5 Electronic Journal of Alied Statistical Analysis 133 and for i = 2, 3,, m 1, P (R i = r i R i 1 = r i 1, R 1 = r 1 = ( n m i 1 l= r l r i (1 n m i 1 l= r l (8 r i We further assume that R i is indeendent of X i for all i. Hence the likelihood function takes the following form where L (α,, ; x = L (α, ; x R = r P (R = r, (9 P (R = r = P (R 1 = r 1 P (R 2 = r 2 R 1 = r 1 P (R 3 = r 3 R 2 = r 2, R 1 = r 1 P (R m 1 = r m 1 R m 2 = r m 2, R 1 = r 1. (1 Substituting (7 and (8 into (1, we get P (R = r = (n m! m 1 r i (1 (m 1(n m m 1 (m ir i ( n m i 1 l=1 r i! m 1 r i! (11 Using (6,(9 and (11, we can write the likelihood function in the following form : where A = and c (n m! (n m i 1 l=1 r i! m 1 r i! L 1 (α; = L (α,, ; x = AL 1 (α, L 2 (, (12, m α [ 1 (1 + x i α] { 1 1 [ 1 (1 + x i α] } r i (1 + xi (α+1 (13 L 2 ( = m 1 r i (1 (m 1(n m m 1 (m ir i. (14 3 Classical and Bayesian Estimation of Parameters 3.1 Maximum Likelihood Estimation In this section, we have obtained the MLE of the arameters, α and based on rogressive tye-ii censored data with binomial removals. We observe from (12, (13 and (14 that likelihood function is multilication of three terms, namely, A, L 1 and L 2. Out of these, A does not deendent on the arameters α, and ; thus, it behaves as constant for maximum likelihood estimation. L 1 does not involved and can be treated as function of α and only, where as L 2 involves only. Therefore, the MLE s of α and can be derived by maximizing L 1 with resect to α and. Similarly, the MLE of can be obtained by maximizing L 2.

6 134 Singh et al. Taking log of both sides of (13, we have ln L 1 (α; = m ln ( + m ln (α m + ( 1 ln [ 1 (1 + x i α] + m r i ln [1 [ 1 (1 + x i α] ] (α + 1 m ln (1 + x i Thus, the normal equations can be obtained by differentiating (15 with resect to α and and equating these to zero; i.e., MLE s ˆα and ˆ of α and resectively, can be obtained by simultaneously solving the following normal equations: and m m α + ( 1 m m (1 + x i α ln (1 + x i 1 (1 + x i α m ln (1 + x i r i [ 1 (1 + xi α] 1 (1 + xi α ln (1 + x i 1 [ 1 (1 + x i α] = m m + ln [ 1 (1 + x i α] r i [ 1 (1 + xi α] ln [ 1 (1 + xi α] 1 [ 1 (1 + x i α] = It may be noted that (16 and (17 can not be solved simultaneously to rovide a nice closed form for the estimators. Therefore, we roose to use fixed oint iteration method for solving these equations numerically. For details about the roosed method readers may refer Jain, Iyengar and Jain (1984. When rocedure for obtaining the iteration function and the choice of initial guesses, based on maximum absolute row sum norms, have been discussed. The log of L 2 ( takes the following form [ m 1 ln L 2 ( = ln r i + ln (1 (m 1 (n m m 1 r i (m i ] (15 (16 (17. (18 The first order derivative of ln L 2 ( with resect to is ln L 2 ( m 1 = r i (m 1 (n m m 1 r i (m i. (19 1

7 Electronic Journal of Alied Statistical Analysis 135 =, we get the normal equation for. Solving this equation for, we get the MLE of as m 1 ˆ M = r i (m 1 (n m m 1 r i (m i 1. (2 Setting ln L 2( 3.2 Bayes rocedure In this section, we have obtained the Bayes estimators of the arameters α, and based on rogressively tye-ii censored data with Binomial removals. In order to obtain the Bayes estimator, we must assume that the arameters α, and are random variables. We further assume that these are indeendently distributed. The random variables α and have non-informative rior distribution with resective rior dfs g 1 (α = 1 c ; < α < c (21 and g 2 ( = 1 ; > (22 where as has Beta distribution of first kind with known arameters a, b. The rior df of is given by g 3 ( = 1 B (a, b a 1 (1 b 1 ; < < 1, a >, b > (23 Based on the assumtions stated above, the joint rior df of α, and is g (α,, = g 1 (α g 2 ( g 3 ( ; < α < c, >, < < 1 (24 Combining the riors given by (21, (22 and (23 with likelihood given by (12, we can easily obtain joint osterior df of (α,, as π (α,, x, r = J 1 J where J 1 = α m m 1 m [ 1 (1 + xi α] { 1 1 [ 1 (1 + x i α] } r i (1 + x i (α+1 m 1 ri+a 1 (1 (m 1(n m m 1 (m ir i+b 1 and J = 1 c J 1dαdd Hence, the resective marginal osterior df s of α, and are given by π 1 (α x, r = π 2 ( x, r = 1 1 c J 1 J dd (25 J 1 J dαd (26

8 136 Singh et al. and π 3 ( x, r = c J 1 J dαd (27 Usually the Bayes estimators are obtained under Square Error Loss Function (SELF l 1 (φ, ˆφ = 1 (φ ˆφ 2 ; 1 > (28 Where ˆφ is the estimate of the arameter φ and the Bayes estimator ˆφ S of φ comes out to be E φ [φ], where E φ denotes the osterior exectation. This loss function is a symmetric loss function and can only be justified, if over estimation and under estimation of equal magnitude are of equal seriousness. A number of asymmetric loss functions are also available in statistical literature. Let us consider the General Entroy Loss Function (GELF, roosed by Calabria and Pulcini (1996, defined as follows : ( δ ( l 2 (φ, ˆφ = 2 ˆφ ˆφ δ ln 1 ; 2 > (29 φ φ The constant δ, involved in (29, is its shae arameter. It reflects dearture from symmetry. When δ >, it considers over estimation (i.e., ositive error to be more serious than under estimation (i.e., negative error and converse for δ <. The Bayes estimator ˆφ E of φ under GELF is given by, ˆφ E = [E φ (φ δ] ( 1 δ (3 rovided the osterior exectation exits. It may be noted here that for δ = 1, the Bayes estimator under loss (28 coincides with the Bayes estimator under SELF l 1. Exressions for the Bayes estimators ˆα E, ˆ E and ˆ E for α, and resectively under GELF can be given as [ c ] 1 ( αˆ E = α δ δ π 1 (α x, r dα (31 [ ] 1 ( ˆ E = δ δ π 1 ( x, r d (32 and [ 1 ]( 1 ˆ E = δ δ π 3 ( x, r d (33 Substituting the osterior dfs from (25, (26 and (27 in (31, (32 and (33 resectively and then simlifying, we get the Bayes Estimators ˆα E, ˆ E and ˆ E of α, and as follows : ˆα E = m 1 c αm δ m m 1 c αm m [ 1 (1 + xi α] { 1 1 [ 1 (1 + x i α] } r i (1 + xi (α+1 dαd [ 1 (1 + xi α] { 1 1 [ 1 (1 + x i α] } r i (1 + xi (α+1 dαd 1/δ (34,

9 Electronic Journal of Alied Statistical Analysis 137 ˆ E = and m δ 1 c αm m m 1 c αm m [ 1 (1 + xi α] { 1 1 [ 1 (1 + x i α] } r i (1 + xi (α+1 dαd [ 1 (1 + xi α] { 1 1 [ 1 (1 + x i α] } r i (1 + xi (α+1 dαd [ B( m 1 ˆ E = r i + a δ, (m 1(n m ] m 1 1/δ (m ir i + b B( m 1 r i + a, (m 1(n m m 1 (m ir. (36 i + b It may be noted that the integrals involved in the exressions for the Bayes estimators ˆα E and ˆ E are not reducible in nice closed form. Therefore, we roose the use of numerical integration method for obtaining the estimates. We have used Gauss-quadrature formula for evaluating the inner integral which has finite range and for evaluating the outer integral having range (,, we used Gauss-Laguerre formula. The evaluations were done by develoing rogram in R-language. 4 Algorithm to Simulate Progressive Tye-II Censored Samle With Binomial Removal For the study of behavior of the estimators obtained in revious sections, we need to simulate rogressive tye-ii censored samles with Binomial removals from secified EPD. To get such a samle, we roose the use of following algorithm : I. Secify the value of n. II. Secify the value of m. III. Secify the value of arameters α, and. IV. Generate a random samle (S r of size n from EP (α,. V. Generate random number r i ( from B n m i 1 l= r l,, at i th stage for i = 1, 2, 3,, m 1. (r = VI. Get ordered samle S o from S r to choose the minimum which will be first observation in desired rogressive tye-ii censored samle S. VII. Dro the observation selected at VI from S r to have a random samle S r of size n (less 1 than that of S r. VIII. Generate r i integers (at i th stage between 1 to n and observations corresonding to these numbers are droed from S r to have a random samle S r of size n = n r i and re-designate the random samle in hand as new random samle S r. 1/δ (35,

10 138 Singh et al. IX. Reeat stes V to VIII (m 1 times. X. Set r m according to the following relation. { n m m 1 l=1 r m = r l if n m m 1 l=1 r l > otherwise and discard all the remaining r m observations. 5 A Numerical Illustration For illustration of the roosed estimation rocedures given in sections 2 and 3, we have considered below a simulated rogressive tye-ii censored samle with Binomial removals from exonentiated Pareto distribution with =.3, m = 14, α = 2, and =.5, using the above algorithm. The samle thus obtained is given below : ( ,2, ( ,1, ( ,, ( ,, ( ,, ( ,, ( ,1, ( ,1, ( ,, ( ,, ( ,, ( ,1, ( ,, ( ,. Substituting value of (x i, r i i = 1,, 14 in (16 and (17, we obtain the normal equations which are solved simultaneously using fixed oint iteration rocedure to get the MLE s of α and. The MLE of is obtained from (2. The Bayes Estimates are obtained from (36, using the rocedure exlained in section 3. The results are summarized in table 1. It may be noted from the table that for the samle in hand, all the estimates are larger than the corresonding true values. However the Bayes estimates ˆα E, ˆ E and ˆ E with δ = 1.5 is nearer to the true value as comared to other estimates. For δ = 1.5, the Bayes estimate ˆα S, ˆ S and MLE ˆ M are nearer to the true as comared to other estimates. But on the basis of this articular samle, it will be illogical to conclude that under GELF, we should use ˆα E, ˆ E and ˆ E if over estimation is more serious than under estimation for their long run use. Therefore, we study the behavior of risks of the estimators in the next section. Table 1: MLE and Bayes estimate of α, and under rogressive tye-ii censored data with random removal. δ α ˆα M ˆα S ˆα E ˆM ˆS ˆE ˆ M ˆ S ˆ E

11 Electronic Journal of Alied Statistical Analysis Simulation Studies The estimators ˆα M, ˆ M and ˆ M denote the MLE s of the arameters α, and resectively while ˆα S, ˆ S and ˆ S are corresonding Bayes estimators under SELF and ˆα E, ˆ E and ˆ E are the corresonding Bayes estimators under GELF. We comare the estimators obtained under GELF with corresonding MLE s and Bayes estimators under SELF. The comarisons are based on the simulated risks (average loss over samle sace under GELF. It may be mentioned here that the exact exressions for the risks can not be obtained because estimators are not in nice closed form. Therefore, the risks of the estimators are estimated on the basis of Monte-carlo simulation study of 1 samles. It may be noted that the risks of the estimators will deend on values of n, m,, α,, c, δ, a and b. In order to consider variation in the values of these, we have obtained the simulated risks for n = 2 [1] 4, m = 12 [2] 18,c = 4 [1] 8, a = 2 [2] 8, b = 4 [2] 1, =.5 [.5] 2, α =.5 [.5] 2.5, δ = ±1.5 and =.3 [.1].6. Generating the rogressive samle as mentioned in section 4, the simulated risks under GELF have been obtained for selected values of n, m,, α,, b, a, δ and c. The results are summarized in (tables ( Discussion of the results It is interesting to note that as n increases, keeing the effective samle size m fixed, the risks of all the estimators increase, in general, excet for ˆα M and ˆ M (see tables 2 and 8. On other hand, if the effective samle size m increases for given samle size n, the risks of the estimators of α and decrease as m increases, in general, excet for the estimators ˆα M and ˆα S (when δ = 1.5. We have noticed that the risks of all the estimators of increase, often, as m increases (see tables 3 and 9. Due to the increase in the value of, the risks of the estimators of α and, in general, first decrease then increase. But the risks of the estimators of do not follow a trend (see table 4 and 1. Further, for the variation in the value of arameter α (when δ = 1.5, we observe that risks of all estimators of α decrease as α increases but risk of estimators of and do not follow a definite trend (see table 5. On the other hand, for δ = 1.5, risks of the estimator ˆα S and ˆα E decrease as α increases but for the rest of the estimators no definite trend is observed (see table 11. The risk of estimators of (when δ = 1.5 often decrease as increases and the risk of estimators ˆα M and ˆ M increases as increases but no definite trend for the risks of rest of the estimators of α and is noted as increases (see table 6. For δ = 1.5, risk of estimators of and increases as increases, but the risks of estimators of α first decrease then increase when increases (see table 12. For the variation of the rior arameter c, when δ = 1.5, we observe that no definite trend is followed by the risks of the estimators of α, and excet for ˆα S which increases as c increases (see table 7. But, for δ = 1.5, the risk of estimators of α increase as c increases and for rest of estimators no definite trend is observed (see table 13.

12 14 Singh et al. It is worth while to remark here that a variation in the values of hyer arameters a and b do not effect the risks of the estimators of α and. It is further seen that for variation in hyer arameter a, the risks of estimators of do not have a definite trend, excet for the risk of ˆ S (when δ = 1.5 and ˆ E (when δ = 1.5 which increase as a increases (see table 14. For the variation of arameter b, the risk of estimators of, excet for ˆ M (for δ = 1.5 and ˆ S (for δ = 1.5 decrease for increase in the value of b in the beginning and then increase for further increase in the value of b (see table 15. It may further be noted from the tables that for δ = 1.5, the risk of ˆ E is smaller than the risks of ˆ M and ˆ S. However, for δ = 1.5, the risk of ˆ M is smaller than those of ˆ E and ˆ S. Among the estimator of α, ˆα E has smaller risks than those of others for δ = 1.5, but for δ = 1.5, none of the estimators ˆα M, ˆα S and ˆα E outerforms the other although there is little variation in the risk. The risk of ˆ E when δ = 1.5 whereas ˆ S when δ = 1.5 is smaller than the risks of other estimators of. 7 Conclusions It is exected that the estimator obtained under a articular loss function shall in general erform better than the estimators obtained under other loss functions, but this could not be established in the resent study. We have seen above that risk under GELF for the estimators ˆα E, ˆ E and ˆ E is not always less than those of ˆα M, ˆ M, ˆ M, ˆα S, ˆ S, ˆ S. Although the risks associated with ˆ E is smaller than the risk associated with other estimators in most of the cases. The risks associated with ˆα E and ˆ E are smaller than those of other estimators when δ = 1.5. On other hand if δ = 1.5 risk associated with ˆα M and ˆ S are noted to be smaller than other estimators. Therefore, for the use of the roosed estimators, we may recommend the following: 1. ˆE may be used as an estimator of when over estimation is more serious than under estimation. On the other hand, if under estimation is more serious than over estimation, ˆ M may be used. 2. ˆα E may be used as an estimator of α when under estimation is more serious than over estimation. Otherwise, ˆα M may be used. 3. ˆ E may be used an estimator of when under estimation is more serious than over estimation; Otherwise, ˆ S may be used. Acknowledgements The authors are grateful to the editor and the learned referee for making useful comments which resulted in much imrovement in the aer. The second author is also grateful to the UGC-BHU, Varanasi-2215, India for roviding financial assistance.

13 Electronic Journal of Alied Statistical Analysis 141 Table 2: Risk of Estimators of α, and under GELF for fixed, m = 14, α = 2, =.5, =.3, c = 5, a = 2, b = 6, δ = 1.5. RE (ˆαM RE (ˆαS RE (ˆαE RE RE RE RE (ˆM RE (ˆS RE (ˆE n α (ˆM (ˆS (ˆE Table 3: Risk of Estimators of α, and under GELF for fixed, n = 2, α = 2, =.5, =.3, c = 5, a = 2, b = 6, δ = 1.5. RE (ˆαM RE (ˆαS RE (ˆαE RE RE RE RE (ˆM RE (ˆS RE (ˆE m α (ˆM (ˆS (ˆE

14 142 Singh et al. Table 4: Risk of Estimators of α, and under GELF for fixed, m = 14, α = 2, n = 2, =.3, c = 5, a = 2, b = 6, δ = 1.5. RE (ˆαM RE (ˆαS RE (ˆαE RE RE RE RE (ˆM RE (ˆS RE (ˆE α (ˆM (ˆS (ˆE Table 5: Risk of Estimators of α, and under GELF for fixed, n = 2, m = 14, =.5, =.3, c = 5, a = 2, b = 6, δ = 1.5. α α RE (ˆαM RE (ˆαS RE (ˆαE RE (ˆM RE (ˆS RE (ˆE RE (ˆM RE (ˆS RE (ˆE

15 Electronic Journal of Alied Statistical Analysis 143 Table 6: Risk of Estimators of α, and under GELF for fixed, n = 2, m = 14, α = 2, =.5, c = 5, a = 2, b = 6, δ = 1.5. RE (ˆαM RE (ˆαS RE (ˆαE RE RE RE RE (ˆM RE (ˆS RE (ˆE α (ˆM (ˆS (ˆE Table 7: Risk of Estimators of α, and under GELF for fixed, n = 2, m = 14, α = 2, =.5, =.3, a = 2, b = 6, δ = 1.5. c α RE (ˆαM RE (ˆαS RE (ˆαE RE (ˆM RE (ˆS RE (ˆE RE (ˆM RE (ˆS RE (ˆE

16 144 Singh et al. Table 8: Risk of Estimators of α, and under GELF for fixed, m = 13, α = 2, =.5, =.3, c = 5, a = 2, b = 6, δ = 1.5. n α RE (ˆαM RE (ˆαS RE (ˆαE RE (ˆM RE (ˆS RE (ˆE RE (ˆM RE (ˆS RE (ˆE Table 9: Risk of Estimators of α, and under GELF for fixed, n = 2, α = 2, =.5, =.3, c = 5, a = 2, b = 6, δ = 1.5. m α RE (ˆαM RE (ˆαS RE (ˆαE RE (ˆM RE (ˆS RE (ˆE RE (ˆM RE (ˆS RE (ˆE

17 Electronic Journal of Alied Statistical Analysis 145 Table 1: Risk of Estimators of α, and under GELF for fixed, m = 14, α = 2, n = 2, =.3, c = 5, a = 2, b = 6, δ = 1.5. α RE (ˆαM RE (ˆαS RE (ˆαE RE (ˆM RE (ˆS RE (ˆE RE (ˆM RE (ˆS RE (ˆE Table 11: Risk of Estimators of α, and under GELF for fixed, n = 2, m = 14, =.5, =.3, c = 5, a = 2, b = 6, δ = 1.5. α α RE (ˆαM RE (ˆαS RE (ˆαE RE (ˆM RE (ˆS RE (ˆE RE (ˆM RE (ˆS RE (ˆE

18 146 Singh et al. Table 12: Risk of Estimators of α, and under GELF for fixed, n = 2, m = 14, α = 2, =.5, c = 5, a = 2, b = 6, δ = 1.5. α RE (ˆαM RE (ˆαS RE (ˆαE RE (ˆM RE (ˆS RE (ˆE RE (ˆM RE (ˆS RE (ˆE Table 13: Risk of Estimators of α, and under GELF for fixed, n = 2, m = 14, α = 2, =.5, =.3, a = 2, b = 6, δ = 1.5. c α RE (ˆαM RE (ˆαS RE (ˆαE RE (ˆM RE (ˆS RE (ˆE RE (ˆM RE (ˆS RE (ˆE

19 Electronic Journal of Alied Statistical Analysis 147 Table 14: Risk of Estimators of α, and under GELF for fixed, n = 2, m = 14, α = 2, =.5, =.3, c = 5, b = 6, δ = ±1.5. a δ = 1.5 δ = 1.5 R E (ˆ M R E (ˆ S R E (ˆ E R E (ˆ M R E (ˆ S R E (ˆ E Table 15: Risk of Estimators of α, and under GELF for fixed, n = 2, m = 14, α = 2, =.5, =.3, c = 5, a = 2, δ = ±1.5. b δ = 1.5 δ = 1.5 R E (ˆ M R E (ˆ S R E (ˆ E R E (ˆ M R E (ˆ S R E (ˆ E References Afify, W. M (21. On estimation of the exonentiated areto distribution under different samle scheme. Alied Mathematical Sciences, 4(8, Ashour, S. K. and Afify, W. M. (27. Statistical analysis of exonentiated Weibull Family under tye-i Progressive Interval Censoring with random removals. Journal of Alied Sciences Research, 3(12, Ashour, S. K. and Afify, W. M. (28. Estimations of The Parameters of Exonentiated Weibull Family with tye-ii Progressive Interval Censoring with random removals. Journal of Alied Sciences Research, 4(11, Balakrishnan, N (27. Progressive methodology : An araisal (with discussion. Test, 16 (2, Balakrishnan, N. and Kannan, N. (21. Point and Interval Estimation for Parameters of the Logistic Distribution Based on Progressively tye-ii Censored Samles, in Handbook of Statisticsm N. Balakrishnan and C. R. Rao, volume 2. Eds. Amsterdam, North- Holand. Calabria, R. and Pulcini, G. (1996. Point estimation under - asymmetric loss functions for life - truncated exonential samles. Commun. statist. Theory meth., 25(3,

20 148 Singh et al. Childs, A. and Balakrishnan, N. (2. Conditional inference rocedures for the lalace distribution when the observed samles are rogressively censored. Metrika, 52, Cohen, A. C (1963. Progressively censored samles in life testing. Technometrics, ages Guta, R. C., Guta, R. D., and Guta, P. L. (1998. Modeling failure time data by lehman alternatives. Commun. Statist. - Theory Meth., 27(4, Jain, M. K., Iyengar, S. R. K., and Jain, R. K., (1984. Numerical Methods for Scientific and Engineering Comutation. New Age International (P Limited, Publishers, New Delhi, fifth edition Johanson, N. L., Kotz, S., and Balakrishnan, N. (1994. Continuous Univariate Distributions, volume 1. Wiley, New York, 2 edition. Mousa, M. and Jaheen,Z. (22. Statistical inference for the burr model based on rogressively censored data. An International Comuters and Mathematics with Alications, 43, Ng, K., Chan, P. S., and Balakrishnan, N. (22. Estimation of arameters from rogressively censored data using an algorithm. Comutational Statistics and Data Analysis, 39, Sarhan, A. M. and Abuamooh, A. (28. (28. Statistical inference using rogressively tye-ii censored data with random scheme. International Mathematical Forum, 35, Shawky, A. I., Abu-Zinadah, and Hanna H. (29. Exonentiated areto distribution : Different method of estimations. Int. J.Contem. Math. Sciences, 4(14, Tse, S. K., Yang, C., and Yuen, H. K. (2. Statistical analysis of Weibull distributed life time data under tye-ii rogressive censoring with binomial removals. Journal of Alied Statistics, 27, Wu, S. J. and Chang, C. T. (22. Parameter estimations based on exonential rogressive tye-ii censored with binomial removals. International Journal of Information and Management Sciences, 13, Yuen, H. K. and Tse, S. K. (1996. Parameters estimation for Weibull distribution under rogressive censoring with random removal. Journal Statis. Comut. Simul, 55, Wu, S. J. and Chang, C. T. (23. Inference in the areto distribution based on rogressive tye-ii censoring with random removales. Journal of Alied Statistics, 3,

BAYESIAN ESTIMATION OF THE EXPONENTI- ATED GAMMA PARAMETER AND RELIABILITY FUNCTION UNDER ASYMMETRIC LOSS FUNC- TION

BAYESIAN ESTIMATION OF THE EXPONENTI- ATED GAMMA PARAMETER AND RELIABILITY FUNCTION UNDER ASYMMETRIC LOSS FUNC- TION REVSTAT Statistical Journal Volume 9, Number 3, November 211, 247 26 BAYESIAN ESTIMATION OF THE EXPONENTI- ATED GAMMA PARAMETER AND RELIABILITY FUNCTION UNDER ASYMMETRIC LOSS FUNC- TION Authors: Sanjay

More information

On Wrapping of Exponentiated Inverted Weibull Distribution

On Wrapping of Exponentiated Inverted Weibull Distribution IJIRST International Journal for Innovative Research in Science & Technology Volume 3 Issue 11 Aril 217 ISSN (online): 2349-61 On Wraing of Exonentiated Inverted Weibull Distribution P.Srinivasa Subrahmanyam

More information

Statistical Inference Using Progressively Type-II Censored Data with Random Scheme

Statistical Inference Using Progressively Type-II Censored Data with Random Scheme International Mathematical Forum, 3, 28, no. 35, 1713-1725 Statistical Inference Using Progressively Type-II Censored Data with Random Scheme Ammar M. Sarhan 1 and A. Abuammoh Department of Statistics

More information

A Comparison between Biased and Unbiased Estimators in Ordinary Least Squares Regression

A 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 information

Estimation of the large covariance matrix with two-step monotone missing data

Estimation 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 information

An Improved Generalized Estimation Procedure of Current Population Mean in Two-Occasion Successive Sampling

An 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 information

Research Note REGRESSION ANALYSIS IN MARKOV CHAIN * A. Y. ALAMUTI AND M. R. MESHKANI **

Research Note REGRESSION ANALYSIS IN MARKOV CHAIN * A. Y. ALAMUTI AND M. R. MESHKANI ** Iranian Journal of Science & Technology, Transaction A, Vol 3, No A3 Printed in The Islamic Reublic of Iran, 26 Shiraz University Research Note REGRESSION ANALYSIS IN MARKOV HAIN * A Y ALAMUTI AND M R

More information

ESTIMATION OF THE RECIPROCAL OF THE MEAN OF THE INVERSE GAUSSIAN DISTRIBUTION WITH PRIOR INFORMATION

ESTIMATION OF THE RECIPROCAL OF THE MEAN OF THE INVERSE GAUSSIAN DISTRIBUTION WITH PRIOR INFORMATION STATISTICA, anno LXVIII, n., 008 ESTIMATION OF THE RECIPROCAL OF THE MEAN OF THE INVERSE GAUSSIAN DISTRIBUTION WITH PRIOR INFORMATION 1. INTRODUCTION The Inverse Gaussian distribution was first introduced

More information

Introduction to Probability and Statistics

Introduction to Probability and Statistics Introduction to Probability and Statistics Chater 8 Ammar M. Sarhan, asarhan@mathstat.dal.ca Deartment of Mathematics and Statistics, Dalhousie University Fall Semester 28 Chater 8 Tests of Hyotheses Based

More information

A MIXED CONTROL CHART ADAPTED TO THE TRUNCATED LIFE TEST BASED ON THE WEIBULL DISTRIBUTION

A 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 information

Slash Distributions and Applications

Slash Distributions and Applications CHAPTER 2 Slash Distributions and Alications 2.1 Introduction The concet of slash distributions was introduced by Kafadar (1988) as a heavy tailed alternative to the normal distribution. Further literature

More information

STA 250: Statistics. Notes 7. Bayesian Approach to Statistics. Book chapters: 7.2

STA 250: Statistics. Notes 7. Bayesian Approach to Statistics. Book chapters: 7.2 STA 25: Statistics Notes 7. Bayesian Aroach to Statistics Book chaters: 7.2 1 From calibrating a rocedure to quantifying uncertainty We saw that the central idea of classical testing is to rovide a rigorous

More information

Information collection on a graph

Information 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 information

Using the Divergence Information Criterion for the Determination of the Order of an Autoregressive Process

Using the Divergence Information Criterion for the Determination of the Order of an Autoregressive Process Using the Divergence Information Criterion for the Determination of the Order of an Autoregressive Process P. Mantalos a1, K. Mattheou b, A. Karagrigoriou b a.deartment of Statistics University of Lund

More information

Information collection on a graph

Information 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 information

Lower Confidence Bound for Process-Yield Index S pk with Autocorrelated Process Data

Lower 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 information

Lecture 23 Maximum Likelihood Estimation and Bayesian Inference

Lecture 23 Maximum Likelihood Estimation and Bayesian Inference Lecture 23 Maximum Likelihood Estimation and Bayesian Inference Thais Paiva STA 111 - Summer 2013 Term II August 7, 2013 1 / 31 Thais Paiva STA 111 - Summer 2013 Term II Lecture 23, 08/07/2013 Lecture

More information

Towards 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. 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 information

arxiv: v1 [physics.data-an] 26 Oct 2012

arxiv: 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 information

Combining Logistic Regression with Kriging for Mapping the Risk of Occurrence of Unexploded Ordnance (UXO)

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 information

On split sample and randomized confidence intervals for binomial proportions

On split sample and randomized confidence intervals for binomial proportions On slit samle and randomized confidence intervals for binomial roortions Måns Thulin Deartment of Mathematics, Usala University arxiv:1402.6536v1 [stat.me] 26 Feb 2014 Abstract Slit samle methods have

More information

Deriving Indicator Direct and Cross Variograms from a Normal Scores Variogram Model (bigaus-full) David F. Machuca Mory and Clayton V.

Deriving 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 information

On the Estimation Of Population Mean Under Systematic Sampling Using Auxiliary Attributes

On the Estimation Of Population Mean Under Systematic Sampling Using Auxiliary Attributes Oriental Journal of Physical Sciences Vol 1 (1 & ) 17 (016) On the Estimation Of Poulation Mean Under Systematic Samling Using Auxiliary Attributes Usman Shahzad Deartment of Mathematics Statistics PMAS

More information

A New Asymmetric Interaction Ridge (AIR) Regression Method

A New Asymmetric Interaction Ridge (AIR) Regression Method A New Asymmetric Interaction Ridge (AIR) Regression Method by Kristofer Månsson, Ghazi Shukur, and Pär Sölander The Swedish Retail Institute, HUI Research, Stockholm, Sweden. Deartment of Economics and

More information

MATHEMATICAL MODELLING OF THE WIRELESS COMMUNICATION NETWORK

MATHEMATICAL 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 information

Bayesian Spatially Varying Coefficient Models in the Presence of Collinearity

Bayesian 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 information

4. Score normalization technical details We now discuss the technical details of the score normalization method.

4. 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 information

Yixi Shi. Jose Blanchet. IEOR Department Columbia University New York, NY 10027, USA. IEOR Department Columbia University New York, NY 10027, USA

Yixi Shi. Jose Blanchet. IEOR Department Columbia University New York, NY 10027, USA. IEOR Department Columbia University New York, NY 10027, USA Proceedings of the 2011 Winter Simulation Conference S. Jain, R. R. Creasey, J. Himmelsach, K. P. White, and M. Fu, eds. EFFICIENT RARE EVENT SIMULATION FOR HEAVY-TAILED SYSTEMS VIA CROSS ENTROPY Jose

More information

Bayesian 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 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 information

CHAPTER-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 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 information

#A64 INTEGERS 18 (2018) APPLYING MODULAR ARITHMETIC TO DIOPHANTINE EQUATIONS

#A64 INTEGERS 18 (2018) APPLYING MODULAR ARITHMETIC TO DIOPHANTINE EQUATIONS #A64 INTEGERS 18 (2018) APPLYING MODULAR ARITHMETIC TO DIOPHANTINE EQUATIONS Ramy F. Taki ElDin Physics and Engineering Mathematics Deartment, Faculty of Engineering, Ain Shams University, Cairo, Egyt

More information

System Reliability Estimation and Confidence Regions from Subsystem and Full System Tests

System 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 information

John Weatherwax. Analysis of Parallel Depth First Search Algorithms

John 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 information

Bayesian Networks Practice

Bayesian Networks Practice Bayesian Networks Practice Part 2 2016-03-17 Byoung-Hee Kim, Seong-Ho Son Biointelligence Lab, CSE, Seoul National University Agenda Probabilistic Inference in Bayesian networks Probability basics D-searation

More information

Availability and Maintainability. Piero Baraldi

Availability and Maintainability. Piero Baraldi Availability and Maintainability 1 Introduction: reliability and availability System tyes Non maintained systems: they cannot be reaired after a failure (a telecommunication satellite, a F1 engine, a vessel

More information

General Linear Model Introduction, Classes of Linear models and Estimation

General 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 information

Hotelling s Two- Sample T 2

Hotelling 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 information

New Information Measures for the Generalized Normal Distribution

New Information Measures for the Generalized Normal Distribution Information,, 3-7; doi:.339/info3 OPEN ACCESS information ISSN 75-7 www.mdi.com/journal/information Article New Information Measures for the Generalized Normal Distribution Christos P. Kitsos * and Thomas

More information

State Estimation with ARMarkov Models

State 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 information

Approximating min-max k-clustering

Approximating 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 information

CHAPTER 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

CHAPTER 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 information

Probability Estimates for Multi-class Classification by Pairwise Coupling

Probability 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

DETC2003/DAC AN EFFICIENT ALGORITHM FOR CONSTRUCTING OPTIMAL DESIGN OF COMPUTER EXPERIMENTS

DETC2003/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 information

Performance of lag length selection criteria in three different situations

Performance of lag length selection criteria in three different situations MPRA Munich Personal RePEc Archive Performance of lag length selection criteria in three different situations Zahid Asghar and Irum Abid Quaid-i-Azam University, Islamabad Aril 2007 Online at htts://mra.ub.uni-muenchen.de/40042/

More information

DEPARTMENT OF ECONOMICS ISSN DISCUSSION PAPER 20/07 TWO NEW EXPONENTIAL FAMILIES OF LORENZ CURVES

DEPARTMENT OF ECONOMICS ISSN DISCUSSION PAPER 20/07 TWO NEW EXPONENTIAL FAMILIES OF LORENZ CURVES DEPARTMENT OF ECONOMICS ISSN 1441-549 DISCUSSION PAPER /7 TWO NEW EXPONENTIAL FAMILIES OF LORENZ CURVES ZuXiang Wang * & Russell Smyth ABSTRACT We resent two new Lorenz curve families by using the basic

More information

MODELING THE RELIABILITY OF C4ISR SYSTEMS HARDWARE/SOFTWARE COMPONENTS USING AN IMPROVED MARKOV MODEL

MODELING 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 information

1. INTRODUCTION. Fn 2 = F j F j+1 (1.1)

1. INTRODUCTION. Fn 2 = F j F j+1 (1.1) CERTAIN CLASSES OF FINITE SUMS THAT INVOLVE GENERALIZED FIBONACCI AND LUCAS NUMBERS The beautiful identity R.S. Melham Deartment of Mathematical Sciences, University of Technology, Sydney PO Box 23, Broadway,

More information

On parameter estimation in deformable models

On parameter estimation in deformable models Downloaded from orbitdtudk on: Dec 7, 07 On arameter estimation in deformable models Fisker, Rune; Carstensen, Jens Michael Published in: Proceedings of the 4th International Conference on Pattern Recognition

More information

Optimal Learning Policies for the Newsvendor Problem with Censored Demand and Unobservable Lost Sales

Optimal Learning Policies for the Newsvendor Problem with Censored Demand and Unobservable Lost Sales Otimal Learning Policies for the Newsvendor Problem with Censored Demand and Unobservable Lost Sales Diana Negoescu Peter Frazier Warren Powell Abstract In this aer, we consider a version of the newsvendor

More information

Tests for Two Proportions in a Stratified Design (Cochran/Mantel-Haenszel Test)

Tests for Two Proportions in a Stratified Design (Cochran/Mantel-Haenszel Test) Chater 225 Tests for Two Proortions in a Stratified Design (Cochran/Mantel-Haenszel Test) Introduction In a stratified design, the subects are selected from two or more strata which are formed from imortant

More information

Monte Carlo Studies. Monte Carlo Studies. Sampling Distribution

Monte Carlo Studies. Monte Carlo Studies. Sampling Distribution Monte Carlo Studies Do not let yourself be intimidated by the material in this lecture This lecture involves more theory but is meant to imrove your understanding of: Samling distributions and tests of

More information

Morten Frydenberg Section for Biostatistics Version :Friday, 05 September 2014

Morten 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 information

Plotting the Wilson distribution

Plotting 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 information

substantial literature on emirical likelihood indicating that it is widely viewed as a desirable and natural aroach to statistical inference in a vari

substantial literature on emirical likelihood indicating that it is widely viewed as a desirable and natural aroach to statistical inference in a vari Condence tubes for multile quantile lots via emirical likelihood John H.J. Einmahl Eindhoven University of Technology Ian W. McKeague Florida State University May 7, 998 Abstract The nonarametric emirical

More information

Casimir Force Between the Two Moving Conductive Plates.

Casimir Force Between the Two Moving Conductive Plates. Casimir Force Between the Two Moving Conductive Plates. Jaroslav Hynecek 1 Isetex, Inc., 95 Pama Drive, Allen, TX 751 ABSTRACT This article resents the derivation of the Casimir force for the two moving

More information

ASYMPTOTIC RESULTS OF A HIGH DIMENSIONAL MANOVA TEST AND POWER COMPARISON WHEN THE DIMENSION IS LARGE COMPARED TO THE SAMPLE SIZE

ASYMPTOTIC RESULTS OF A HIGH DIMENSIONAL MANOVA TEST AND POWER COMPARISON WHEN THE DIMENSION IS LARGE COMPARED TO THE SAMPLE SIZE J Jaan Statist Soc Vol 34 No 2004 9 26 ASYMPTOTIC RESULTS OF A HIGH DIMENSIONAL MANOVA TEST AND POWER COMPARISON WHEN THE DIMENSION IS LARGE COMPARED TO THE SAMPLE SIZE Yasunori Fujikoshi*, Tetsuto Himeno

More information

AI*IA 2003 Fusion of Multiple Pattern Classifiers PART III

AI*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 information

Statistics II Logistic Regression. So far... Two-way repeated measures ANOVA: an example. RM-ANOVA example: the data after log transform

Statistics II Logistic Regression. So far... Two-way repeated measures ANOVA: an example. RM-ANOVA example: the data after log transform Statistics II Logistic Regression Çağrı Çöltekin Exam date & time: June 21, 10:00 13:00 (The same day/time lanned at the beginning of the semester) University of Groningen, Det of Information Science May

More information

Shadow Computing: An Energy-Aware Fault Tolerant Computing Model

Shadow 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 information

Analysis of some entrance probabilities for killed birth-death processes

Analysis of some entrance probabilities for killed birth-death processes Analysis of some entrance robabilities for killed birth-death rocesses Master s Thesis O.J.G. van der Velde Suervisor: Dr. F.M. Sieksma July 5, 207 Mathematical Institute, Leiden University Contents Introduction

More information

Lower bound solutions for bearing capacity of jointed rock

Lower 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 information

Estimating function analysis for a class of Tweedie regression models

Estimating 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 information

Topic: Lower Bounds on Randomized Algorithms Date: September 22, 2004 Scribe: Srinath Sridhar

Topic: Lower Bounds on Randomized Algorithms Date: September 22, 2004 Scribe: Srinath Sridhar 15-859(M): Randomized Algorithms Lecturer: Anuam Guta Toic: Lower Bounds on Randomized Algorithms Date: Setember 22, 2004 Scribe: Srinath Sridhar 4.1 Introduction In this lecture, we will first consider

More information

ECON 4130 Supplementary Exercises 1-4

ECON 4130 Supplementary Exercises 1-4 HG Set. 0 ECON 430 Sulementary Exercises - 4 Exercise Quantiles (ercentiles). Let X be a continuous random variable (rv.) with df f( x ) and cdf F( x ). For 0< < we define -th quantile (or 00-th ercentile),

More information

Research Article A Note on the Modified q-bernoulli Numbers and Polynomials with Weight α

Research Article A Note on the Modified q-bernoulli Numbers and Polynomials with Weight α Abstract and Alied Analysis Volume 20, Article ID 54534, 8 ages doi:0.55/20/54534 Research Article A Note on the Modified -Bernoulli Numbers and Polynomials with Weight α T. Kim, D. V. Dolgy, 2 S. H. Lee,

More information

Scaling Multiple Point Statistics for Non-Stationary Geostatistical Modeling

Scaling Multiple Point Statistics for Non-Stationary Geostatistical Modeling 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

More information

Machine Learning: Homework 4

Machine Learning: Homework 4 10-601 Machine Learning: Homework 4 Due 5.m. Monday, February 16, 2015 Instructions Late homework olicy: Homework is worth full credit if submitted before the due date, half credit during the next 48 hours,

More information

Analysis of M/M/n/K Queue with Multiple Priorities

Analysis 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 information

On Wald-Type Optimal Stopping for Brownian Motion

On Wald-Type Optimal Stopping for Brownian Motion J Al Probab Vol 34, No 1, 1997, (66-73) Prerint Ser No 1, 1994, Math Inst Aarhus On Wald-Tye Otimal Stoing for Brownian Motion S RAVRSN and PSKIR The solution is resented to all otimal stoing roblems of

More information

Optimal Design of Truss Structures Using a Neutrosophic Number Optimization Model under an Indeterminate Environment

Optimal 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 information

arxiv: v3 [physics.data-an] 23 May 2011

arxiv: v3 [physics.data-an] 23 May 2011 Date: October, 8 arxiv:.7v [hysics.data-an] May -values for Model Evaluation F. Beaujean, A. Caldwell, D. Kollár, K. Kröninger Max-Planck-Institut für Physik, München, Germany CERN, Geneva, Switzerland

More information

Collaborative Place Models Supplement 1

Collaborative Place Models Supplement 1 Collaborative Place Models Sulement Ber Kaicioglu Foursquare Labs ber.aicioglu@gmail.com Robert E. Schaire Princeton University schaire@cs.rinceton.edu David S. Rosenberg P Mobile Labs david.davidr@gmail.com

More information

Named Entity Recognition using Maximum Entropy Model SEEM5680

Named Entity Recognition using Maximum Entropy Model SEEM5680 Named Entity Recognition using Maximum Entroy Model SEEM5680 Named Entity Recognition System Named Entity Recognition (NER): Identifying certain hrases/word sequences in a free text. Generally it involves

More information

Pairwise active appearance model and its application to echocardiography tracking

Pairwise 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 information

Infinite Number of Twin Primes

Infinite Number of Twin Primes dvances in Pure Mathematics, 06, 6, 95-97 htt://wwwscirorg/journal/am ISSN Online: 60-08 ISSN Print: 60-068 Infinite Number of Twin Primes S N Baibeov, Durmagambetov LN Gumilyov Eurasian National University,

More information

Preconditioning techniques for Newton s method for the incompressible Navier Stokes equations

Preconditioning techniques for Newton s method for the incompressible Navier Stokes equations Preconditioning techniques for Newton s method for the incomressible Navier Stokes equations H. C. ELMAN 1, D. LOGHIN 2 and A. J. WATHEN 3 1 Deartment of Comuter Science, University of Maryland, College

More information

A New Perspective on Learning Linear Separators with Large L q L p Margins

A New Perspective on Learning Linear Separators with Large L q L p Margins A New Persective on Learning Linear Searators with Large L q L Margins Maria-Florina Balcan Georgia Institute of Technology Christoher Berlind Georgia Institute of Technology Abstract We give theoretical

More information

Maximum Entropy and the Stress Distribution in Soft Disk Packings Above Jamming

Maximum Entropy and the Stress Distribution in Soft Disk Packings Above Jamming Maximum Entroy and the Stress Distribution in Soft Disk Packings Above Jamming Yegang Wu and S. Teitel Deartment of Physics and Astronomy, University of ochester, ochester, New York 467, USA (Dated: August

More information

The Mathematics of Winning Streaks

The Mathematics of Winning Streaks The Mathematics of Winning Streaks Erik Leffler lefflererik@gmail.com under the direction of Prof. Henrik Eriksson Deartment of Comuter Science and Communications Royal Institute of Technology Research

More information

arxiv: v1 [hep-ex] 8 Jun 2017

arxiv: v1 [hep-ex] 8 Jun 2017 UCHEP 17 05 6 Aril 017 Prosects for time-deendent mixing and CP-violation measurements at Belle II arxiv:1706.0363v1 [he-ex] 8 Jun 017 Physics Deartment, University of Cincinnati, Cincinnati, Ohio 51 E-mail:

More information

Research of PMU Optimal Placement in Power Systems

Research of PMU Optimal Placement in Power Systems Proceedings of the 5th WSEAS/IASME Int. Conf. on SYSTEMS THEORY and SCIENTIFIC COMPUTATION, Malta, Setember 15-17, 2005 (38-43) Research of PMU Otimal Placement in Power Systems TIAN-TIAN CAI, QIAN AI

More information

Adaptive estimation with change detection for streaming data

Adaptive 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 information

Review of Probability Theory II

Review of Probability Theory II Review of Probability Theory II January 9-3, 008 Exectation If the samle sace Ω = {ω, ω,...} is countable and g is a real-valued function, then we define the exected value or the exectation of a function

More information

SIMULATED ANNEALING AND JOINT MANUFACTURING BATCH-SIZING. Ruhul SARKER. Xin YAO

SIMULATED ANNEALING AND JOINT MANUFACTURING BATCH-SIZING. Ruhul SARKER. Xin YAO Yugoslav Journal of Oerations Research 13 (003), Number, 45-59 SIMULATED ANNEALING AND JOINT MANUFACTURING BATCH-SIZING Ruhul SARKER School of Comuter Science, The University of New South Wales, ADFA,

More information

Covariance Matrix Estimation for Reinforcement Learning

Covariance Matrix Estimation for Reinforcement Learning Covariance Matrix Estimation for Reinforcement Learning Tomer Lancewicki Deartment of Electrical Engineering and Comuter Science University of Tennessee Knoxville, TN 37996 tlancewi@utk.edu Itamar Arel

More information

An Investigation on the Numerical Ill-conditioning of Hybrid State Estimators

An Investigation on the Numerical Ill-conditioning of Hybrid State Estimators An Investigation on the Numerical Ill-conditioning of Hybrid State Estimators S. K. Mallik, Student Member, IEEE, S. Chakrabarti, Senior Member, IEEE, S. N. Singh, Senior Member, IEEE Deartment of Electrical

More information

José Alberto Mauricio. Instituto Complutense de Análisis Económico. Universidad Complutense de Madrid. Campus de Somosaguas Madrid - SPAIN

José Alberto Mauricio. Instituto Complutense de Análisis Económico. Universidad Complutense de Madrid. Campus de Somosaguas Madrid - SPAIN A CORRECTED ALGORITHM FOR COMPUTING THE THEORETICAL AUTOCOVARIANCE MATRICES OF A VECTOR ARMA MODEL José Alberto Mauricio Instituto Comlutense de Análisis Económico Universidad Comlutense de Madrid Camus

More information

Lecture: Condorcet s Theorem

Lecture: Condorcet s Theorem Social Networs and Social Choice Lecture Date: August 3, 00 Lecture: Condorcet s Theorem Lecturer: Elchanan Mossel Scribes: J. Neeman, N. Truong, and S. Troxler Condorcet s theorem, the most basic jury

More information

STK4900/ Lecture 7. Program

STK4900/ Lecture 7. Program STK4900/9900 - Lecture 7 Program 1. Logistic regression with one redictor 2. Maximum likelihood estimation 3. Logistic regression with several redictors 4. Deviance and likelihood ratio tests 5. A comment

More information

Modeling and Estimation of Full-Chip Leakage Current Considering Within-Die Correlation

Modeling 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 information

A Qualitative Event-based Approach to Multiple Fault Diagnosis in Continuous Systems using Structural Model Decomposition

A 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 information

Elementary Analysis in Q p

Elementary Analysis in Q p Elementary Analysis in Q Hannah Hutter, May Szedlák, Phili Wirth November 17, 2011 This reort follows very closely the book of Svetlana Katok 1. 1 Sequences and Series In this section we will see some

More information

Location of solutions for quasi-linear elliptic equations with general gradient dependence

Location of solutions for quasi-linear elliptic equations with general gradient dependence Electronic Journal of Qualitative Theory of Differential Equations 217, No. 87, 1 1; htts://doi.org/1.14232/ejqtde.217.1.87 www.math.u-szeged.hu/ejqtde/ Location of solutions for quasi-linear ellitic equations

More information

One-way ANOVA Inference for one-way ANOVA

One-way ANOVA Inference for one-way ANOVA One-way ANOVA Inference for one-way ANOVA IPS Chater 12.1 2009 W.H. Freeman and Comany Objectives (IPS Chater 12.1) Inference for one-way ANOVA Comaring means The two-samle t statistic An overview of ANOVA

More information

Discrete alpha-skew-laplace distribution

Discrete alpha-skew-laplace distribution Statistics & Oerations Research Transactions SORT 39 () January-June 25, 7-84 ISSN: 696-228 eissn: 23-883 www.idescat.cat/sort/ Statistics & Oerations Research c Institut d Estadstica de Catalunya Transactions

More information

Asymptotically Optimal Simulation Allocation under Dependent Sampling

Asymptotically 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 information

The Binomial Approach for Probability of Detection

The Binomial Approach for Probability of Detection Vol. No. (Mar 5) - The e-journal of Nondestructive Testing - ISSN 45-494 www.ndt.net/?id=7498 The Binomial Aroach for of Detection Carlos Correia Gruo Endalloy C.A. - Caracas - Venezuela www.endalloy.net

More information

Notes on Instrumental Variables Methods

Notes on Instrumental Variables Methods Notes on Instrumental Variables Methods Michele Pellizzari IGIER-Bocconi, IZA and frdb 1 The Instrumental Variable Estimator Instrumental variable estimation is the classical solution to the roblem of

More information

Chapter 13 Variable Selection and Model Building

Chapter 13 Variable Selection and Model Building Chater 3 Variable Selection and Model Building The comlete regsion analysis deends on the exlanatory variables ent in the model. It is understood in the regsion analysis that only correct and imortant

More information

Use of Transformations and the Repeated Statement in PROC GLM in SAS Ed Stanek

Use 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 information