What is the PDF of $X_1X_2-X_3X_4$ conditionally on $sumlimits_{k=1}^4X_k=1$, for $(X_k)$ i.i.d. uniform on...












1












$begingroup$


Consider the following expression



$$D=X-(X+Y)(X+Z) text{under the constraint} X+Y+Z leq 1 tag{1}$$



(being understood that $X, Y$ and $Z$ are $geq 0$).



This expression can also be written under the form :



$$D=XT-YZ=begin{vmatrix}X&Y\Z&Tend{vmatrix} text{given that} X+Y+Z+T=1.tag{2}$$



where $(X,Y,Z,T)$ take their values in $[0,1]^4$ ($[0,1]$ being the set of real values between $0$ and $1$).



(all these constraints are due to the fact that we deal with a complete system of events ; see "Origin of the question" below with definitions (4)).



My question : what is the pdf (probability density function) of random variable D ?


My work :



I have attempted different computations, but they have failed to be conclusive ; I hit on the difficulty to assign a distribution to the random vector $(X,Y,Z,T)$ taking into account constraint $X+Y+Z+T=1$. I have in particular attempted to replace uniformly distributed $X,Y,Z,T$ on $mathbb{R}$ by random variables taking discrete values ${0/n,1/n,cdots (n-1)/n,1}$ without success.



I have made simulations displaying a rather "smooth" symmetrical behavior (see Fig. below displaying an histogram for a large number of simulations) with a relative spike around value $d=0$.



enter image description here



The origin of this question : A recent question about extreme values of what could be called a "degree of dependence':



$$D=P(A cap B) - P(A)P(B) tag{3}$$



(that can be transformed into (1)) and the answer I gave : https://math.stackexchange.com/q/3090278
has triggered (as often is the case) a different question : instead of restricting our attention to the extreme values of (1) (or the equivalent expression (2)), I asked myself the more general question of the pdf of the values taken by (1) with the following definitions :



$$X=P(A cap B), Y=P(A cap B^c), Z=P(A^c cap B), T=P(A^c cap B^c).tag{4}$$



as one can see globally on the following diagram :



enter image description here



Besides, if one drops constraint $X+Y+Z+T=1$ (yielding a very different issue), what would become the pdf of $D$ ?(I give below a histogram gathering results of simulation with this hypothesis).



enter image description here










share|cite|improve this question











$endgroup$








  • 1




    $begingroup$
    First, note that you can eliminate the condition by doing a few operations on the determinant to transform it into: $D=X-(X+Y)(X+Z)$. (Just add the first column to the last column and then the first row to the last row.) So that's already one less constraint.
    $endgroup$
    – Raskolnikov
    Jan 28 at 16:47






  • 1




    $begingroup$
    Then the condition $Dleq d$ determines a section of space bounded by a hyperboloid. Adding the conditions that $X,Y,Z in [0,1]$ means that you have to find the volume of that section of space within the unit cube to find the probability.
    $endgroup$
    – Raskolnikov
    Jan 28 at 16:52












  • $begingroup$
    @Raskolnikov In fact expression $D=X-(X+Y)(X+Z)=P(A cap B)-P(A)P(B)$ is the initial expression ... and you are right that, although less symmetrical, this expression avoids to use constraint (2).
    $endgroup$
    – Jean Marie
    Jan 28 at 17:24










  • $begingroup$
    @Raskolnikov ... but it will remain the constraint $X+Y+Z leq 1$.
    $endgroup$
    – Jean Marie
    Jan 28 at 17:34










  • $begingroup$
    @Raskolnikov I have slightly modified my first sentences in order to take into account your remarks.
    $endgroup$
    – Jean Marie
    Jan 28 at 18:50
















1












$begingroup$


Consider the following expression



$$D=X-(X+Y)(X+Z) text{under the constraint} X+Y+Z leq 1 tag{1}$$



(being understood that $X, Y$ and $Z$ are $geq 0$).



This expression can also be written under the form :



$$D=XT-YZ=begin{vmatrix}X&Y\Z&Tend{vmatrix} text{given that} X+Y+Z+T=1.tag{2}$$



where $(X,Y,Z,T)$ take their values in $[0,1]^4$ ($[0,1]$ being the set of real values between $0$ and $1$).



(all these constraints are due to the fact that we deal with a complete system of events ; see "Origin of the question" below with definitions (4)).



