/*! 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} Q28E Let \({{\rm{X}}_{\rm{1}}}{\rm{,}... [FREE SOLUTION] | 91影视

91影视

Let \({{\rm{X}}_{\rm{1}}}{\rm{,}}{{\rm{X}}_{\rm{2}}}.....{\rm{,}}{{\rm{X}}_{\rm{n}}}\) represent a random sample from the Rayleigh distribution with density function given in Exercise \({\rm{15}}\). Determine a. The maximum likelihood estimator of \({\rm{\theta }}\), and then calculate the estimate for the vibratory stress data given in that exercise. Is this estimator the same as the unbiased estimator suggested in Exercise \({\rm{15}}\)? b. The mle of the median of the vibratory stress distribution. (Hint: First express the median in terms of \({\rm{\theta }}\).)

Short Answer

Expert verified

(a) Maximum likelihood estimator is \({\rm{\hat \theta = }}\frac{{\rm{1}}}{{{\rm{2n}}}}\sum\limits_{{\rm{i = 1}}}^{\rm{n}} {{\rm{X}}_{\rm{i}}^{\rm{2}}} \).

(b) The mle of median is \({\rm{\hat \sim x = }}\sqrt {{\rm{1}}{\rm{.3863 \times \hat \theta }}} \).

Step by step solution

01

Define equations

A mathematical language that asserts that two algebraic expressions must be equal in nature is known as an equation.

02

Explanation

(a) The pdf of Rayleigh distribution is,

\({\rm{f(x;\theta ) = }}\frac{{\rm{x}}}{{\rm{\theta }}}{{\rm{e}}^{{\rm{ - }}{{\rm{x}}^{\rm{2}}}{\rm{/(2\theta )}}}}{\rm{, x > 0}}\)

Allow joint pdf or pmf for random variables\({{\rm{X}}_{\rm{1}}}{\rm{,}}{{\rm{X}}_{\rm{2}}}{\rm{, \ldots ,}}{{\rm{X}}_{\rm{n}}}\).

\({\rm{f}}\left( {{{\rm{x}}_{\rm{1}}}{\rm{,}}{{\rm{x}}_{\rm{2}}}{\rm{, \ldots ,}}{{\rm{x}}_{\rm{n}}}{\rm{;}}{{\rm{\theta }}_{\rm{1}}}{\rm{,}}{{\rm{\theta }}_{\rm{2}}}{\rm{, \ldots ,}}{{\rm{\theta }}_{\rm{m}}}} \right){\rm{, n,m}} \in {\rm{N}}\)

where\({{\rm{\theta }}_{\rm{i}}}{\rm{,i = 1,2, \ldots ,m}}\)are unknown parameters. The likelihood function is defined as a function of parameters\({{\rm{\theta }}_{\rm{i}}}{\rm{,i = 1,2, \ldots ,m}}\)where function f is a function of parameter. The maximum likelihood estimates (mle's), or values\(\widehat {{{\rm{\theta }}_{\rm{i}}}}\)for which the likelihood function is maximised, are the maximum likelihood estimates,

\({\rm{f}}\left( {{{\rm{x}}_{\rm{1}}}{\rm{,}}{{\rm{x}}_{\rm{2}}}{\rm{, \ldots ,}}{{\rm{x}}_{\rm{n}}}{\rm{;}}{{{\rm{\hat \theta }}}_{\rm{1}}}{\rm{,}}{{{\rm{\hat \theta }}}_{\rm{2}}}{\rm{, \ldots ,}}{{{\rm{\hat \theta }}}_{\rm{m}}}} \right) \ge {\rm{f}}\left( {{{\rm{x}}_{\rm{1}}}{\rm{,}}{{\rm{x}}_{\rm{2}}}{\rm{, \ldots ,}}{{\rm{x}}_{\rm{n}}}{\rm{;}}{{\rm{\theta }}_{\rm{1}}}{\rm{,}}{{\rm{\theta }}_{\rm{2}}}{\rm{, \ldots ,}}{{\rm{\theta }}_{\rm{m}}}} \right)\)

As,\({\rm{i = 1,2, \ldots ,m}}\)for every\({{\rm{\theta }}_{\rm{i}}}\). Maximum likelihood estimators are derived by replacing\({{\rm{X}}_{\rm{i}}}\)with\({{\rm{x}}_{\rm{i}}}\).

Because of the independence, the likelihood function becomes,

\(\begin{array}{c}{\rm{f}}\left( {{{\rm{x}}_{\rm{1}}}{\rm{,}}{{\rm{x}}_{\rm{x}}}{\rm{, \ldots ,}}{{\rm{x}}_{\rm{n}}}{\rm{;\theta }}} \right){\rm{ = }}\frac{{{{\rm{x}}_{\rm{1}}}}}{{\rm{\theta }}}{{\rm{e}}^{{\rm{ - x}}_{\rm{1}}^{\rm{2}}{\rm{/(2\theta )}}}}{\rm{ \times }}\frac{{{{\rm{x}}_{\rm{2}}}}}{{\rm{\theta }}}{{\rm{e}}^{{\rm{ - x}}_{\rm{2}}^{\rm{2}}{\rm{/(2\theta )}}}}{\rm{ \times \ldots \times }}\frac{{{{\rm{x}}_{\rm{n}}}}}{{\rm{\theta }}}{{\rm{e}}^{{\rm{ - x}}_{\rm{n}}^{\rm{2}}{\rm{/(2\theta )}}}}\\{\rm{ = }}\left( {{{\rm{x}}_{\rm{1}}}{\rm{ \times }}{{\rm{x}}_{\rm{2}}}{\rm{ \times \ldots \times }}{{\rm{x}}_{\rm{n}}}} \right){\rm{ \times }}\frac{{\rm{1}}}{{{{\rm{\theta }}^{\rm{n}}}}}{\rm{ \times exp}}\left\{ {{\rm{ - }}\frac{{\sum\limits_{{\rm{i = 1}}}^{\rm{n}} {{\rm{x}}_{\rm{i}}^{\rm{2}}} }}{{{\rm{2\theta }}}}} \right\}\end{array}\)

03

Evaluating the maximum likelihood estimators

Look at the log likelihood function to determine the maximum.

\(\begin{array}{c}{\rm{lnf}}\left( {{{\rm{x}}_{\rm{1}}}{\rm{,}}{{\rm{x}}_{\rm{x}}}{\rm{, \ldots ,}}{{\rm{x}}_{\rm{n}}}{\rm{;\theta }}} \right){\rm{ = ln}}\left( {\left( {{{\rm{x}}_{\rm{1}}}{\rm{ \times }}{{\rm{x}}_{\rm{2}}}{\rm{ \times \ldots \times }}{{\rm{x}}_{\rm{n}}}} \right){\rm{ \times }}\frac{{\rm{1}}}{{{{\rm{\theta }}^{\rm{n}}}}}{\rm{ \times exp}}\left\{ {{\rm{ - }}\frac{{\sum\limits_{{\rm{i = 1}}}^{\rm{n}} {{\rm{x}}_{\rm{i}}^{\rm{2}}} }}{{{\rm{2\theta }}}}} \right\}} \right)\\{\rm{ = ln}}\left( {{{\rm{x}}_{\rm{1}}}{\rm{ \times }}{{\rm{x}}_{\rm{2}}}{\rm{ \times \ldots \times }}{{\rm{x}}_{\rm{n}}}} \right){\rm{ - nln\theta - }}\frac{{\sum\limits_{{\rm{i = 1}}}^{\rm{n}} {{\rm{x}}_{\rm{i}}^{\rm{2}}} }}{{{\rm{2\theta }}}}\end{array}\)

The maximum likelihood estimator is generated by taking the derivative of the log likelihood function in regard to\({\rm{\theta }}\)and equating it to\({\rm{0}}\).

As a result, the derivative,

\(\begin{array}{c}\frac{{\rm{d}}}{{{\rm{d\theta }}}}{\rm{f}}\left( {{{\rm{x}}_{\rm{1}}}{\rm{,}}{{\rm{x}}_{\rm{x}}}{\rm{, \ldots ,}}{{\rm{x}}_{\rm{n}}}{\rm{;\theta }}} \right){\rm{ = }}\frac{{\rm{d}}}{{{\rm{d\theta }}}}\left( {{\rm{ln}}\left( {{{\rm{x}}_{\rm{1}}}{\rm{ \times }}{{\rm{x}}_{\rm{2}}}{\rm{ \times \ldots \times }}{{\rm{x}}_{\rm{n}}}} \right){\rm{ - nln\theta - }}\frac{{\sum\limits_{{\rm{i = 1}}}^{\rm{n}} {{\rm{x}}_{\rm{i}}^{\rm{2}}} }}{{{\rm{2\theta }}}}} \right)\\{\rm{ = 0 - n \times }}\frac{{\rm{1}}}{{\rm{\theta }}}{\rm{ + }}\frac{{\rm{1}}}{{{\rm{2}}{{\rm{\theta }}^{\rm{2}}}}}\sum\limits_{{\rm{i = 1}}}^{\rm{n}} {{\rm{x}}_{\rm{i}}^{\rm{2}}} \\{\rm{ = }}\frac{{\rm{1}}}{{{\rm{2}}{{\rm{\theta }}^{\rm{2}}}}}\sum\limits_{{\rm{i = 1}}}^{\rm{n}} {{\rm{x}}_{\rm{i}}^{\rm{2}}} {\rm{ - n \times }}\frac{{\rm{1}}}{{\rm{\theta }}}\end{array}\)

As a result, solving equation provides the maximum likelihood estimator.

\(\begin{array}{c}\frac{{\rm{1}}}{{{\rm{2}}{{{\rm{\hat \theta }}}^{\rm{2}}}}}\sum\limits_{{\rm{i = 1}}}^{\rm{n}} {{\rm{x}}_{\rm{i}}^{\rm{2}}} {\rm{ - n \times }}\frac{{\rm{1}}}{{{\rm{\hat \theta }}}}{\rm{ = 0}}\\{\rm{n\hat \theta = }}\frac{{\rm{1}}}{{\rm{2}}}\sum\limits_{{\rm{i = 1}}}^{\rm{n}} {{\rm{x}}_{\rm{i}}^{\rm{2}}} \end{array}\)

For\({\rm{\theta }}\). Hence, the maximum likelihood estimator is,

\({\rm{\hat \theta = }}\frac{{\rm{1}}}{{{\rm{2n}}}}\sum\limits_{{\rm{i = 1}}}^{\rm{n}} {{\rm{X}}_{\rm{i}}^{\rm{2}}} \)

04

Explanation

(b) The median can be calculated by equating the cdf of the Rayleigh distribution with \({\rm{0}}{\rm{.5}}\), or equally in terms of \({\rm{\theta }}\),

\({\rm{F(x;\theta ) = 0}}{\rm{.5}}\)

For\({\rm{x > 0}}\), the cdf is,

\(\begin{array}{c}{\rm{ F(x;\theta ) = P(X}} \le {\rm{x)}}\\{\rm{ = }}\int_{\rm{0}}^{\rm{x}} {\rm{f}} {\rm{(t;\theta )dt}}\\{\rm{ = 1 - }}{{\rm{e}}^{{\rm{ - }}{{\rm{x}}^{\rm{2}}}{\rm{/(2\theta )}}}}\end{array}\)

where this cdf can be found (you can solve integral and obtain it as well). Hence,

\(\begin{array}{c}{\rm{F(x;\theta ) = 0}}{\rm{.5}}\\{\rm{1 - }}{{\rm{e}}^{{\rm{ - }}{{\rm{x}}^{\rm{2}}}{\rm{/(2\theta )}}}}{\rm{ = 0}}{\rm{.5}}\\{{\rm{e}}^{{\rm{ - }}{{\rm{x}}^{\rm{2}}}{\rm{/(2\theta )}}}}{\rm{ = 0}}{\rm{.5}}\\\frac{{{\rm{ - }}{{\rm{x}}^{\rm{2}}}}}{{{\rm{2\theta }}}}{\rm{ = ln0}}{\rm{.5}}\\{{\rm{x}}^{\rm{2}}}{\rm{ = - 2\theta \times - 0}}{\rm{.6931}}\\{\rm{x = }}\sqrt {{\rm{1}}{\rm{.3863 \times \theta }}} \end{array}\)

Let\({\rm{i = 1,2,}}......{\rm{,n}}\)be the maximum likelihood estimates of the parameters\(\widehat {{{\rm{\theta }}_{\rm{i}}}}{\rm{,i = 1,2,}}......{\rm{,n}}\). The function of the mle's\(\widehat {{{\rm{\theta }}_{\rm{i}}}}\); is the mle of any function of parameters\(\widehat {{{\rm{\theta }}_{\rm{i}}}}\).

Because of the invariance principle and the fact that the median is a function of\({\rm{\theta }}\), the maximum likelihood estimator of the median is a function\(\widehat {\rm{\theta }}\)of or equally is,

\({\rm{\hat \sim x = }}\sqrt {{\rm{1}}{\rm{.3863 \times \hat \theta }}} \).

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

The first four deviations from the mean in a sample of n=5 reaction times were .3, .9, 1.0, and 1.3. What is the fifth deviation from the mean? Give a sample for which these are the five deviations from the mean.

The accompanying summary data on CeO2 particlesizes (nm) under certain experimental conditions wasread from a graph in the article 鈥淣anoceria鈥擡nergetics of Surfaces, Interfaces and WaterAdsorption鈥 (J. of the Amer. Ceramic Soc., 2011:3992鈥3999):

3.0鈭<3.5 3.5鈭<4.0 4.0鈭<4.5 4.5鈭<5.0 5.0鈭<5.5

5 15 27 34 22

5.5鈭<6.0 6.0鈭<6.5 6.5鈭<7.0 7.0鈭<7.5 7.5鈭<8.0

14 7 2 4 1

a. What proportion of the observations are less than 5?

b. What proportion of the observations are at least 6?

c. Construct a histogram with relative frequency on the vertical axis and comment on interesting features. In particular, does the distribution of particle sizes appear to be reasonably symmetric or somewhat skewed? (Note:The investigators fit lognormaldistribution to the data; this is discussed in Chapter 4.)

d. Construct a histogram with density on the vertical axis and compare to the histogram in (c).

Once an individual has been infected with a certain disease, let \({\rm{X}}\) represent the time (days) that elapses before the individual becomes infectious. The article proposes a Weibull distribution with \({\rm{\alpha = 2}}{\rm{.2}}\), \({\rm{\beta = 1}}{\rm{.1}}\), and \({\rm{\gamma = 0}}{\rm{.5}}\).

a. Calculate \({\rm{P(1 < X < 2)}}\).

b. Calculate \({\rm{P(X > 1}}{\rm{.5)}}\).

c. What is the \({\rm{90th}}\) percentile of the distribution?

d. What are the mean and standard deviation of \({\rm{X}}\)?

Observations on burst strength (lb/in 2) were obtained both for test nozzle closure welds and for production cannister nozzle welds (鈥淧roper Procedures Are the Key to Welding Radioactive Waste Cannisters,鈥
Welding J., Aug. 1997: 61鈥67).
Test 7200 6100 7300 7300 8000 7400
7300 7300 8000 6700 8300
Cannister 5250 5625 5900 5900 5700 6050
5800 6000 5875 6100 5850 6600


Construct a comparative boxplot and comment on interesting features (the cited article did not include such a picture, but the authors commented that they had looked
at one).

The article cited in Example 1.2 also gave the accompanying strength observations for cylinders:

6.1

5.8

7.8

7.1

7.2

9.2

6.6

8.3

7.0

8.3

7.8

8.1

7.4

8.5

8.9

9.8

9.7

14.1

12.6

11.2


a. Construct a comparative stem-and-leaf display(see the previous exercise) of the beam and cylinder data, and then answer the questions in parts(b)鈥(d) of Exercise 10 for the observations oncylinders.

b. In what ways are the two sides of the display similar? Are there any obvious differences between the beam observations and the cylinder observations?
c. Construct a dotplot of the cylinder data.

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.