Ahh. For two conditions F should be the abs(T**2). You can just use a t-test for indpendent samples here instead as statfun.<br><div class="gmail_quote"><div dir="ltr">On Mon, 12 Jun 2017 at 10:13, Talitha Ford <<a href="mailto:tcford@swin.edu.au">tcford@swin.edu.au</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div style="word-wrap:break-word">
Hi Dennis,<span><br>
Thank you, this is the script I’ve been working from. The problem I am having though, is that as f-stats are >0, they do not indicate which group is larger than the other, which is what I would like to know. I have tried to use scipy.stats.ttest_ind but I get
this error:<br>
</span><span>ValueError: could not broadcast input array from shape (1000) into shape (614520)<br>
</span><span><br>
</span>
<div>The command is:</div>
<div> T_obs, clusters, cluster_p_values, H0 = clu =\</div>
<div><span class="m_1290456277243168132Apple-tab-span" style="white-space:pre-wrap"></span>spatio_temporal_cluster_test(all_data, connectivity=connectivity, n_jobs=2, stat_fun= scipy.stats.ttest_ind)</div>
<div><br>
</div>
<div>all_data is 2 lists (2 groups) of 16 and 19 participants, with 30 time points of 20484 vertices for each participant.</div>
<div><br>
</div>
<div>I hope that makes sense (and there is a possible work around!). Thanks again for you help,</div>
<div><br>
<div>
<div style="color:rgb(0,0,0);letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;word-wrap:break-word">
<div style="color:rgb(0,0,0);letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;word-wrap:break-word">
<div>
<div style="word-wrap:break-word">
Talitha </div></div></div></div></div></div></div><div style="word-wrap:break-word"><div><div><div style="color:rgb(0,0,0);letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;word-wrap:break-word"><div style="color:rgb(0,0,0);letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;word-wrap:break-word"><div>
<div style="word-wrap:break-word">
<br>
</div>
<div style="word-wrap:break-word">
<br>
On 12 Jun 2017, at 05:33, Denis-Alexander Engemann <<a href="mailto:denis.engemann@gmail.com" target="_blank">denis.engemann@gmail.com</a>> wrote:</div>
</div></div></div></div></div></div><div style="word-wrap:break-word"><div>
<div><br class="m_1290456277243168132Apple-interchange-newline">
<div>
<div dir="ltr">Hi Talitha,
<div><br>
</div>
<div>you can run the permutation clustering with a wide array of contrasts. This might be what you are looking for:</div>
<div><br>
</div>
<div><a href="http://martinos.org/mne/dev/auto_tutorials/plot_stats_cluster_spatio_temporal_2samp.html" target="_blank">http://martinos.org/mne/dev/auto_tutorials/plot_stats_cluster_spatio_temporal_2samp.html</a><br>
</div>
<div><br>
</div>
<div>I hope this helps,</div>
<div>Denis</div>
<div><br>
<div class="gmail_quote">
<div dir="ltr">On Sun, Jun 11, 2017 at 12:00 PM Talitha Ford <<a href="mailto:tcford@swin.edu.au" target="_blank">tcford@swin.edu.au</a>> wrote:<br>
</div>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Dear all,<br>
<br>
Similar to conducing a pair samples t-test on source data using spatio-temporal clustering, allowing the visualisation of clusters where condition A > condition B and vice versa, is it possible to run an independent samples t-test to visualise differences between
two groups? The 2 samples permutation tests currently available are limited in that they plot F statistics that don’t give a direction of the difference between the groups. Basically, is it possible to identify clusters that differ significantly between the
groups, as well as identify/visualise the direction in which they differ?<br>
<br>
I am currently attempting to extract the cluster values for each participant to get an overall mean for each vertex within the cluster to compare between groups, but this seems very inefficient.<br>
<br>
Cheers,<br>
Talitha<br>
<br>
_______________________________________________<br>
Mne_analysis mailing list<br>
<a href="mailto:Mne_analysis@nmr.mgh.harvard.edu" target="_blank">Mne_analysis@nmr.mgh.harvard.edu</a><br>
<a href="https://mail.nmr.mgh.harvard.edu/mailman/listinfo/mne_analysis" rel="noreferrer" target="_blank">https://mail.nmr.mgh.harvard.edu/mailman/listinfo/mne_analysis</a><br>
<br>
<br>
The information in this e-mail is intended only for the person to whom it is<br>
addressed. If you believe this e-mail was sent to you in error and the e-mail<br>
contains patient information, please contact the Partners Compliance HelpLine at<br>
<a href="http://www.partners.org/complianceline" rel="noreferrer" target="_blank">http://www.partners.org/complianceline</a> . If the e-mail was sent to you in error<br>
but does not contain patient information, please contact the sender and properly<br>
dispose of the e-mail.<br>
</blockquote>
</div>
</div>
</div>
_______________________________________________<br>
Mne_analysis mailing list<br>
<a href="mailto:Mne_analysis@nmr.mgh.harvard.edu" target="_blank">Mne_analysis@nmr.mgh.harvard.edu</a><br>
<a href="https://mail.nmr.mgh.harvard.edu/mailman/listinfo/mne_analysis" target="_blank">https://mail.nmr.mgh.harvard.edu/mailman/listinfo/mne_analysis</a><br>
<br>
<br>
The information in this e-mail is intended only for the person to whom it is<br>
addressed. If you believe this e-mail was sent to you in error and the e-mail<br>
contains patient information, please contact the Partners Compliance HelpLine at<br>
<a href="http://www.partners.org/complianceline" target="_blank">http://www.partners.org/complianceline</a> . If the e-mail was sent to you in error<br>
but does not contain patient information, please contact the sender and properly<br>
dispose of the e-mail.<br>
</div>
</div>
<br>
</div></div>
_______________________________________________<br>
Mne_analysis mailing list<br>
<a href="mailto:Mne_analysis@nmr.mgh.harvard.edu" target="_blank">Mne_analysis@nmr.mgh.harvard.edu</a><br>
<a href="https://mail.nmr.mgh.harvard.edu/mailman/listinfo/mne_analysis" rel="noreferrer" target="_blank">https://mail.nmr.mgh.harvard.edu/mailman/listinfo/mne_analysis</a><br>
<br>
<br>
The information in this e-mail is intended only for the person to whom it is<br>
addressed. If you believe this e-mail was sent to you in error and the e-mail<br>
contains patient information, please contact the Partners Compliance HelpLine at<br>
<a href="http://www.partners.org/complianceline" rel="noreferrer" target="_blank">http://www.partners.org/complianceline</a> . If the e-mail was sent to you in error<br>
but does not contain patient information, please contact the sender and properly<br>
dispose of the e-mail.<br>
</blockquote></div>