Expectation of Product of Subset of Jointly Gaussian Random Vector

Multi tool use
Multi tool use

The name of the pictureThe name of the pictureThe name of the pictureClash Royale CLAN TAG#URR8PPP











up vote
0
down vote

favorite












I have a scenario, where I know that $boldsymbolz sim mathcalN( boldsymbolmu, boldsymbolSigma ) $ and I need to compute $ mathbbE( z_i^2 z_j ) $, where $z_i$ is the i-th component of vector $boldsymbolz$.



Is there any way of computing this, assuming that $boldsymbolz$ is a random vector of dimension n>2, say e.g. n=5?



If there is not a general way of computing this, is there anything for the special case of $boldsymbolmu=boldsymbol0$?



The closest thing I found so far is Isserlis' theorem, but it still seems to handle a slightly different case.







share|cite|improve this question
























    up vote
    0
    down vote

    favorite












    I have a scenario, where I know that $boldsymbolz sim mathcalN( boldsymbolmu, boldsymbolSigma ) $ and I need to compute $ mathbbE( z_i^2 z_j ) $, where $z_i$ is the i-th component of vector $boldsymbolz$.



    Is there any way of computing this, assuming that $boldsymbolz$ is a random vector of dimension n>2, say e.g. n=5?



    If there is not a general way of computing this, is there anything for the special case of $boldsymbolmu=boldsymbol0$?



    The closest thing I found so far is Isserlis' theorem, but it still seems to handle a slightly different case.







    share|cite|improve this question






















      up vote
      0
      down vote

      favorite









      up vote
      0
      down vote

      favorite











      I have a scenario, where I know that $boldsymbolz sim mathcalN( boldsymbolmu, boldsymbolSigma ) $ and I need to compute $ mathbbE( z_i^2 z_j ) $, where $z_i$ is the i-th component of vector $boldsymbolz$.



      Is there any way of computing this, assuming that $boldsymbolz$ is a random vector of dimension n>2, say e.g. n=5?



      If there is not a general way of computing this, is there anything for the special case of $boldsymbolmu=boldsymbol0$?



      The closest thing I found so far is Isserlis' theorem, but it still seems to handle a slightly different case.







      share|cite|improve this question












      I have a scenario, where I know that $boldsymbolz sim mathcalN( boldsymbolmu, boldsymbolSigma ) $ and I need to compute $ mathbbE( z_i^2 z_j ) $, where $z_i$ is the i-th component of vector $boldsymbolz$.



      Is there any way of computing this, assuming that $boldsymbolz$ is a random vector of dimension n>2, say e.g. n=5?



      If there is not a general way of computing this, is there anything for the special case of $boldsymbolmu=boldsymbol0$?



      The closest thing I found so far is Isserlis' theorem, but it still seems to handle a slightly different case.









      share|cite|improve this question











      share|cite|improve this question




      share|cite|improve this question










      asked Aug 22 at 9:26









      bonanza

      196213




      196213




















          1 Answer
          1






          active

          oldest

          votes

















          up vote
          0
          down vote













          We can find $a,b$ such that $z_i$ and $az_i+bz_j$ are independent: just make the covariance $0$. [Since they are jointly normal this is enough to say they are independent]. Now $Ez_i^2(az_i+bz_j)=Ez_i^2E(az_i+bz_j)$. The right side can be computed and the left side is $Ez_i^2(az_i+bz_j)=aEz_i^3+bEz_i^2 z_j$. From this you can compute $Ez_i^2z_j$.






          share|cite|improve this answer






















          • Thanks! But how to find corresponding a and b?
            – bonanza
            Aug 22 at 10:25










          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',
          convertImagesToLinks: true,
          noModals: false,
          showLowRepImageUploadWarning: true,
          reputationToPostImages: 10,
          bindNavPrevention: true,
          postfix: "",
          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%2f2890843%2fexpectation-of-product-of-subset-of-jointly-gaussian-random-vector%23new-answer', 'question_page');

          );

          Post as a guest






























          1 Answer
          1






          active

          oldest

          votes








          1 Answer
          1






          active

          oldest

          votes









          active

          oldest

          votes






          active

          oldest

          votes








          up vote
          0
          down vote













          We can find $a,b$ such that $z_i$ and $az_i+bz_j$ are independent: just make the covariance $0$. [Since they are jointly normal this is enough to say they are independent]. Now $Ez_i^2(az_i+bz_j)=Ez_i^2E(az_i+bz_j)$. The right side can be computed and the left side is $Ez_i^2(az_i+bz_j)=aEz_i^3+bEz_i^2 z_j$. From this you can compute $Ez_i^2z_j$.






          share|cite|improve this answer






















          • Thanks! But how to find corresponding a and b?
            – bonanza
            Aug 22 at 10:25














          up vote
          0
          down vote













          We can find $a,b$ such that $z_i$ and $az_i+bz_j$ are independent: just make the covariance $0$. [Since they are jointly normal this is enough to say they are independent]. Now $Ez_i^2(az_i+bz_j)=Ez_i^2E(az_i+bz_j)$. The right side can be computed and the left side is $Ez_i^2(az_i+bz_j)=aEz_i^3+bEz_i^2 z_j$. From this you can compute $Ez_i^2z_j$.






          share|cite|improve this answer






















          • Thanks! But how to find corresponding a and b?
            – bonanza
            Aug 22 at 10:25












          up vote
          0
          down vote










          up vote
          0
          down vote









          We can find $a,b$ such that $z_i$ and $az_i+bz_j$ are independent: just make the covariance $0$. [Since they are jointly normal this is enough to say they are independent]. Now $Ez_i^2(az_i+bz_j)=Ez_i^2E(az_i+bz_j)$. The right side can be computed and the left side is $Ez_i^2(az_i+bz_j)=aEz_i^3+bEz_i^2 z_j$. From this you can compute $Ez_i^2z_j$.






          share|cite|improve this answer














          We can find $a,b$ such that $z_i$ and $az_i+bz_j$ are independent: just make the covariance $0$. [Since they are jointly normal this is enough to say they are independent]. Now $Ez_i^2(az_i+bz_j)=Ez_i^2E(az_i+bz_j)$. The right side can be computed and the left side is $Ez_i^2(az_i+bz_j)=aEz_i^3+bEz_i^2 z_j$. From this you can compute $Ez_i^2z_j$.







          share|cite|improve this answer














          share|cite|improve this answer



          share|cite|improve this answer








          edited Aug 22 at 9:53

























          answered Aug 22 at 9:42









          Kavi Rama Murthy

          23.5k2933




          23.5k2933











          • Thanks! But how to find corresponding a and b?
            – bonanza
            Aug 22 at 10:25
















          • Thanks! But how to find corresponding a and b?
            – bonanza
            Aug 22 at 10:25















          Thanks! But how to find corresponding a and b?
          – bonanza
          Aug 22 at 10:25




          Thanks! But how to find corresponding a and b?
          – bonanza
          Aug 22 at 10:25












           

          draft saved


          draft discarded


























           


          draft saved


          draft discarded














          StackExchange.ready(
          function ()
          StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fmath.stackexchange.com%2fquestions%2f2890843%2fexpectation-of-product-of-subset-of-jointly-gaussian-random-vector%23new-answer', 'question_page');

          );

          Post as a guest













































































          aW,41IxbaVXUi9gXOWCgfZ7upT,opyeWkMlzY1ow f,3xjnPgkt53KW,M,OTlCX
          S88GEr,hgSOf JTMrqVXf8DMJc3pQ2dMaDOpiPCi MvYKWYfpXuBlNXZO kk02K4LcrN,7z

          這個網誌中的熱門文章

          How to combine Bézier curves to a surface?

          Propositional logic and tautologies

          Distribution of Stopped Wiener Process with Stochastic Volatility