Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- {
- "cells": [
- {
- "cell_type": "code",
- "execution_count": 8,
- "metadata": {
- "collapsed": false
- },
- "outputs": [
- {
- "name": "stdout",
- "output_type": "stream",
- "text": [
- "Populating the interactive namespace from numpy and matplotlib\n"
- ]
- }
- ],
- "source": [
- "%pylab inline\n",
- "import seaborn as sns\n",
- "import pandas as pd\n",
- "import statsmodels.formula.api as sf"
- ]
- },
- {
- "cell_type": "code",
- "execution_count": 3,
- "metadata": {
- "collapsed": true
- },
- "outputs": [],
- "source": [
- "df = pd.read_csv(\"combined_demographics_imaging.csv\",index_col=0)"
- ]
- },
- {
- "cell_type": "code",
- "execution_count": 28,
- "metadata": {
- "collapsed": false
- },
- "outputs": [
- {
- "name": "stdout",
- "output_type": "stream",
- "text": [
- "Index([u'connection', u'r00_correlation', u'r00_partial correlation', u'r00_tangent', u'r00_fa', u'r00_md', u'r00_ad', u'r00_rd', u'r00_mt_map', u'r00_mtr_map', u'r00_r1_map', u'mse', u'MS Group Subject ID', u'MS Group Exam ID', u'type', u'DCM_SeriesDescription', u'DCM_StudyDate', u'mseid', u'msid', u'age_at_visit', u'edss', u'date', u'disease_course', u'timed_walk_trial1_time', u'nhpt_dominant_time', u'fssc_bowel', u'fssc_brainstem', u'fssc_cerebellar', u'fssc_mental', u'fssc_pyramidal', u'fssc_sensory', u'fssc_visual', u'epicid', u'gender', u'msgb', u'handedness', u'smoker_status', u'msgb_non_mhc_snps', u'year_converted_rr_to_sp', u'date_reformat'], dtype='object')\n"
- ]
- }
- ],
- "source": [
- "print df.columns"
- ]
- },
- {
- "cell_type": "code",
- "execution_count": 45,
- "metadata": {
- "collapsed": false
- },
- "outputs": [],
- "source": [
- "results = pd.DataFrame()\n",
- "for conn, x in df.groupby(\"connection\"):\n",
- " model = sf.ols(\"edss ~ r00_correlation + r00_mtr_map +r00_md+age_at_visit +gender\", x).fit()\n",
- " foo = model.pvalues\n",
- " foo[\"r2\"] = model.rsquared\n",
- " foo[\"measure\"] = conn\n",
- " results = results.append(foo,ignore_index=True)\n"
- ]
- },
- {
- "cell_type": "code",
- "execution_count": 46,
- "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>Intercept</th>\n",
- " <th>age_at_visit</th>\n",
- " <th>gender[T.M]</th>\n",
- " <th>measure</th>\n",
- " <th>r00_correlation</th>\n",
- " <th>r00_md</th>\n",
- " <th>r00_mtr_map</th>\n",
- " <th>r2</th>\n",
- " </tr>\n",
- " </thead>\n",
- " <tbody>\n",
- " <tr>\n",
- " <th>7 </th>\n",
- " <td> 0.139689</td>\n",
- " <td> 0.000889</td>\n",
- " <td> 0.784138</td>\n",
- " <td> L_S-M1_L_PFC</td>\n",
- " <td> 0.000316</td>\n",
- " <td> 0.831087</td>\n",
- " <td> 0.033602</td>\n",
- " <td> 0.238137</td>\n",
- " </tr>\n",
- " <tr>\n",
- " <th>18</th>\n",
- " <td> 0.503964</td>\n",
- " <td> 0.003005</td>\n",
- " <td> 0.725810</td>\n",
- " <td> R_S-M1_R_PFC</td>\n",
- " <td> 0.000253</td>\n",
- " <td> 0.162028</td>\n",
- " <td> 0.755075</td>\n",
- " <td> 0.216037</td>\n",
- " </tr>\n",
- " <tr>\n",
- " <th>4 </th>\n",
- " <td> 0.723769</td>\n",
- " <td> 0.000205</td>\n",
- " <td> 0.762192</td>\n",
- " <td> L_AS-Sens-C_R_deepgray</td>\n",
- " <td> 0.024560</td>\n",
- " <td> 0.304218</td>\n",
- " <td> 0.259391</td>\n",
- " <td> 0.206945</td>\n",
- " </tr>\n",
- " <tr>\n",
- " <th>19</th>\n",
- " <td> 0.181365</td>\n",
- " <td> 0.002724</td>\n",
- " <td> 0.699412</td>\n",
- " <td> R_S-M1_R_S2</td>\n",
- " <td> 0.007781</td>\n",
- " <td> 0.016981</td>\n",
- " <td> 0.534828</td>\n",
- " <td> 0.203086</td>\n",
- " </tr>\n",
- " <tr>\n",
- " <th>15</th>\n",
- " <td> 0.128126</td>\n",
- " <td> 0.002279</td>\n",
- " <td> 0.772100</td>\n",
- " <td> R_PFC_L_deepgray</td>\n",
- " <td> 0.004811</td>\n",
- " <td> 0.029533</td>\n",
- " <td> 0.320453</td>\n",
- " <td> 0.200328</td>\n",
- " </tr>\n",
- " <tr>\n",
- " <th>14</th>\n",
- " <td> 0.905604</td>\n",
- " <td> 0.000916</td>\n",
- " <td> 0.991047</td>\n",
- " <td> R_AS-Sens-C_R_S2</td>\n",
- " <td> 0.058617</td>\n",
- " <td> 0.070504</td>\n",
- " <td> 0.500667</td>\n",
- " <td> 0.199615</td>\n",
- " </tr>\n",
- " <tr>\n",
- " <th>20</th>\n",
- " <td> 0.538895</td>\n",
- " <td> 0.000267</td>\n",
- " <td> 0.961521</td>\n",
- " <td> R_S-M1_R_deepgray</td>\n",
- " <td> 0.013669</td>\n",
- " <td> 0.088691</td>\n",
- " <td> 0.918744</td>\n",
- " <td> 0.197836</td>\n",
- " </tr>\n",
- " <tr>\n",
- " <th>13</th>\n",
- " <td> 0.484858</td>\n",
- " <td> 0.000346</td>\n",
- " <td> 0.837863</td>\n",
- " <td> R_AS-Sens-C_L_deepgray</td>\n",
- " <td> 0.009018</td>\n",
- " <td> 0.627967</td>\n",
- " <td> 0.154479</td>\n",
- " <td> 0.193647</td>\n",
- " </tr>\n",
- " <tr>\n",
- " <th>6 </th>\n",
- " <td> 0.440572</td>\n",
- " <td> 0.001456</td>\n",
- " <td> 0.843126</td>\n",
- " <td> L_PFC_R_PFC</td>\n",
- " <td> 0.008558</td>\n",
- " <td> 0.139267</td>\n",
- " <td> 0.571633</td>\n",
- " <td> 0.187000</td>\n",
- " </tr>\n",
- " <tr>\n",
- " <th>3 </th>\n",
- " <td> 0.595859</td>\n",
- " <td> 0.000703</td>\n",
- " <td> 0.898476</td>\n",
- " <td> L_AS-Sens-C_R_S2</td>\n",
- " <td> 0.028592</td>\n",
- " <td> 0.570917</td>\n",
- " <td> 0.191679</td>\n",
- " <td> 0.186707</td>\n",
- " </tr>\n",
- " <tr>\n",
- " <th>8 </th>\n",
- " <td> 0.479252</td>\n",
- " <td> 0.000214</td>\n",
- " <td> 0.831458</td>\n",
- " <td> L_S-M1_L_deepgray</td>\n",
- " <td> 0.047262</td>\n",
- " <td> 0.461983</td>\n",
- " <td> 0.142359</td>\n",
- " <td> 0.184781</td>\n",
- " </tr>\n",
- " <tr>\n",
- " <th>23</th>\n",
- " <td> 0.447515</td>\n",
- " <td> 0.000180</td>\n",
- " <td> 0.744982</td>\n",
- " <td> R_deepgray_R_PFC</td>\n",
- " <td> 0.006440</td>\n",
- " <td> 0.229144</td>\n",
- " <td> 0.728392</td>\n",
- " <td> 0.184582</td>\n",
- " </tr>\n",
- " <tr>\n",
- " <th>0 </th>\n",
- " <td> 0.630970</td>\n",
- " <td> 0.000807</td>\n",
- " <td> 0.946772</td>\n",
- " <td> L_AS-Sens-C_L_S-M1</td>\n",
- " <td> 0.184438</td>\n",
- " <td> 0.201691</td>\n",
- " <td> 0.095884</td>\n",
- " <td> 0.181971</td>\n",
- " </tr>\n",
- " <tr>\n",
- " <th>22</th>\n",
- " <td> 0.648730</td>\n",
- " <td> 0.000526</td>\n",
- " <td> 0.688840</td>\n",
- " <td> R_deepgray_R_AS-Sens-C</td>\n",
- " <td> 0.010884</td>\n",
- " <td> 0.279351</td>\n",
- " <td> 0.991240</td>\n",
- " <td> 0.181893</td>\n",
- " </tr>\n",
- " <tr>\n",
- " <th>5 </th>\n",
- " <td> 0.490455</td>\n",
- " <td> 0.000912</td>\n",
- " <td> 0.552587</td>\n",
- " <td> L_PFC_L_deepgray</td>\n",
- " <td> 0.043813</td>\n",
- " <td> 0.374006</td>\n",
- " <td> 0.187795</td>\n",
- " <td> 0.181429</td>\n",
- " </tr>\n",
- " <tr>\n",
- " <th>1 </th>\n",
- " <td> 0.461557</td>\n",
- " <td> 0.000675</td>\n",
- " <td> 0.722178</td>\n",
- " <td> L_AS-Sens-C_L_deepgray</td>\n",
- " <td> 0.009938</td>\n",
- " <td> 0.141772</td>\n",
- " <td> 0.828873</td>\n",
- " <td> 0.180617</td>\n",
- " </tr>\n",
- " <tr>\n",
- " <th>10</th>\n",
- " <td> 0.395775</td>\n",
- " <td> 0.000173</td>\n",
- " <td> 0.812032</td>\n",
- " <td> L_S2_L_S-M1</td>\n",
- " <td> 0.093203</td>\n",
- " <td> 0.561145</td>\n",
- " <td> 0.105617</td>\n",
- " <td> 0.178244</td>\n",
- " </tr>\n",
- " <tr>\n",
- " <th>16</th>\n",
- " <td> 0.615387</td>\n",
- " <td> 0.000078</td>\n",
- " <td> 0.934255</td>\n",
- " <td> R_S-M1_L_S-M1</td>\n",
- " <td> 0.074787</td>\n",
- " <td> 0.099295</td>\n",
- " <td> 0.683167</td>\n",
- " <td> 0.171543</td>\n",
- " </tr>\n",
- " <tr>\n",
- " <th>21</th>\n",
- " <td> 0.848833</td>\n",
- " <td> 0.001316</td>\n",
- " <td> 0.674588</td>\n",
- " <td> R_deepgray_L_PFC</td>\n",
- " <td> 0.041422</td>\n",
- " <td> 0.208550</td>\n",
- " <td> 0.446654</td>\n",
- " <td> 0.170309</td>\n",
- " </tr>\n",
- " <tr>\n",
- " <th>24</th>\n",
- " <td> 0.254644</td>\n",
- " <td> 0.000211</td>\n",
- " <td> 0.924624</td>\n",
- " <td> R_deepgray_R_S2</td>\n",
- " <td> 0.064068</td>\n",
- " <td> 0.093760</td>\n",
- " <td> 0.501135</td>\n",
- " <td> 0.167555</td>\n",
- " </tr>\n",
- " <tr>\n",
- " <th>2 </th>\n",
- " <td> 0.633750</td>\n",
- " <td> 0.000357</td>\n",
- " <td> 0.905607</td>\n",
- " <td> L_AS-Sens-C_R_AS-Sens-C</td>\n",
- " <td> 0.109588</td>\n",
- " <td> 0.696519</td>\n",
- " <td> 0.212974</td>\n",
- " <td> 0.164318</td>\n",
- " </tr>\n",
- " <tr>\n",
- " <th>17</th>\n",
- " <td> 0.752484</td>\n",
- " <td> 0.000689</td>\n",
- " <td> 0.952216</td>\n",
- " <td> R_S-M1_R_AS-Sens-C</td>\n",
- " <td> 0.059903</td>\n",
- " <td> 0.345851</td>\n",
- " <td> 0.648331</td>\n",
- " <td> 0.160446</td>\n",
- " </tr>\n",
- " <tr>\n",
- " <th>11</th>\n",
- " <td> 0.265150</td>\n",
- " <td> 0.000309</td>\n",
- " <td> 0.608655</td>\n",
- " <td> L_S2_L_deepgray</td>\n",
- " <td> 0.096973</td>\n",
- " <td> 0.811663</td>\n",
- " <td> 0.147194</td>\n",
- " <td> 0.160357</td>\n",
- " </tr>\n",
- " <tr>\n",
- " <th>9 </th>\n",
- " <td> 0.695413</td>\n",
- " <td> 0.000511</td>\n",
- " <td> 0.823294</td>\n",
- " <td> L_S2_L_AS-Sens-C</td>\n",
- " <td> 0.098898</td>\n",
- " <td> 0.284123</td>\n",
- " <td> 0.963899</td>\n",
- " <td> 0.155628</td>\n",
- " </tr>\n",
- " <tr>\n",
- " <th>12</th>\n",
- " <td> 0.433166</td>\n",
- " <td> 0.000247</td>\n",
- " <td> 0.875090</td>\n",
- " <td> L_S2_R_AS-Sens-C</td>\n",
- " <td> 0.083617</td>\n",
- " <td> 0.602741</td>\n",
- " <td> 0.513080</td>\n",
- " <td> 0.143404</td>\n",
- " </tr>\n",
- " </tbody>\n",
- "</table>\n",
- "</div>"
- ],
- "text/plain": [
- " Intercept age_at_visit gender[T.M] measure \\\n",
- "7 0.139689 0.000889 0.784138 L_S-M1_L_PFC \n",
- "18 0.503964 0.003005 0.725810 R_S-M1_R_PFC \n",
- "4 0.723769 0.000205 0.762192 L_AS-Sens-C_R_deepgray \n",
- "19 0.181365 0.002724 0.699412 R_S-M1_R_S2 \n",
- "15 0.128126 0.002279 0.772100 R_PFC_L_deepgray \n",
- "14 0.905604 0.000916 0.991047 R_AS-Sens-C_R_S2 \n",
- "20 0.538895 0.000267 0.961521 R_S-M1_R_deepgray \n",
- "13 0.484858 0.000346 0.837863 R_AS-Sens-C_L_deepgray \n",
- "6 0.440572 0.001456 0.843126 L_PFC_R_PFC \n",
- "3 0.595859 0.000703 0.898476 L_AS-Sens-C_R_S2 \n",
- "8 0.479252 0.000214 0.831458 L_S-M1_L_deepgray \n",
- "23 0.447515 0.000180 0.744982 R_deepgray_R_PFC \n",
- "0 0.630970 0.000807 0.946772 L_AS-Sens-C_L_S-M1 \n",
- "22 0.648730 0.000526 0.688840 R_deepgray_R_AS-Sens-C \n",
- "5 0.490455 0.000912 0.552587 L_PFC_L_deepgray \n",
- "1 0.461557 0.000675 0.722178 L_AS-Sens-C_L_deepgray \n",
- "10 0.395775 0.000173 0.812032 L_S2_L_S-M1 \n",
- "16 0.615387 0.000078 0.934255 R_S-M1_L_S-M1 \n",
- "21 0.848833 0.001316 0.674588 R_deepgray_L_PFC \n",
- "24 0.254644 0.000211 0.924624 R_deepgray_R_S2 \n",
- "2 0.633750 0.000357 0.905607 L_AS-Sens-C_R_AS-Sens-C \n",
- "17 0.752484 0.000689 0.952216 R_S-M1_R_AS-Sens-C \n",
- "11 0.265150 0.000309 0.608655 L_S2_L_deepgray \n",
- "9 0.695413 0.000511 0.823294 L_S2_L_AS-Sens-C \n",
- "12 0.433166 0.000247 0.875090 L_S2_R_AS-Sens-C \n",
- "\n",
- " r00_correlation r00_md r00_mtr_map r2 \n",
- "7 0.000316 0.831087 0.033602 0.238137 \n",
- "18 0.000253 0.162028 0.755075 0.216037 \n",
- "4 0.024560 0.304218 0.259391 0.206945 \n",
- "19 0.007781 0.016981 0.534828 0.203086 \n",
- "15 0.004811 0.029533 0.320453 0.200328 \n",
- "14 0.058617 0.070504 0.500667 0.199615 \n",
- "20 0.013669 0.088691 0.918744 0.197836 \n",
- "13 0.009018 0.627967 0.154479 0.193647 \n",
- "6 0.008558 0.139267 0.571633 0.187000 \n",
- "3 0.028592 0.570917 0.191679 0.186707 \n",
- "8 0.047262 0.461983 0.142359 0.184781 \n",
- "23 0.006440 0.229144 0.728392 0.184582 \n",
- "0 0.184438 0.201691 0.095884 0.181971 \n",
- "22 0.010884 0.279351 0.991240 0.181893 \n",
- "5 0.043813 0.374006 0.187795 0.181429 \n",
- "1 0.009938 0.141772 0.828873 0.180617 \n",
- "10 0.093203 0.561145 0.105617 0.178244 \n",
- "16 0.074787 0.099295 0.683167 0.171543 \n",
- "21 0.041422 0.208550 0.446654 0.170309 \n",
- "24 0.064068 0.093760 0.501135 0.167555 \n",
- "2 0.109588 0.696519 0.212974 0.164318 \n",
- "17 0.059903 0.345851 0.648331 0.160446 \n",
- "11 0.096973 0.811663 0.147194 0.160357 \n",
- "9 0.098898 0.284123 0.963899 0.155628 \n",
- "12 0.083617 0.602741 0.513080 0.143404 "
- ]
- },
- "execution_count": 46,
- "metadata": {},
- "output_type": "execute_result"
- }
- ],
- "source": [
- "results.sort(\"r2\")[::-1]"
- ]
- },
- {
- "cell_type": "code",
- "execution_count": 30,
- "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>Intercept</th>\n",
- " <th>age_at_visit</th>\n",
- " <th>gender[T.M]</th>\n",
- " <th>measure</th>\n",
- " <th>r00_correlation</th>\n",
- " <th>r00_mt_map</th>\n",
- " <th>r2</th>\n",
- " </tr>\n",
- " </thead>\n",
- " <tbody>\n",
- " <tr>\n",
- " <th>8 </th>\n",
- " <td> 0.102274</td>\n",
- " <td> 0.000004</td>\n",
- " <td> 0.260861</td>\n",
- " <td> L_S-M1_L_deepgray</td>\n",
- " <td> 0.110876</td>\n",
- " <td> 0.124282</td>\n",
- " <td> 0.175911</td>\n",
- " </tr>\n",
- " <tr>\n",
- " <th>6 </th>\n",
- " <td> 0.367417</td>\n",
- " <td> 0.000111</td>\n",
- " <td> 0.136767</td>\n",
- " <td> L_PFC_R_PFC</td>\n",
- " <td> 0.119911</td>\n",
- " <td> 0.211216</td>\n",
- " <td> 0.174641</td>\n",
- " </tr>\n",
- " <tr>\n",
- " <th>11</th>\n",
- " <td> 0.315755</td>\n",
- " <td> 0.000005</td>\n",
- " <td> 0.299819</td>\n",
- " <td> L_S2_L_deepgray</td>\n",
- " <td> 0.130187</td>\n",
- " <td> 0.336098</td>\n",
- " <td> 0.168614</td>\n",
- " </tr>\n",
- " <tr>\n",
- " <th>18</th>\n",
- " <td> 0.528601</td>\n",
- " <td> 0.000094</td>\n",
- " <td> 0.143693</td>\n",
- " <td> R_S-M1_R_PFC</td>\n",
- " <td> 0.199807</td>\n",
- " <td> 0.744565</td>\n",
- " <td> 0.163194</td>\n",
- " </tr>\n",
- " <tr>\n",
- " <th>20</th>\n",
- " <td> 0.391619</td>\n",
- " <td> 0.000003</td>\n",
- " <td> 0.244441</td>\n",
- " <td> R_S-M1_R_deepgray</td>\n",
- " <td> 0.224516</td>\n",
- " <td> 0.496301</td>\n",
- " <td> 0.163075</td>\n",
- " </tr>\n",
- " <tr>\n",
- " <th>4 </th>\n",
- " <td> 0.066396</td>\n",
- " <td> 0.000003</td>\n",
- " <td> 0.185651</td>\n",
- " <td> L_AS-Sens-C_R_deepgray</td>\n",
- " <td> 0.262431</td>\n",
- " <td> 0.039983</td>\n",
- " <td> 0.179182</td>\n",
- " </tr>\n",
- " <tr>\n",
- " <th>19</th>\n",
- " <td> 0.330611</td>\n",
- " <td> 0.000116</td>\n",
- " <td> 0.144397</td>\n",
- " <td> R_S-M1_R_S2</td>\n",
- " <td> 0.321852</td>\n",
- " <td> 0.361948</td>\n",
- " <td> 0.164666</td>\n",
- " </tr>\n",
- " <tr>\n",
- " <th>1 </th>\n",
- " <td> 0.635464</td>\n",
- " <td> 0.000004</td>\n",
- " <td> 0.242858</td>\n",
- " <td> L_AS-Sens-C_L_deepgray</td>\n",
- " <td> 0.429835</td>\n",
- " <td> 0.720494</td>\n",
- " <td> 0.157303</td>\n",
- " </tr>\n",
- " <tr>\n",
- " <th>13</th>\n",
- " <td> 0.187570</td>\n",
- " <td> 0.000004</td>\n",
- " <td> 0.173831</td>\n",
- " <td> R_AS-Sens-C_L_deepgray</td>\n",
- " <td> 0.488205</td>\n",
- " <td> 0.163867</td>\n",
- " <td> 0.165688</td>\n",
- " </tr>\n",
- " <tr>\n",
- " <th>22</th>\n",
- " <td> 0.453454</td>\n",
- " <td> 0.000006</td>\n",
- " <td> 0.240347</td>\n",
- " <td> R_deepgray_R_AS-Sens-C</td>\n",
- " <td> 0.512192</td>\n",
- " <td> 0.513965</td>\n",
- " <td> 0.157535</td>\n",
- " </tr>\n",
- " <tr>\n",
- " <th>24</th>\n",
- " <td> 0.578129</td>\n",
- " <td> 0.000005</td>\n",
- " <td> 0.265644</td>\n",
- " <td> R_deepgray_R_S2</td>\n",
- " <td> 0.581158</td>\n",
- " <td> 0.671305</td>\n",
- " <td> 0.155675</td>\n",
- " </tr>\n",
- " <tr>\n",
- " <th>10</th>\n",
- " <td> 0.086713</td>\n",
- " <td> 0.000008</td>\n",
- " <td> 0.283359</td>\n",
- " <td> L_S2_L_S-M1</td>\n",
- " <td> 0.609827</td>\n",
- " <td> 0.090642</td>\n",
- " <td> 0.169303</td>\n",
- " </tr>\n",
- " <tr>\n",
- " <th>23</th>\n",
- " <td> 0.957143</td>\n",
- " <td> 0.000008</td>\n",
- " <td> 0.238453</td>\n",
- " <td> R_deepgray_R_PFC</td>\n",
- " <td> 0.614807</td>\n",
- " <td> 0.849846</td>\n",
- " <td> 0.154907</td>\n",
- " </tr>\n",
- " <tr>\n",
- " <th>7 </th>\n",
- " <td> 0.098529</td>\n",
- " <td> 0.000023</td>\n",
- " <td> 0.156098</td>\n",
- " <td> L_S-M1_L_PFC</td>\n",
- " <td> 0.627938</td>\n",
- " <td> 0.103787</td>\n",
- " <td> 0.169381</td>\n",
- " </tr>\n",
- " <tr>\n",
- " <th>17</th>\n",
- " <td> 0.414644</td>\n",
- " <td> 0.000026</td>\n",
- " <td> 0.183286</td>\n",
- " <td> R_S-M1_R_AS-Sens-C</td>\n",
- " <td> 0.631614</td>\n",
- " <td> 0.489520</td>\n",
- " <td> 0.157911</td>\n",
- " </tr>\n",
- " <tr>\n",
- " <th>15</th>\n",
- " <td> 0.662763</td>\n",
- " <td> 0.000012</td>\n",
- " <td> 0.230427</td>\n",
- " <td> R_PFC_L_deepgray</td>\n",
- " <td> 0.665886</td>\n",
- " <td> 0.824459</td>\n",
- " <td> 0.154979</td>\n",
- " </tr>\n",
- " <tr>\n",
- " <th>9 </th>\n",
- " <td> 0.653553</td>\n",
- " <td> 0.000005</td>\n",
- " <td> 0.284665</td>\n",
- " <td> L_S2_L_AS-Sens-C</td>\n",
- " <td> 0.707206</td>\n",
- " <td> 0.795643</td>\n",
- " <td> 0.154416</td>\n",
- " </tr>\n",
- " <tr>\n",
- " <th>12</th>\n",
- " <td> 0.845060</td>\n",
- " <td> 0.000002</td>\n",
- " <td> 0.273083</td>\n",
- " <td> L_S2_R_AS-Sens-C</td>\n",
- " <td> 0.741972</td>\n",
- " <td> 0.626063</td>\n",
- " <td> 0.155587</td>\n",
- " </tr>\n",
- " <tr>\n",
- " <th>16</th>\n",
- " <td> 0.425286</td>\n",
- " <td> 0.000019</td>\n",
- " <td> 0.209314</td>\n",
- " <td> R_S-M1_L_S-M1</td>\n",
- " <td> 0.812435</td>\n",
- " <td> 0.438653</td>\n",
- " <td> 0.157009</td>\n",
- " </tr>\n",
- " <tr>\n",
- " <th>0 </th>\n",
- " <td> 0.134851</td>\n",
- " <td> 0.000014</td>\n",
- " <td> 0.218120</td>\n",
- " <td> L_AS-Sens-C_L_S-M1</td>\n",
- " <td> 0.814273</td>\n",
- " <td> 0.143198</td>\n",
- " <td> 0.165081</td>\n",
- " </tr>\n",
- " <tr>\n",
- " <th>21</th>\n",
- " <td> 0.312527</td>\n",
- " <td> 0.000015</td>\n",
- " <td> 0.201064</td>\n",
- " <td> R_deepgray_L_PFC</td>\n",
- " <td> 0.881631</td>\n",
- " <td> 0.368024</td>\n",
- " <td> 0.157931</td>\n",
- " </tr>\n",
- " <tr>\n",
- " <th>5 </th>\n",
- " <td> 0.413435</td>\n",
- " <td> 0.000022</td>\n",
- " <td> 0.222746</td>\n",
- " <td> L_PFC_L_deepgray</td>\n",
- " <td> 0.889638</td>\n",
- " <td> 0.494838</td>\n",
- " <td> 0.156065</td>\n",
- " </tr>\n",
- " <tr>\n",
- " <th>3 </th>\n",
- " <td> 0.304998</td>\n",
- " <td> 0.000016</td>\n",
- " <td> 0.189790</td>\n",
- " <td> L_AS-Sens-C_R_S2</td>\n",
- " <td> 0.938010</td>\n",
- " <td> 0.351768</td>\n",
- " <td> 0.158349</td>\n",
- " </tr>\n",
- " <tr>\n",
- " <th>14</th>\n",
- " <td> 0.169929</td>\n",
- " <td> 0.000037</td>\n",
- " <td> 0.190827</td>\n",
- " <td> R_AS-Sens-C_R_S2</td>\n",
- " <td> 0.954645</td>\n",
- " <td> 0.185670</td>\n",
- " <td> 0.163750</td>\n",
- " </tr>\n",
- " <tr>\n",
- " <th>2 </th>\n",
- " <td> 0.335984</td>\n",
- " <td> 0.000013</td>\n",
- " <td> 0.202978</td>\n",
- " <td> L_AS-Sens-C_R_AS-Sens-C</td>\n",
- " <td> 0.998496</td>\n",
- " <td> 0.389903</td>\n",
- " <td> 0.157566</td>\n",
- " </tr>\n",
- " </tbody>\n",
- "</table>\n",
- "</div>"
- ],
- "text/plain": [
- " Intercept age_at_visit gender[T.M] measure \\\n",
- "8 0.102274 0.000004 0.260861 L_S-M1_L_deepgray \n",
- "6 0.367417 0.000111 0.136767 L_PFC_R_PFC \n",
- "11 0.315755 0.000005 0.299819 L_S2_L_deepgray \n",
- "18 0.528601 0.000094 0.143693 R_S-M1_R_PFC \n",
- "20 0.391619 0.000003 0.244441 R_S-M1_R_deepgray \n",
- "4 0.066396 0.000003 0.185651 L_AS-Sens-C_R_deepgray \n",
- "19 0.330611 0.000116 0.144397 R_S-M1_R_S2 \n",
- "1 0.635464 0.000004 0.242858 L_AS-Sens-C_L_deepgray \n",
- "13 0.187570 0.000004 0.173831 R_AS-Sens-C_L_deepgray \n",
- "22 0.453454 0.000006 0.240347 R_deepgray_R_AS-Sens-C \n",
- "24 0.578129 0.000005 0.265644 R_deepgray_R_S2 \n",
- "10 0.086713 0.000008 0.283359 L_S2_L_S-M1 \n",
- "23 0.957143 0.000008 0.238453 R_deepgray_R_PFC \n",
- "7 0.098529 0.000023 0.156098 L_S-M1_L_PFC \n",
- "17 0.414644 0.000026 0.183286 R_S-M1_R_AS-Sens-C \n",
- "15 0.662763 0.000012 0.230427 R_PFC_L_deepgray \n",
- "9 0.653553 0.000005 0.284665 L_S2_L_AS-Sens-C \n",
- "12 0.845060 0.000002 0.273083 L_S2_R_AS-Sens-C \n",
- "16 0.425286 0.000019 0.209314 R_S-M1_L_S-M1 \n",
- "0 0.134851 0.000014 0.218120 L_AS-Sens-C_L_S-M1 \n",
- "21 0.312527 0.000015 0.201064 R_deepgray_L_PFC \n",
- "5 0.413435 0.000022 0.222746 L_PFC_L_deepgray \n",
- "3 0.304998 0.000016 0.189790 L_AS-Sens-C_R_S2 \n",
- "14 0.169929 0.000037 0.190827 R_AS-Sens-C_R_S2 \n",
- "2 0.335984 0.000013 0.202978 L_AS-Sens-C_R_AS-Sens-C \n",
- "\n",
- " r00_correlation r00_mt_map r2 \n",
- "8 0.110876 0.124282 0.175911 \n",
- "6 0.119911 0.211216 0.174641 \n",
- "11 0.130187 0.336098 0.168614 \n",
- "18 0.199807 0.744565 0.163194 \n",
- "20 0.224516 0.496301 0.163075 \n",
- "4 0.262431 0.039983 0.179182 \n",
- "19 0.321852 0.361948 0.164666 \n",
- "1 0.429835 0.720494 0.157303 \n",
- "13 0.488205 0.163867 0.165688 \n",
- "22 0.512192 0.513965 0.157535 \n",
- "24 0.581158 0.671305 0.155675 \n",
- "10 0.609827 0.090642 0.169303 \n",
- "23 0.614807 0.849846 0.154907 \n",
- "7 0.627938 0.103787 0.169381 \n",
- "17 0.631614 0.489520 0.157911 \n",
- "15 0.665886 0.824459 0.154979 \n",
- "9 0.707206 0.795643 0.154416 \n",
- "12 0.741972 0.626063 0.155587 \n",
- "16 0.812435 0.438653 0.157009 \n",
- "0 0.814273 0.143198 0.165081 \n",
- "21 0.881631 0.368024 0.157931 \n",
- "5 0.889638 0.494838 0.156065 \n",
- "3 0.938010 0.351768 0.158349 \n",
- "14 0.954645 0.185670 0.163750 \n",
- "2 0.998496 0.389903 0.157566 "
- ]
- },
- "execution_count": 30,
- "metadata": {},
- "output_type": "execute_result"
- }
- ],
- "source": [
- "results = pd.DataFrame()\n",
- "for conn, x in df.groupby(\"connection\"):\n",
- " model = sf.ols(\"timed_walk_trial1_time ~ r00_correlation + r00_mt_map +age_at_visit +gender\", x).fit()\n",
- " foo = model.pvalues\n",
- " foo[\"r2\"] = model.rsquared\n",
- " foo[\"measure\"] = conn\n",
- " results = results.append(foo,ignore_index=True)\n",
- "results.sort(\"r00_correlation\")"
- ]
- },
- {
- "cell_type": "code",
- "execution_count": 31,
- "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>Intercept</th>\n",
- " <th>age_at_visit</th>\n",
- " <th>gender[T.M]</th>\n",
- " <th>measure</th>\n",
- " <th>r00_correlation</th>\n",
- " <th>r00_mt_map</th>\n",
- " <th>r2</th>\n",
- " </tr>\n",
- " </thead>\n",
- " <tbody>\n",
- " <tr>\n",
- " <th>18</th>\n",
- " <td> 0.001050</td>\n",
- " <td> 0.027799</td>\n",
- " <td> 0.054800</td>\n",
- " <td> R_S-M1_R_PFC</td>\n",
- " <td> 0.042501</td>\n",
- " <td> 0.219691</td>\n",
- " <td> 0.136168</td>\n",
- " </tr>\n",
- " <tr>\n",
- " <th>6 </th>\n",
- " <td> 0.000064</td>\n",
- " <td> 0.023402</td>\n",
- " <td> 0.043774</td>\n",
- " <td> L_PFC_R_PFC</td>\n",
- " <td> 0.056876</td>\n",
- " <td> 0.027625</td>\n",
- " <td> 0.149855</td>\n",
- " </tr>\n",
- " <tr>\n",
- " <th>16</th>\n",
- " <td> 0.003051</td>\n",
- " <td> 0.009895</td>\n",
- " <td> 0.049720</td>\n",
- " <td> R_S-M1_L_S-M1</td>\n",
- " <td> 0.229823</td>\n",
- " <td> 0.076170</td>\n",
- " <td> 0.128553</td>\n",
- " </tr>\n",
- " <tr>\n",
- " <th>19</th>\n",
- " <td> 0.002014</td>\n",
- " <td> 0.015761</td>\n",
- " <td> 0.050947</td>\n",
- " <td> R_S-M1_R_S2</td>\n",
- " <td> 0.232018</td>\n",
- " <td> 0.163115</td>\n",
- " <td> 0.124931</td>\n",
- " </tr>\n",
- " <tr>\n",
- " <th>17</th>\n",
- " <td> 0.000062</td>\n",
- " <td> 0.014243</td>\n",
- " <td> 0.051321</td>\n",
- " <td> R_S-M1_R_AS-Sens-C</td>\n",
- " <td> 0.237646</td>\n",
- " <td> 0.067674</td>\n",
- " <td> 0.133613</td>\n",
- " </tr>\n",
- " <tr>\n",
- " <th>15</th>\n",
- " <td> 0.012735</td>\n",
- " <td> 0.002615</td>\n",
- " <td> 0.018871</td>\n",
- " <td> R_PFC_L_deepgray</td>\n",
- " <td> 0.376680</td>\n",
- " <td> 0.820765</td>\n",
- " <td> 0.108716</td>\n",
- " </tr>\n",
- " <tr>\n",
- " <th>3 </th>\n",
- " <td> 0.000012</td>\n",
- " <td> 0.008224</td>\n",
- " <td> 0.054157</td>\n",
- " <td> L_AS-Sens-C_R_S2</td>\n",
- " <td> 0.412252</td>\n",
- " <td> 0.061450</td>\n",
- " <td> 0.128726</td>\n",
- " </tr>\n",
- " <tr>\n",
- " <th>7 </th>\n",
- " <td> 0.000000</td>\n",
- " <td> 0.005829</td>\n",
- " <td> 0.044702</td>\n",
- " <td> L_S-M1_L_PFC</td>\n",
- " <td> 0.451307</td>\n",
- " <td> 0.006729</td>\n",
- " <td> 0.148246</td>\n",
- " </tr>\n",
- " <tr>\n",
- " <th>2 </th>\n",
- " <td> 0.000037</td>\n",
- " <td> 0.005917</td>\n",
- " <td> 0.044495</td>\n",
- " <td> L_AS-Sens-C_R_AS-Sens-C</td>\n",
- " <td> 0.462218</td>\n",
- " <td> 0.106427</td>\n",
- " <td> 0.123201</td>\n",
- " </tr>\n",
- " <tr>\n",
- " <th>23</th>\n",
- " <td> 0.000828</td>\n",
- " <td> 0.006420</td>\n",
- " <td> 0.021273</td>\n",
- " <td> R_deepgray_R_PFC</td>\n",
- " <td> 0.521677</td>\n",
- " <td> 0.171136</td>\n",
- " <td> 0.118901</td>\n",
- " </tr>\n",
- " <tr>\n",
- " <th>14</th>\n",
- " <td> 0.000064</td>\n",
- " <td> 0.011467</td>\n",
- " <td> 0.038808</td>\n",
- " <td> R_AS-Sens-C_R_S2</td>\n",
- " <td> 0.588705</td>\n",
- " <td> 0.041388</td>\n",
- " <td> 0.132721</td>\n",
- " </tr>\n",
- " <tr>\n",
- " <th>12</th>\n",
- " <td> 0.016344</td>\n",
- " <td> 0.002152</td>\n",
- " <td> 0.029329</td>\n",
- " <td> L_S2_R_AS-Sens-C</td>\n",
- " <td> 0.595380</td>\n",
- " <td> 0.963862</td>\n",
- " <td> 0.105163</td>\n",
- " </tr>\n",
- " <tr>\n",
- " <th>8 </th>\n",
- " <td> 0.000037</td>\n",
- " <td> 0.004482</td>\n",
- " <td> 0.021598</td>\n",
- " <td> L_S-M1_L_deepgray</td>\n",
- " <td> 0.603762</td>\n",
- " <td> 0.015219</td>\n",
- " <td> 0.136280</td>\n",
- " </tr>\n",
- " <tr>\n",
- " <th>5 </th>\n",
- " <td> 0.000197</td>\n",
- " <td> 0.012670</td>\n",
- " <td> 0.019629</td>\n",
- " <td> L_PFC_L_deepgray</td>\n",
- " <td> 0.657948</td>\n",
- " <td> 0.063127</td>\n",
- " <td> 0.124463</td>\n",
- " </tr>\n",
- " <tr>\n",
- " <th>21</th>\n",
- " <td> 0.000007</td>\n",
- " <td> 0.009961</td>\n",
- " <td> 0.028479</td>\n",
- " <td> R_deepgray_L_PFC</td>\n",
- " <td> 0.681705</td>\n",
- " <td> 0.023235</td>\n",
- " <td> 0.133790</td>\n",
- " </tr>\n",
- " <tr>\n",
- " <th>9 </th>\n",
- " <td> 0.002011</td>\n",
- " <td> 0.004132</td>\n",
- " <td> 0.025987</td>\n",
- " <td> L_S2_L_AS-Sens-C</td>\n",
- " <td> 0.716290</td>\n",
- " <td> 0.289003</td>\n",
- " <td> 0.112135</td>\n",
- " </tr>\n",
- " <tr>\n",
- " <th>13</th>\n",
- " <td> 0.000002</td>\n",
- " <td> 0.003724</td>\n",
- " <td> 0.038902</td>\n",
- " <td> R_AS-Sens-C_L_deepgray</td>\n",
- " <td> 0.730586</td>\n",
- " <td> 0.020673</td>\n",
- " <td> 0.134931</td>\n",
- " </tr>\n",
- " <tr>\n",
- " <th>22</th>\n",
- " <td> 0.000032</td>\n",
- " <td> 0.007084</td>\n",
- " <td> 0.017682</td>\n",
- " <td> R_deepgray_R_AS-Sens-C</td>\n",
- " <td> 0.780205</td>\n",
- " <td> 0.066916</td>\n",
- " <td> 0.123963</td>\n",
- " </tr>\n",
- " <tr>\n",
- " <th>11</th>\n",
- " <td> 0.000081</td>\n",
- " <td> 0.009850</td>\n",
- " <td> 0.015995</td>\n",
- " <td> L_S2_L_deepgray</td>\n",
- " <td> 0.796521</td>\n",
- " <td> 0.031031</td>\n",
- " <td> 0.129467</td>\n",
- " </tr>\n",
- " <tr>\n",
- " <th>1 </th>\n",
- " <td> 0.001578</td>\n",
- " <td> 0.004778</td>\n",
- " <td> 0.017942</td>\n",
- " <td> L_AS-Sens-C_L_deepgray</td>\n",
- " <td> 0.813290</td>\n",
- " <td> 0.215427</td>\n",
- " <td> 0.112836</td>\n",
- " </tr>\n",
- " <tr>\n",
- " <th>24</th>\n",
- " <td> 0.002798</td>\n",
- " <td> 0.002842</td>\n",
- " <td> 0.017301</td>\n",
- " <td> R_deepgray_R_S2</td>\n",
- " <td> 0.825346</td>\n",
- " <td> 0.203149</td>\n",
- " <td> 0.112634</td>\n",
- " </tr>\n",
- " <tr>\n",
- " <th>0 </th>\n",
- " <td> 0.000000</td>\n",
- " <td> 0.008066</td>\n",
- " <td> 0.036150</td>\n",
- " <td> L_AS-Sens-C_L_S-M1</td>\n",
- " <td> 0.828228</td>\n",
- " <td> 0.001947</td>\n",
- " <td> 0.159483</td>\n",
- " </tr>\n",
- " <tr>\n",
- " <th>4 </th>\n",
- " <td> 0.000001</td>\n",
- " <td> 0.002742</td>\n",
- " <td> 0.025568</td>\n",
- " <td> L_AS-Sens-C_R_deepgray</td>\n",
- " <td> 0.922103</td>\n",
- " <td> 0.009649</td>\n",
- " <td> 0.141273</td>\n",
- " </tr>\n",
- " <tr>\n",
- " <th>10</th>\n",
- " <td> 0.000002</td>\n",
- " <td> 0.003864</td>\n",
- " <td> 0.021712</td>\n",
- " <td> L_S2_L_S-M1</td>\n",
- " <td> 0.988556</td>\n",
- " <td> 0.005241</td>\n",
- " <td> 0.147528</td>\n",
- " </tr>\n",
- " <tr>\n",
- " <th>20</th>\n",
- " <td> 0.000133</td>\n",
- " <td> 0.004122</td>\n",
- " <td> 0.033329</td>\n",
- " <td> R_S-M1_R_deepgray</td>\n",
- " <td> 0.998675</td>\n",
- " <td> 0.084061</td>\n",
- " <td> 0.120437</td>\n",
- " </tr>\n",
- " </tbody>\n",
- "</table>\n",
- "</div>"
- ],
- "text/plain": [
- " Intercept age_at_visit gender[T.M] measure \\\n",
- "18 0.001050 0.027799 0.054800 R_S-M1_R_PFC \n",
- "6 0.000064 0.023402 0.043774 L_PFC_R_PFC \n",
- "16 0.003051 0.009895 0.049720 R_S-M1_L_S-M1 \n",
- "19 0.002014 0.015761 0.050947 R_S-M1_R_S2 \n",
- "17 0.000062 0.014243 0.051321 R_S-M1_R_AS-Sens-C \n",
- "15 0.012735 0.002615 0.018871 R_PFC_L_deepgray \n",
- "3 0.000012 0.008224 0.054157 L_AS-Sens-C_R_S2 \n",
- "7 0.000000 0.005829 0.044702 L_S-M1_L_PFC \n",
- "2 0.000037 0.005917 0.044495 L_AS-Sens-C_R_AS-Sens-C \n",
- "23 0.000828 0.006420 0.021273 R_deepgray_R_PFC \n",
- "14 0.000064 0.011467 0.038808 R_AS-Sens-C_R_S2 \n",
- "12 0.016344 0.002152 0.029329 L_S2_R_AS-Sens-C \n",
- "8 0.000037 0.004482 0.021598 L_S-M1_L_deepgray \n",
- "5 0.000197 0.012670 0.019629 L_PFC_L_deepgray \n",
- "21 0.000007 0.009961 0.028479 R_deepgray_L_PFC \n",
- "9 0.002011 0.004132 0.025987 L_S2_L_AS-Sens-C \n",
- "13 0.000002 0.003724 0.038902 R_AS-Sens-C_L_deepgray \n",
- "22 0.000032 0.007084 0.017682 R_deepgray_R_AS-Sens-C \n",
- "11 0.000081 0.009850 0.015995 L_S2_L_deepgray \n",
- "1 0.001578 0.004778 0.017942 L_AS-Sens-C_L_deepgray \n",
- "24 0.002798 0.002842 0.017301 R_deepgray_R_S2 \n",
- "0 0.000000 0.008066 0.036150 L_AS-Sens-C_L_S-M1 \n",
- "4 0.000001 0.002742 0.025568 L_AS-Sens-C_R_deepgray \n",
- "10 0.000002 0.003864 0.021712 L_S2_L_S-M1 \n",
- "20 0.000133 0.004122 0.033329 R_S-M1_R_deepgray \n",
- "\n",
- " r00_correlation r00_mt_map r2 \n",
- "18 0.042501 0.219691 0.136168 \n",
- "6 0.056876 0.027625 0.149855 \n",
- "16 0.229823 0.076170 0.128553 \n",
- "19 0.232018 0.163115 0.124931 \n",
- "17 0.237646 0.067674 0.133613 \n",
- "15 0.376680 0.820765 0.108716 \n",
- "3 0.412252 0.061450 0.128726 \n",
- "7 0.451307 0.006729 0.148246 \n",
- "2 0.462218 0.106427 0.123201 \n",
- "23 0.521677 0.171136 0.118901 \n",
- "14 0.588705 0.041388 0.132721 \n",
- "12 0.595380 0.963862 0.105163 \n",
- "8 0.603762 0.015219 0.136280 \n",
- "5 0.657948 0.063127 0.124463 \n",
- "21 0.681705 0.023235 0.133790 \n",
- "9 0.716290 0.289003 0.112135 \n",
- "13 0.730586 0.020673 0.134931 \n",
- "22 0.780205 0.066916 0.123963 \n",
- "11 0.796521 0.031031 0.129467 \n",
- "1 0.813290 0.215427 0.112836 \n",
- "24 0.825346 0.203149 0.112634 \n",
- "0 0.828228 0.001947 0.159483 \n",
- "4 0.922103 0.009649 0.141273 \n",
- "10 0.988556 0.005241 0.147528 \n",
- "20 0.998675 0.084061 0.120437 "
- ]
- },
- "execution_count": 31,
- "metadata": {},
- "output_type": "execute_result"
- }
- ],
- "source": [
- "results = pd.DataFrame()\n",
- "for conn, x in df.groupby(\"connection\"):\n",
- " model = sf.ols(\"nhpt_dominant_time ~ r00_correlation + r00_mt_map +age_at_visit +gender\", x).fit()\n",
- " foo = model.pvalues\n",
- " foo[\"r2\"] = model.rsquared\n",
- " foo[\"measure\"] = conn\n",
- " results = results.append(foo,ignore_index=True)\n",
- "results.sort(\"r00_correlation\")"
- ]
- },
- {
- "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
- }
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement