/*! This file is auto-generated */ .wp-block-button__link{color:#fff;background-color:#32373c;border-radius:9999px;box-shadow:none;text-decoration:none;padding:calc(.667em + 2px) calc(1.333em + 2px);font-size:1.125em}.wp-block-file__button{background:#32373c;color:#fff;text-decoration:none} Q21E Let聽\({\rm{X}}\)have a Weibull ... [FREE SOLUTION] | 91影视

91影视

Let\({\rm{X}}\)have a Weibull distribution with parameters\({\rm{\alpha }}\)and\({\rm{\beta }}\), so

\(\begin{array}{l}{\rm{E(X) = \beta \times \Gamma (1 + 1/\alpha )V(X)}}\\{\rm{ = }}{{\rm{\beta }}^{\rm{2}}}\left\{ {{\rm{\Gamma (1 + 2/\alpha ) - (\Gamma (1 + 1/\alpha )}}{{\rm{)}}^{\rm{2}}}} \right\}\end{array}\)

a. Based on a random sample\({{\rm{X}}_{\rm{1}}}{\rm{, \ldots ,}}{{\rm{X}}_{\rm{n}}}\), write equations for the method of moments estimators of\({\rm{\beta }}\)and\({\rm{\alpha }}\). Show that, once the estimate of\({\rm{\alpha }}\)has been obtained, the estimate of\({\rm{\beta }}\)can be found from a table of the gamma function and that the estimate of\({\rm{\alpha }}\)is the solution to a complicated equation involving the gamma function.

b. If\({\rm{n = 20,\bar x = 28}}{\rm{.0}}\), and\({\rm{\Sigma x}}_{\rm{i}}^{\rm{2}}{\rm{ = 16,500}}\), compute the estimates. (Hint:\(\left. {{{{\rm{(\Gamma (1}}{\rm{.2))}}}^{\rm{2}}}{\rm{/\Gamma (1}}{\rm{.4) = }}{\rm{.95}}{\rm{.}}} \right)\)

Short Answer

Expert verified

a) The moment estimator \({\rm{\hat \beta }}\)can be obtained once the equation is solved. The goal was to figure out how to get the method estimators alpha and\({\rm{\hat \beta }}\).

b) The estimates are\(\hat \alpha = 5;\hat \beta = \frac{{28}}{{\Gamma (1.2)}}\)

Step by step solution

01

Introduction

An estimator is a rule for computing an estimate of a given quantity based on observable data: the rule (estimator), the quantity of interest (estimate), and the output (estimate) are all distinct.

02

Explanation

a)

Let random variables \({{\rm{X}}_{\rm{1}}}{\rm{,}}{{\rm{X}}_{\rm{2}}}{\rm{, \ldots ,}}{{\rm{X}}_{\rm{n}}}\) have the same distribution as pmf or pdf \(f\left( {x;{\theta _1},{\theta _2}, \ldots ,{\theta _m}} \right),m \^I N\) with unknown parameters \({{\rm{\theta }}_{\rm{i}}}{\rm{,i = 1,2, \ldots ,m}}\).

By equating sample moments to corresponding population moments and solving for unknown parameters\(\widehat {{{\rm{\theta }}_{\rm{1}}}}{\rm{,}}\widehat {{{\rm{\theta }}_{\rm{2}}}}{\rm{, \ldots ,}}\widehat {{{\rm{\theta }}_{\rm{m}}}}\), the moment estimators\({{\rm{\theta }}_{\rm{i}}}{\rm{,i = 1,2, \ldots ,m}}\)may be obtained.

The specified distribution in this exercise is Weibull's distribution with parameters\({\rm{\alpha }}\)and\({\rm{\beta }}\),forwhichmomentestimatorsmustbedeveloped.

The sample moment of first order is

\({\rm{\bar X = }}\frac{{\rm{1}}}{{\rm{n}}}\left( {{{\rm{X}}_{\rm{1}}}{\rm{ + }}{{\rm{X}}_{\rm{2}}}{\rm{ + \ldots + }}{{\rm{X}}_{\rm{n}}}} \right)\)

and the population moment of first order is

\({\rm{E(X) = \beta \times \Gamma }}\left( {{\rm{1 + }}\frac{{\rm{1}}}{{\rm{\alpha }}}} \right)\)

The first equation in the system of equations from which the moment estimators are obtained is

\({\rm{\bar X = E(X)}}\)

The sample moment of second order is:

\(\frac{{\rm{1}}}{{\rm{n}}}\left( {{\rm{X}}_{\rm{1}}^{\rm{2}}{\rm{ + X}}_{\rm{2}}^{\rm{2}}{\rm{ + \ldots + X}}_{\rm{n}}^{\rm{2}}} \right)\)

And the population moment of second order is:

\({\rm{E}}\left( {{{\rm{X}}^{\rm{2}}}} \right){\rm{ = V(X) + (E(X)}}{{\rm{)}}^{\rm{2}}}{\rm{ - }}{{\rm{\beta }}^{\rm{2}}}\left\{ {{\rm{\Gamma }}\left( {{\rm{1 + }}\frac{{\rm{2}}}{{\rm{\alpha }}}} \right){\rm{ - }}{{\left( {{\rm{\Gamma }}\left( {{\rm{1 + }}\frac{{\rm{1}}}{{\rm{\alpha }}}} \right)} \right)}^{\rm{2}}}} \right\}{\rm{ - }}{\left\{ {{\rm{\beta \times \Gamma }}\left( {{\rm{1 + }}\frac{{\rm{1}}}{{\rm{\alpha }}}} \right)} \right\}^{\rm{2}}}{\rm{ - }}{{\rm{\beta }}^{\rm{2}}}{\rm{\Gamma }}\left( {{\rm{1 + }}\frac{{\rm{2}}}{{\rm{\alpha }}}} \right)\)

The moment estimators are derived from the second equation in the system of equations.

\(\frac{{\rm{1}}}{{\rm{n}}}\sum\limits_{{\rm{i = 1}}}^{\rm{n}} {{\rm{X}}_{\rm{i}}^{\rm{2}}} {\rm{ - E}}\left( {{{\rm{X}}^{\rm{2}}}} \right)\)

As a result, the system of equations that must be solved for\({\rm{\hat \alpha }}\) and \({\rm{\hat \beta }}\) is

\(\begin{array}{*{20}{r}}{{\rm{\bar X - \hat \beta \times \Gamma }}\left( {{\rm{1 + }}\frac{{\rm{1}}}{{{\rm{\hat \alpha }}}}} \right){\rm{,}}}\\{\frac{{\rm{1}}}{{\rm{n}}}\sum\limits_{{\rm{i = 1}}}^{\rm{n}} {{\rm{X}}_{\rm{i}}^{\rm{2}}} {\rm{ - }}{{{\rm{\hat \beta }}}^{\rm{2}}}{\rm{\Gamma }}\left( {{\rm{1 + }}\frac{{\rm{2}}}{{{\rm{\hat \alpha }}}}} \right){\rm{.}}}\end{array}\)

Hence, \(\beta \)can be computed from the first equation as

\({\rm{\hat \beta = }}\frac{{{\rm{\bar X}}}}{{{\rm{\Gamma }}\left( {{\rm{1 + }}\frac{{\rm{1}}}{{{\rm{\dot \alpha }}}}} \right)}}{\rm{.}}\)

In order to compute\({\rm{\hat \beta }}\), first \({\rm{\hat \alpha }}\) needs to be determined and gamma functioned evaluated.

From first equation, by squaring both sides, the following stands

\({{\rm{\bar X}}^{\rm{2}}}{\rm{ = }}{{\rm{\hat \beta }}^{\rm{2}}}{{\rm{\Gamma }}^{\rm{2}}}\left( {{\rm{1 + }}\frac{{\rm{1}}}{{{\rm{\hat \alpha }}}}} \right){\rm{.}}\)

or equally

\({{\rm{\hat \beta }}^{\rm{2}}}{\rm{ = }}\frac{{{{{\rm{\bar X}}}^{\rm{2}}}}}{{{{\rm{\Gamma }}^{\rm{2}}}\left( {{\rm{1 + }}\frac{{\rm{1}}}{{{\rm{\hat \alpha }}}}} \right)}}{\rm{.}}\)

Now plug in in the second equation:

