Merge pull request #196 from mingujo/min-multi

deleted the two files that are requested by Soazig

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',\

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',\

code/utils/scripts/multi_comparison_script.py

@@ -10,21 +10,23 @@
 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
 #from convolution_high_res_script import X_matrix_high_res
 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
+