[med-svn] [python-mne] 199/376: misc in tf cluster example

[Yaroslav Halchenko]

This is an automated email from the git hooks/post-receive script.

yoh pushed a commit to annotated tag v0.1
in repository python-mne.

commit 945bf97a16eb3b56d537120f466d0c9a9b424f8e
Author: Alexandre Gramfort <alexandre.gramfort at inria.fr>
Date:   Thu Apr 14 22:18:06 2011 -0400

    misc in tf cluster example
---
 examples/stats/plot_cluster_1samp_test_time_frequency.py | 7 +++----
 1 file changed, 3 insertions(+), 4 deletions(-)

diff --git a/examples/stats/plot_cluster_1samp_test_time_frequency.py b/examples/stats/plot_cluster_1samp_test_time_frequency.py
index 6807948..20ddeb4 100755
--- a/examples/stats/plot_cluster_1samp_test_time_frequency.py
+++ b/examples/stats/plot_cluster_1samp_test_time_frequency.py
@@ -33,14 +33,13 @@ from mne.datasets import sample
 # Set parameters
 data_path = sample.data_path('..')
 raw_fname = data_path + '/MEG/sample/sample_audvis_raw.fif'
-event_fname = data_path + '/MEG/sample/sample_audvis_raw-eve.fif'
 event_id = 1
 tmin = -0.2
 tmax = 0.5
 
 # Setup for reading the raw data
 raw = fiff.Raw(raw_fname)
-events = mne.read_events(event_fname)
+events = mne.find_events(raw)
 
 include = []
 exclude = raw.info['bads'] + ['MEG 2443', 'EEG 053'] # bads + 2 more
@@ -81,12 +80,12 @@ evoked_data = evoked_data[:,time_mask]
 times = times[time_mask]
 
 epochs_power = epochs_power[:,0,:,:]
-epochs_power = np.log(epochs_power) # take log of ratio
+epochs_power = np.log10(epochs_power) # take log of ratio
 # under the null hypothesis epochs_power should be now be 0
 
 ###############################################################################
 # Compute statistic
-threshold = 3
+threshold = 2
 T_obs, clusters, cluster_p_values, H0 = \
                    permutation_cluster_t_test(epochs_power,
                                n_permutations=100, threshold=threshold, tail=0)

-- 
Alioth's /usr/local/bin/git-commit-notice on /srv/git.debian.org/git/debian-med/python-mne.git