\(\frac{{\rm{1}}}{{\rm{n}}}\sum\limits_{{\rm{i = 1}}}^{\rm{n}} {{\rm{X}}_{\rm{i}}^{\rm{2}}} {\rm{ = }}\frac{{{{{\rm{\bar X}}}^{\rm{2}}}}}{{{{\rm{\Gamma }}^{\rm{2}}}\left( {{\rm{1 + }}\frac{{\rm{1}}}{{\rm{\alpha }}}} \right)}}{\rm{ \times \Gamma }}\left( {{\rm{1 + }}\frac{{\rm{2}}}{{{\rm{\hat \alpha }}}}} \right)\)

or equally

\(\frac{{\rm{1}}}{{\rm{n}}}\frac{{\sum\limits_{{\rm{i = 1}}}^{\rm{n}} {{\rm{X}}_{\rm{i}}^{\rm{2}}} }}{{{{{\rm{\bar X}}}^{\rm{2}}}}}{\rm{ = }}\frac{{{\rm{\Gamma }}\left( {{\rm{1 + }}\frac{{\rm{2}}}{{\rm{\alpha }}}} \right)}}{{{{\rm{\Gamma }}^{\rm{2}}}\left( {{\rm{1 + }}\frac{{\rm{1}}}{{\rm{\alpha }}}} \right)}}{\rm{.}}\)

The only unknown variable in the last equation is \(\hat \alpha \), and the moment estimator alpha may be obtained by solving this equation. The moment estimator \({\rm{\hat \beta }}\)can be obtained once the equation is solved.

Because, as stated in the exercise, solving this equation is difficult, there is no need to solve it. The goal was to figure out how to get the method estimators alpha and \({\rm{\hat \beta }}\).

03

Explanation

b)

Consider the given information,

\(\frac{{\rm{1}}}{{\rm{n}}}\frac{{\sum\limits_{{\rm{i = 1}}}^{\rm{n}} {{\rm{X}}_{\rm{i}}^{\rm{2}}} }}{{{{{\rm{\bar X}}}^{\rm{2}}}}}{\rm{ = }}\frac{{{\rm{\Gamma }}\left( {{\rm{1 + }}\frac{{\rm{2}}}{{{\rm{\hat \alpha }}}}} \right)}}{{{{\rm{\Gamma }}^{\rm{2}}}\left( {{\rm{1 + }}\frac{{\rm{1}}}{{\rm{\alpha }}}} \right)}}\)

For given\({\rm{n = 20,\bar x = 28}}\), and \(\sum {{\rm{x}}_{\rm{i}}^{\rm{2}}} {\rm{ = 16,500,\hat \alpha }}\) needs to be found.

The following is true:

\(\frac{{\rm{1}}}{{{\rm{20}}}}{\rm{ \times }}\left( {\frac{{{\rm{16,500}}}}{{{\rm{2}}{{\rm{8}}^{\rm{2}}}}}} \right){\rm{ = 1}}{\rm{.05}}\)

therefore,

\(\frac{{{\rm{\Gamma }}\left( {{\rm{1 + }}\frac{{\rm{2}}}{{\frac{{\rm{\alpha }}}{{\rm{\alpha }}}}}} \right)}}{{{{\rm{\Gamma }}^{\rm{2}}}\left( {{\rm{1 + }}\frac{{\rm{1}}}{{\rm{\alpha }}}} \right)}}{\rm{ = 1}}{\rm{.05}}{\rm{.}}\)

From the hint

\(\frac{{{{{\rm{(\Gamma (1}}{\rm{.2))}}}^{\rm{2}}}}}{{{\rm{\Gamma (1}}{\rm{.4)}}}}{\rm{ = 0}}{\rm{.95}}\)

or equally

\(\begin{array}{l}\frac{{{\rm{\Gamma (1 + 0}}{\rm{.4)}}}}{{{{\rm{\Gamma }}^{\rm{2}}}{\rm{(1 + 0}}{\rm{.2)}}}}{\rm{ = }}\frac{{\rm{1}}}{{{\rm{0}}{\rm{.95}}}}\\\frac{{{\rm{\Gamma (1 + 0}}{\rm{.4)}}}}{{{{\rm{\Gamma }}^{\rm{2}}}{\rm{(1 + 0}}{\rm{.2)}}}}{\rm{ = 1}}{\rm{.05}}\end{array}\)

