Menu

Why the design matrix needs to be orthogonalized?

Help
qwer1304
2013-03-20
2013-03-20
  • qwer1304

    qwer1304 - 2013-03-20

    Hi,

    Why the design matrix needs to be orthogonalized? What if it weren't?

    Thx,
    David

     
  • Dominik Bach

    Dominik Bach - 2013-03-20

    The design matrix is NOT orthogonalised.

    Only pmods are orthogonalised with respect to the canonical response.

    The reason is that any event elicits a response. A pmod allows you to assess how this general response is modulated by a continuous variable.

     
  • qwer1304

    qwer1304 - 2013-03-20

    Well, in the code, under 'create design matrix', apparently it does:

    % orthogonalize after convolution
    foo=spm_orth(tmp.Xc{cond});
    

    That's done for all conditions.

    Thx

     
  • Dominik Bach

    Dominik Bach - 2013-03-20

    No, the design matrix is not orthogonalised

    As you write, the code snippet you copied works PER condition, i. e. it orthogonolises the basis functions and all pmods WITHIN one condition - this is not the full design matrix.

    Best,
    Dom

     

    Last edit: Dominik Bach 2013-03-20
  • qwer1304

    qwer1304 - 2013-03-20

    I see, thx for clarifying this.

    So why is orthogonalization within a condition required?

    David

     
  • Dominik Bach

    Dominik Bach - 2013-03-20

    Because you want to separate the overall response from the impact of the pmod.

     

Log in to post a comment.