Some tips to use wavelet toolbox

55 sec read

Wavelet toolbox is a useful tool to study hyperscanning data. Many recent publications on NIRS hyperscanning use wavelet coherence to quantify the relationship between two interacting brains (e.g. Baker et al 2016, Nozawa et al 2016). You can see more information about wavelet coherence at https://www.alivelearn.net/?p=1169

wavelet
wavelet

Here are some tips to use the toolbox:

1. It often takes a long time to run Monte Carlo simulation. You may use ‘mcc’=0 to disable it.

figure;wtc(signal(:,jj),signal(:,jj+22),'mcc',0);

2. If you need to get the values of the result (instead of the graphic), you may specify the return value

[Rsq,period,scale,coi,sig95] = wtc(signal1,signal2,'mcc',0); %Rsq is a complex matrix

3. If you are only interested in a certain band, you can specify the MaxScale (i.e. ms) parameter. More information at https://www.alivelearn.net/?p=1518

[Rsq,period,scale,coi,sig95] = wtc(signal1,signal2,'mcc',0, 'ms', 128);

4. If you are interested in finding the “phase” information (visualized by the arrows), you may use xwt function. The returned value is a complex matrix and you can calculate the phase.

[Rsq,period,scale,coi,sig95] = xwt(signal(:,jj),signal(:,jj+22));
figure;imagesc(angle(Rsq));

5. To visualize the power of a single signal, you may use wt, which I personally feel much better than FFT.

figure;wt(signal(:,jj));

6. To change the density of the arrows, you may specify the ArrowDensity parameter

figure;wtc(signal(:,jj),signal(:,jj+22),'mcc',0,'ArrowDensity',[30 30]);

Do you have any tips? Please let me know.


Receive email notification via email 博客有新内容通知我
文献鸟,永远不错过新文献!

Don't want to miss new papers in your field? Check out Stork we developed:


第二十六期 fNIRS Journal Club 视频

来自上海精神卫生中心的胡志善博士为大家介绍了一篇他们去年发表的综述(静息态近红外在研究发育的应用)。视频如下: Youtube:https://youtu.be/C7xIRrBUhhg Youku: https://v.youku.com/v_show/id_XNTgzNTQ5ODg3Mg==.html
Xu Cui
4 sec read

第二十六期 fNIRS Journal Club 通知 2022/1/15,10:00am

来自上海精神卫生中心的胡志善博士将为大家介绍具体一篇他们去年发表的综述(静息态近红外在研究发育的应用)。热烈欢迎大家参与讨论。 时间: 北京时间2022年1月15日上午10:00地点: https://zoom.com房间号: 897 8179 5073密码: 333858 胡博士要讲解的文章如下: Hu, Liu, Dong, Niu (2020) Applications of Resting-State fNIRS in the Developing Brain: A Review...
Xu Cui
6 sec read

第二十五期 fNIRS Journal Club 视频

来自华东师范大学的胡银莹博士为大家介绍了具体如何计算和验证超扫描实验中脑间同步性,并分享做fNIRS实验相关方面的应用经验。视频如下: Youtube:https://youtu.be/xI_6wKd06io Youku:上传后被youku屏蔽,原因不详。
Xu Cui
3 sec read

6 Replies to “Some tips to use wavelet toolbox”

  1. Dear Cui,

    Thanks a lot for your post.
    Nevertheless it is not clear to me why we must extract the phases on the cross-correlation and not on the coherence whereas this one is present in the figure. Could you explain ? And then, il web want to extract not only the mean but the dynamic (fluctuation) over time, Is there a solution other than compute a loop?
    Thanks a lot !

  2. Hi Dr. Xu,

    Thanks so much for sharing all this info on your site! I was hoping to get your insight on an analysis decision with WTC: in your 2012 paper, you analyzed the coherence band at period=64 because there was good reason to think that the activity pattern would be at that frequency (task timing = 7s). What would you suggest for hyperscanning studies that aren’t event related, and are instead one long scanning period over the entirety of the subjects’ interaction? Is there any theory about which frequency band is most important in that case, or is the best option just to average the coherence value over a large frequency band?

    Thanks so much!

  3. @Shannon
    It’s a good question, and a difficult question! Let me share briefly with you the thought behind the band choice:

    When I was doing data analysis and saw the coherence figures, there was an obvious red band in the plot. So this is something interesting, I thought. Then I examined its frequency, and found it is at ~7s, which turned out to be our task frequency. So in our paper we wrote we use the coherence value around the task frequency band.

    As you see, we discovered the pattern (red band) first, then try to rationalize this pattern (find a reason choose this band). In the paper it sounds we already have a good reason to choose the band in the first place, but the reality is reverse.

    So while I can’t answer your question directly, I do encourage you to see if there are any patterns emerging in your figure.

  4. Dear Dr.Cui:
    Thanks so much for sharing these ideas,It’s very useful.
    Btw. I want to ask you whether there wiil be some differences between the reult figures that one I did Monte Carlo simulation and one did not(‘mcc’=0).

    sincerely ,it really cost so much time to do Monte Carlo! =_=

Leave a Reply

Your email address will not be published. Required fields are marked *