My question : what is the pdf (probability density function) of random variable D ?


My work :



I have attempted different computations, but they have failed to be conclusive ; I hit on the difficulty to assign a distribution to the random vector $(X,Y,Z,T)$ taking into account constraint $X+Y+Z+T=1$. I have in particular attempted to replace uniformly distributed $X,Y,Z,T$ on $mathbb{R}$ by random variables taking discrete values ${0/n,1/n,cdots (n-1)/n,1}$ without success.



I have made simulations displaying a rather "smooth" symmetrical behavior (see Fig. below displaying an histogram for a large number of simulations) with a relative spike around value $d=0$.



enter image description here



The origin of this question : A recent question about extreme values of what could be called a "degree of dependence':



$$D=P(A cap B) - P(A)P(B) tag{3}$$



(that can be transformed into (1)) and the answer I gave : https://math.stackexchange.com/q/3090278
has triggered (as often is the case) a different question : instead of restricting our attention to the extreme values of (1) (or the equivalent expression (2)), I asked myself the more general question of the pdf of the values taken by (1) with the following definitions :



$$X=P(A cap B), Y=P(A cap B^c), Z=P(A^c cap B), T=P(A^c cap B^c).tag{4}$$



as one can see globally on the following diagram :



enter image description here



Besides, if one drops constraint $X+Y+Z+T=1$ (yielding a very different issue), what would become the pdf of $D$ ?(I give below a histogram gathering results of simulation with this hypothesis).



enter image description here










share|cite|improve this question











$endgroup$








  • 1




    $begingroup$
    First, note that you can eliminate the condition by doing a few operations on the determinant to transform it into: $D=X-(X+Y)(X+Z)$. (Just add the first column to the last column and then the first row to the last row.) So that's already one less constraint.
    $endgroup$
    – Raskolnikov
    Jan 28 at 16:47






  • 1




    $begingroup$
    Then the condition $Dleq d$ determines a section of space bounded by a hyperboloid. Adding the conditions that $X,Y,Z in [0,1]$ means that you have to find the volume of that section of space within the unit cube to find the probability.
    $endgroup$
    – Raskolnikov
    Jan 28 at 16:52












  • $begingroup$
    @Raskolnikov In fact expression $D=X-(X+Y)(X+Z)=P(A cap B)-P(A)P(B)$ is the initial expression ... and you are right that, although less symmetrical, this expression avoids to use constraint (2).
    $endgroup$
    – Jean Marie
    Jan 28 at 17:24










  • $begingroup$
    @Raskolnikov ... but it will remain the constraint $X+Y+Z leq 1$.
    $endgroup$
    – Jean Marie
    Jan 28 at 17:34










  • $begingroup$
    @Raskolnikov I have slightly modified my first sentences in order to take into account your remarks.
    $endgroup$
    – Jean Marie
    Jan 28 at 18:50














1












1








1


1



$begingroup$


Consider the following expression



$$D=X-(X+Y)(X+Z) text{under the constraint} X+Y+Z leq 1 tag{1}$$



(being understood that $X, Y$ and $Z$ are $geq 0$).



This expression can also be written under the form :



$$D=XT-YZ=begin{vmatrix}X&Y\Z&Tend{vmatrix} text{given that} X+Y+Z+T=1.tag{2}$$



where $(X,Y,Z,T)$ take their values in $[0,1]^4$ ($[0,1]$ being the set of real values between $0$ and $1$).



(all these constraints are due to the fact that we deal with a complete system of events ; see "Origin of the question" below with definitions (4)).



My question : what is the pdf (probability density function) of random variable D ?


My work :



I have attempted different computations, but they have failed to be conclusive ; I hit on the difficulty to assign a distribution to the random vector $(X,Y,Z,T)$ taking into account constraint $X+Y+Z+T=1$. I have in particular attempted to replace uniformly distributed $X,Y,Z,T$ on $mathbb{R}$ by random variables taking discrete values ${0/n,1/n,cdots (n-1)/n,1}$ without success.



I have made simulations displaying a rather "smooth" symmetrical behavior (see Fig. below displaying an histogram for a large number of simulations) with a relative spike around value $d=0$.



enter image description here



The origin of this question : A recent question about extreme values of what could be called a "degree of dependence':



$$D=P(A cap B) - P(A)P(B) tag{3}$$



