<div dir="ltr">Hi Alex, <div><br></div><div>Thank you for your reply. </div><div><br></div><div>This is my code, with two versions of a diagonal connectivity matrix (using version 14.1 of MNE): </div><div><div># Implement cluster permutation</div><div>numChan = Data.shape[1] # Data is shaped: Trials (115) x channels (217) x timepoints (30) numpy array</div><div>#connectivity = np.diagflat( range(0,numChan), k = 0)</div><div>connectivity = np.diagflat( np.ones((1,numChan)), k = 0) # N x N diagonal matrix with N being # of electrodes</div><div>TestData = np.swapaxes(Data, 1,2) # trials x time x space (electrodes)</div><div>T_obs, clusters, p_values, H0 = mne.stats.permutation_cluster_1samp_test(TestData, n_permutations = 500, tail = 0, connectivity = connectivity)</div></div><div><br></div><div>Both versions of the connectivity matrix now give me the following error: </div><div>'numpy.ndarray' object has no attribute 'tocsr'<br></div><div><br></div><div>Thanks,</div><div>Yvonne</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Sep 19, 2017 at 1:02 AM, Alexandre Gramfort <span dir="ltr"><<a href="mailto:alexandre.gramfort@inria.fr" target="_blank">alexandre.gramfort@inria.fr</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi Yvonne<br>
<br>
I would have said that a sparse diagonal matrix should work.<br>
<br>
can you share the code snippet you used?<br>
<br>
Alex<br>
<div><div class="h5"><br>
<br>
On Tue, Sep 19, 2017 at 3:17 AM, Yvonne Fonken <<a href="mailto:yfonken@berkeley.edu">yfonken@berkeley.edu</a>> wrote:<br>
> Hi all,<br>
><br>
> I am trying to apply the cluster permutation function in MNE on intracranial<br>
> data (electrodes by time points), either by using<br>
> mne.stats.permutation_cluster_<wbr>test, or<br>
> mne.stats.permutation_cluster_<wbr>1samp_test. In this case, I do not want to<br>
> cluster across electrodes, but I want to include this dimension in order to<br>
> correct for multiple comparisons.<br>
><br>
> My question is: how do I specify to only cluster over time, and not<br>
> electrodes? Can I do this using the connectivity matrix, and if yes, what<br>
> would a connectivity matrix look like when all electrodes do not have<br>
> neighbors? I've tried using diagonal matrices already, but I still get<br>
> clusters spanning multiple electrodes.<br>
><br>
> Thanks,<br>
> Yvonne<br>
><br>
> --<br>
> Yvonne Fonken<br>
> PhD Candidate<br>
> Helen Wills Neuroscience Institute<br>
> University of California, Berkeley<br>
><br>
</div></div>> ______________________________<wbr>_________________<br>
> Mne_analysis mailing list<br>
> <a href="mailto:Mne_analysis@nmr.mgh.harvard.edu">Mne_analysis@nmr.mgh.harvard.<wbr>edu</a><br>
> <a href="https://mail.nmr.mgh.harvard.edu/mailman/listinfo/mne_analysis" rel="noreferrer" target="_blank">https://mail.nmr.mgh.harvard.<wbr>edu/mailman/listinfo/mne_<wbr>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<br>
> e-mail<br>
> contains patient information, please contact the Partners Compliance<br>
> HelpLine at<br>
> <a href="http://www.partners.org/complianceline" rel="noreferrer" target="_blank">http://www.partners.org/<wbr>complianceline</a> . If the e-mail was sent to you in<br>
> error<br>
> but does not contain patient information, please contact the sender and<br>
> properly<br>
> dispose of the e-mail.<br>
><br>
______________________________<wbr>_________________<br>
Mne_analysis mailing list<br>
<a href="mailto:Mne_analysis@nmr.mgh.harvard.edu">Mne_analysis@nmr.mgh.harvard.<wbr>edu</a><br>
<a href="https://mail.nmr.mgh.harvard.edu/mailman/listinfo/mne_analysis" rel="noreferrer" target="_blank">https://mail.nmr.mgh.harvard.<wbr>edu/mailman/listinfo/mne_<wbr>analysis</a><br>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature">Yvonne Fonken<br>Graduate student<br>Helen Wills Neuroscience Institute<br>University of California, Berkeley<br></div>
</div>