/*! 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} Q33E At time \({\rm{t = 0}}\), there ... [FREE SOLUTION] | 91影视

91影视

At time \({\rm{t = 0}}\), there is one individual alive in a certain population. A pure birth process then unfolds as follows. The time until the first birth is exponentially distributed with parameter \({\rm{\lambda }}\). After the first birth, there are two individuals alive. The time until the first gives birth again is exponential with parameter \({\rm{\lambda }}\), and similarly for the second individual. Therefore, the time until the next birth is the minimum of two exponential (\({\rm{\lambda }}\)) variables, which is exponential with parameter \({\rm{2\lambda }}\). Similarly, once the second birth has occurred, there are three individuals alive, so the time until the next birth is an exponential \({\rm{rv}}\) with parameter \({\rm{3\lambda }}\), and so on (the memoryless property of the exponential distribution is being used here). Suppose the process is observed until the sixth birth has occurred and the successive birth times are \({\rm{25}}{\rm{.2,41}}{\rm{.7,51}}{\rm{.2,55}}{\rm{.5,59}}{\rm{.5,61}}{\rm{.8}}\) (from which you should calculate the times between successive births). Derive the mle of l. (Hint: The likelihood is a product of exponential terms.)

Short Answer

Expert verified

The value is \({\rm{\hat \lambda = 0}}{\rm{.0436}}\).

Step by step solution

01

Define exponential function

A function that increases or decays at a rate proportional to its present value is called an exponential function.

02

Explanation

As, X is a random variable that can be used with pdf,

\({{\rm{f}}_{\rm{X}}}{\rm{(x) = }}\left\{ {\begin{array}{*{20}{l}}{{\rm{\lambda }}{{\rm{e}}^{{\rm{ - \lambda x}}}}}&{{\rm{,x}} \ge {\rm{0}}}\\{\rm{0}}&{{\rm{,x < 0}}}\end{array}} \right.\)

With parameter\({\rm{\lambda }}\), it is said to have an exponential distribution.

Denote with

\({{\rm{X}}_{\rm{1}}}\)= time till the first birth,

\({{\rm{X}}_{\rm{2}}}{\rm{ = }}\) time between the first and second births,

.

.

\({{\rm{X}}_{\rm{n}}}{\rm{ = }}\) time between the \({{\rm{n}}^{{\rm{th }}}}\) and \({{\rm{(n - 1)}}^{{\rm{th}}}}\) births,

As, \({{\rm{X}}_{\rm{i}}}{\rm{,i = 1,2, \ldots ,n}}\) are independent and exponentially distributed random variables with appropriate parameters \({\rm{ - }}{{\rm{X}}_{\rm{1}}}\) with parameter \({\rm{\lambda }}\), \({{\rm{X}}_{\rm{2}}}\) with parameter \({\rm{2\lambda }}\),..., \({{\rm{X}}_{\rm{n}}}\) with parameter \({\rm{n\lambda }}\).

Allow joint pdf or pmb 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}}}\).

As, \({{\rm{X}}_{\rm{i}}}{\rm{,i = 1,2, \ldots ,n}}\) are independent and exponentially distributed random variables with appropriate parameters \({\rm{ - }}{{\rm{X}}_{\rm{1}}}\) with parameter \({\rm{\lambda }}\), \({{\rm{X}}_{\rm{2}}}\)

03

Evaluating the maximum likelihood estimate

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{;\lambda }}} \right){\rm{ = \lambda }}{{\rm{e}}^{{\rm{ - \lambda }}{{\rm{x}}_{\rm{1}}}}}{\rm{ \times (2\lambda )}}{{\rm{e}}^{{\rm{ - 2\lambda }}{{\rm{x}}_{\rm{2}}}}}{\rm{ \times \ldots \times (n\lambda )}}{{\rm{e}}^{{\rm{ - n\lambda }}\left( {{{\rm{x}}_{\rm{n}}}{\rm{ - \theta }}} \right)}}\\{\rm{ = n \times (n - 1) \times \ldots \times 1 \times }}{{\rm{\lambda }}^{\rm{n}}}{\rm{ \times exp}}\left\{ {{\rm{\lambda }}\sum\limits_{{\rm{i = 1}}}^{\rm{n}} {\rm{i}} {{\rm{x}}_{\rm{i}}}} \right\}\\{\rm{ = n! \times }}{{\rm{\lambda }}^{\rm{n}}}{\rm{ \times exp}}\left\{ {{\rm{\lambda }}\sum\limits_{{\rm{i = 1}}}^{\rm{n}} {\rm{i}} {{\rm{x}}_{\rm{i}}}} \right\}\end{array}\)

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{;\lambda ,}}} \right){\rm{ = ln}}\left( {{\rm{n! \times}}{{\rm{\lambda }}^{\rm{n}}}{\rm{ \times exp}}\left\{ {{\rm{\lambda }}\sum\limits_{{\rm{i = 1}}}^{\rm{n}} {\rm{i}} {{\rm{x}}_{\rm{i}}}} \right\}} \right)\\{\rm{ = lnn! + nln\lambda - \lambda }}\sum\limits_{{\rm{i = 1}}}^{\rm{n}} {\rm{i}} {{\rm{x}}_{\rm{i}}}\end{array}\)

The maximum likelihood estimator is generated by taking the derivative of the log likelihood function in regard to \({\rm{\lambda }}\) 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{;\lambda }}} \right){\rm{ = }}\frac{{\rm{d}}}{{{\rm{d\lambda }}}}\left( {{\rm{lnn! + nln\lambda - \lambda }}\sum\limits_{{\rm{i = 1}}}^{\rm{n}} {\rm{i}} {{\rm{x}}_{\rm{i}}}} \right)\\{\rm{ = }}\frac{{\rm{n}}}{{\rm{\lambda }}}{\rm{ - }}\sum\limits_{{\rm{i = 1}}}^{\rm{n}} {\rm{i}} {{\rm{x}}_{\rm{i}}}\end{array}\)

As a result, solving equation provides the maximum likelihood estimator\(\widehat {\rm{\lambda }}\).

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

For\({\rm{\hat \lambda }}\). Hence, the maximum likelihood estimator of parameter\({\rm{\lambda }}\)is,

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

04

Evaluating the value

We need the values x, which are the times between successive births, based on the serial birth times. The following formula is used to calculate the values:

\({{\rm{x}}_{\rm{1}}}{\rm{ = 25}}{\rm{.2}}\) the time of the first birth;

\({{\rm{x}}_{\rm{2}}}{\rm{ = 41}}{\rm{.7 - 25}}{\rm{.2 = 16}}{\rm{.5}}\) the period of time between the first and second births;

\({{\rm{x}}_{\rm{3}}}{\rm{ = 51}}{\rm{.2 - 41}}{\rm{.7 = 9}}{\rm{.5}}\) the time between the second and third births;

\({{\rm{x}}_{\rm{4}}}{\rm{ = 55}}{\rm{.5 - 51 = 4}}{\rm{.3}}\) the time between the third and fourth births;

\({{\rm{x}}_{\rm{5}}}{\rm{ = 59}}{\rm{.5 - 55}}{\rm{.5 = 4}}\) the time between the fourth and fifth births;

\({{\rm{x}}_{\rm{6}}}{\rm{ = 61}}{\rm{.8 - 59}}{\rm{.5 = 2}}{\rm{.3}}\) the time between the fifth and sixth births;

As a result, the total can now be calculated as follows:

\(\begin{array}{c}\sum\limits_{{\rm{i = 1}}}^{\rm{n}} {\rm{i}} {{\rm{x}}_{\rm{i}}}{\rm{ = 1 \times }}{{\rm{x}}_{\rm{1}}}{\rm{ + 2 \times }}{{\rm{x}}_{\rm{2}}}{\rm{ + \ldots + 6 \times }}{{\rm{x}}_{\rm{6}}}\\{\rm{ = 1 \times 25}}{\rm{.2 + 2 \times 16}}{\rm{.5 + \ldots + 6 \times 2}}{\rm{.3}}\\{\rm{ = 137}}{\rm{.7}}\end{array}\)

Last but not least, the maximum likelihood estimates of\({\rm{\lambda }}\)is,

\(\begin{array}{c}{\rm{\hat \lambda = }}\frac{{\rm{n}}}{{\sum\limits_{{\rm{i = 1}}}^{\rm{n}} {\rm{i}} {{\rm{x}}_{\rm{i}}}}}\\{\rm{ = }}\frac{{\rm{6}}}{{{\rm{137}}{\rm{.7}}}}\\{\rm{ = 0}}{\rm{.0436}}\end{array}\)

Therefore, \({\rm{\hat \lambda = 0}}{\rm{.0436}}\).

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

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}}\))

The shear strength of each of ten test spot welds is determined, yielding the following data (psi):

\(\begin{array}{*{20}{l}}{{\rm{392}}}&{{\rm{376}}}&{{\rm{401}}}&{{\rm{367}}}&{{\rm{389}}}&{{\rm{362}}}&{{\rm{409}}}&{{\rm{415}}}&{{\rm{358}}}&{{\rm{375}}}\end{array}\)

a. Assuming that shear strength is normally distributed, estimate the true average shear strength and standard deviation of shear strength using the method of maximum likelihood.

b. Again assuming a normal distribution, estimate the strength value below which\({\rm{95\% }}\)of all welds will have their strengths. (Hint: What is the\({\rm{95 th}}\)percentile in terms of\({\rm{\mu }}\)and\({\rm{\sigma }}\)? Now use the invariance principle.)

c. Suppose we decide to examine another test spot weld. Let\({\rm{X = }}\)shear strength of the weld. Use the given data to obtain the mle of\({\rm{P(X拢400)}}{\rm{.(Hint:P(X拢400) = \Phi ((400 - \mu )/\sigma )}}{\rm{.)}}\)

We defined a negative binomial\({\rm{rv}}\)as the number of failures that occur before the\({\rm{rth}}\)success in a sequence of independent and identical success/failure trials. The probability mass function (\({\rm{pmf}}\)) of\({\rm{X}}\)is\({\rm{nb(x,r,p) = }}\)\(\left( {\begin{array}{*{20}{c}}{{\rm{x + r - 1}}}\\{\rm{x}}\end{array}} \right){{\rm{p}}^{\rm{r}}}{{\rm{(1 - p)}}^{\rm{x}}}\quad {\rm{x = 0,1,2, \ldots }}\)

a. Suppose that. Show that\({\rm{\hat p = (r - 1)/(X + r - 1)}}\)is an unbiased estimator for\({\rm{p}}\). (Hint: Write out\({\rm{E(\hat p)}}\)and cancel\({\rm{x + r - 1}}\)inside the sum.)

b. A reporter wishing to interview five individuals who support a certain candidate begins asking people whether\({\rm{(S)}}\)or not\({\rm{(F)}}\)they support the candidate. If the sequence of responses is SFFSFFFSSS, estimate\({\rm{p = }}\)the true proportion who support the candidate.

\({{\rm{X}}_{\rm{1}}}{\rm{,}}.....{\rm{,}}{{\rm{X}}_{\rm{n}}}\)be a random sample from a gamma distribution with parameters \({\rm{\alpha }}\) and \({\rm{\beta }}\). a. Derive the equations whose solutions yield the maximum likelihood estimators of \({\rm{\alpha }}\) and \({\rm{\beta }}\). Do you think they can be solved explicitly? b. Show that the mle of \({\rm{\mu = \alpha \beta }}\) is \(\widehat {\rm{\mu }}{\rm{ = }}\overline {\rm{X}} \).

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)}}\)?)

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.