(that can be transformed into (1)) and the answer I gave : https://math.stackexchange.com/q/3090278
has triggered (as often is the case) a different question : instead of restricting our attention to the extreme values of (1) (or the equivalent expression (2)), I asked myself the more general question of the pdf of the values taken by (1) with the following definitions :



$$X=P(A cap B), Y=P(A cap B^c), Z=P(A^c cap B), T=P(A^c cap B^c).tag{4}$$



as one can see globally on the following diagram :



enter image description here



Besides, if one drops constraint $X+Y+Z+T=1$ (yielding a very different issue), what would become the pdf of $D$ ?(I give below a histogram gathering results of simulation with this hypothesis).



enter image description here










share|cite|improve this question











$endgroup$




Consider the following expression



$$D=X-(X+Y)(X+Z) text{under the constraint} X+Y+Z leq 1 tag{1}$$



(being understood that $X, Y$ and $Z$ are $geq 0$).



This expression can also be written under the form :



$$D=XT-YZ=begin{vmatrix}X&Y\Z&Tend{vmatrix} text{given that} X+Y+Z+T=1.tag{2}$$



where $(X,Y,Z,T)$ take their values in $[0,1]^4$ ($[0,1]$ being the set of real values between $0$ and $1$).



(all these constraints are due to the fact that we deal with a complete system of events ; see "Origin of the question" below with definitions (4)).



My question : what is the pdf (probability density function) of random variable D ?


My work :



I have attempted different computations, but they have failed to be conclusive ; I hit on the difficulty to assign a distribution to the random vector $(X,Y,Z,T)$ taking into account constraint $X+Y+Z+T=1$. I have in particular attempted to replace uniformly distributed $X,Y,Z,T$ on $mathbb{R}$ by random variables taking discrete values ${0/n,1/n,cdots (n-1)/n,1}$ without success.



I have made simulations displaying a rather "smooth" symmetrical behavior (see Fig. below displaying an histogram for a large number of simulations) with a relative spike around value $d=0$.



enter image description here



The origin of this question : A recent question about extreme values of what could be called a "degree of dependence':



$$D=P(A cap B) - P(A)P(B) tag{3}$$



(that can be transformed into (1)) and the answer I gave : https://math.stackexchange.com/q/3090278
has triggered (as often is the case) a different question : instead of restricting our attention to the extreme values of (1) (or the equivalent expression (2)), I asked myself the more general question of the pdf of the values taken by (1) with the following definitions :



$$X=P(A cap B), Y=P(A cap B^c), Z=P(A^c cap B), T=P(A^c cap B^c).tag{4}$$



as one can see globally on the following diagram :



enter image description here



Besides, if one drops constraint $X+Y+Z+T=1$ (yielding a very different issue), what would become the pdf of $D$ ?(I give below a histogram gathering results of simulation with this hypothesis).



enter image description here







probability probability-distributions






share|cite|improve this question















share|cite|improve this question













share|cite|improve this question




share|cite|improve this question








edited Jan 28 at 19:00









Did

249k23226466




249k23226466










asked Jan 28 at 10:41









Jean MarieJean Marie

31k42255