which means that

\(\begin{array}{l}\frac{{{\rm{\Gamma }}\left( {{\rm{1 + }}\frac{{\rm{2}}}{{{\rm{\bar \alpha }}}}} \right)}}{{{{\rm{\Gamma }}^{\rm{2}}}\left( {{\rm{1 + }}\frac{{\rm{1}}}{{{\rm{\dot \alpha }}}}} \right)}}{\rm{ = 1}}{\rm{.05}}\\{\rm{ = }}\frac{{{\rm{\Gamma (1 + 0}}{\rm{.4)}}}}{{{{\rm{\Gamma }}^{\rm{2}}}{\rm{(1 + 0}}{\rm{.2)}}}}\end{array}\)

Therefore, because of this equality, the following must hold

\(\frac{{\rm{2}}}{{{\rm{\hat \alpha }}}}{\rm{ = 0}}{\rm{.4}}\)

hence,

\({\rm{\hat \alpha = 5}}{\rm{.}}\)

From the estimator

\({\rm{\hat \beta = }}\frac{{{\rm{\bar X}}}}{{{\rm{\Gamma }}\left( {{\rm{1 + }}\frac{{\rm{1}}}{{{\rm{\bar \alpha }}}}} \right)}}\)

The estimate is computed as follows:\({\rm{\hat \beta = }}\frac{{{\rm{28}}}}{{{\rm{\Gamma (1}}{\rm{.2)}}}}\)

04

Additional

The average value of a function \(f\left( x,y \right)\) over a rectangle \(R\) is defined to be \({{f}_{ave}}=\frac{1}{A\left( R \right)}\iint\limits_{R}{f\left( x,y \right)dA}\).

Find the average value of \(f\) over the given rectangle, \(f\left( x,y \right)={{e}^{y}}\sqrt{x+{{e}^{y}}}\), \(R=\left( 0,4 \right)\times \left( 0,1 \right)\).

Given: \(f\left( x,y \right)={{e}^{y}}\sqrt{x+{{e}^{y}}}\)

\(R=\left( 0,4 \right)\times \left( 0,1 \right)\)

To find: average value of f.

Unlock Step-by-Step Solutions & Ace Your Exams!

  • Full Textbook Solutions

    Get detailed explanations and key concepts

  • Unlimited Al creation

    Al flashcards, explanations, exams and more...

  • Ads-free access

    To over 500 millions flashcards

  • Money-back guarantee

    We refund you if you fail your exam.

Over 30 million students worldwide already upgrade their learning with 91影视!

One App. One Place for Learning.

All the tools & learning materials you need for study success - in one app.

Get started for free

Most popular questions from this chapter

A vehicle with a particular defect in its emission control system is taken to a succession of randomly selected mechanics until\({\rm{r = 3}}\)of them have correctly diagnosed the problem. Suppose that this requires diagnoses by\({\rm{20}}\)different mechanics (so there were\({\rm{17}}\)incorrect diagnoses). Let\({\rm{p = P}}\)(correct diagnosis), so\({\rm{p}}\)is the proportion of all mechanics who would correctly diagnose the problem. What is the mle of\({\rm{p}}\)? Is it the same as the mle if a random sample of\({\rm{20}}\)mechanics results in\({\rm{3}}\)correct diagnoses? Explain. How does the mle compare to the estimate resulting from the use of the unbiased estimator?

Let\({{\rm{X}}_{\rm{1}}}{\rm{,}}{{\rm{X}}_{\rm{2}}}{\rm{, \ldots ,}}{{\rm{X}}_{\rm{n}}}\)represent a random sample from a Rayleigh distribution with pdf

\({\rm{f(x,\theta ) = }}\frac{{\rm{x}}}{{\rm{\theta }}}{{\rm{e}}^{{\rm{ - }}{{\rm{x}}^{\rm{2}}}{\rm{/(2\theta )}}}}\quad {\rm{x > 0}}\)a. It can be shown that\({\rm{E}}\left( {{{\rm{X}}^{\rm{2}}}} \right){\rm{ = 2\theta }}\). Use this fact to construct an unbiased estimator of\({\rm{\theta }}\)based on\({\rm{\Sigma X}}_{\rm{i}}^{\rm{2}}\)(and use rules of expected value to show that it is unbiased).

