$alpha$-Quantile for Gamma distribution $Gamma(n,frac{1}{2})$












1












$begingroup$


Let $XsimGamma(n,frac{1}{2})$ and $alphain(0,1)$. I know that



begin{align}
P(Xle s)&=int_0^ sfrac{1}{2^nGamma(n)}x^{n-1}e^{-frac{x}{2}}dx\
&=frac{1}{2(n-1)!}int_0^sBig(frac{x}{2}Big)^{n-1}e^{-frac{x}{2}}dx\
&=frac{1}{(n-1)!}int_0^{s/2} t^{n-1} e^{-t}dt\
end{align}



Now I wonder how can I find $u,v$ such that



$$frac{alpha}{2}ge frac{1}{(n-1)!}int_0^{u/2} t^{n-1} e^{-t}dt=P(Xle u),$$
$$1-frac{alpha}{2} le frac{1}{(n-1)!}int_0^{v/2} t^{n-1} e^{-t}dt= P(Xle v)text{ ?}$$
Maybe there is an easy argument... Thanks in advance for any help!










share|cite|improve this question











$endgroup$












  • $begingroup$
    Of course, the integral in the end is the same as in the definition of the Gammafunction.
    $endgroup$
    – user408858
    Dec 18 '18 at 0:12










  • $begingroup$
    If I understand your notation and if you need the answers instead of hints how to do the calculus, then here is R code that will give you the answer for $n = 1, 2, dots, 50$ and $alpha = 0.05:$ n = 1:50; L = qgamma(.025, n, 2); U = qgamma(.975, n ,2); cbind(n, L, U). For $n = 1,$ you have the exponential distribution, for which the calculus is trivial. For $n > 200,$ you should get good answers with normal approximation.
    $endgroup$
    – BruceET
    Dec 18 '18 at 0:28












  • $begingroup$
    Oh! Do I get you right, that there are these tables in backmatter of statistic books for exactly these kind of questions? :-)
    $endgroup$
    – user408858
    Dec 18 '18 at 0:30












  • $begingroup$
    You might be able to use widely available printed tables of chi-squared distributions (with appriat modifications; halve DF to get gamma shape). Some printed tables have only the upper tail. Look around. If you're allergic to R, I suppose MatLab and such will do much the same computations. // But R is free, devotedly bug free for basic stuff like this, and trivial to install. See Wikipedia. // Code qchisq(.95, 4) and qgamma(.95, 2, 1/2) both return 9.487729
    $endgroup$
    – BruceET
    Dec 18 '18 at 1:08












  • $begingroup$
    Actually one can just take the $1-frac{alpha}{2}$-quantile by definition(?) or just calculate it like @BruceET did. :-)
    $endgroup$
    – user408858
    Dec 20 '18 at 0:46
















1












$begingroup$


Let $XsimGamma(n,frac{1}{2})$ and $alphain(0,1)$. I know that



begin{align}
P(Xle s)&=int_0^ sfrac{1}{2^nGamma(n)}x^{n-1}e^{-frac{x}{2}}dx\
&=frac{1}{2(n-1)!}int_0^sBig(frac{x}{2}Big)^{n-1}e^{-frac{x}{2}}dx\
&=frac{1}{(n-1)!}int_0^{s/2} t^{n-1} e^{-t}dt\
end{align}



Now I wonder how can I find $u,v$ such that



$$frac{alpha}{2}ge frac{1}{(n-1)!}int_0^{u/2} t^{n-1} e^{-t}dt=P(Xle u),$$
$$1-frac{alpha}{2} le frac{1}{(n-1)!}int_0^{v/2} t^{n-1} e^{-t}dt= P(Xle v)text{ ?}$$
Maybe there is an easy argument... Thanks in advance for any help!










share|cite|improve this question











$endgroup$












  • $begingroup$
    Of course, the integral in the end is the same as in the definition of the Gammafunction.
    $endgroup$
    – user408858
    Dec 18 '18 at 0:12










  • $begingroup$
    If I understand your notation and if you need the answers instead of hints how to do the calculus, then here is R code that will give you the answer for $n = 1, 2, dots, 50$ and $alpha = 0.05:$ n = 1:50; L = qgamma(.025, n, 2); U = qgamma(.975, n ,2); cbind(n, L, U). For $n = 1,$ you have the exponential distribution, for which the calculus is trivial. For $n > 200,$ you should get good answers with normal approximation.
    $endgroup$
    – BruceET
    Dec 18 '18 at 0:28












  • $begingroup$
    Oh! Do I get you right, that there are these tables in backmatter of statistic books for exactly these kind of questions? :-)
    $endgroup$
    – user408858
    Dec 18 '18 at 0:30












  • $begingroup$
    You might be able to use widely available printed tables of chi-squared distributions (with appriat modifications; halve DF to get gamma shape). Some printed tables have only the upper tail. Look around. If you're allergic to R, I suppose MatLab and such will do much the same computations. // But R is free, devotedly bug free for basic stuff like this, and trivial to install. See Wikipedia. // Code qchisq(.95, 4) and qgamma(.95, 2, 1/2) both return 9.487729
    $endgroup$
    – BruceET
    Dec 18 '18 at 1:08












  • $begingroup$
    Actually one can just take the $1-frac{alpha}{2}$-quantile by definition(?) or just calculate it like @BruceET did. :-)
    $endgroup$
    – user408858
    Dec 20 '18 at 0:46














1












1








1





$begingroup$


Let $XsimGamma(n,frac{1}{2})$ and $alphain(0,1)$. I know that



begin{align}
P(Xle s)&=int_0^ sfrac{1}{2^nGamma(n)}x^{n-1}e^{-frac{x}{2}}dx\
&=frac{1}{2(n-1)!}int_0^sBig(frac{x}{2}Big)^{n-1}e^{-frac{x}{2}}dx\
&=frac{1}{(n-1)!}int_0^{s/2} t^{n-1} e^{-t}dt\
end{align}



Now I wonder how can I find $u,v$ such that



$$frac{alpha}{2}ge frac{1}{(n-1)!}int_0^{u/2} t^{n-1} e^{-t}dt=P(Xle u),$$
$$1-frac{alpha}{2} le frac{1}{(n-1)!}int_0^{v/2} t^{n-1} e^{-t}dt= P(Xle v)text{ ?}$$
Maybe there is an easy argument... Thanks in advance for any help!










share|cite|improve this question











$endgroup$




Let $XsimGamma(n,frac{1}{2})$ and $alphain(0,1)$. I know that



begin{align}
P(Xle s)&=int_0^ sfrac{1}{2^nGamma(n)}x^{n-1}e^{-frac{x}{2}}dx\
&=frac{1}{2(n-1)!}int_0^sBig(frac{x}{2}Big)^{n-1}e^{-frac{x}{2}}dx\
&=frac{1}{(n-1)!}int_0^{s/2} t^{n-1} e^{-t}dt\
end{align}



Now I wonder how can I find $u,v$ such that



$$frac{alpha}{2}ge frac{1}{(n-1)!}int_0^{u/2} t^{n-1} e^{-t}dt=P(Xle u),$$
$$1-frac{alpha}{2} le frac{1}{(n-1)!}int_0^{v/2} t^{n-1} e^{-t}dt= P(Xle v)text{ ?}$$
Maybe there is an easy argument... Thanks in advance for any help!







real-analysis probability analysis probability-theory statistics






share|cite|improve this question















share|cite|improve this question













share|cite|improve this question




share|cite|improve this question








edited Dec 18 '18 at 0:21







user408858

















asked Dec 18 '18 at 0:11









user408858user408858

482213




482213












  • $begingroup$
    Of course, the integral in the end is the same as in the definition of the Gammafunction.
    $endgroup$
    – user408858
    Dec 18 '18 at 0:12










  • $begingroup$
    If I understand your notation and if you need the answers instead of hints how to do the calculus, then here is R code that will give you the answer for $n = 1, 2, dots, 50$ and $alpha = 0.05:$ n = 1:50; L = qgamma(.025, n, 2); U = qgamma(.975, n ,2); cbind(n, L, U). For $n = 1,$ you have the exponential distribution, for which the calculus is trivial. For $n > 200,$ you should get good answers with normal approximation.
    $endgroup$
    – BruceET
    Dec 18 '18 at 0:28












  • $begingroup$
    Oh! Do I get you right, that there are these tables in backmatter of statistic books for exactly these kind of questions? :-)
    $endgroup$
    – user408858
    Dec 18 '18 at 0:30












  • $begingroup$
    You might be able to use widely available printed tables of chi-squared distributions (with appriat modifications; halve DF to get gamma shape). Some printed tables have only the upper tail. Look around. If you're allergic to R, I suppose MatLab and such will do much the same computations. // But R is free, devotedly bug free for basic stuff like this, and trivial to install. See Wikipedia. // Code qchisq(.95, 4) and qgamma(.95, 2, 1/2) both return 9.487729
    $endgroup$
    – BruceET
    Dec 18 '18 at 1:08












  • $begingroup$
    Actually one can just take the $1-frac{alpha}{2}$-quantile by definition(?) or just calculate it like @BruceET did. :-)
    $endgroup$
    – user408858
    Dec 20 '18 at 0:46


















  • $begingroup$
    Of course, the integral in the end is the same as in the definition of the Gammafunction.
    $endgroup$
    – user408858
    Dec 18 '18 at 0:12










  • $begingroup$
    If I understand your notation and if you need the answers instead of hints how to do the calculus, then here is R code that will give you the answer for $n = 1, 2, dots, 50$ and $alpha = 0.05:$ n = 1:50; L = qgamma(.025, n, 2); U = qgamma(.975, n ,2); cbind(n, L, U). For $n = 1,$ you have the exponential distribution, for which the calculus is trivial. For $n > 200,$ you should get good answers with normal approximation.
    $endgroup$
    – BruceET
    Dec 18 '18 at 0:28












  • $begingroup$
    Oh! Do I get you right, that there are these tables in backmatter of statistic books for exactly these kind of questions? :-)
    $endgroup$
    – user408858
    Dec 18 '18 at 0:30












  • $begingroup$
    You might be able to use widely available printed tables of chi-squared distributions (with appriat modifications; halve DF to get gamma shape). Some printed tables have only the upper tail. Look around. If you're allergic to R, I suppose MatLab and such will do much the same computations. // But R is free, devotedly bug free for basic stuff like this, and trivial to install. See Wikipedia. // Code qchisq(.95, 4) and qgamma(.95, 2, 1/2) both return 9.487729
    $endgroup$
    – BruceET
    Dec 18 '18 at 1:08












  • $begingroup$
    Actually one can just take the $1-frac{alpha}{2}$-quantile by definition(?) or just calculate it like @BruceET did. :-)
    $endgroup$
    – user408858
    Dec 20 '18 at 0:46
















$begingroup$
Of course, the integral in the end is the same as in the definition of the Gammafunction.
$endgroup$
– user408858
Dec 18 '18 at 0:12




$begingroup$
Of course, the integral in the end is the same as in the definition of the Gammafunction.
$endgroup$
– user408858
Dec 18 '18 at 0:12












$begingroup$
If I understand your notation and if you need the answers instead of hints how to do the calculus, then here is R code that will give you the answer for $n = 1, 2, dots, 50$ and $alpha = 0.05:$ n = 1:50; L = qgamma(.025, n, 2); U = qgamma(.975, n ,2); cbind(n, L, U). For $n = 1,$ you have the exponential distribution, for which the calculus is trivial. For $n > 200,$ you should get good answers with normal approximation.
$endgroup$
– BruceET
Dec 18 '18 at 0:28






$begingroup$
If I understand your notation and if you need the answers instead of hints how to do the calculus, then here is R code that will give you the answer for $n = 1, 2, dots, 50$ and $alpha = 0.05:$ n = 1:50; L = qgamma(.025, n, 2); U = qgamma(.975, n ,2); cbind(n, L, U). For $n = 1,$ you have the exponential distribution, for which the calculus is trivial. For $n > 200,$ you should get good answers with normal approximation.
$endgroup$
– BruceET
Dec 18 '18 at 0:28














$begingroup$
Oh! Do I get you right, that there are these tables in backmatter of statistic books for exactly these kind of questions? :-)
$endgroup$
– user408858
Dec 18 '18 at 0:30






$begingroup$
Oh! Do I get you right, that there are these tables in backmatter of statistic books for exactly these kind of questions? :-)
$endgroup$
– user408858
Dec 18 '18 at 0:30














$begingroup$
You might be able to use widely available printed tables of chi-squared distributions (with appriat modifications; halve DF to get gamma shape). Some printed tables have only the upper tail. Look around. If you're allergic to R, I suppose MatLab and such will do much the same computations. // But R is free, devotedly bug free for basic stuff like this, and trivial to install. See Wikipedia. // Code qchisq(.95, 4) and qgamma(.95, 2, 1/2) both return 9.487729
$endgroup$
– BruceET
Dec 18 '18 at 1:08






$begingroup$
You might be able to use widely available printed tables of chi-squared distributions (with appriat modifications; halve DF to get gamma shape). Some printed tables have only the upper tail. Look around. If you're allergic to R, I suppose MatLab and such will do much the same computations. // But R is free, devotedly bug free for basic stuff like this, and trivial to install. See Wikipedia. // Code qchisq(.95, 4) and qgamma(.95, 2, 1/2) both return 9.487729
$endgroup$
– BruceET
Dec 18 '18 at 1:08














$begingroup$
Actually one can just take the $1-frac{alpha}{2}$-quantile by definition(?) or just calculate it like @BruceET did. :-)
$endgroup$
– user408858
Dec 20 '18 at 0:46




$begingroup$
Actually one can just take the $1-frac{alpha}{2}$-quantile by definition(?) or just calculate it like @BruceET did. :-)
$endgroup$
– user408858
Dec 20 '18 at 0:46










0






active

oldest

votes











Your Answer





StackExchange.ifUsing("editor", function () {
return StackExchange.using("mathjaxEditing", function () {
StackExchange.MarkdownEditor.creationCallbacks.add(function (editor, postfix) {
StackExchange.mathjaxEditing.prepareWmdForMathJax(editor, postfix, [["$", "$"], ["\\(","\\)"]]);
});
});
}, "mathjax-editing");

StackExchange.ready(function() {
var channelOptions = {
tags: "".split(" "),
id: "69"
};
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: true,
noModals: true,
showLowRepImageUploadWarning: true,
reputationToPostImages: 10,
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
},
noCode: true, onDemand: true,
discardSelector: ".discard-answer"
,immediatelyShowMarkdownHelp:true
});


}
});














draft saved

draft discarded


















StackExchange.ready(
function () {
StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fmath.stackexchange.com%2fquestions%2f3044625%2falpha-quantile-for-gamma-distribution-gamman-frac12%23new-answer', 'question_page');
}
);

Post as a guest















Required, but never shown

























0






active

oldest

votes








0






active

oldest

votes









active

oldest

votes






active

oldest

votes
















draft saved

draft discarded




















































Thanks for contributing an answer to Mathematics Stack Exchange!


  • 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%2fmath.stackexchange.com%2fquestions%2f3044625%2falpha-quantile-for-gamma-distribution-gamman-frac12%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

Plaza Victoria

Puebla de Zaragoza

Musa