31k42255








  • 1




    $begingroup$
    First, note that you can eliminate the condition by doing a few operations on the determinant to transform it into: $D=X-(X+Y)(X+Z)$. (Just add the first column to the last column and then the first row to the last row.) So that's already one less constraint.
    $endgroup$
    – Raskolnikov
    Jan 28 at 16:47






  • 1




    $begingroup$
    Then the condition $Dleq d$ determines a section of space bounded by a hyperboloid. Adding the conditions that $X,Y,Z in [0,1]$ means that you have to find the volume of that section of space within the unit cube to find the probability.
    $endgroup$
    – Raskolnikov
    Jan 28 at 16:52












  • $begingroup$
    @Raskolnikov In fact expression $D=X-(X+Y)(X+Z)=P(A cap B)-P(A)P(B)$ is the initial expression ... and you are right that, although less symmetrical, this expression avoids to use constraint (2).
    $endgroup$
    – Jean Marie
    Jan 28 at 17:24










  • $begingroup$
    @Raskolnikov ... but it will remain the constraint $X+Y+Z leq 1$.
    $endgroup$
    – Jean Marie
    Jan 28 at 17:34










  • $begingroup$
    @Raskolnikov I have slightly modified my first sentences in order to take into account your remarks.
    $endgroup$
    – Jean Marie
    Jan 28 at 18:50














  • 1




    $begingroup$
    First, note that you can eliminate the condition by doing a few operations on the determinant to transform it into: $D=X-(X+Y)(X+Z)$. (Just add the first column to the last column and then the first row to the last row.) So that's already one less constraint.
    $endgroup$
    – Raskolnikov
    Jan 28 at 16:47






  • 1




    $begingroup$
    Then the condition $Dleq d$ determines a section of space bounded by a hyperboloid. Adding the conditions that $X,Y,Z in [0,1]$ means that you have to find the volume of that section of space within the unit cube to find the probability.
    $endgroup$
    – Raskolnikov
    Jan 28 at 16:52












  • $begingroup$
    @Raskolnikov In fact expression $D=X-(X+Y)(X+Z)=P(A cap B)-P(A)P(B)$ is the initial expression ... and you are right that, although less symmetrical, this expression avoids to use constraint (2).
    $endgroup$
    – Jean Marie
    Jan 28 at 17:24










  • $begingroup$
    @Raskolnikov ... but it will remain the constraint $X+Y+Z leq 1$.
    $endgroup$
    – Jean Marie
    Jan 28 at 17:34










  • $begingroup$
    @Raskolnikov I have slightly modified my first sentences in order to take into account your remarks.
    $endgroup$
    – Jean Marie
    Jan 28 at 18:50








1




1




$begingroup$
First, note that you can eliminate the condition by doing a few operations on the determinant to transform it into: $D=X-(X+Y)(X+Z)$. (Just add the first column to the last column and then the first row to the last row.) So that's already one less constraint.
$endgroup$
– Raskolnikov
Jan 28 at 16:47




$begingroup$
First, note that you can eliminate the condition by doing a few operations on the determinant to transform it into: $D=X-(X+Y)(X+Z)$. (Just add the first column to the last column and then the first row to the last row.) So that's already one less constraint.
$endgroup$
– Raskolnikov
Jan 28 at 16:47




1




1




$begingroup$
Then the condition $Dleq d$ determines a section of space bounded by a hyperboloid. Adding the conditions that $X,Y,Z in [0,1]$ means that you have to find the volume of that section of space within the unit cube to find the probability.
$endgroup$
– Raskolnikov
Jan 28 at 16:52






$begingroup$
Then the condition $Dleq d$ determines a section of space bounded by a hyperboloid. Adding the conditions that $X,Y,Z in [0,1]$ means that you have to find the volume of that section of space within the unit cube to find the probability.
$endgroup$
– Raskolnikov
Jan 28 at 16:52














$begingroup$
@Raskolnikov In fact expression $D=X-(X+Y)(X+Z)=P(A cap B)-P(A)P(B)$ is the initial expression ... and you are right that, although less symmetrical, this expression avoids to use constraint (2).
$endgroup$
– Jean Marie
Jan 28 at 17:24




$begingroup$
@Raskolnikov In fact expression $D=X-(X+Y)(X+Z)=P(A cap B)-P(A)P(B)$ is the initial expression ... and you are right that, although less symmetrical, this expression avoids to use constraint (2).
$endgroup$
– Jean Marie
Jan 28 at 17:24












$begingroup$
@Raskolnikov ... but it will remain the constraint $X+Y+Z leq 1$.
$endgroup$
– Jean Marie
Jan 28 at 17:34




$begingroup$
@Raskolnikov ... but it will remain the constraint $X+Y+Z leq 1$.
$endgroup$
– Jean Marie
Jan 28 at 17:34












$begingroup$
@Raskolnikov I have slightly modified my first sentences in order to take into account your remarks.
$endgroup$
– Jean Marie
Jan 28 at 18:50




$begingroup$
@Raskolnikov I have slightly modified my first sentences in order to take into account your remarks.
$endgroup$
– Jean Marie
Jan 28 at 18:50










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%2f3090704%2fwhat-is-the-pdf-of-x-1x-2-x-3x-4-conditionally-on-sum-limits-k-14x-k-1%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%2f3090704%2fwhat-is-the-pdf-of-x-1x-2-x-3x-4-conditionally-on-sum-limits-k-14x-k-1%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

Can a sorcerer learn a 5th-level spell early by creating spell slots using the Font of Magic feature?

ts Property 'filter' does not exist on type '{}'

Notepad++ export/extract a list of installed plugins