How do we build a confidence interval for the parameter of the exponential distribution? Announcing the arrival of Valued Associate #679: Cesar Manara Planned maintenance scheduled April 17/18, 2019 at 00:00UTC (8:00pm US/Eastern)How to compute confidence interval from a confidence distributionparameter and prediction confidence intervalsConfidence interval for known non-normal estimation?Confidence interval for exponential distributionA confidence area for an Archimedean's copula familyIs the canonical parameter (and therefore the canonical link function) for a Gamma not unique?UMAU confidence interval for $theta$ in a shifted exponential distributionCalculate the constants and the MSE from two estimators related to a uniform distributionBinomial distributed random sample: find the least variance from the set of all unbiased estimators of $theta$Build an approximated confidence interval for $sigma$ based on its maximum likelihood estimator

The following signatures were invalid: EXPKEYSIG 1397BC53640DB551

Antler Helmet: Can it work?

Can a non-EU citizen traveling with me come with me through the EU passport line?

How should I respond to a player wanting to catch a sword between their hands?

Slither Like a Snake

Area of a 2D convex hull

Why does tar appear to skip file contents when output file is /dev/null?

How do we build a confidence interval for the parameter of the exponential distribution?

Direct Experience of Meditation

Is dark matter really a meaningful hypothesis?

Passing functions in C++

How do I keep my slimes from escaping their pens?

Single author papers against my advisor's will?

If I can make up priors, why can't I make up posteriors?

How to say 'striped' in Latin

Blender game recording at the wrong time

What LEGO pieces have "real-world" functionality?

Why does this iterative way of solving of equation work?

Estimate capacitor parameters

Mortgage adviser recommends a longer term than necessary combined with overpayments

If A makes B more likely then B makes A more likely"

What is the order of Mitzvot in Rambam's Sefer Hamitzvot?

Did the new image of black hole confirm the general theory of relativity?

What do you call a plan that's an alternative plan in case your initial plan fails?



How do we build a confidence interval for the parameter of the exponential distribution?



Announcing the arrival of Valued Associate #679: Cesar Manara
Planned maintenance scheduled April 17/18, 2019 at 00:00UTC (8:00pm US/Eastern)How to compute confidence interval from a confidence distributionparameter and prediction confidence intervalsConfidence interval for known non-normal estimation?Confidence interval for exponential distributionA confidence area for an Archimedean's copula familyIs the canonical parameter (and therefore the canonical link function) for a Gamma not unique?UMAU confidence interval for $theta$ in a shifted exponential distributionCalculate the constants and the MSE from two estimators related to a uniform distributionBinomial distributed random sample: find the least variance from the set of all unbiased estimators of $theta$Build an approximated confidence interval for $sigma$ based on its maximum likelihood estimator



.everyoneloves__top-leaderboard:empty,.everyoneloves__mid-leaderboard:empty,.everyoneloves__bot-mid-leaderboard:empty margin-bottom:0;








1












$begingroup$


EDIT



Let $X_1,X_2,ldots,X_n$ be a random sample whose distribution is given by $textExp(theta)$, where $theta$ is not known. Precisely, $f(x|theta) = (1/thetaexp)(-x/theta)$ Describe a method to build a confidence interval with confidence coefficient $1 - alpha$ for $theta$.



MY ATTEMPT



