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 <[email protected]> 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 _______________________________________________ debian-med-commit mailing list [email protected] http://lists.alioth.debian.org/cgi-bin/mailman/listinfo/debian-med-commit
