$alpha$-Quantile for Gamma distribution $Gamma(n,frac{1}{2})$
$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!
real-analysis probability analysis probability-theory statistics
$endgroup$
add a comment |
$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!
real-analysis probability analysis probability-theory statistics
$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. // Codeqchisq(.95, 4)
andqgamma(.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
add a comment |
$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!
real-analysis probability analysis probability-theory statistics
$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
real-analysis probability analysis probability-theory statistics
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. // Codeqchisq(.95, 4)
andqgamma(.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
add a comment |
$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. // Codeqchisq(.95, 4)
andqgamma(.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
add a comment |
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
});
}
});
Sign up or log in
StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
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
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.
Sign up or log in
StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
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
Sign up or log in
StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
Sign up or log in
StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
Sign up or log in
StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
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
$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)
andqgamma(.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