roi_analysis.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 24,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "import os\n",
    "import os.path as op\n",
    "import glob\n",
    "import numpy as np\n",
    "import pandas as pd\n",
    "from pandas import DataFrame\n",
    "import nibabel as nib\n",
    "import scipy\n",
    "import matplotlib\n",
    "import seaborn as sns\n",
    "import statsmodels\n",
    "import matplotlib.pyplot as plt\n",
    "import statsmodels.formula.api as smf\n",
    "from statsmodels.stats.multitest import multipletests, fdrcorrection\n",
    "%matplotlib inline"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 25,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "#preliminary housekeeping\n",
    "home_dir = '/data/home/iballard/mag/'\n",
    "subj_file = home_dir + 'subjects.txt'\n",
    "subs = list(np.loadtxt(subj_file,'string'))\n",
    "os.chdir(home_dir)\n",
    "runs = map(str,range(1,5))"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Define experimental info and ROIs"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 26,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "smoothing = 'smoothed'\n",
    "regspace = 'epi'\n",
    "\n",
    "exp = 'mag-SVtotaldiff'\n",
    "\n",
    "rois = ['_amfg',\n",
    "        '_pmfg',\n",
    "        '_fpc',\n",
    "        '_ifg',\n",
    "        '_psfs',\n",
    "        '_ifs',\n",
    "       '_ains',\n",
    "       '.vmpfc_']"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Main functions for extracting and analysing ROI data"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 27,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "def extract_copes(exp,regspace,smoothing,contrast):\n",
    "    df = pd.DataFrame(columns=('sub','contrast','mask_vox','hemi','label','value'))\n",
    "\n",
    "    for sub in subs:\n",
    "        cope = home_dir + '/analysis/' + exp + '/' + sub + '/ffx/' + regspace + \\\n",
    "            '/' + smoothing + '/' + contrast + '/cope1.nii.gz'\n",
    "        cope = nib.load(cope).get_data()\n",
    "\n",
    "        masks = glob.glob(home_dir + '/data/' + sub + '/masks/*.*nii.gz')\n",
    "        masks = [m for m in masks if any(substring in m for substring in rois)]\n",
    "\n",
    "        for m in masks:\n",
    "            if len(m.split('/')[-1].split('.'))==4:\n",
    "                hemi = m.split('/')[-1].split('.')[0]\n",
    "                label = m.split('/')[-1].split('.')[1]\n",
    "            else:\n",
    "                hemi = 'na'\n",
    "                label = m.split('/')[-1].split('.')[0]\n",
    "\n",
    "            m = nib.load(m).get_data().astype(bool)\n",
    "            num_voxels = m.sum()\n",
    "            func_masked = cope[m]\n",
    "            row = pd.DataFrame([dict(sub = sub, \n",
    "                                     contrast =contrast, \n",
    "                                     mask_vox = num_voxels, \n",
    "                                     hemi = hemi,\n",
    "                                     label = label, \n",
    "                                     value = np.nanmean(func_masked)), ])\n",
    "            df = df.append(row, ignore_index = True)\n",
    "    return df"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 28,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "def plot_fig(res,contrast):\n",
    "    plt.figure(figsize=(10, 6))\n",
    "\n",
    "    #relabel\n",
    "    res = res.copy()\n",
    "    res['ROI'] = [label_map[x] for x in res['label'].values]\n",
    "    \n",
    "    #sort dataframe\n",
    "    res['ROI'] = pd.Categorical(res['ROI'],roi_plotting_order)\n",
    "    res = res.sort('ROI')\n",
    "\n",
    "    #plot\n",
    "    sns.barplot(x='ROI',y='val',data = res,\n",
    "                **{'yerr':res['err'],\n",
    "                   'error_kw':dict(ecolor='grey', lw=2.5)}\n",
    "                )\n",
    "    sns.despine()\n",
    "    plt.ylabel('Beta (AU)')\n",
    "    plt.title(title_map[contrast])\n",
    "    plt.savefig('/home/iballard/Dropbox/mag_figs/' + contrast + '.svg')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 29,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "def run_stats(df):\n",
    "    results_df =pd.DataFrame(columns=('contrast','label','nvox',\n",
    "                                      'std_vox','p_unc','val','err','std_err','dz','t'))\n",
    "    for roi in set(df['label']):\n",
    "\n",
    "        new = df[df['label'] == roi]\n",
    "\n",
    "        #run mixed effects model\n",
    "        model_str = 'value ~ 1'\n",
    "        mixedfx = smf.mixedlm(model_str, new, groups = new['sub'])    \n",
    "        mixedfx = mixedfx.fit()\n",
    "        \n",
    "        #extract p value and errors\n",
    "        p = mixedfx.pvalues['Intercept']\n",
    "        val = mixedfx.params[0]\n",
    "        err = mixedfx.bse_fe[0] #standard error\n",
    "        t = mixedfx.tvalues[0]\n",
    "        dz = val / (err * np.sqrt(len(subs))) #cohen's dz\n",
    "        conf_int = mixedfx.conf_int()\n",
    "        conf_int = (conf_int.loc['Intercept',1] - conf_int.loc['Intercept',0]) / 2.0 #plots 95% conf\n",
    "        \n",
    "        #update dataframe\n",
    "        row = pd.DataFrame([dict(contrast = title_map[contrast], \n",
    "                                 label = roi,\n",
    "                                 nvox = np.mean(new['mask_vox']),\n",
    "                                 val = val,\n",
    "                                 dz = dz,\n",
    "                                 t = t,\n",
    "                                 std_err = err,\n",
    "                                 err = conf_int,\n",
    "                                p_unc = p), ])\n",
    "        results_df = results_df.append(row, ignore_index = True)\n",
    "        \n",
    "    #correct for multiple comparisons\n",
    "    sig,pcorr = fdrcorrection(results_df['p_unc'],.05)\n",
    "    results_df['p_corr'] = pcorr\n",
    "    results_df['sig'] = sig\n",
    "    return results_df\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 30,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "def run_contrast(exp,regspace,smoothing,contrast):\n",
    "    df = extract_copes(exp,regspace,smoothing,contrast)\n",
    "    df['ROI'] = [label_map[x] for x in df['label'].values]\n",
    "    \n",
    "    df.to_csv('/home/iballard/Dropbox/Magnitude_fMRI/roi_analysis/' + exp  + '_' + contrast + '.csv',\n",
    "              index = False)\n",
    "    \n",
    "    results = run_stats(df)\n",
    "    results['ROI'] = [label_map[x] for x in results['label'].values]\n",
    "\n",
    "    plot_fig(results,contrast)\n",
    "\n",
    "    return results"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Set up figure aesthetics"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 31,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "label_map = {'yeo17_amfg':'aMFG','yeo17_pmfg':'pMFG','yeo17_ains':'aINS','yeo17_fpc':'FPC',\n",
    "            'yeo17_ifs':'IFS','yeo17_ifg':'IFG','yeo17_psfs':'FEF','vmpfc_bartre':'vmPFC'}\n",
    "roi_plotting_order = ['IFS','FPC','pMFG','FEF','aINS','IFG','aMFG','vmPFC']"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 40,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAcwAAABGCAYAAABBh6SMAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAApVJREFUeJzt2jFPE3Ecx+EfFEoHAw5eDPEFsHTsa/AtOPRdOLr7Blxd\nXDr4BnwTsJjc4uDMUERJMBTIWQeD0oTKb/n3D83zjHfX5JvcpZ+kvY35fD4PAOC/NmsPAIDHQDAB\nIEEwASBBMAEgQTABIGFr2YnZbBZt20bTNNHr9Va5CQCq6LouptNpDIfDGAwGC+eWBrNt2xiPx8XH\nAcBDM5lMYjQaLRxbGsymaSIi4t3BQTT9ftlllXx9dVp7QlEfX7ysPaGoN4dHtScUtXsxuv+iR+z9\n9/V+Pg/jS+0JRb3+8Kn2hCJONy/j7e7nvw28bWkwb36Gbfr92N/ZKbeuorNn27UnFLXzfK/2hKL2\nn67nc3ljr7/e9+/J9X7tCUVtx0ntCUU1vwb3X/SI3fVXpJd+ACBBMAEgQTABIEEwASBBMAEgQTAB\nIEEwASBBMAEgQTABIEEwASBBMAEgQTABIEEwASBBMAEgQTABIEEwASBBMAEgQTABIEEwASBBMAEg\nQTABIEEwASBBMAEgQTABIEEwASBBMAEgQTABIEEwASBBMAEgQTABIEEwASBBMAEgQTABIEEwASBB\nMAEgQTABIEEwASBBMAEgQTABIEEwASBBMAEgQTABIEEwASBBMAEgQTABIEEwASBBMAEgQTABIEEw\nASBBMAEgQTABIEEwASBha9mJrusiImJ6dbWyMat2cnJde0JRlztntScUdfzjsvaEon5erPf9Oz8/\nrj2hqOv4VntCUdPNWe0JRZxu/vleuWngbRvz+Xx+14eOjo5iPB6XXQYAD9BkMonRaLRwbGkwZ7NZ\ntG0bTdNEr9dbyUAAqKnruphOpzEcDmMwGCycWxpMAOAfL/0AQIJgAkCCYAJAgmACQMJvCo1l7Seu\n4MUAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7fe241122310>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "#get colors from ROI figure and match them\n",
    "ctab = '/home/iballard/mag/data/fsaverage/label/pfc.ctab'\n",
    "color = np.loadtxt(ctab,object)\n",
    "custom_palette = []\n",
    "for row in color:\n",
    "    rgb = (int(row[2]),int(row[3]),int(row[4]))\n",
    "    hex_rgb = '#%02x%02x%02x' % rgb\n",
    "    roi_id = label_map[row[1]]\n",
    "    custom_palette.append((hex_rgb, roi_id))\n",
    "\n",
    "custom_palette.sort(key=lambda x: roi_plotting_order.index(x[1]))\n",
    "custom_palette = [x[0] for x in custom_palette]\n",
    "\n",
    "sns.set_palette(custom_palette)\n",
    "sns.palplot(sns.color_palette())"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 38,
   "metadata": {
    "collapsed": false,
    "scrolled": true
   },
   "outputs": [],
   "source": [
    "sns.set_context('paper',font_scale = 2.5)\n",
    "sns.set(font=\"Times New Roman\")\n",
    "sns.set_style(\"white\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 39,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "#for figure plotting\n",
    "title_map = {'decision': 'Decision Period', \n",
    "             'SV_total': 'Subjective Value',\n",
    "             'neg_SV_diff': 'Difficulty',\n",
    "             'low_minus_high_SV_diff': 'Difficulty (High minus Low Magnitude)',\n",
    "             'high_minus_low_SV': 'Subjective Value (High minus Low Magnitude)',\n",
    "            'high_minus_low': 'High minus Low Magnitude'}\n",
    "contrasts = title_map.keys()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Run analysis"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 42,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "all_results = []\n",
    "\n",
    "for contrast in contrasts:\n",
    "    results = run_contrast(exp,regspace,smoothing,contrast)\n",
    "    all_results.append(results.copy())\n",
    "    \n",
    "all_results = pd.concat(all_results)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 36,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div style=\"max-height:1000px;max-width:1500px;overflow:auto;\">\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>ROI</th>\n",
       "      <th>val</th>\n",
       "      <th>std_err</th>\n",
       "      <th>t</th>\n",
       "      <th>p_unc</th>\n",
       "      <th>p_corr</th>\n",
       "      <th>sig</th>\n",
       "      <th>dz</th>\n",
       "      <th>95_CI_low</th>\n",
       "      <th>95_CI_high</th>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>contrast</th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>Subjective Value</th>\n",
       "      <td>IFG</td>\n",
       "      <td>-0.227902</td>\n",
       "      <td>0.343442</td>\n",
       "      <td>-0.663584</td>\n",
       "      <td>5.069563e-01</td>\n",
       "      <td>7.039552e-01</td>\n",
       "      <td>False</td>\n",
       "      <td>-0.152237</td>\n",
       "      <td>-0.901036</td>\n",
       "      <td>0.445231</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Subjective Value</th>\n",
       "      <td>FPC</td>\n",
       "      <td>-0.075737</td>\n",
       "      <td>0.383605</td>\n",
       "      <td>-0.197434</td>\n",
       "      <td>8.434879e-01</td>\n",
       "      <td>8.434879e-01</td>\n",
       "      <td>False</td>\n",
       "      <td>-0.045294</td>\n",
       "      <td>-0.827588</td>\n",
       "      <td>0.676115</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Subjective Value</th>\n",
       "      <td>aMFG</td>\n",
       "      <td>-0.295708</td>\n",
       "      <td>0.364751</td>\n",
       "      <td>-0.810714</td>\n",
       "      <td>4.175301e-01</td>\n",
       "      <td>7.039552e-01</td>\n",
       "      <td>False</td>\n",
       "      <td>-0.185990</td>\n",
       "      <td>-1.010606</td>\n",
       "      <td>0.419190</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Subjective Value</th>\n",
       "      <td>aINS</td>\n",
       "      <td>-0.077649</td>\n",
       "      <td>0.320861</td>\n",
       "      <td>-0.242003</td>\n",
       "      <td>8.087778e-01</td>\n",
       "      <td>8.434879e-01</td>\n",
       "      <td>False</td>\n",
       "      <td>-0.055519</td>\n",
       "      <td>-0.706525</td>\n",
       "      <td>0.551226</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Subjective Value</th>\n",
       "      <td>IFS</td>\n",
       "      <td>0.877997</td>\n",
       "      <td>0.425268</td>\n",
       "      <td>2.064572</td>\n",
       "      <td>3.896351e-02</td>\n",
       "      <td>1.558540e-01</td>\n",
       "      <td>False</td>\n",
       "      <td>0.473645</td>\n",
       "      <td>0.044487</td>\n",
       "      <td>1.711508</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Subjective Value</th>\n",
       "      <td>pMFG</td>\n",
       "      <td>-0.392609</td>\n",
       "      <td>0.355032</td>\n",
       "      <td>-1.105841</td>\n",
       "      <td>2.687953e-01</td>\n",
       "      <td>7.039552e-01</td>\n",
       "      <td>False</td>\n",
       "      <td>-0.253697</td>\n",
       "      <td>-1.088458</td>\n",
       "      <td>0.303241</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Subjective Value</th>\n",
       "      <td>vmPFC</td>\n",
       "      <td>1.727025</td>\n",
       "      <td>0.551990</td>\n",
       "      <td>3.128725</td>\n",
       "      <td>1.755664e-03</td>\n",
       "      <td>1.404531e-02</td>\n",
       "      <td>True</td>\n",
       "      <td>0.717779</td>\n",
       "      <td>0.645144</td>\n",
       "      <td>2.808905</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Subjective Value</th>\n",
       "      <td>FEF</td>\n",
       "      <td>-0.361484</td>\n",
       "      <td>0.572772</td>\n",
       "      <td>-0.631113</td>\n",
       "      <td>5.279664e-01</td>\n",
       "      <td>7.039552e-01</td>\n",
       "      <td>False</td>\n",
       "      <td>-0.144787</td>\n",
       "      <td>-1.484097</td>\n",
       "      <td>0.761129</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Decision Period</th>\n",
       "      <td>IFG</td>\n",
       "      <td>-8.018357</td>\n",
       "      <td>1.609433</td>\n",
       "      <td>-4.982100</td>\n",
       "      <td>6.289789e-07</td>\n",
       "      <td>1.257958e-06</td>\n",
       "      <td>True</td>\n",
       "      <td>-1.142972</td>\n",
       "      <td>-11.172789</td>\n",
       "      <td>-4.863926</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Decision Period</th>\n",
       "      <td>FPC</td>\n",
       "      <td>-4.997577</td>\n",
       "      <td>2.017170</td>\n",
       "      <td>-2.477519</td>\n",
       "      <td>1.322992e-02</td>\n",
       "      <td>2.116787e-02</td>\n",
       "      <td>True</td>\n",
       "      <td>-0.568382</td>\n",
       "      <td>-8.951157</td>\n",
       "      <td>-1.043997</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Decision Period</th>\n",
       "      <td>aMFG</td>\n",
       "      <td>-3.004705</td>\n",
       "      <td>1.558327</td>\n",
       "      <td>-1.928161</td>\n",
       "      <td>5.383516e-02</td>\n",
       "      <td>7.178021e-02</td>\n",
       "      <td>False</td>\n",
       "      <td>-0.442350</td>\n",
       "      <td>-6.058971</td>\n",
       "      <td>0.049560</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Decision Period</th>\n",
       "      <td>aINS</td>\n",
       "      <td>9.606990</td>\n",
       "      <td>1.621377</td>\n",
       "      <td>5.925203</td>\n",
       "      <td>3.119109e-09</td>\n",
       "      <td>1.247644e-08</td>\n",
       "      <td>True</td>\n",
       "      <td>1.359335</td>\n",
       "      <td>6.429149</td>\n",
       "      <td>12.784831</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Decision Period</th>\n",
       "      <td>IFS</td>\n",
       "      <td>10.328838</td>\n",
       "      <td>1.913305</td>\n",
       "      <td>5.398427</td>\n",
       "      <td>6.722772e-08</td>\n",
       "      <td>1.792739e-07</td>\n",
       "      <td>True</td>\n",
       "      <td>1.238484</td>\n",
       "      <td>6.578829</td>\n",
       "      <td>14.078847</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Decision Period</th>\n",
       "      <td>pMFG</td>\n",
       "      <td>1.606482</td>\n",
       "      <td>1.820506</td>\n",
       "      <td>0.882437</td>\n",
       "      <td>3.775405e-01</td>\n",
       "      <td>3.775405e-01</td>\n",
       "      <td>False</td>\n",
       "      <td>0.202445</td>\n",
       "      <td>-1.961644</td>\n",
       "      <td>5.174607</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Decision Period</th>\n",
       "      <td>vmPFC</td>\n",
       "      <td>-33.996927</td>\n",
       "      <td>2.624094</td>\n",
       "      <td>-12.955682</td>\n",
       "      <td>2.181896e-38</td>\n",
       "      <td>1.745517e-37</td>\n",
       "      <td>True</td>\n",
       "      <td>-2.972237</td>\n",
       "      <td>-39.140056</td>\n",
       "      <td>-28.853797</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Decision Period</th>\n",
       "      <td>FEF</td>\n",
       "      <td>-2.140544</td>\n",
       "      <td>2.099276</td>\n",
       "      <td>-1.019658</td>\n",
       "      <td>3.078905e-01</td>\n",
       "      <td>3.518749e-01</td>\n",
       "      <td>False</td>\n",
       "      <td>-0.233926</td>\n",
       "      <td>-6.255050</td>\n",
       "      <td>1.973961</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Difficulty</th>\n",
       "      <td>IFG</td>\n",
       "      <td>-0.231936</td>\n",
       "      <td>0.653274</td>\n",
       "      <td>-0.355036</td>\n",
       "      <td>7.225626e-01</td>\n",
       "      <td>7.225626e-01</td>\n",
       "      <td>False</td>\n",
       "      <td>-0.081451</td>\n",
       "      <td>-1.512329</td>\n",
       "      <td>1.048458</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Difficulty</th>\n",
       "      <td>FPC</td>\n",
       "      <td>1.912566</td>\n",
       "      <td>0.710861</td>\n",
       "      <td>2.690492</td>\n",
       "      <td>7.134664e-03</td>\n",
       "      <td>1.286546e-02</td>\n",
       "      <td>True</td>\n",
       "      <td>0.617241</td>\n",
       "      <td>0.519304</td>\n",
       "      <td>3.305828</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Difficulty</th>\n",
       "      <td>aMFG</td>\n",
       "      <td>0.745508</td>\n",
       "      <td>0.447358</td>\n",
       "      <td>1.666469</td>\n",
       "      <td>9.562007e-02</td>\n",
       "      <td>1.274934e-01</td>\n",
       "      <td>False</td>\n",
       "      <td>0.382314</td>\n",
       "      <td>-0.131297</td>\n",
       "      <td>1.622313</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Difficulty</th>\n",
       "      <td>aINS</td>\n",
       "      <td>1.154436</td>\n",
       "      <td>0.356968</td>\n",
       "      <td>3.234007</td>\n",
       "      <td>1.220667e-03</td>\n",
       "      <td>3.255111e-03</td>\n",
       "      <td>True</td>\n",
       "      <td>0.741932</td>\n",
       "      <td>0.454792</td>\n",
       "      <td>1.854081</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Difficulty</th>\n",
       "      <td>IFS</td>\n",
       "      <td>2.450334</td>\n",
       "      <td>0.481019</td>\n",
       "      <td>5.094045</td>\n",
       "      <td>3.505033e-07</td>\n",
       "      <td>2.804027e-06</td>\n",
       "      <td>True</td>\n",
       "      <td>1.168654</td>\n",
       "      <td>1.507554</td>\n",
       "      <td>3.393115</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Difficulty</th>\n",
       "      <td>pMFG</td>\n",
       "      <td>2.078334</td>\n",
       "      <td>0.632626</td>\n",
       "      <td>3.285247</td>\n",
       "      <td>1.018930e-03</td>\n",
       "      <td>3.255111e-03</td>\n",
       "      <td>True</td>\n",
       "      <td>0.753687</td>\n",
       "      <td>0.838409</td>\n",
       "      <td>3.318259</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Difficulty</th>\n",
       "      <td>vmPFC</td>\n",
       "      <td>-1.919153</td>\n",
       "      <td>0.724114</td>\n",
       "      <td>-2.650347</td>\n",
       "      <td>8.040915e-03</td>\n",
       "      <td>1.286546e-02</td>\n",
       "      <td>True</td>\n",
       "      <td>-0.608031</td>\n",
       "      <td>-3.338390</td>\n",
       "      <td>-0.499916</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Difficulty</th>\n",
       "      <td>FEF</td>\n",
       "      <td>0.390121</td>\n",
       "      <td>0.623569</td>\n",
       "      <td>0.625627</td>\n",
       "      <td>5.315599e-01</td>\n",
       "      <td>6.074971e-01</td>\n",
       "      <td>False</td>\n",
       "      <td>0.143529</td>\n",
       "      <td>-0.832051</td>\n",
       "      <td>1.612293</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>High minus Low Magnitude</th>\n",
       "      <td>IFG</td>\n",
       "      <td>-4.058899</td>\n",
       "      <td>0.753600</td>\n",
       "      <td>-5.386009</td>\n",
       "      <td>7.203916e-08</td>\n",
       "      <td>5.763133e-07</td>\n",
       "      <td>True</td>\n",
       "      <td>-1.235635</td>\n",
       "      <td>-5.535928</td>\n",
       "      <td>-2.581869</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>High minus Low Magnitude</th>\n",
       "      <td>FPC</td>\n",
       "      <td>0.258278</td>\n",
       "      <td>1.016781</td>\n",
       "      <td>0.254015</td>\n",
       "      <td>7.994840e-01</td>\n",
       "      <td>7.994840e-01</td>\n",
       "      <td>False</td>\n",
       "      <td>0.058275</td>\n",
       "      <td>-1.734576</td>\n",
       "      <td>2.251131</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>High minus Low Magnitude</th>\n",
       "      <td>aMFG</td>\n",
       "      <td>-0.804398</td>\n",
       "      <td>0.778557</td>\n",
       "      <td>-1.033190</td>\n",
       "      <td>3.015150e-01</td>\n",
       "      <td>6.030300e-01</td>\n",
       "      <td>False</td>\n",
       "      <td>-0.237030</td>\n",
       "      <td>-2.330342</td>\n",
       "      <td>0.721547</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>High minus Low Magnitude</th>\n",
       "      <td>aINS</td>\n",
       "      <td>-1.406107</td>\n",
       "      <td>0.847746</td>\n",
       "      <td>-1.658641</td>\n",
       "      <td>9.718815e-02</td>\n",
       "      <td>3.887526e-01</td>\n",
       "      <td>False</td>\n",
       "      <td>-0.380518</td>\n",
       "      <td>-3.067660</td>\n",
       "      <td>0.255446</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>High minus Low Magnitude</th>\n",
       "      <td>IFS</td>\n",
       "      <td>0.306222</td>\n",
       "      <td>1.025211</td>\n",
       "      <td>0.298691</td>\n",
       "      <td>7.651756e-01</td>\n",
       "      <td>7.994840e-01</td>\n",
       "      <td>False</td>\n",
       "      <td>0.068524</td>\n",
       "      <td>-1.703155</td>\n",
       "      <td>2.315598</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>High minus Low Magnitude</th>\n",
       "      <td>pMFG</td>\n",
       "      <td>-1.233085</td>\n",
       "      <td>0.988291</td>\n",
       "      <td>-1.247695</td>\n",
       "      <td>2.121429e-01</td>\n",
       "      <td>5.657144e-01</td>\n",
       "      <td>False</td>\n",
       "      <td>-0.286241</td>\n",
       "      <td>-3.170099</td>\n",
       "      <td>0.703929</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>High minus Low Magnitude</th>\n",
       "      <td>vmPFC</td>\n",
       "      <td>-1.035973</td>\n",
       "      <td>1.841250</td>\n",
       "      <td>-0.562646</td>\n",
       "      <td>5.736757e-01</td>\n",
       "      <td>7.649009e-01</td>\n",
       "      <td>False</td>\n",
       "      <td>-0.129080</td>\n",
       "      <td>-4.644757</td>\n",
       "      <td>2.572812</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>High minus Low Magnitude</th>\n",
       "      <td>FEF</td>\n",
       "      <td>-0.659646</td>\n",
       "      <td>0.818984</td>\n",
       "      <td>-0.805445</td>\n",
       "      <td>4.205629e-01</td>\n",
       "      <td>6.729006e-01</td>\n",
       "      <td>False</td>\n",
       "      <td>-0.184782</td>\n",
       "      <td>-2.264825</td>\n",
       "      <td>0.945532</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Difficulty (High minus Low Magnitude)</th>\n",
       "      <td>IFG</td>\n",
       "      <td>1.407193</td>\n",
       "      <td>0.539457</td>\n",
       "      <td>2.608535</td>\n",
       "      <td>9.093070e-03</td>\n",
       "      <td>3.956508e-02</td>\n",
       "      <td>True</td>\n",
       "      <td>0.598439</td>\n",
       "      <td>0.349876</td>\n",
       "      <td>2.464510</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Difficulty (High minus Low Magnitude)</th>\n",
       "      <td>FPC</td>\n",
       "      <td>0.371296</td>\n",
       "      <td>0.626634</td>\n",
       "      <td>0.592524</td>\n",
       "      <td>5.535000e-01</td>\n",
       "      <td>6.143256e-01</td>\n",
       "      <td>False</td>\n",
       "      <td>0.135934</td>\n",
       "      <td>-0.856885</td>\n",
       "      <td>1.599477</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Difficulty (High minus Low Magnitude)</th>\n",
       "      <td>aMFG</td>\n",
       "      <td>0.251362</td>\n",
       "      <td>0.498825</td>\n",
       "      <td>0.503909</td>\n",
       "      <td>6.143256e-01</td>\n",
       "      <td>6.143256e-01</td>\n",
       "      <td>False</td>\n",
       "      <td>0.115605</td>\n",
       "      <td>-0.726317</td>\n",
       "      <td>1.229042</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Difficulty (High minus Low Magnitude)</th>\n",
       "      <td>aINS</td>\n",
       "      <td>0.717854</td>\n",
       "      <td>0.515254</td>\n",
       "      <td>1.393206</td>\n",
       "      <td>1.635576e-01</td>\n",
       "      <td>3.044696e-01</td>\n",
       "      <td>False</td>\n",
       "      <td>0.319623</td>\n",
       "      <td>-0.292024</td>\n",
       "      <td>1.727733</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Difficulty (High minus Low Magnitude)</th>\n",
       "      <td>IFS</td>\n",
       "      <td>1.131639</td>\n",
       "      <td>0.438687</td>\n",
       "      <td>2.579607</td>\n",
       "      <td>9.891270e-03</td>\n",
       "      <td>3.956508e-02</td>\n",
       "      <td>True</td>\n",
       "      <td>0.591803</td>\n",
       "      <td>0.271829</td>\n",
       "      <td>1.991450</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Difficulty (High minus Low Magnitude)</th>\n",
       "      <td>pMFG</td>\n",
       "      <td>1.147320</td>\n",
       "      <td>0.526411</td>\n",
       "      <td>2.179511</td>\n",
       "      <td>2.929371e-02</td>\n",
       "      <td>7.811655e-02</td>\n",
       "      <td>False</td>\n",
       "      <td>0.500014</td>\n",
       "      <td>0.115572</td>\n",
       "      <td>2.179067</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Difficulty (High minus Low Magnitude)</th>\n",
       "      <td>vmPFC</td>\n",
       "      <td>-0.786301</td>\n",
       "      <td>0.600362</td>\n",
       "      <td>-1.309711</td>\n",
       "      <td>1.902935e-01</td>\n",
       "      <td>3.044696e-01</td>\n",
       "      <td>False</td>\n",
       "      <td>-0.300468</td>\n",
       "      <td>-1.962988</td>\n",
       "      <td>0.390387</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Difficulty (High minus Low Magnitude)</th>\n",
       "      <td>FEF</td>\n",
       "      <td>0.338785</td>\n",
       "      <td>0.315980</td>\n",
       "      <td>1.072172</td>\n",
       "      <td>2.836426e-01</td>\n",
       "      <td>3.781902e-01</td>\n",
       "      <td>False</td>\n",
       "      <td>0.245973</td>\n",
       "      <td>-0.280524</td>\n",
       "      <td>0.958094</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Subjective Value (High minus Low Magnitude)</th>\n",
       "      <td>IFG</td>\n",
       "      <td>0.508003</td>\n",
       "      <td>0.623231</td>\n",
       "      <td>0.815112</td>\n",
       "      <td>4.150083e-01</td>\n",
       "      <td>7.460049e-01</td>\n",
       "      <td>False</td>\n",
       "      <td>0.186999</td>\n",
       "      <td>-0.713507</td>\n",
       "      <td>1.729514</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Subjective Value (High minus Low Magnitude)</th>\n",
       "      <td>FPC</td>\n",
       "      <td>-0.250288</td>\n",
       "      <td>0.458286</td>\n",
       "      <td>-0.546140</td>\n",
       "      <td>5.849694e-01</td>\n",
       "      <td>7.460049e-01</td>\n",
       "      <td>False</td>\n",
       "      <td>-0.125293</td>\n",
       "      <td>-1.148512</td>\n",
       "      <td>0.647935</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Subjective Value (High minus Low Magnitude)</th>\n",
       "      <td>aMFG</td>\n",
       "      <td>0.280754</td>\n",
       "      <td>0.461147</td>\n",
       "      <td>0.608817</td>\n",
       "      <td>5.426458e-01</td>\n",
       "      <td>7.460049e-01</td>\n",
       "      <td>False</td>\n",
       "      <td>0.139672</td>\n",
       "      <td>-0.623078</td>\n",
       "      <td>1.184586</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Subjective Value (High minus Low Magnitude)</th>\n",
       "      <td>aINS</td>\n",
       "      <td>0.453992</td>\n",
       "      <td>0.532264</td>\n",
       "      <td>0.852944</td>\n",
       "      <td>3.936905e-01</td>\n",
       "      <td>7.460049e-01</td>\n",
       "      <td>False</td>\n",
       "      <td>0.195679</td>\n",
       "      <td>-0.589227</td>\n",
       "      <td>1.497211</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Subjective Value (High minus Low Magnitude)</th>\n",
       "      <td>IFS</td>\n",
       "      <td>0.272213</td>\n",
       "      <td>0.605001</td>\n",
       "      <td>0.449939</td>\n",
       "      <td>6.527543e-01</td>\n",
       "      <td>7.460049e-01</td>\n",
       "      <td>False</td>\n",
       "      <td>0.103223</td>\n",
       "      <td>-0.913566</td>\n",
       "      <td>1.457993</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Subjective Value (High minus Low Magnitude)</th>\n",
       "      <td>pMFG</td>\n",
       "      <td>-0.447035</td>\n",
       "      <td>0.547123</td>\n",
       "      <td>-0.817065</td>\n",
       "      <td>4.138914e-01</td>\n",
       "      <td>7.460049e-01</td>\n",
       "      <td>False</td>\n",
       "      <td>-0.187448</td>\n",
       "      <td>-1.519375</td>\n",
       "      <td>0.625306</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Subjective Value (High minus Low Magnitude)</th>\n",
       "      <td>vmPFC</td>\n",
       "      <td>0.037648</td>\n",
       "      <td>0.715305</td>\n",
       "      <td>0.052632</td>\n",
       "      <td>9.580254e-01</td>\n",
       "      <td>9.580254e-01</td>\n",
       "      <td>False</td>\n",
       "      <td>0.012075</td>\n",
       "      <td>-1.364323</td>\n",
       "      <td>1.439619</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Subjective Value (High minus Low Magnitude)</th>\n",
       "      <td>FEF</td>\n",
       "      <td>0.900952</td>\n",
       "      <td>0.319484</td>\n",
       "      <td>2.820019</td>\n",
       "      <td>4.802083e-03</td>\n",
       "      <td>3.841667e-02</td>\n",
       "      <td>True</td>\n",
       "      <td>0.646957</td>\n",
       "      <td>0.274774</td>\n",
       "      <td>1.527130</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "                                               ROI        val   std_err  \\\n",
       "contrast                                                                  \n",
       "Subjective Value                               IFG  -0.227902  0.343442   \n",
       "Subjective Value                               FPC  -0.075737  0.383605   \n",
       "Subjective Value                              aMFG  -0.295708  0.364751   \n",
       "Subjective Value                              aINS  -0.077649  0.320861   \n",
       "Subjective Value                               IFS   0.877997  0.425268   \n",
       "Subjective Value                              pMFG  -0.392609  0.355032   \n",
       "Subjective Value                             vmPFC   1.727025  0.551990   \n",
       "Subjective Value                               FEF  -0.361484  0.572772   \n",
       "Decision Period                                IFG  -8.018357  1.609433   \n",
       "Decision Period                                FPC  -4.997577  2.017170   \n",
       "Decision Period                               aMFG  -3.004705  1.558327   \n",
       "Decision Period                               aINS   9.606990  1.621377   \n",
       "Decision Period                                IFS  10.328838  1.913305   \n",
       "Decision Period                               pMFG   1.606482  1.820506   \n",
       "Decision Period                              vmPFC -33.996927  2.624094   \n",
       "Decision Period                                FEF  -2.140544  2.099276   \n",
       "Difficulty                                     IFG  -0.231936  0.653274   \n",
       "Difficulty                                     FPC   1.912566  0.710861   \n",
       "Difficulty                                    aMFG   0.745508  0.447358   \n",
       "Difficulty                                    aINS   1.154436  0.356968   \n",
       "Difficulty                                     IFS   2.450334  0.481019   \n",
       "Difficulty                                    pMFG   2.078334  0.632626   \n",
       "Difficulty                                   vmPFC  -1.919153  0.724114   \n",
       "Difficulty                                     FEF   0.390121  0.623569   \n",
       "High minus Low Magnitude                       IFG  -4.058899  0.753600   \n",
       "High minus Low Magnitude                       FPC   0.258278  1.016781   \n",
       "High minus Low Magnitude                      aMFG  -0.804398  0.778557   \n",
       "High minus Low Magnitude                      aINS  -1.406107  0.847746   \n",
       "High minus Low Magnitude                       IFS   0.306222  1.025211   \n",
       "High minus Low Magnitude                      pMFG  -1.233085  0.988291   \n",
       "High minus Low Magnitude                     vmPFC  -1.035973  1.841250   \n",
       "High minus Low Magnitude                       FEF  -0.659646  0.818984   \n",
       "Difficulty (High minus Low Magnitude)          IFG   1.407193  0.539457   \n",
       "Difficulty (High minus Low Magnitude)          FPC   0.371296  0.626634   \n",
       "Difficulty (High minus Low Magnitude)         aMFG   0.251362  0.498825   \n",
       "Difficulty (High minus Low Magnitude)         aINS   0.717854  0.515254   \n",
       "Difficulty (High minus Low Magnitude)          IFS   1.131639  0.438687   \n",
       "Difficulty (High minus Low Magnitude)         pMFG   1.147320  0.526411   \n",
       "Difficulty (High minus Low Magnitude)        vmPFC  -0.786301  0.600362   \n",
       "Difficulty (High minus Low Magnitude)          FEF   0.338785  0.315980   \n",
       "Subjective Value (High minus Low Magnitude)    IFG   0.508003  0.623231   \n",
       "Subjective Value (High minus Low Magnitude)    FPC  -0.250288  0.458286   \n",
       "Subjective Value (High minus Low Magnitude)   aMFG   0.280754  0.461147   \n",
       "Subjective Value (High minus Low Magnitude)   aINS   0.453992  0.532264   \n",
       "Subjective Value (High minus Low Magnitude)    IFS   0.272213  0.605001   \n",
       "Subjective Value (High minus Low Magnitude)   pMFG  -0.447035  0.547123   \n",
       "Subjective Value (High minus Low Magnitude)  vmPFC   0.037648  0.715305   \n",
       "Subjective Value (High minus Low Magnitude)    FEF   0.900952  0.319484   \n",
       "\n",
       "                                                     t         p_unc  \\\n",
       "contrast                                                               \n",
       "Subjective Value                             -0.663584  5.069563e-01   \n",
       "Subjective Value                             -0.197434  8.434879e-01   \n",
       "Subjective Value                             -0.810714  4.175301e-01   \n",
       "Subjective Value                             -0.242003  8.087778e-01   \n",
       "Subjective Value                              2.064572  3.896351e-02   \n",
       "Subjective Value                             -1.105841  2.687953e-01   \n",
       "Subjective Value                              3.128725  1.755664e-03   \n",
       "Subjective Value                             -0.631113  5.279664e-01   \n",
       "Decision Period                              -4.982100  6.289789e-07   \n",
       "Decision Period                              -2.477519  1.322992e-02   \n",
       "Decision Period                              -1.928161  5.383516e-02   \n",
       "Decision Period                               5.925203  3.119109e-09   \n",
       "Decision Period                               5.398427  6.722772e-08   \n",
       "Decision Period                               0.882437  3.775405e-01   \n",
       "Decision Period                             -12.955682  2.181896e-38   \n",
       "Decision Period                              -1.019658  3.078905e-01   \n",
       "Difficulty                                   -0.355036  7.225626e-01   \n",
       "Difficulty                                    2.690492  7.134664e-03   \n",
       "Difficulty                                    1.666469  9.562007e-02   \n",
       "Difficulty                                    3.234007  1.220667e-03   \n",
       "Difficulty                                    5.094045  3.505033e-07   \n",
       "Difficulty                                    3.285247  1.018930e-03   \n",
       "Difficulty                                   -2.650347  8.040915e-03   \n",
       "Difficulty                                    0.625627  5.315599e-01   \n",
       "High minus Low Magnitude                     -5.386009  7.203916e-08   \n",
       "High minus Low Magnitude                      0.254015  7.994840e-01   \n",
       "High minus Low Magnitude                     -1.033190  3.015150e-01   \n",
       "High minus Low Magnitude                     -1.658641  9.718815e-02   \n",
       "High minus Low Magnitude                      0.298691  7.651756e-01   \n",
       "High minus Low Magnitude                     -1.247695  2.121429e-01   \n",
       "High minus Low Magnitude                     -0.562646  5.736757e-01   \n",
       "High minus Low Magnitude                     -0.805445  4.205629e-01   \n",
       "Difficulty (High minus Low Magnitude)         2.608535  9.093070e-03   \n",
       "Difficulty (High minus Low Magnitude)         0.592524  5.535000e-01   \n",
       "Difficulty (High minus Low Magnitude)         0.503909  6.143256e-01   \n",
       "Difficulty (High minus Low Magnitude)         1.393206  1.635576e-01   \n",
       "Difficulty (High minus Low Magnitude)         2.579607  9.891270e-03   \n",
       "Difficulty (High minus Low Magnitude)         2.179511  2.929371e-02   \n",
       "Difficulty (High minus Low Magnitude)        -1.309711  1.902935e-01   \n",
       "Difficulty (High minus Low Magnitude)         1.072172  2.836426e-01   \n",
       "Subjective Value (High minus Low Magnitude)   0.815112  4.150083e-01   \n",
       "Subjective Value (High minus Low Magnitude)  -0.546140  5.849694e-01   \n",
       "Subjective Value (High minus Low Magnitude)   0.608817  5.426458e-01   \n",
       "Subjective Value (High minus Low Magnitude)   0.852944  3.936905e-01   \n",
       "Subjective Value (High minus Low Magnitude)   0.449939  6.527543e-01   \n",
       "Subjective Value (High minus Low Magnitude)  -0.817065  4.138914e-01   \n",
       "Subjective Value (High minus Low Magnitude)   0.052632  9.580254e-01   \n",
       "Subjective Value (High minus Low Magnitude)   2.820019  4.802083e-03   \n",
       "\n",
       "                                                   p_corr    sig        dz  \\\n",
       "contrast                                                                     \n",
       "Subjective Value                             7.039552e-01  False -0.152237   \n",
       "Subjective Value                             8.434879e-01  False -0.045294   \n",
       "Subjective Value                             7.039552e-01  False -0.185990   \n",
       "Subjective Value                             8.434879e-01  False -0.055519   \n",
       "Subjective Value                             1.558540e-01  False  0.473645   \n",
       "Subjective Value                             7.039552e-01  False -0.253697   \n",
       "Subjective Value                             1.404531e-02   True  0.717779   \n",
       "Subjective Value                             7.039552e-01  False -0.144787   \n",
       "Decision Period                              1.257958e-06   True -1.142972   \n",
       "Decision Period                              2.116787e-02   True -0.568382   \n",
       "Decision Period                              7.178021e-02  False -0.442350   \n",
       "Decision Period                              1.247644e-08   True  1.359335   \n",
       "Decision Period                              1.792739e-07   True  1.238484   \n",
       "Decision Period                              3.775405e-01  False  0.202445   \n",
       "Decision Period                              1.745517e-37   True -2.972237   \n",
       "Decision Period                              3.518749e-01  False -0.233926   \n",
       "Difficulty                                   7.225626e-01  False -0.081451   \n",
       "Difficulty                                   1.286546e-02   True  0.617241   \n",
       "Difficulty                                   1.274934e-01  False  0.382314   \n",
       "Difficulty                                   3.255111e-03   True  0.741932   \n",
       "Difficulty                                   2.804027e-06   True  1.168654   \n",
       "Difficulty                                   3.255111e-03   True  0.753687   \n",
       "Difficulty                                   1.286546e-02   True -0.608031   \n",
       "Difficulty                                   6.074971e-01  False  0.143529   \n",
       "High minus Low Magnitude                     5.763133e-07   True -1.235635   \n",
       "High minus Low Magnitude                     7.994840e-01  False  0.058275   \n",
       "High minus Low Magnitude                     6.030300e-01  False -0.237030   \n",
       "High minus Low Magnitude                     3.887526e-01  False -0.380518   \n",
       "High minus Low Magnitude                     7.994840e-01  False  0.068524   \n",
       "High minus Low Magnitude                     5.657144e-01  False -0.286241   \n",
       "High minus Low Magnitude                     7.649009e-01  False -0.129080   \n",
       "High minus Low Magnitude                     6.729006e-01  False -0.184782   \n",
       "Difficulty (High minus Low Magnitude)        3.956508e-02   True  0.598439   \n",
       "Difficulty (High minus Low Magnitude)        6.143256e-01  False  0.135934   \n",
       "Difficulty (High minus Low Magnitude)        6.143256e-01  False  0.115605   \n",
       "Difficulty (High minus Low Magnitude)        3.044696e-01  False  0.319623   \n",
       "Difficulty (High minus Low Magnitude)        3.956508e-02   True  0.591803   \n",
       "Difficulty (High minus Low Magnitude)        7.811655e-02  False  0.500014   \n",
       "Difficulty (High minus Low Magnitude)        3.044696e-01  False -0.300468   \n",
       "Difficulty (High minus Low Magnitude)        3.781902e-01  False  0.245973   \n",
       "Subjective Value (High minus Low Magnitude)  7.460049e-01  False  0.186999   \n",
       "Subjective Value (High minus Low Magnitude)  7.460049e-01  False -0.125293   \n",
       "Subjective Value (High minus Low Magnitude)  7.460049e-01  False  0.139672   \n",
       "Subjective Value (High minus Low Magnitude)  7.460049e-01  False  0.195679   \n",
       "Subjective Value (High minus Low Magnitude)  7.460049e-01  False  0.103223   \n",
       "Subjective Value (High minus Low Magnitude)  7.460049e-01  False -0.187448   \n",
       "Subjective Value (High minus Low Magnitude)  9.580254e-01  False  0.012075   \n",
       "Subjective Value (High minus Low Magnitude)  3.841667e-02   True  0.646957   \n",
       "\n",
       "                                             95_CI_low  95_CI_high  \n",
       "contrast                                                            \n",
       "Subjective Value                             -0.901036    0.445231  \n",
       "Subjective Value                             -0.827588    0.676115  \n",
       "Subjective Value                             -1.010606    0.419190  \n",
       "Subjective Value                             -0.706525    0.551226  \n",
       "Subjective Value                              0.044487    1.711508  \n",
       "Subjective Value                             -1.088458    0.303241  \n",
       "Subjective Value                              0.645144    2.808905  \n",
       "Subjective Value                             -1.484097    0.761129  \n",
       "Decision Period                             -11.172789   -4.863926  \n",
       "Decision Period                              -8.951157   -1.043997  \n",
       "Decision Period                              -6.058971    0.049560  \n",
       "Decision Period                               6.429149   12.784831  \n",
       "Decision Period                               6.578829   14.078847  \n",
       "Decision Period                              -1.961644    5.174607  \n",
       "Decision Period                             -39.140056  -28.853797  \n",
       "Decision Period                              -6.255050    1.973961  \n",
       "Difficulty                                   -1.512329    1.048458  \n",
       "Difficulty                                    0.519304    3.305828  \n",
       "Difficulty                                   -0.131297    1.622313  \n",
       "Difficulty                                    0.454792    1.854081  \n",
       "Difficulty                                    1.507554    3.393115  \n",
       "Difficulty                                    0.838409    3.318259  \n",
       "Difficulty                                   -3.338390   -0.499916  \n",
       "Difficulty                                   -0.832051    1.612293  \n",
       "High minus Low Magnitude                     -5.535928   -2.581869  \n",
       "High minus Low Magnitude                     -1.734576    2.251131  \n",
       "High minus Low Magnitude                     -2.330342    0.721547  \n",
       "High minus Low Magnitude                     -3.067660    0.255446  \n",
       "High minus Low Magnitude                     -1.703155    2.315598  \n",
       "High minus Low Magnitude                     -3.170099    0.703929  \n",
       "High minus Low Magnitude                     -4.644757    2.572812  \n",
       "High minus Low Magnitude                     -2.264825    0.945532  \n",
       "Difficulty (High minus Low Magnitude)         0.349876    2.464510  \n",
       "Difficulty (High minus Low Magnitude)        -0.856885    1.599477  \n",
       "Difficulty (High minus Low Magnitude)        -0.726317    1.229042  \n",
       "Difficulty (High minus Low Magnitude)        -0.292024    1.727733  \n",
       "Difficulty (High minus Low Magnitude)         0.271829    1.991450  \n",
       "Difficulty (High minus Low Magnitude)         0.115572    2.179067  \n",
       "Difficulty (High minus Low Magnitude)        -1.962988    0.390387  \n",
       "Difficulty (High minus Low Magnitude)        -0.280524    0.958094  \n",
       "Subjective Value (High minus Low Magnitude)  -0.713507    1.729514  \n",
       "Subjective Value (High minus Low Magnitude)  -1.148512    0.647935  \n",
       "Subjective Value (High minus Low Magnitude)  -0.623078    1.184586  \n",
       "Subjective Value (High minus Low Magnitude)  -0.589227    1.497211  \n",
       "Subjective Value (High minus Low Magnitude)  -0.913566    1.457993  \n",
       "Subjective Value (High minus Low Magnitude)  -1.519375    0.625306  \n",
       "Subjective Value (High minus Low Magnitude)  -1.364323    1.439619  \n",
       "Subjective Value (High minus Low Magnitude)   0.274774    1.527130  "
      ]
     },
     "execution_count": 36,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "#compute 95% confidence interval\n",
    "all_results['95_CI_high'] = all_results['val'] + all_results['err']\n",
    "all_results['95_CI_low'] = all_results['val'] - all_results['err']\n",
    "\n",
    "all_results = all_results[['contrast','ROI','val','std_err','t','p_unc','p_corr','sig','dz','95_CI_low','95_CI_high']]\n",
    "all_results.to_csv('/home/iballard/Dropbox/mag_figs/roi_results.csv',\n",
    "                   float_format='%.3f',\n",
    "                   index = False)\n",
    "\n",
    "all_results = all_results.set_index('contrast')\n",
    "all_results"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "collapsed": false
   },
   "source": [
    "# Run brain-behavior correlation"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 138,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "#load data\n",
    "params = pd.read_csv('/home/iballard/Dropbox/mag_figs/model_fits.csv')\n",
    "params = params[params['sub'] != 'asu09'] #subject excluded for strong anti-magnitude effect\n",
    "params = params.reset_index()\n",
    "\n",
    "#fix subject id coding\n",
    "for n,s in enumerate(params['sub']):\n",
    "    if s == 'asu04':\n",
    "        params.loc[n,'sub'] = 'asu4'\n",
    "    if s == 'asu08':\n",
    "        params.loc[n,'sub'] = 'asu8'"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 139,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "params['mag_effect'] = [(np.log(low) - np.log(high)) for low,high in zip(params['k_low'],params['k_high'])]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 190,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "contrast = 'low_minus_high_SV_diff'\n",
    "ROI = 'IFS'"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 191,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "#extract data\n",
    "df = extract_copes(exp,regspace,smoothing,contrast)\n",
    "df['ROI'] = [label_map[x] for x in df['label'].values]\n",
    "df.set_index('ROI').loc[ROI]\n",
    "df = df.groupby('sub').mean().reset_index()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 192,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "df = pd.merge(df,\n",
    "         params,\n",
    "        how = 'left',\n",
    "        on = 'sub')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 193,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "                    Robust linear Model Regression Results                    \n",
      "==============================================================================\n",
      "Dep. Variable:                  value   No. Observations:                   19\n",
      "Model:                            RLM   Df Residuals:                       17\n",
      "Method:                          IRLS   Df Model:                            1\n",
      "Norm:                          HuberT                                         \n",
      "Scale Est.:                       mad                                         \n",
      "Cov Type:                          H1                                         \n",
      "Date:                Wed, 05 Apr 2017                                         \n",
      "Time:                        14:38:53                                         \n",
      "No. Iterations:                     9                                         \n",
      "==============================================================================\n",
      "                 coef    std err          z      P>|z|      [95.0% Conf. Int.]\n",
      "------------------------------------------------------------------------------\n",
      "const         -0.6385      0.934     -0.683      0.494        -2.470     1.193\n",
      "mag_effect     0.3817      0.285      1.340      0.180        -0.177     0.940\n",
      "==============================================================================\n",
      "\n",
      "If the model instance has been used for another fit with different fit\n",
      "parameters, then the fit options might not be the correct ones anymore .\n"
     ]
    }
   ],
   "source": [
    "exog = statsmodels.tools.add_constant(df['mag_effect'])\n",
    "huber_t = smf.RLM(df['value'], exog, M=sm.robust.norms.HuberT())\n",
    "hub_results = huber_t.fit()\n",
    "print hub_results.summary()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 194,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "                            OLS Regression Results                            \n",
      "==============================================================================\n",
      "Dep. Variable:                  value   R-squared:                       0.108\n",
      "Model:                            OLS   Adj. R-squared:                  0.056\n",
      "Method:                 Least Squares   F-statistic:                     2.058\n",
      "Date:                Wed, 05 Apr 2017   Prob (F-statistic):              0.170\n",
      "Time:                        14:38:53   Log-Likelihood:                -28.944\n",
      "No. Observations:                  19   AIC:                             61.89\n",
      "Df Residuals:                      17   BIC:                             63.78\n",
      "Df Model:                           1                                         \n",
      "Covariance Type:            nonrobust                                         \n",
      "==============================================================================\n",
      "                 coef    std err          t      P>|t|      [95.0% Conf. Int.]\n",
      "------------------------------------------------------------------------------\n",
      "const         -0.6285      0.879     -0.715      0.484        -2.484     1.227\n",
      "mag_effect     0.3847      0.268      1.435      0.170        -0.181     0.950\n",
      "==============================================================================\n",
      "Omnibus:                        0.072   Durbin-Watson:                   1.715\n",
      "Prob(Omnibus):                  0.965   Jarque-Bera (JB):                0.162\n",
      "Skew:                           0.113   Prob(JB):                        0.922\n",
      "Kurtosis:                       2.608   Cond. No.                         11.6\n",
      "==============================================================================\n",
      "\n",
      "Warnings:\n",
      "[1] Standard Errors assume that the covariance matrix of the errors is correctly specified.\n"
     ]
    }
   ],
   "source": [
    "res = smf.OLS(df['value'], exog)\n",
    "res = res.fit()\n",
    "print res.summary()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 2",
   "language": "python",
   "name": "python2"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 2
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython2",
   "version": "2.7.11"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 0
}