b. Estimate\({\rm{\theta }}\)from the following\({\rm{n = 10}}\)observations on vibratory stress of a turbine blade under specified conditions:

\(\begin{array}{*{20}{l}}{{\rm{16}}{\rm{.88}}}&{{\rm{10}}{\rm{.23}}}&{{\rm{4}}{\rm{.59}}}&{{\rm{6}}{\rm{.66}}}&{{\rm{13}}{\rm{.68}}}\\{{\rm{14}}{\rm{.23}}}&{{\rm{19}}{\rm{.87}}}&{{\rm{9}}{\rm{.40}}}&{{\rm{6}}{\rm{.51}}}&{{\rm{10}}{\rm{.95}}}\end{array}\)

Urinary angiotensinogen (AGT) level is one quantitative indicator of kidney function. The article 鈥淯rinary Angiotensinogen as a Potential Biomarker of Chronic Kidney Diseases鈥 (J. of the Amer. Society of Hypertension, \({\rm{2008: 349 - 354}}\)) describes a study in which urinary AGT level \({\rm{(\mu g)}}\) was determined for a sample of adults with chronic kidney disease. Here is representative data (consistent with summary quantities and descriptions in the cited article):

An appropriate probability plot supports the use of the lognormal distribution (see Section \({\rm{4}}{\rm{.5}}\)) as a reasonable model for urinary AGT level (this is what the investigators did).

a. Estimate the parameters of the distribution. (Hint: Rem ember that \({\rm{X}}\) has a lognormal distribution with parameters \({\rm{\mu }}\) and \({{\rm{\sigma }}^{\rm{2}}}\) if \({\rm{ln(X)}}\) is normally distributed with mean \({\rm{\mu }}\) and variance \({{\rm{\sigma }}^{\rm{2}}}\).)

b. Use the estimates of part (a) to calculate an estimate of the expected value of AGT level. (Hint: What is \({\rm{E(X)}}\)?)

Each of 150 newly manufactured items is examined and the number of scratches per item is recorded (the items are supposed to be free of scratches), yielding the following data:

Assume that X has a Poisson distribution with parameter \({\bf{\mu }}.\)and that X represents the number of scratches on a randomly picked item.

a. Calculate the estimate for the data using an unbiased \({\bf{\mu }}.\)estimator. (Hint: for X Poisson, \({\rm{E(X) = \mu }}\) ,therefore \({\rm{E(\bar X) = ?)}}\)

c. What is your estimator's standard deviation (standard error)? Calculate the standard error estimate. (Hint: \({\rm{\sigma }}_{\rm{X}}^{\rm{2}}{\rm{ = \mu }}\), \({\rm{X}}\))

Consider randomly selecting \({\rm{n}}\) segments of pipe and determining the corrosion loss (mm) in the wall thickness for each one. Denote these corrosion losses by \({{\rm{Y}}_{\rm{1}}}{\rm{,}}.....{\rm{,}}{{\rm{Y}}_{\rm{n}}}\). The article 鈥淎 Probabilistic Model for a Gas Explosion Due to Leakages in the Grey Cast Iron Gas Mains鈥 (Reliability Engr. and System Safety (\({\rm{(2013:270 - 279)}}\)) proposes a linear corrosion model: \({{\rm{Y}}_{\rm{i}}}{\rm{ = }}{{\rm{t}}_{\rm{i}}}{\rm{R}}\), where \({{\rm{t}}_{\rm{i}}}\) is the age of the pipe and \({\rm{R}}\), the corrosion rate, is exponentially distributed with parameter \({\rm{\lambda }}\). Obtain the maximum likelihood estimator of the exponential parameter (the resulting mle appears in the cited article). (Hint: If \({\rm{c > 0}}\) and \({\rm{X}}\) has an exponential distribution, so does \({\rm{cX}}\).)

See all solutions

Recommended explanations on Math Textbooks

View all explanations

What do you think about this solution?

We value your feedback to improve our textbook solutions.

Study anywhere. Anytime. Across all devices.