Since the distribution in discussion is not normal and I do not know the size of the sample, I think we cannot apply the central limit theorem. One possible approach is to consider the maximum likelihood estimator of $theta$, whose distribution is approximately $mathcalN(theta,(nI_F(theta)^-1)$. Another possible approach consists in using the score function, whose distribution is approximately $mathcalN(0,nI_F(theta))$. However, in both cases, it is assumed the CLT is applicable.



The exercise also provides the following hint: find $c_1$ and $c_2$ such that
beginalign*
textbfPleft(c_1 < frac1thetasum_i=1^n X_i < c_2right) = 1 -alpha
endalign*



Can someone help me out? Thanks in advance!










share|cite|improve this question











$endgroup$







  • 1




    $begingroup$
    You should clarify which parameterization of the exponential distribution you're using. From the later parts of your post it looks like you're using the scale parameterization rather than the rate parameterization but you should be explicit, not leave it to people to guess.
    $endgroup$
    – Glen_b
    4 hours ago










  • $begingroup$
    Thanks for the comment and sorry for the inconvenience. I edited the question.
    $endgroup$
    – user1337
    4 hours ago






  • 1




    $begingroup$
    Okay, you've defined it as the rate parameterization, which is fine, but then the hint at the end is wrong.
    $endgroup$
    – Glen_b
    4 hours ago











  • $begingroup$
    For rather large $n$ an approach using the CLT might provide a useful approximation. My answer gives an exact CI that works even for small $n.$
    $endgroup$
    – BruceET
    3 hours ago










  • $begingroup$
    There are so many options here because there are different choices of pivots. A C.I. could also be found using $min X_i$ which also has an exp distribution, but this won't be as 'good' as the one based on $sum X_i$.
    $endgroup$
    – StubbornAtom
    27 mins ago

















1












$begingroup$


EDIT



Let $X_1,X_2,ldots,X_n$ be a random sample whose distribution is given by $textExp(theta)$, where $theta$ is not known. Precisely, $f(x|theta) = (1/thetaexp)(-x/theta)$ Describe a method to build a confidence interval with confidence coefficient $1 - alpha$ for $theta$.



MY ATTEMPT



Since the distribution in discussion is not normal and I do not know the size of the sample, I think we cannot apply the central limit theorem. One possible approach is to consider the maximum likelihood estimator of $theta$, whose distribution is approximately $mathcalN(theta,(nI_F(theta)^-1)$. Another possible approach consists in using the score function, whose distribution is approximately $mathcalN(0,nI_F(theta))$. However, in both cases, it is assumed the CLT is applicable.



The exercise also provides the following hint: find $c_1$ and $c_2$ such that
beginalign*
textbfPleft(c_1 < frac1thetasum_i=1^n X_i < c_2right) = 1 -alpha
endalign*



Can someone help me out? Thanks in advance!










share|cite|improve this question











$endgroup$







  • 1




    $begingroup$
    You should clarify which parameterization of the exponential distribution you're using. From the later parts of your post it looks like you're using the scale parameterization rather than the rate parameterization but you should be explicit, not leave it to people to guess.
    $endgroup$
    – Glen_b
    4 hours ago










  • $begingroup$
    Thanks for the comment and sorry for the inconvenience. I edited the question.
    $endgroup$
    – user1337
    4 hours ago






  • 1




    $begingroup$
    Okay, you've defined it as the rate parameterization, which is fine, but then the hint at the end is wrong.
    $endgroup$
    – Glen_b
    4 hours ago











  • $begingroup$
    For rather large $n$ an approach using the CLT might provide a useful approximation. My answer gives an exact CI that works even for small $n.$
    $endgroup$
    – BruceET
    3 hours ago










  • $begingroup$
    There are so many options here because there are different choices of pivots. A C.I. could also be found using $min X_i$ which also has an exp distribution, but this won't be as 'good' as the one based on $sum X_i$.
    $endgroup$
    – StubbornAtom
    27 mins ago













1












1








1


0



$begingroup$


EDIT



Let $X_1,X_2,ldots,X_n$ be a random sample whose distribution is given by $textExp(theta)$, where $theta$ is not known. Precisely, $f(x|theta) = (1/thetaexp)(-x/theta)$ Describe a method to build a confidence interval with confidence coefficient $1 - alpha$ for $theta$.



MY ATTEMPT



Since the distribution in discussion is not normal and I do not know the size of the sample, I think we cannot apply the central limit theorem. One possible approach is to consider the maximum likelihood estimator of $theta$, whose distribution is approximately $mathcalN(theta,(nI_F(theta)^-1)$. Another possible approach consists in using the score function, whose distribution is approximately $mathcalN(0,nI_F(theta))$. However, in both cases, it is assumed the CLT is applicable.



The exercise also provides the following hint: find $c_1$ and $c_2$ such that
beginalign*
textbfPleft(c_1 < frac1thetasum_i=1^n X_i < c_2right) = 1 -alpha
endalign*



Can someone help me out? Thanks in advance!










share|cite|improve this question











$endgroup$




EDIT



Let $X_1,X_2,ldots,X_n$ be a random sample whose distribution is given by $textExp(theta)$, where $theta$ is not known. Precisely, $f(x|theta) = (1/thetaexp)(-x/theta)$ Describe a method to build a confidence interval with confidence coefficient $1 - alpha$ for $theta$.



MY ATTEMPT



Since the distribution in discussion is not normal and I do not know the size of the sample, I think we cannot apply the central limit theorem. One possible approach is to consider the maximum likelihood estimator of $theta$, whose distribution is approximately $mathcalN(theta,(nI_F(theta)^-1)$. Another possible approach consists in using the score function, whose distribution is approximately $mathcalN(0,nI_F(theta))$. However, in both cases, it is assumed the CLT is applicable.



The exercise also provides the following hint: find $c_1$ and $c_2$ such that
beginalign*
textbfPleft(c_1 < frac1thetasum_i=1^n X_i < c_2right) = 1 -alpha
endalign*



Can someone help me out? Thanks in advance!







self-study confidence-interval exponential-distribution






share|cite|improve this question















share|cite|improve this question













share|cite|improve this question




share|cite|improve this question








edited 1 hour ago







user1337

















asked 7 hours ago









user1337user1337

1795




1795







  • 1




    $begingroup$
    You should clarify which parameterization of the exponential distribution you're using. From the later parts of your post it looks like you're using the scale parameterization rather than the rate parameterization but you should be explicit, not leave it to people to guess.
    $endgroup$
    – Glen_b
    4 hours ago










  • $begingroup$
    Thanks for the comment and sorry for the inconvenience. I edited the question.
    $endgroup$
    – user1337
    4 hours ago






  • 1




    $begingroup$
    Okay, you've defined it as the rate parameterization, which is fine, but then the hint at the end is wrong.
    $endgroup$
    – Glen_b
    4 hours ago











  • $begingroup$
    For rather large $n$ an approach using the CLT might provide a useful approximation. My answer gives an exact CI that works even for small $n.$
    $endgroup$
    – BruceET
    3 hours ago










  • $begingroup$
    There are so many options here because there are different choices of pivots. A C.I. could also be found using $min X_i$ which also has an exp distribution, but this won't be as 'good' as the one based on $sum X_i$.
    $endgroup$
    – StubbornAtom
    27 mins ago












  • 1




    $begingroup$
    You should clarify which parameterization of the exponential distribution you're using. From the later parts of your post it looks like you're using the scale parameterization rather than the rate parameterization but you should be explicit, not leave it to people to guess.
    $endgroup$
    – Glen_b
    4 hours ago










  • $begingroup$
    Thanks for the comment and sorry for the inconvenience. I edited the question.
    $endgroup$
    – user1337
    4 hours ago






  • 1




    $begingroup$
    Okay, you've defined it as the rate parameterization, which is fine, but then the hint at the end is wrong.
    $endgroup$
    – Glen_b
    4 hours ago











  • $begingroup$
    For rather large $n$ an approach using the CLT might provide a useful approximation. My answer gives an exact CI that works even for small $n.$
    $endgroup$
    – BruceET
    3 hours ago










  • $begingroup$
    There are so many options here because there are different choices of pivots. A C.I. could also be found using $min X_i$ which also has an exp distribution, but this won't be as 'good' as the one based on $sum X_i$.
    $endgroup$
    – StubbornAtom
    27 mins ago







1




1




$begingroup$
You should clarify which parameterization of the exponential distribution you're using. From the later parts of your post it looks like you're using the scale parameterization rather than the rate parameterization but you should be explicit, not leave it to people to guess.
$endgroup$
– Glen_b
4 hours ago




$begingroup$
You should clarify which parameterization of the exponential distribution you're using. From the later parts of your post it looks like you're using the scale parameterization rather than the rate parameterization but you should be explicit, not leave it to people to guess.
$endgroup$
– Glen_b
4 hours ago












$begingroup$
Thanks for the comment and sorry for the inconvenience. I edited the question.
$endgroup$
– user1337
4 hours ago




$begingroup$
Thanks for the comment and sorry for the inconvenience. I edited the question.
$endgroup$
– user1337
4 hours ago




1




1




$begingroup$
Okay, you've defined it as the rate parameterization, which is fine, but then the hint at the end is wrong.
$endgroup$
– Glen_b
4 hours ago





$begingroup$
Okay, you've defined it as the rate parameterization, which is fine, but then the hint at the end is wrong.
$endgroup$
– Glen_b
4 hours ago













$begingroup$
For rather large $n$ an approach using the CLT might provide a useful approximation. My answer gives an exact CI that works even for small $n.$
$endgroup$
– BruceET
3 hours ago




$begingroup$
For rather large $n$ an approach using the CLT might provide a useful approximation. My answer gives an exact CI that works even for small $n.$
$endgroup$
– BruceET
3 hours ago












$begingroup$
There are so many options here because there are different choices of pivots. A C.I. could also be found using $min X_i$ which also has an exp distribution, but this won't be as 'good' as the one based on $sum X_i$.
$endgroup$
– StubbornAtom
27 mins ago




$begingroup$
There are so many options here because there are different choices of pivots. A C.I. could also be found using $min X_i$ which also has an exp distribution, but this won't be as 'good' as the one based on $sum X_i$.
$endgroup$
– StubbornAtom
27 mins ago










2 Answers
2






active

oldest

votes


















2












$begingroup$

You don't say how the exponential distribution is
parameterized. Two parameterizations are in common use--mean and rate.



Let $E(X_i) = mu.$ Then one
can show that $$frac 1 mu sum_i=1^n X_i sim
mathsfGamma(textshape = n, textrate=scale = 1).$$



In R statistical software the exponential distribution is parameterized according rate $lambda = 1/mu.$ Let $n = 10$ and $lambda = 1/5,$ so that $mu = 5.$ The following program simulates $m = 10^6$ samples of size $n = 10$ from $mathsfExp(textrate = lambda = 1/5),$ finds $$Q = frac 1 mu sum_i=1^n X_i =
lambda sum_i=1^n X_i$$
for each sample, and plots the histogram of the one million $Q$'s, The figure
illustrates that $Q sim mathsfGamma(10, 1).$
(Use MGFs for a formal proof.)



set.seed(414) # for reproducibility
q = replicate(10^5, sum(rexp(10, 1/5))/5)
lbl = "Simulated Dist'n of Q with Density of GAMMA(10, 1)"
hist(q, prob=T, br=30, col="skyblue2", main=lbl)
curve(dgamma(x,10,1), col="red", add=T)


enter image description here



Thus, for $n = 10$ the constants $c_1 = 4.975$ and
$c_2 = 17.084$ for
a 95% confidence interval are quantiles 0.025 and 0.975, respectively, of $Q sim mathsfGamma(10, 1).$



qgamma(c(.025, .975), 10, 1)
[1] 4.795389 17.084803


In particular, for the exponential sample shown below (second row),
a 95% confidence interval is $(2.224, 7.922).$ Notice the reversal of the quantiles in 'pivoting' $Q,$ which
has $mu$ in the denominator.



set.seed(1234); x = sort(round(rexp(10, 1/5), 2)); x
[1] 0.03 0.45 1.01 1.23 1.94 3.80 4.12 4.19 8.71 12.51
t = sum(x); t
[1] 37.99
t/qgamma(c(.975, .025), 10, 1)
[1] 2.223614 7.922194


Note: Because the chi-squared distribution is a member of the gamma family, it is possible to find endpoints for such a confidence interval in terms of a chi-squared distribution.



See Wikipedia on exponential distributions under 'confidence intervals'. (That discussion uses rate parameter $lambda$ for the exponential distribution, instead of $mu.)$






share|cite|improve this answer











$endgroup$




















    1












    $begingroup$

    Taking $theta$ as the scale parameter, it can be shown that:



    $$fracn barXtheta sim textGa(n,1).$$



    To form a confidence interval we choose any critical points $c_1 < c_2$ from the $textGa(n,1)$ distribution such that these points contain probability $1-alpha$ of the distribution. Using the above pivotal quantity we then have:



    $$mathbbP Bigg( c_1 leqslant fracn barXtheta leqslant c_2 Bigg) = 1-alpha
    quad quad quad quad quad
    int limits_c_1^c_2 textGa(r|n,1) dr = 1 - alpha.$$



    Re-arranging the inequality in this probability statement and substituting the observed sample mean gives the confidence interval:



    $$textCI_theta(1-alpha) = Bigg[ fracn barxc_2 , fracn barxc_1 Bigg].$$



    This confidence interval is valid for any choice of $c_1<c_2$ so long as it obeys the required integral condition. For simplicity, many analysts use the symmetric critical points. However, it is possible to optimise the confidence interval by minimising its length, which we show below.




    Optimising the confidence interval: The length of this confidence interval is proportional to $1/c_1-1/c_2$, and so we minimise the length of the interval by choosing the critical points to minimise this distance. This can be done using the nlm function in R. In the following code we give a function for the minimum-length confidence interval for this problem, which we apply to some simulated data.



    #Set the objective function for minimisation
    OBJECTIVE <- function(c1, n, alpha)
    pp <- pgamma(c1, n, 1, lower.tail = TRUE);
    c2 <- qgamma(1 - alpha + pp, n, 1, lower.tail = TRUE);
    1/c1 - 1/c2;

    #Find the minimum-length confidence interval
    CONF_INT <- function(n, alpha, xbar)
    START_c1 <- qgamma(alpha/2, n, 1, lower.tail = TRUE);
    MINIMISE <- nlm(f = OBJECTIVE, p = START_c1, n = n, alpha = alpha);
    c1 <- MINIMISE$estimate;
    pp <- pgamma(c1, n, 1, lower.tail = TRUE);
    c2 <- qgamma(1 - alpha + pp, n, 1, lower.tail = TRUE);
    c(n*xbar/c2, n*xbar/c1);

    #Generate simulation data
    set.seed(921730198);
    n <- 300;
    scale <- 25.4;
    DATA <- rexp(n, rate = 1/scale);

    #Application of confidence interval to simulated data
    n <- length(DATA);
    xbar <- mean(DATA);
    alpha <- 0.05;

    CONF_INT(n, alpha, xbar);

    [1] 23.32040 29.24858





    share|cite|improve this answer











    $endgroup$













      Your Answer








      StackExchange.ready(function()
      var channelOptions =
      tags: "".split(" "),
      id: "65"
      ;
      initTagRenderer("".split(" "), "".split(" "), channelOptions);

      StackExchange.using("externalEditor", function()
      // Have to fire editor after snippets, if snippets enabled
      if (StackExchange.settings.snippets.snippetsEnabled)
      StackExchange.using("snippets", function()
      createEditor();
      );

      else
      createEditor();

      );

      function createEditor()
      StackExchange.prepareEditor(
      heartbeatType: 'answer',
      autoActivateHeartbeat: false,
      convertImagesToLinks: false,
      noModals: true,
      showLowRepImageUploadWarning: true,
      reputationToPostImages: null,
      bindNavPrevention: true,
      postfix: "",
      imageUploader:
      brandingHtml: "Powered by u003ca class="icon-imgur-white" href="https://imgur.com/"u003eu003c/au003e",
      contentPolicyHtml: "User contributions licensed under u003ca href="https://creativecommons.org/licenses/by-sa/3.0/"u003ecc by-sa 3.0 with attribution requiredu003c/au003e u003ca href="https://stackoverflow.com/legal/content-policy"u003e(content policy)u003c/au003e",
      allowUrls: true
      ,
      onDemand: true,
      discardSelector: ".discard-answer"
      ,immediatelyShowMarkdownHelp:true
      );



      );













      draft saved

      draft discarded


















      StackExchange.ready(
      function ()
      StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fstats.stackexchange.com%2fquestions%2f403059%2fhow-do-we-build-a-confidence-interval-for-the-parameter-of-the-exponential-distr%23new-answer', 'question_page');

      );

      Post as a guest















      Required, but never shown

























      2 Answers
      2






      active

      oldest

      votes








      2 Answers
      2






      active

      oldest

      votes









      active

      oldest

      votes






      active

      oldest

      votes









      2












      $begingroup$

      You don't say how the exponential distribution is
      parameterized. Two parameterizations are in common use--mean and rate.



      Let $E(X_i) = mu.$ Then one
      can show that $$frac 1 mu sum_i=1^n X_i sim
      mathsfGamma(textshape = n, textrate=scale = 1).$$



      In R statistical software the exponential distribution is parameterized according rate $lambda = 1/mu.$ Let $n = 10$ and $lambda = 1/5,$ so that $mu = 5.$ The following program simulates $m = 10^6$ samples of size $n = 10$ from $mathsfExp(textrate = lambda = 1/5),$ finds $$Q = frac 1 mu sum_i=1^n X_i =
      lambda sum_i=1^n X_i$$
      for each sample, and plots the histogram of the one million $Q$'s, The figure
      illustrates that $Q sim mathsfGamma(10, 1).$
      (Use MGFs for a formal proof.)



      set.seed(414) # for reproducibility
      q = replicate(10^5, sum(rexp(10, 1/5))/5)
      lbl = "Simulated Dist'n of Q with Density of GAMMA(10, 1)"
      hist(q, prob=T, br=30, col="skyblue2", main=lbl)
      curve(dgamma(x,10,1), col="red", add=T)


      enter image description here



      Thus, for $n = 10$ the constants $c_1 = 4.975$ and
      $c_2 = 17.084$ for
      a 95% confidence interval are quantiles 0.025 and 0.975, respectively, of $Q sim mathsfGamma(10, 1).$



      qgamma(c(.025, .975), 10, 1)
      [1] 4.795389 17.084803


      In particular, for the exponential sample shown below (second row),
      a 95% confidence interval is $(2.224, 7.922).$ Notice the reversal of the quantiles in 'pivoting' $Q,$ which
      has $mu$ in the denominator.



      set.seed(1234); x = sort(round(rexp(10, 1/5), 2)); x
      [1] 0.03 0.45 1.01 1.23 1.94 3.80 4.12 4.19 8.71 12.51
      t = sum(x); t
      [1] 37.99
      t/qgamma(c(.975, .025), 10, 1)
      [1] 2.223614 7.922194


      Note: Because the chi-squared distribution is a member of the gamma family, it is possible to find endpoints for such a confidence interval in terms of a chi-squared distribution.



      See Wikipedia on exponential distributions under 'confidence intervals'. (That discussion uses rate parameter $lambda$ for the exponential distribution, instead of $mu.)$






      share|cite|improve this answer











      $endgroup$

















        2












        $begingroup$

        You don't say how the exponential distribution is
        parameterized. Two parameterizations are in common use--mean and rate.



        Let $E(X_i) = mu.$ Then one
        can show that $$frac 1 mu sum_i=1^n X_i sim
        mathsfGamma(textshape = n, textrate=scale = 1).$$



        In R statistical software the exponential distribution is parameterized according rate $lambda = 1/mu.$ Let $n = 10$ and $lambda = 1/5,$ so that $mu = 5.$ The following program simulates $m = 10^6$ samples of size $n = 10$ from $mathsfExp(textrate = lambda = 1/5),$ finds $$Q = frac 1 mu sum_i=1^n X_i =
        lambda sum_i=1^n X_i$$
        for each sample, and plots the histogram of the one million $Q$'s, The figure
        illustrates that $Q sim mathsfGamma(10, 1).$
        (Use MGFs for a formal proof.)



        set.seed(414) # for reproducibility
        q = replicate(10^5, sum(rexp(10, 1/5))/5)
        lbl = "Simulated Dist'n of Q with Density of GAMMA(10, 1)"
        hist(q, prob=T, br=30, col="skyblue2", main=lbl)
        curve(dgamma(x,10,1), col="red", add=T)


        enter image description here



        Thus, for $n = 10$ the constants $c_1 = 4.975$ and
        $c_2 = 17.084$ for
        a 95% confidence interval are quantiles 0.025 and 0.975, respectively, of $Q sim mathsfGamma(10, 1).$



        qgamma(c(.025, .975), 10, 1)
        [1] 4.795389 17.084803


        In particular, for the exponential sample shown below (second row),
        a 95% confidence interval is $(2.224, 7.922).$ Notice the reversal of the quantiles in 'pivoting' $Q,$ which
        has $mu$ in the denominator.



        set.seed(1234); x = sort(round(rexp(10, 1/5), 2)); x
        [1] 0.03 0.45 1.01 1.23 1.94 3.80 4.12 4.19 8.71 12.51
        t = sum(x); t
        [1] 37.99
        t/qgamma(c(.975, .025), 10, 1)
        [1] 2.223614 7.922194


        Note: Because the chi-squared distribution is a member of the gamma family, it is possible to find endpoints for such a confidence interval in terms of a chi-squared distribution.



        See Wikipedia on exponential distributions under 'confidence intervals'. (That discussion uses rate parameter $lambda$ for the exponential distribution, instead of $mu.)$






        share|cite|improve this answer











        $endgroup$















          2












          2








          2





          $begingroup$

          You don't say how the exponential distribution is
          parameterized. Two parameterizations are in common use--mean and rate.



          Let $E(X_i) = mu.$ Then one
          can show that $$frac 1 mu sum_i=1^n X_i sim
          mathsfGamma(textshape = n, textrate=scale = 1).$$



          In R statistical software the exponential distribution is parameterized according rate $lambda = 1/mu.$ Let $n = 10$ and $lambda = 1/5,$ so that $mu = 5.$ The following program simulates $m = 10^6$ samples of size $n = 10$ from $mathsfExp(textrate = lambda = 1/5),$ finds $$Q = frac 1 mu sum_i=1^n X_i =
          lambda sum_i=1^n X_i$$
          for each sample, and plots the histogram of the one million $Q$'s, The figure
          illustrates that $Q sim mathsfGamma(10, 1).$
          (Use MGFs for a formal proof.)



          set.seed(414) # for reproducibility
          q = replicate(10^5, sum(rexp(10, 1/5))/5)
          lbl = "Simulated Dist'n of Q with Density of GAMMA(10, 1)"
          hist(q, prob=T, br=30, col="skyblue2", main=lbl)
          curve(dgamma(x,10,1), col="red", add=T)


          enter image description here



          Thus, for $n = 10$ the constants $c_1 = 4.975$ and
          $c_2 = 17.084$ for
          a 95% confidence interval are quantiles 0.025 and 0.975, respectively, of $Q sim mathsfGamma(10, 1).$



          qgamma(c(.025, .975), 10, 1)
          [1] 4.795389 17.084803


          In particular, for the exponential sample shown below (second row),
          a 95% confidence interval is $(2.224, 7.922).$ Notice the reversal of the quantiles in 'pivoting' $Q,$ which
          has $mu$ in the denominator.



          set.seed(1234); x = sort(round(rexp(10, 1/5), 2)); x
          [1] 0.03 0.45 1.01 1.23 1.94 3.80 4.12 4.19 8.71 12.51
          t = sum(x); t
          [1] 37.99
          t/qgamma(c(.975, .025), 10, 1)
          [1] 2.223614 7.922194


          Note: Because the chi-squared distribution is a member of the gamma family, it is possible to find endpoints for such a confidence interval in terms of a chi-squared distribution.



          See Wikipedia on exponential distributions under 'confidence intervals'. (That discussion uses rate parameter $lambda$ for the exponential distribution, instead of $mu.)$






          share|cite|improve this answer











          $endgroup$



          You don't say how the exponential distribution is
          parameterized. Two parameterizations are in common use--mean and rate.



          Let $E(X_i) = mu.$ Then one
          can show that $$frac 1 mu sum_i=1^n X_i sim
          mathsfGamma(textshape = n, textrate=scale = 1).$$



          In R statistical software the exponential distribution is parameterized according rate $lambda = 1/mu.$ Let $n = 10$ and $lambda = 1/5,$ so that $mu = 5.$ The following program simulates $m = 10^6$ samples of size $n = 10$ from $mathsfExp(textrate = lambda = 1/5),$ finds $$Q = frac 1 mu sum_i=1^n X_i =
          lambda sum_i=1^n X_i$$
          for each sample, and plots the histogram of the one million $Q$'s, The figure
          illustrates that $Q sim mathsfGamma(10, 1).$
          (Use MGFs for a formal proof.)



          set.seed(414) # for reproducibility
          q = replicate(10^5, sum(rexp(10, 1/5))/5)
          lbl = "Simulated Dist'n of Q with Density of GAMMA(10, 1)"
          hist(q, prob=T, br=30, col="skyblue2", main=lbl)
          curve(dgamma(x,10,1), col="red", add=T)


          enter image description here



          Thus, for $n = 10$ the constants $c_1 = 4.975$ and
          $c_2 = 17.084$ for
          a 95% confidence interval are quantiles 0.025 and 0.975, respectively, of $Q sim mathsfGamma(10, 1).$



          qgamma(c(.025, .975), 10, 1)
          [1] 4.795389 17.084803


          In particular, for the exponential sample shown below (second row),
          a 95% confidence interval is $(2.224, 7.922).$ Notice the reversal of the quantiles in 'pivoting' $Q,$ which
          has $mu$ in the denominator.



          set.seed(1234); x = sort(round(rexp(10, 1/5), 2)); x
          [1] 0.03 0.45 1.01 1.23 1.94 3.80 4.12 4.19 8.71 12.51
          t = sum(x); t
          [1] 37.99
          t/qgamma(c(.975, .025), 10, 1)
          [1] 2.223614 7.922194


          Note: Because the chi-squared distribution is a member of the gamma family, it is possible to find endpoints for such a confidence interval in terms of a chi-squared distribution.



          See Wikipedia on exponential distributions under 'confidence intervals'. (That discussion uses rate parameter $lambda$ for the exponential distribution, instead of $mu.)$







          share|cite|improve this answer














          share|cite|improve this answer



          share|cite|improve this answer








          edited 4 hours ago

























          answered 4 hours ago









          BruceETBruceET

          6,6331721




          6,6331721























              1












              $begingroup$

              Taking $theta$ as the scale parameter, it can be shown that:



              $$fracn barXtheta sim textGa(n,1).$$



              To form a confidence interval we choose any critical points $c_1 < c_2$ from the $textGa(n,1)$ distribution such that these points contain probability $1-alpha$ of the distribution. Using the above pivotal quantity we then have:



              $$mathbbP Bigg( c_1 leqslant fracn barXtheta leqslant c_2 Bigg) = 1-alpha
              quad quad quad quad quad
              int limits_c_1^c_2 textGa(r|n,1) dr = 1 - alpha.$$



              Re-arranging the inequality in this probability statement and substituting the observed sample mean gives the confidence interval:



              $$textCI_theta(1-alpha) = Bigg[ fracn barxc_2 , fracn barxc_1 Bigg].$$



              This confidence interval is valid for any choice of $c_1<c_2$ so long as it obeys the required integral condition. For simplicity, many analysts use the symmetric critical points. However, it is possible to optimise the confidence interval by minimising its length, which we show below.




              Optimising the confidence interval: The length of this confidence interval is proportional to $1/c_1-1/c_2$, and so we minimise the length of the interval by choosing the critical points to minimise this distance. This can be done using the nlm function in R. In the following code we give a function for the minimum-length confidence interval for this problem, which we apply to some simulated data.



              #Set the objective function for minimisation
              OBJECTIVE <- function(c1, n, alpha)
              pp <- pgamma(c1, n, 1, lower.tail = TRUE);
              c2 <- qgamma(1 - alpha + pp, n, 1, lower.tail = TRUE);
              1/c1 - 1/c2;

              #Find the minimum-length confidence interval
              CONF_INT <- function(n, alpha, xbar)
              START_c1 <- qgamma(alpha/2, n, 1, lower.tail = TRUE);
              MINIMISE <- nlm(f = OBJECTIVE, p = START_c1, n = n, alpha = alpha);
              c1 <- MINIMISE$estimate;
              pp <- pgamma(c1, n, 1, lower.tail = TRUE);
              c2 <- qgamma(1 - alpha + pp, n, 1, lower.tail = TRUE);
              c(n*xbar/c2, n*xbar/c1);

              #Generate simulation data
              set.seed(921730198);
              n <- 300;
              scale <- 25.4;
              DATA <- rexp(n, rate = 1/scale);

              #Application of confidence interval to simulated data
              n <- length(DATA);
              xbar <- mean(DATA);
              alpha <- 0.05;

              CONF_INT(n, alpha, xbar);

              [1] 23.32040 29.24858





              share|cite|improve this answer











              $endgroup$

















                1












                $begingroup$

                Taking $theta$ as the scale parameter, it can be shown that:



                $$fracn barXtheta sim textGa(n,1).$$



                To form a confidence interval we choose any critical points $c_1 < c_2$ from the $textGa(n,1)$ distribution such that these points contain probability $1-alpha$ of the distribution. Using the above pivotal quantity we then have:



                $$mathbbP Bigg( c_1 leqslant fracn barXtheta leqslant c_2 Bigg) = 1-alpha
                quad quad quad quad quad
                int limits_c_1^c_2 textGa(r|n,1) dr = 1 - alpha.$$



                Re-arranging the inequality in this probability statement and substituting the observed sample mean gives the confidence interval:



                $$textCI_theta(1-alpha) = Bigg[ fracn barxc_2 , fracn barxc_1 Bigg].$$



                This confidence interval is valid for any choice of $c_1<c_2$ so long as it obeys the required integral condition. For simplicity, many analysts use the symmetric critical points. However, it is possible to optimise the confidence interval by minimising its length, which we show below.




                Optimising the confidence interval: The length of this confidence interval is proportional to $1/c_1-1/c_2$, and so we minimise the length of the interval by choosing the critical points to minimise this distance. This can be done using the nlm function in R. In the following code we give a function for the minimum-length confidence interval for this problem, which we apply to some simulated data.



                #Set the objective function for minimisation
                OBJECTIVE <- function(c1, n, alpha)
                pp <- pgamma(c1, n, 1, lower.tail = TRUE);
                c2 <- qgamma(1 - alpha + pp, n, 1, lower.tail = TRUE);
                1/c1 - 1/c2;

                #Find the minimum-length confidence interval
                CONF_INT <- function(n, alpha, xbar)
                START_c1 <- qgamma(alpha/2, n, 1, lower.tail = TRUE);
                MINIMISE <- nlm(f = OBJECTIVE, p = START_c1, n = n, alpha = alpha);
                c1 <- MINIMISE$estimate;
                pp <- pgamma(c1, n, 1, lower.tail = TRUE);
                c2 <- qgamma(1 - alpha + pp, n, 1, lower.tail = TRUE);
                c(n*xbar/c2, n*xbar/c1);

                #Generate simulation data
                set.seed(921730198);
                n <- 300;
                scale <- 25.4;
                DATA <- rexp(n, rate = 1/scale);

                #Application of confidence interval to simulated data
                n <- length(DATA);
                xbar <- mean(DATA);
                alpha <- 0.05;

                CONF_INT(n, alpha, xbar);

                [1] 23.32040 29.24858





                share|cite|improve this answer











                $endgroup$















                  1












                  1








                  1





                  $begingroup$

                  Taking $theta$ as the scale parameter, it can be shown that:



                  $$fracn barXtheta sim textGa(n,1).$$



                  To form a confidence interval we choose any critical points $c_1 < c_2$ from the $textGa(n,1)$ distribution such that these points contain probability $1-alpha$ of the distribution. Using the above pivotal quantity we then have:



                  $$mathbbP Bigg( c_1 leqslant fracn barXtheta leqslant c_2 Bigg) = 1-alpha
                  quad quad quad quad quad
                  int limits_c_1^c_2 textGa(r|n,1) dr = 1 - alpha.$$



                  Re-arranging the inequality in this probability statement and substituting the observed sample mean gives the confidence interval:



                  $$textCI_theta(1-alpha) = Bigg[ fracn barxc_2 , fracn barxc_1 Bigg].$$



                  This confidence interval is valid for any choice of $c_1<c_2$ so long as it obeys the required integral condition. For simplicity, many analysts use the symmetric critical points. However, it is possible to optimise the confidence interval by minimising its length, which we show below.




                  Optimising the confidence interval: The length of this confidence interval is proportional to $1/c_1-1/c_2$, and so we minimise the length of the interval by choosing the critical points to minimise this distance. This can be done using the nlm function in R. In the following code we give a function for the minimum-length confidence interval for this problem, which we apply to some simulated data.



                  #Set the objective function for minimisation
                  OBJECTIVE <- function(c1, n, alpha)
                  pp <- pgamma(c1, n, 1, lower.tail = TRUE);
                  c2 <- qgamma(1 - alpha + pp, n, 1, lower.tail = TRUE);
                  1/c1 - 1/c2;

                  #Find the minimum-length confidence interval
                  CONF_INT <- function(n, alpha, xbar)
                  START_c1 <- qgamma(alpha/2, n, 1, lower.tail = TRUE);
                  MINIMISE <- nlm(f = OBJECTIVE, p = START_c1, n = n, alpha = alpha);
                  c1 <- MINIMISE$estimate;
                  pp <- pgamma(c1, n, 1, lower.tail = TRUE);
                  c2 <- qgamma(1 - alpha + pp, n, 1, lower.tail = TRUE);
                  c(n*xbar/c2, n*xbar/c1);

                  #Generate simulation data
                  set.seed(921730198);
                  n <- 300;
                  scale <- 25.4;
                  DATA <- rexp(n, rate = 1/scale);

                  #Application of confidence interval to simulated data
                  n <- length(DATA);
                  xbar <- mean(DATA);
                  alpha <- 0.05;

                  CONF_INT(n, alpha, xbar);

                  [1] 23.32040 29.24858





                  share|cite|improve this answer











                  $endgroup$



                  Taking $theta$ as the scale parameter, it can be shown that:



                  $$fracn barXtheta sim textGa(n,1).$$



                  To form a confidence interval we choose any critical points $c_1 < c_2$ from the $textGa(n,1)$ distribution such that these points contain probability $1-alpha$ of the distribution. Using the above pivotal quantity we then have:



                  $$mathbbP Bigg( c_1 leqslant fracn barXtheta leqslant c_2 Bigg) = 1-alpha
                  quad quad quad quad quad
                  int limits_c_1^c_2 textGa(r|n,1) dr = 1 - alpha.$$



                  Re-arranging the inequality in this probability statement and substituting the observed sample mean gives the confidence interval:



                  $$textCI_theta(1-alpha) = Bigg[ fracn barxc_2 , fracn barxc_1 Bigg].$$



                  This confidence interval is valid for any choice of $c_1<c_2$ so long as it obeys the required integral condition. For simplicity, many analysts use the symmetric critical points. However, it is possible to optimise the confidence interval by minimising its length, which we show below.




                  Optimising the confidence interval: The length of this confidence interval is proportional to $1/c_1-1/c_2$, and so we minimise the length of the interval by choosing the critical points to minimise this distance. This can be done using the nlm function in R. In the following code we give a function for the minimum-length confidence interval for this problem, which we apply to some simulated data.



                  #Set the objective function for minimisation
                  OBJECTIVE <- function(c1, n, alpha)
                  pp <- pgamma(c1, n, 1, lower.tail = TRUE);
                  c2 <- qgamma(1 - alpha + pp, n, 1, lower.tail = TRUE);
                  1/c1 - 1/c2;

                  #Find the minimum-length confidence interval
                  CONF_INT <- function(n, alpha, xbar)
                  START_c1 <- qgamma(alpha/2, n, 1, lower.tail = TRUE);
                  MINIMISE <- nlm(f = OBJECTIVE, p = START_c1, n = n, alpha = alpha);
                  c1 <- MINIMISE$estimate;
                  pp <- pgamma(c1, n, 1, lower.tail = TRUE);
                  c2 <- qgamma(1 - alpha + pp, n, 1, lower.tail = TRUE);
                  c(n*xbar/c2, n*xbar/c1);

                  #Generate simulation data
                  set.seed(921730198);
                  n <- 300;
                  scale <- 25.4;
                  DATA <- rexp(n, rate = 1/scale);

                  #Application of confidence interval to simulated data
                  n <- length(DATA);
                  xbar <- mean(DATA);
                  alpha <- 0.05;

                  CONF_INT(n, alpha, xbar);

                  [1] 23.32040 29.24858






                  share|cite|improve this answer














                  share|cite|improve this answer



                  share|cite|improve this answer








                  edited 1 hour ago

























                  answered 2 hours ago









                  BenBen

                  28.3k233128




                  28.3k233128



























                      draft saved

                      draft discarded
















































                      Thanks for contributing an answer to Cross Validated!


                      • Please be sure to answer the question. Provide details and share your research!

                      But avoid


                      • Asking for help, clarification, or responding to other answers.

                      • Making statements based on opinion; back them up with references or personal experience.

                      Use MathJax to format equations. MathJax reference.


                      To learn more, see our tips on writing great answers.




                      draft saved


                      draft discarded














                      StackExchange.ready(
                      function ()
                      StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fstats.stackexchange.com%2fquestions%2f403059%2fhow-do-we-build-a-confidence-interval-for-the-parameter-of-the-exponential-distr%23new-answer', 'question_page');

                      );

                      Post as a guest















                      Required, but never shown





















































                      Required, but never shown














                      Required, but never shown












                      Required, but never shown







                      Required, but never shown

































                      Required, but never shown














                      Required, but never shown












                      Required, but never shown







                      Required, but never shown







                      Popular posts from this blog

                      Are there any AGPL-style licences that require source code modifications to be public? Planned maintenance scheduled April 23, 2019 at 23:30 UTC (7:30pm US/Eastern) Announcing the arrival of Valued Associate #679: Cesar Manara Unicorn Meta Zoo #1: Why another podcast?Force derivative works to be publicAre there any GPL like licenses for Apple App Store?Do you violate the GPL if you provide source code that cannot be compiled?GPL - is it distribution to use libraries in an appliance loaned to customers?Distributing App for free which uses GPL'ed codeModifications of server software under GPL, with web/CLI interfaceDoes using an AGPLv3-licensed library prevent me from dual-licensing my own source code?Can I publish only select code under GPLv3 from a private project?Is there published precedent regarding the scope of covered work that uses AGPL software?If MIT licensed code links to GPL licensed code what should be the license of the resulting binary program?If I use a public API endpoint that has its source code licensed under AGPL in my app, do I need to disclose my source?

                      2013 GY136 Descoberta | Órbita | Referências Menu de navegação«List Of Centaurs and Scattered-Disk Objects»«List of Known Trans-Neptunian Objects»

                      Mortes em março de 2019 Referências Menu de navegação«Zhores Alferov, Nobel de Física bielorrusso, morre aos 88 anos - Ciência»«Fallece Rafael Torija, o bispo emérito de Ciudad Real»«Peter Hurford dies at 88»«Keith Flint, vocalista do The Prodigy, morre aos 49 anos»«Luke Perry, ator de 'Barrados no baile' e 'Riverdale', morre aos 52 anos»«Former Rangers and Scotland captain Eric Caldow dies, aged 84»«Morreu, aos 61 anos, a antiga lenda do wrestling King Kong Bundy»«Fallece el actor y director teatral Abraham Stavans»«In Memoriam Guillaume Faye»«Sidney Sheinberg, a Force Behind Universal and Spielberg, Is Dead at 84»«Carmine Persico, Colombo Crime Family Boss, Is Dead at 85»«Dirigent Michael Gielen gestorben»«Ciclista tricampeã mundial e prata na Rio 2016 é encontrada morta em casa aos 23 anos»«Pagan Community Notes: Raven Grimassi dies, Indianapolis pop-up event cancelled, Circle Sanctuary announces new podcast, and more!»«Hal Blaine, Wrecking Crew Drummer, Dies at 90»«Morre Coutinho, que editou dupla lendária com Pelé no Santos»«Cantor Demétrius, ídolo da Jovem Guarda, morre em SP»«Ex-presidente do Vasco, Eurico Miranda morre no Rio de Janeiro»«Bronze no Mundial de basquete de 1971, Laís Elena morre aos 76 anos»«Diretor de Corridas da F1, Charlie Whiting morre aos 66 anos às vésperas do GP da Austrália»«Morreu o cardeal Danneels, da Bélgica»«Morreu o cartoonista Augusto Cid»«Morreu a atriz Maria Isabel de Lizandra, de "Vale Tudo" e novelas da Tupi»«WS Merwin, prize-winning poet of nature, dies at 91»«Atriz Márcia Real morre em São Paulo aos 88 anos»«Mauritanie: décès de l'ancien président Mohamed Mahmoud ould Louly»«Morreu Dick Dale, o rei da surf guitar e de "Pulp Fiction"»«Falleció Víctor Genes»«João Carlos Marinho, autor de 'O Gênio do Crime', morre em SP»«Legendary Horror Director and SFX Artist John Carl Buechler Dies at 66»«Morre em Salvador a religiosa Makota Valdina»«مرگ بازیکن‌ سابق نساجی بر اثر سقوط سنگ در مازندران»«Domingos Oliveira morre no Rio»«Morre Airton Ravagniani, ex-São Paulo, Fla, Vasco, Grêmio e Sport - Notícias»«Morre o escritor Flavio Moreira da Costa»«Larry Cohen, Writer-Director of 'It's Alive' and 'Hell Up in Harlem,' Dies at 77»«Scott Walker, experimental singer-songwriter, dead at 76»«Joseph Pilato, Day of the Dead Star and Horror Favorite, Dies at 70»«Sheffield United set to pay tribute to legendary goalkeeper Ted Burgin who has died at 91»«Morre Rafael Henzel, sobrevivente de acidente aéreo da Chapecoense»«Morre Valery Bykovsky, um dos primeiros cosmonautas da União Soviética»«Agnès Varda, cineasta da Nouvelle Vague, morre aos 90 anos»«Agnès Varda, cineasta francesa, morre aos 90 anos»«Tania Mallet, James Bond Actress and Helen Mirren's Cousin, Dies at 77»e