From 31a5231b184d48be6732a6eb500984bdb78dd807 Mon Sep 17 00:00:00 2001 From: Min Gu Jo Date: Fri, 11 Dec 2015 20:20:49 -0800 Subject: [PATCH] task plot on progress (not uploading beta values anymore --- code/utils/scripts/convolution_high_res_script.py | 3 +-- code/utils/scripts/convolution_normal_script.py | 1 - code/utils/scripts/multi_comparison_script.py | 13 ++++++++++--- 3 files changed, 11 insertions(+), 6 deletions(-) diff --git a/code/utils/scripts/convolution_high_res_script.py b/code/utils/scripts/convolution_high_res_script.py index 3750281..616141f 100644 --- a/code/utils/scripts/convolution_high_res_script.py +++ b/code/utils/scripts/convolution_high_res_script.py @@ -20,14 +20,13 @@ from __future__ import absolute_import, division, print_function import sys, os #TODO : i'm gonna fix the location when i finish the test for stimuli.py -sys.path.append(os.path.join(os.path.dirname(__file__), "../")) +sys.path.append(os.path.join(os.path.dirname('__file__'), "../")) import numpy as np import matplotlib.pyplot as plt import nibabel as nib from stimuli import * from scipy.stats import gamma from organize_behavior_data import * -from load_BOLD import * # Create the necessary directories if they do not exist dirs = ['../../../txt_output', '../../../txt_output/conv_high_res',\ diff --git a/code/utils/scripts/convolution_normal_script.py b/code/utils/scripts/convolution_normal_script.py index aa2e48a..a9db87b 100644 --- a/code/utils/scripts/convolution_normal_script.py +++ b/code/utils/scripts/convolution_normal_script.py @@ -26,7 +26,6 @@ from stimuli import * from scipy.stats import gamma from organize_behavior_data import * -from load_BOLD import * # Create the necessary directories if they do not exist dirs = ['../../../txt_output', '../../../txt_output/conv_normal',\ diff --git a/code/utils/scripts/multi_comparison_script.py b/code/utils/scripts/multi_comparison_script.py index 33bea87..9f9b794 100644 --- a/code/utils/scripts/multi_comparison_script.py +++ b/code/utils/scripts/multi_comparison_script.py @@ -10,13 +10,13 @@ 1. calculate the mean of each single beta values across subject and plot them 2. calculate the variance of each single beta values across subject and plot them 3. calculate the t-stat of each single beta values across subject and plot them -4. ????? +4. calculate the p-value of each single betav values across subject and plot them """ import sys, os ##sys.path.append(os.path.join(os.path.dirname(__file__), "../functions/")) -sys.path.append(os.path.join(os.path.dirname(__file__), "../")) +sys.path.append(os.path.join(os.path.dirname('__file__'), "../")) import numpy as np from glm import * #from convolution_normal_script import X_matrix @@ -24,7 +24,9 @@ from load_BOLD import * import nibabel as nib import matplotlib.pyplot as plt +from scipy.stats import sem from smoothing import * +from visulation import * dirs = ['../../../txt_output/multi_beta'] @@ -53,14 +55,19 @@ task_sum +=task[x] task_mean = task_sum/16 +beta_plot = present_3d(task_mean) ##this does not work... +plt.imshow(beta_plot,interpolation='nearest', cmap='seismic') -#calculate variance and plot +#calculate variance and plot +#will use sem here (standard error across 16 subjects) #calculate t-test and plot +#calculate p-value and plot +