Advertisement
Guest User

Untitled

a guest
May 3rd, 2016
61
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 43.34 KB | None | 0 0
  1. {
  2. "cells": [
  3. {
  4. "cell_type": "code",
  5. "execution_count": 8,
  6. "metadata": {
  7. "collapsed": false
  8. },
  9. "outputs": [
  10. {
  11. "name": "stdout",
  12. "output_type": "stream",
  13. "text": [
  14. "Populating the interactive namespace from numpy and matplotlib\n"
  15. ]
  16. }
  17. ],
  18. "source": [
  19. "%pylab inline\n",
  20. "import seaborn as sns\n",
  21. "import pandas as pd\n",
  22. "import statsmodels.formula.api as sf"
  23. ]
  24. },
  25. {
  26. "cell_type": "code",
  27. "execution_count": 3,
  28. "metadata": {
  29. "collapsed": true
  30. },
  31. "outputs": [],
  32. "source": [
  33. "df = pd.read_csv(\"combined_demographics_imaging.csv\",index_col=0)"
  34. ]
  35. },
  36. {
  37. "cell_type": "code",
  38. "execution_count": 28,
  39. "metadata": {
  40. "collapsed": false
  41. },
  42. "outputs": [
  43. {
  44. "name": "stdout",
  45. "output_type": "stream",
  46. "text": [
  47. "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"
  48. ]
  49. }
  50. ],
  51. "source": [
  52. "print df.columns"
  53. ]
  54. },
  55. {
  56. "cell_type": "code",
  57. "execution_count": 45,
  58. "metadata": {
  59. "collapsed": false
  60. },
  61. "outputs": [],
  62. "source": [
  63. "results = pd.DataFrame()\n",
  64. "for conn, x in df.groupby(\"connection\"):\n",
  65. " model = sf.ols(\"edss ~ r00_correlation + r00_mtr_map +r00_md+age_at_visit +gender\", x).fit()\n",
  66. " foo = model.pvalues\n",
  67. " foo[\"r2\"] = model.rsquared\n",
  68. " foo[\"measure\"] = conn\n",
  69. " results = results.append(foo,ignore_index=True)\n"
  70. ]
  71. },
  72. {
  73. "cell_type": "code",
  74. "execution_count": 46,
  75. "metadata": {
  76. "collapsed": false
  77. },
  78. "outputs": [
  79. {
  80. "data": {
  81. "text/html": [
  82. "<div style=\"max-height:1000px;max-width:1500px;overflow:auto;\">\n",
  83. "<table border=\"1\" class=\"dataframe\">\n",
  84. " <thead>\n",
  85. " <tr style=\"text-align: right;\">\n",
  86. " <th></th>\n",
  87. " <th>Intercept</th>\n",
  88. " <th>age_at_visit</th>\n",
  89. " <th>gender[T.M]</th>\n",
  90. " <th>measure</th>\n",
  91. " <th>r00_correlation</th>\n",
  92. " <th>r00_md</th>\n",
  93. " <th>r00_mtr_map</th>\n",
  94. " <th>r2</th>\n",
  95. " </tr>\n",
  96. " </thead>\n",
  97. " <tbody>\n",
  98. " <tr>\n",
  99. " <th>7 </th>\n",
  100. " <td> 0.139689</td>\n",
  101. " <td> 0.000889</td>\n",
  102. " <td> 0.784138</td>\n",
  103. " <td> L_S-M1_L_PFC</td>\n",
  104. " <td> 0.000316</td>\n",
  105. " <td> 0.831087</td>\n",
  106. " <td> 0.033602</td>\n",
  107. " <td> 0.238137</td>\n",
  108. " </tr>\n",
  109. " <tr>\n",
  110. " <th>18</th>\n",
  111. " <td> 0.503964</td>\n",
  112. " <td> 0.003005</td>\n",
  113. " <td> 0.725810</td>\n",
  114. " <td> R_S-M1_R_PFC</td>\n",
  115. " <td> 0.000253</td>\n",
  116. " <td> 0.162028</td>\n",
  117. " <td> 0.755075</td>\n",
  118. " <td> 0.216037</td>\n",
  119. " </tr>\n",
  120. " <tr>\n",
  121. " <th>4 </th>\n",
  122. " <td> 0.723769</td>\n",
  123. " <td> 0.000205</td>\n",
  124. " <td> 0.762192</td>\n",
  125. " <td> L_AS-Sens-C_R_deepgray</td>\n",
  126. " <td> 0.024560</td>\n",
  127. " <td> 0.304218</td>\n",
  128. " <td> 0.259391</td>\n",
  129. " <td> 0.206945</td>\n",
  130. " </tr>\n",
  131. " <tr>\n",
  132. " <th>19</th>\n",
  133. " <td> 0.181365</td>\n",
  134. " <td> 0.002724</td>\n",
  135. " <td> 0.699412</td>\n",
  136. " <td> R_S-M1_R_S2</td>\n",
  137. " <td> 0.007781</td>\n",
  138. " <td> 0.016981</td>\n",
  139. " <td> 0.534828</td>\n",
  140. " <td> 0.203086</td>\n",
  141. " </tr>\n",
  142. " <tr>\n",
  143. " <th>15</th>\n",
  144. " <td> 0.128126</td>\n",
  145. " <td> 0.002279</td>\n",
  146. " <td> 0.772100</td>\n",
  147. " <td> R_PFC_L_deepgray</td>\n",
  148. " <td> 0.004811</td>\n",
  149. " <td> 0.029533</td>\n",
  150. " <td> 0.320453</td>\n",
  151. " <td> 0.200328</td>\n",
  152. " </tr>\n",
  153. " <tr>\n",
  154. " <th>14</th>\n",
  155. " <td> 0.905604</td>\n",
  156. " <td> 0.000916</td>\n",
  157. " <td> 0.991047</td>\n",
  158. " <td> R_AS-Sens-C_R_S2</td>\n",
  159. " <td> 0.058617</td>\n",
  160. " <td> 0.070504</td>\n",
  161. " <td> 0.500667</td>\n",
  162. " <td> 0.199615</td>\n",
  163. " </tr>\n",
  164. " <tr>\n",
  165. " <th>20</th>\n",
  166. " <td> 0.538895</td>\n",
  167. " <td> 0.000267</td>\n",
  168. " <td> 0.961521</td>\n",
  169. " <td> R_S-M1_R_deepgray</td>\n",
  170. " <td> 0.013669</td>\n",
  171. " <td> 0.088691</td>\n",
  172. " <td> 0.918744</td>\n",
  173. " <td> 0.197836</td>\n",
  174. " </tr>\n",
  175. " <tr>\n",
  176. " <th>13</th>\n",
  177. " <td> 0.484858</td>\n",
  178. " <td> 0.000346</td>\n",
  179. " <td> 0.837863</td>\n",
  180. " <td> R_AS-Sens-C_L_deepgray</td>\n",
  181. " <td> 0.009018</td>\n",
  182. " <td> 0.627967</td>\n",
  183. " <td> 0.154479</td>\n",
  184. " <td> 0.193647</td>\n",
  185. " </tr>\n",
  186. " <tr>\n",
  187. " <th>6 </th>\n",
  188. " <td> 0.440572</td>\n",
  189. " <td> 0.001456</td>\n",
  190. " <td> 0.843126</td>\n",
  191. " <td> L_PFC_R_PFC</td>\n",
  192. " <td> 0.008558</td>\n",
  193. " <td> 0.139267</td>\n",
  194. " <td> 0.571633</td>\n",
  195. " <td> 0.187000</td>\n",
  196. " </tr>\n",
  197. " <tr>\n",
  198. " <th>3 </th>\n",
  199. " <td> 0.595859</td>\n",
  200. " <td> 0.000703</td>\n",
  201. " <td> 0.898476</td>\n",
  202. " <td> L_AS-Sens-C_R_S2</td>\n",
  203. " <td> 0.028592</td>\n",
  204. " <td> 0.570917</td>\n",
  205. " <td> 0.191679</td>\n",
  206. " <td> 0.186707</td>\n",
  207. " </tr>\n",
  208. " <tr>\n",
  209. " <th>8 </th>\n",
  210. " <td> 0.479252</td>\n",
  211. " <td> 0.000214</td>\n",
  212. " <td> 0.831458</td>\n",
  213. " <td> L_S-M1_L_deepgray</td>\n",
  214. " <td> 0.047262</td>\n",
  215. " <td> 0.461983</td>\n",
  216. " <td> 0.142359</td>\n",
  217. " <td> 0.184781</td>\n",
  218. " </tr>\n",
  219. " <tr>\n",
  220. " <th>23</th>\n",
  221. " <td> 0.447515</td>\n",
  222. " <td> 0.000180</td>\n",
  223. " <td> 0.744982</td>\n",
  224. " <td> R_deepgray_R_PFC</td>\n",
  225. " <td> 0.006440</td>\n",
  226. " <td> 0.229144</td>\n",
  227. " <td> 0.728392</td>\n",
  228. " <td> 0.184582</td>\n",
  229. " </tr>\n",
  230. " <tr>\n",
  231. " <th>0 </th>\n",
  232. " <td> 0.630970</td>\n",
  233. " <td> 0.000807</td>\n",
  234. " <td> 0.946772</td>\n",
  235. " <td> L_AS-Sens-C_L_S-M1</td>\n",
  236. " <td> 0.184438</td>\n",
  237. " <td> 0.201691</td>\n",
  238. " <td> 0.095884</td>\n",
  239. " <td> 0.181971</td>\n",
  240. " </tr>\n",
  241. " <tr>\n",
  242. " <th>22</th>\n",
  243. " <td> 0.648730</td>\n",
  244. " <td> 0.000526</td>\n",
  245. " <td> 0.688840</td>\n",
  246. " <td> R_deepgray_R_AS-Sens-C</td>\n",
  247. " <td> 0.010884</td>\n",
  248. " <td> 0.279351</td>\n",
  249. " <td> 0.991240</td>\n",
  250. " <td> 0.181893</td>\n",
  251. " </tr>\n",
  252. " <tr>\n",
  253. " <th>5 </th>\n",
  254. " <td> 0.490455</td>\n",
  255. " <td> 0.000912</td>\n",
  256. " <td> 0.552587</td>\n",
  257. " <td> L_PFC_L_deepgray</td>\n",
  258. " <td> 0.043813</td>\n",
  259. " <td> 0.374006</td>\n",
  260. " <td> 0.187795</td>\n",
  261. " <td> 0.181429</td>\n",
  262. " </tr>\n",
  263. " <tr>\n",
  264. " <th>1 </th>\n",
  265. " <td> 0.461557</td>\n",
  266. " <td> 0.000675</td>\n",
  267. " <td> 0.722178</td>\n",
  268. " <td> L_AS-Sens-C_L_deepgray</td>\n",
  269. " <td> 0.009938</td>\n",
  270. " <td> 0.141772</td>\n",
  271. " <td> 0.828873</td>\n",
  272. " <td> 0.180617</td>\n",
  273. " </tr>\n",
  274. " <tr>\n",
  275. " <th>10</th>\n",
  276. " <td> 0.395775</td>\n",
  277. " <td> 0.000173</td>\n",
  278. " <td> 0.812032</td>\n",
  279. " <td> L_S2_L_S-M1</td>\n",
  280. " <td> 0.093203</td>\n",
  281. " <td> 0.561145</td>\n",
  282. " <td> 0.105617</td>\n",
  283. " <td> 0.178244</td>\n",
  284. " </tr>\n",
  285. " <tr>\n",
  286. " <th>16</th>\n",
  287. " <td> 0.615387</td>\n",
  288. " <td> 0.000078</td>\n",
  289. " <td> 0.934255</td>\n",
  290. " <td> R_S-M1_L_S-M1</td>\n",
  291. " <td> 0.074787</td>\n",
  292. " <td> 0.099295</td>\n",
  293. " <td> 0.683167</td>\n",
  294. " <td> 0.171543</td>\n",
  295. " </tr>\n",
  296. " <tr>\n",
  297. " <th>21</th>\n",
  298. " <td> 0.848833</td>\n",
  299. " <td> 0.001316</td>\n",
  300. " <td> 0.674588</td>\n",
  301. " <td> R_deepgray_L_PFC</td>\n",
  302. " <td> 0.041422</td>\n",
  303. " <td> 0.208550</td>\n",
  304. " <td> 0.446654</td>\n",
  305. " <td> 0.170309</td>\n",
  306. " </tr>\n",
  307. " <tr>\n",
  308. " <th>24</th>\n",
  309. " <td> 0.254644</td>\n",
  310. " <td> 0.000211</td>\n",
  311. " <td> 0.924624</td>\n",
  312. " <td> R_deepgray_R_S2</td>\n",
  313. " <td> 0.064068</td>\n",
  314. " <td> 0.093760</td>\n",
  315. " <td> 0.501135</td>\n",
  316. " <td> 0.167555</td>\n",
  317. " </tr>\n",
  318. " <tr>\n",
  319. " <th>2 </th>\n",
  320. " <td> 0.633750</td>\n",
  321. " <td> 0.000357</td>\n",
  322. " <td> 0.905607</td>\n",
  323. " <td> L_AS-Sens-C_R_AS-Sens-C</td>\n",
  324. " <td> 0.109588</td>\n",
  325. " <td> 0.696519</td>\n",
  326. " <td> 0.212974</td>\n",
  327. " <td> 0.164318</td>\n",
  328. " </tr>\n",
  329. " <tr>\n",
  330. " <th>17</th>\n",
  331. " <td> 0.752484</td>\n",
  332. " <td> 0.000689</td>\n",
  333. " <td> 0.952216</td>\n",
  334. " <td> R_S-M1_R_AS-Sens-C</td>\n",
  335. " <td> 0.059903</td>\n",
  336. " <td> 0.345851</td>\n",
  337. " <td> 0.648331</td>\n",
  338. " <td> 0.160446</td>\n",
  339. " </tr>\n",
  340. " <tr>\n",
  341. " <th>11</th>\n",
  342. " <td> 0.265150</td>\n",
  343. " <td> 0.000309</td>\n",
  344. " <td> 0.608655</td>\n",
  345. " <td> L_S2_L_deepgray</td>\n",
  346. " <td> 0.096973</td>\n",
  347. " <td> 0.811663</td>\n",
  348. " <td> 0.147194</td>\n",
  349. " <td> 0.160357</td>\n",
  350. " </tr>\n",
  351. " <tr>\n",
  352. " <th>9 </th>\n",
  353. " <td> 0.695413</td>\n",
  354. " <td> 0.000511</td>\n",
  355. " <td> 0.823294</td>\n",
  356. " <td> L_S2_L_AS-Sens-C</td>\n",
  357. " <td> 0.098898</td>\n",
  358. " <td> 0.284123</td>\n",
  359. " <td> 0.963899</td>\n",
  360. " <td> 0.155628</td>\n",
  361. " </tr>\n",
  362. " <tr>\n",
  363. " <th>12</th>\n",
  364. " <td> 0.433166</td>\n",
  365. " <td> 0.000247</td>\n",
  366. " <td> 0.875090</td>\n",
  367. " <td> L_S2_R_AS-Sens-C</td>\n",
  368. " <td> 0.083617</td>\n",
  369. " <td> 0.602741</td>\n",
  370. " <td> 0.513080</td>\n",
  371. " <td> 0.143404</td>\n",
  372. " </tr>\n",
  373. " </tbody>\n",
  374. "</table>\n",
  375. "</div>"
  376. ],
  377. "text/plain": [
  378. " Intercept age_at_visit gender[T.M] measure \\\n",
  379. "7 0.139689 0.000889 0.784138 L_S-M1_L_PFC \n",
  380. "18 0.503964 0.003005 0.725810 R_S-M1_R_PFC \n",
  381. "4 0.723769 0.000205 0.762192 L_AS-Sens-C_R_deepgray \n",
  382. "19 0.181365 0.002724 0.699412 R_S-M1_R_S2 \n",
  383. "15 0.128126 0.002279 0.772100 R_PFC_L_deepgray \n",
  384. "14 0.905604 0.000916 0.991047 R_AS-Sens-C_R_S2 \n",
  385. "20 0.538895 0.000267 0.961521 R_S-M1_R_deepgray \n",
  386. "13 0.484858 0.000346 0.837863 R_AS-Sens-C_L_deepgray \n",
  387. "6 0.440572 0.001456 0.843126 L_PFC_R_PFC \n",
  388. "3 0.595859 0.000703 0.898476 L_AS-Sens-C_R_S2 \n",
  389. "8 0.479252 0.000214 0.831458 L_S-M1_L_deepgray \n",
  390. "23 0.447515 0.000180 0.744982 R_deepgray_R_PFC \n",
  391. "0 0.630970 0.000807 0.946772 L_AS-Sens-C_L_S-M1 \n",
  392. "22 0.648730 0.000526 0.688840 R_deepgray_R_AS-Sens-C \n",
  393. "5 0.490455 0.000912 0.552587 L_PFC_L_deepgray \n",
  394. "1 0.461557 0.000675 0.722178 L_AS-Sens-C_L_deepgray \n",
  395. "10 0.395775 0.000173 0.812032 L_S2_L_S-M1 \n",
  396. "16 0.615387 0.000078 0.934255 R_S-M1_L_S-M1 \n",
  397. "21 0.848833 0.001316 0.674588 R_deepgray_L_PFC \n",
  398. "24 0.254644 0.000211 0.924624 R_deepgray_R_S2 \n",
  399. "2 0.633750 0.000357 0.905607 L_AS-Sens-C_R_AS-Sens-C \n",
  400. "17 0.752484 0.000689 0.952216 R_S-M1_R_AS-Sens-C \n",
  401. "11 0.265150 0.000309 0.608655 L_S2_L_deepgray \n",
  402. "9 0.695413 0.000511 0.823294 L_S2_L_AS-Sens-C \n",
  403. "12 0.433166 0.000247 0.875090 L_S2_R_AS-Sens-C \n",
  404. "\n",
  405. " r00_correlation r00_md r00_mtr_map r2 \n",
  406. "7 0.000316 0.831087 0.033602 0.238137 \n",
  407. "18 0.000253 0.162028 0.755075 0.216037 \n",
  408. "4 0.024560 0.304218 0.259391 0.206945 \n",
  409. "19 0.007781 0.016981 0.534828 0.203086 \n",
  410. "15 0.004811 0.029533 0.320453 0.200328 \n",
  411. "14 0.058617 0.070504 0.500667 0.199615 \n",
  412. "20 0.013669 0.088691 0.918744 0.197836 \n",
  413. "13 0.009018 0.627967 0.154479 0.193647 \n",
  414. "6 0.008558 0.139267 0.571633 0.187000 \n",
  415. "3 0.028592 0.570917 0.191679 0.186707 \n",
  416. "8 0.047262 0.461983 0.142359 0.184781 \n",
  417. "23 0.006440 0.229144 0.728392 0.184582 \n",
  418. "0 0.184438 0.201691 0.095884 0.181971 \n",
  419. "22 0.010884 0.279351 0.991240 0.181893 \n",
  420. "5 0.043813 0.374006 0.187795 0.181429 \n",
  421. "1 0.009938 0.141772 0.828873 0.180617 \n",
  422. "10 0.093203 0.561145 0.105617 0.178244 \n",
  423. "16 0.074787 0.099295 0.683167 0.171543 \n",
  424. "21 0.041422 0.208550 0.446654 0.170309 \n",
  425. "24 0.064068 0.093760 0.501135 0.167555 \n",
  426. "2 0.109588 0.696519 0.212974 0.164318 \n",
  427. "17 0.059903 0.345851 0.648331 0.160446 \n",
  428. "11 0.096973 0.811663 0.147194 0.160357 \n",
  429. "9 0.098898 0.284123 0.963899 0.155628 \n",
  430. "12 0.083617 0.602741 0.513080 0.143404 "
  431. ]
  432. },
  433. "execution_count": 46,
  434. "metadata": {},
  435. "output_type": "execute_result"
  436. }
  437. ],
  438. "source": [
  439. "results.sort(\"r2\")[::-1]"
  440. ]
  441. },
  442. {
  443. "cell_type": "code",
  444. "execution_count": 30,
  445. "metadata": {
  446. "collapsed": false
  447. },
  448. "outputs": [
  449. {
  450. "data": {
  451. "text/html": [
  452. "<div style=\"max-height:1000px;max-width:1500px;overflow:auto;\">\n",
  453. "<table border=\"1\" class=\"dataframe\">\n",
  454. " <thead>\n",
  455. " <tr style=\"text-align: right;\">\n",
  456. " <th></th>\n",
  457. " <th>Intercept</th>\n",
  458. " <th>age_at_visit</th>\n",
  459. " <th>gender[T.M]</th>\n",
  460. " <th>measure</th>\n",
  461. " <th>r00_correlation</th>\n",
  462. " <th>r00_mt_map</th>\n",
  463. " <th>r2</th>\n",
  464. " </tr>\n",
  465. " </thead>\n",
  466. " <tbody>\n",
  467. " <tr>\n",
  468. " <th>8 </th>\n",
  469. " <td> 0.102274</td>\n",
  470. " <td> 0.000004</td>\n",
  471. " <td> 0.260861</td>\n",
  472. " <td> L_S-M1_L_deepgray</td>\n",
  473. " <td> 0.110876</td>\n",
  474. " <td> 0.124282</td>\n",
  475. " <td> 0.175911</td>\n",
  476. " </tr>\n",
  477. " <tr>\n",
  478. " <th>6 </th>\n",
  479. " <td> 0.367417</td>\n",
  480. " <td> 0.000111</td>\n",
  481. " <td> 0.136767</td>\n",
  482. " <td> L_PFC_R_PFC</td>\n",
  483. " <td> 0.119911</td>\n",
  484. " <td> 0.211216</td>\n",
  485. " <td> 0.174641</td>\n",
  486. " </tr>\n",
  487. " <tr>\n",
  488. " <th>11</th>\n",
  489. " <td> 0.315755</td>\n",
  490. " <td> 0.000005</td>\n",
  491. " <td> 0.299819</td>\n",
  492. " <td> L_S2_L_deepgray</td>\n",
  493. " <td> 0.130187</td>\n",
  494. " <td> 0.336098</td>\n",
  495. " <td> 0.168614</td>\n",
  496. " </tr>\n",
  497. " <tr>\n",
  498. " <th>18</th>\n",
  499. " <td> 0.528601</td>\n",
  500. " <td> 0.000094</td>\n",
  501. " <td> 0.143693</td>\n",
  502. " <td> R_S-M1_R_PFC</td>\n",
  503. " <td> 0.199807</td>\n",
  504. " <td> 0.744565</td>\n",
  505. " <td> 0.163194</td>\n",
  506. " </tr>\n",
  507. " <tr>\n",
  508. " <th>20</th>\n",
  509. " <td> 0.391619</td>\n",
  510. " <td> 0.000003</td>\n",
  511. " <td> 0.244441</td>\n",
  512. " <td> R_S-M1_R_deepgray</td>\n",
  513. " <td> 0.224516</td>\n",
  514. " <td> 0.496301</td>\n",
  515. " <td> 0.163075</td>\n",
  516. " </tr>\n",
  517. " <tr>\n",
  518. " <th>4 </th>\n",
  519. " <td> 0.066396</td>\n",
  520. " <td> 0.000003</td>\n",
  521. " <td> 0.185651</td>\n",
  522. " <td> L_AS-Sens-C_R_deepgray</td>\n",
  523. " <td> 0.262431</td>\n",
  524. " <td> 0.039983</td>\n",
  525. " <td> 0.179182</td>\n",
  526. " </tr>\n",
  527. " <tr>\n",
  528. " <th>19</th>\n",
  529. " <td> 0.330611</td>\n",
  530. " <td> 0.000116</td>\n",
  531. " <td> 0.144397</td>\n",
  532. " <td> R_S-M1_R_S2</td>\n",
  533. " <td> 0.321852</td>\n",
  534. " <td> 0.361948</td>\n",
  535. " <td> 0.164666</td>\n",
  536. " </tr>\n",
  537. " <tr>\n",
  538. " <th>1 </th>\n",
  539. " <td> 0.635464</td>\n",
  540. " <td> 0.000004</td>\n",
  541. " <td> 0.242858</td>\n",
  542. " <td> L_AS-Sens-C_L_deepgray</td>\n",
  543. " <td> 0.429835</td>\n",
  544. " <td> 0.720494</td>\n",
  545. " <td> 0.157303</td>\n",
  546. " </tr>\n",
  547. " <tr>\n",
  548. " <th>13</th>\n",
  549. " <td> 0.187570</td>\n",
  550. " <td> 0.000004</td>\n",
  551. " <td> 0.173831</td>\n",
  552. " <td> R_AS-Sens-C_L_deepgray</td>\n",
  553. " <td> 0.488205</td>\n",
  554. " <td> 0.163867</td>\n",
  555. " <td> 0.165688</td>\n",
  556. " </tr>\n",
  557. " <tr>\n",
  558. " <th>22</th>\n",
  559. " <td> 0.453454</td>\n",
  560. " <td> 0.000006</td>\n",
  561. " <td> 0.240347</td>\n",
  562. " <td> R_deepgray_R_AS-Sens-C</td>\n",
  563. " <td> 0.512192</td>\n",
  564. " <td> 0.513965</td>\n",
  565. " <td> 0.157535</td>\n",
  566. " </tr>\n",
  567. " <tr>\n",
  568. " <th>24</th>\n",
  569. " <td> 0.578129</td>\n",
  570. " <td> 0.000005</td>\n",
  571. " <td> 0.265644</td>\n",
  572. " <td> R_deepgray_R_S2</td>\n",
  573. " <td> 0.581158</td>\n",
  574. " <td> 0.671305</td>\n",
  575. " <td> 0.155675</td>\n",
  576. " </tr>\n",
  577. " <tr>\n",
  578. " <th>10</th>\n",
  579. " <td> 0.086713</td>\n",
  580. " <td> 0.000008</td>\n",
  581. " <td> 0.283359</td>\n",
  582. " <td> L_S2_L_S-M1</td>\n",
  583. " <td> 0.609827</td>\n",
  584. " <td> 0.090642</td>\n",
  585. " <td> 0.169303</td>\n",
  586. " </tr>\n",
  587. " <tr>\n",
  588. " <th>23</th>\n",
  589. " <td> 0.957143</td>\n",
  590. " <td> 0.000008</td>\n",
  591. " <td> 0.238453</td>\n",
  592. " <td> R_deepgray_R_PFC</td>\n",
  593. " <td> 0.614807</td>\n",
  594. " <td> 0.849846</td>\n",
  595. " <td> 0.154907</td>\n",
  596. " </tr>\n",
  597. " <tr>\n",
  598. " <th>7 </th>\n",
  599. " <td> 0.098529</td>\n",
  600. " <td> 0.000023</td>\n",
  601. " <td> 0.156098</td>\n",
  602. " <td> L_S-M1_L_PFC</td>\n",
  603. " <td> 0.627938</td>\n",
  604. " <td> 0.103787</td>\n",
  605. " <td> 0.169381</td>\n",
  606. " </tr>\n",
  607. " <tr>\n",
  608. " <th>17</th>\n",
  609. " <td> 0.414644</td>\n",
  610. " <td> 0.000026</td>\n",
  611. " <td> 0.183286</td>\n",
  612. " <td> R_S-M1_R_AS-Sens-C</td>\n",
  613. " <td> 0.631614</td>\n",
  614. " <td> 0.489520</td>\n",
  615. " <td> 0.157911</td>\n",
  616. " </tr>\n",
  617. " <tr>\n",
  618. " <th>15</th>\n",
  619. " <td> 0.662763</td>\n",
  620. " <td> 0.000012</td>\n",
  621. " <td> 0.230427</td>\n",
  622. " <td> R_PFC_L_deepgray</td>\n",
  623. " <td> 0.665886</td>\n",
  624. " <td> 0.824459</td>\n",
  625. " <td> 0.154979</td>\n",
  626. " </tr>\n",
  627. " <tr>\n",
  628. " <th>9 </th>\n",
  629. " <td> 0.653553</td>\n",
  630. " <td> 0.000005</td>\n",
  631. " <td> 0.284665</td>\n",
  632. " <td> L_S2_L_AS-Sens-C</td>\n",
  633. " <td> 0.707206</td>\n",
  634. " <td> 0.795643</td>\n",
  635. " <td> 0.154416</td>\n",
  636. " </tr>\n",
  637. " <tr>\n",
  638. " <th>12</th>\n",
  639. " <td> 0.845060</td>\n",
  640. " <td> 0.000002</td>\n",
  641. " <td> 0.273083</td>\n",
  642. " <td> L_S2_R_AS-Sens-C</td>\n",
  643. " <td> 0.741972</td>\n",
  644. " <td> 0.626063</td>\n",
  645. " <td> 0.155587</td>\n",
  646. " </tr>\n",
  647. " <tr>\n",
  648. " <th>16</th>\n",
  649. " <td> 0.425286</td>\n",
  650. " <td> 0.000019</td>\n",
  651. " <td> 0.209314</td>\n",
  652. " <td> R_S-M1_L_S-M1</td>\n",
  653. " <td> 0.812435</td>\n",
  654. " <td> 0.438653</td>\n",
  655. " <td> 0.157009</td>\n",
  656. " </tr>\n",
  657. " <tr>\n",
  658. " <th>0 </th>\n",
  659. " <td> 0.134851</td>\n",
  660. " <td> 0.000014</td>\n",
  661. " <td> 0.218120</td>\n",
  662. " <td> L_AS-Sens-C_L_S-M1</td>\n",
  663. " <td> 0.814273</td>\n",
  664. " <td> 0.143198</td>\n",
  665. " <td> 0.165081</td>\n",
  666. " </tr>\n",
  667. " <tr>\n",
  668. " <th>21</th>\n",
  669. " <td> 0.312527</td>\n",
  670. " <td> 0.000015</td>\n",
  671. " <td> 0.201064</td>\n",
  672. " <td> R_deepgray_L_PFC</td>\n",
  673. " <td> 0.881631</td>\n",
  674. " <td> 0.368024</td>\n",
  675. " <td> 0.157931</td>\n",
  676. " </tr>\n",
  677. " <tr>\n",
  678. " <th>5 </th>\n",
  679. " <td> 0.413435</td>\n",
  680. " <td> 0.000022</td>\n",
  681. " <td> 0.222746</td>\n",
  682. " <td> L_PFC_L_deepgray</td>\n",
  683. " <td> 0.889638</td>\n",
  684. " <td> 0.494838</td>\n",
  685. " <td> 0.156065</td>\n",
  686. " </tr>\n",
  687. " <tr>\n",
  688. " <th>3 </th>\n",
  689. " <td> 0.304998</td>\n",
  690. " <td> 0.000016</td>\n",
  691. " <td> 0.189790</td>\n",
  692. " <td> L_AS-Sens-C_R_S2</td>\n",
  693. " <td> 0.938010</td>\n",
  694. " <td> 0.351768</td>\n",
  695. " <td> 0.158349</td>\n",
  696. " </tr>\n",
  697. " <tr>\n",
  698. " <th>14</th>\n",
  699. " <td> 0.169929</td>\n",
  700. " <td> 0.000037</td>\n",
  701. " <td> 0.190827</td>\n",
  702. " <td> R_AS-Sens-C_R_S2</td>\n",
  703. " <td> 0.954645</td>\n",
  704. " <td> 0.185670</td>\n",
  705. " <td> 0.163750</td>\n",
  706. " </tr>\n",
  707. " <tr>\n",
  708. " <th>2 </th>\n",
  709. " <td> 0.335984</td>\n",
  710. " <td> 0.000013</td>\n",
  711. " <td> 0.202978</td>\n",
  712. " <td> L_AS-Sens-C_R_AS-Sens-C</td>\n",
  713. " <td> 0.998496</td>\n",
  714. " <td> 0.389903</td>\n",
  715. " <td> 0.157566</td>\n",
  716. " </tr>\n",
  717. " </tbody>\n",
  718. "</table>\n",
  719. "</div>"
  720. ],
  721. "text/plain": [
  722. " Intercept age_at_visit gender[T.M] measure \\\n",
  723. "8 0.102274 0.000004 0.260861 L_S-M1_L_deepgray \n",
  724. "6 0.367417 0.000111 0.136767 L_PFC_R_PFC \n",
  725. "11 0.315755 0.000005 0.299819 L_S2_L_deepgray \n",
  726. "18 0.528601 0.000094 0.143693 R_S-M1_R_PFC \n",
  727. "20 0.391619 0.000003 0.244441 R_S-M1_R_deepgray \n",
  728. "4 0.066396 0.000003 0.185651 L_AS-Sens-C_R_deepgray \n",
  729. "19 0.330611 0.000116 0.144397 R_S-M1_R_S2 \n",
  730. "1 0.635464 0.000004 0.242858 L_AS-Sens-C_L_deepgray \n",
  731. "13 0.187570 0.000004 0.173831 R_AS-Sens-C_L_deepgray \n",
  732. "22 0.453454 0.000006 0.240347 R_deepgray_R_AS-Sens-C \n",
  733. "24 0.578129 0.000005 0.265644 R_deepgray_R_S2 \n",
  734. "10 0.086713 0.000008 0.283359 L_S2_L_S-M1 \n",
  735. "23 0.957143 0.000008 0.238453 R_deepgray_R_PFC \n",
  736. "7 0.098529 0.000023 0.156098 L_S-M1_L_PFC \n",
  737. "17 0.414644 0.000026 0.183286 R_S-M1_R_AS-Sens-C \n",
  738. "15 0.662763 0.000012 0.230427 R_PFC_L_deepgray \n",
  739. "9 0.653553 0.000005 0.284665 L_S2_L_AS-Sens-C \n",
  740. "12 0.845060 0.000002 0.273083 L_S2_R_AS-Sens-C \n",
  741. "16 0.425286 0.000019 0.209314 R_S-M1_L_S-M1 \n",
  742. "0 0.134851 0.000014 0.218120 L_AS-Sens-C_L_S-M1 \n",
  743. "21 0.312527 0.000015 0.201064 R_deepgray_L_PFC \n",
  744. "5 0.413435 0.000022 0.222746 L_PFC_L_deepgray \n",
  745. "3 0.304998 0.000016 0.189790 L_AS-Sens-C_R_S2 \n",
  746. "14 0.169929 0.000037 0.190827 R_AS-Sens-C_R_S2 \n",
  747. "2 0.335984 0.000013 0.202978 L_AS-Sens-C_R_AS-Sens-C \n",
  748. "\n",
  749. " r00_correlation r00_mt_map r2 \n",
  750. "8 0.110876 0.124282 0.175911 \n",
  751. "6 0.119911 0.211216 0.174641 \n",
  752. "11 0.130187 0.336098 0.168614 \n",
  753. "18 0.199807 0.744565 0.163194 \n",
  754. "20 0.224516 0.496301 0.163075 \n",
  755. "4 0.262431 0.039983 0.179182 \n",
  756. "19 0.321852 0.361948 0.164666 \n",
  757. "1 0.429835 0.720494 0.157303 \n",
  758. "13 0.488205 0.163867 0.165688 \n",
  759. "22 0.512192 0.513965 0.157535 \n",
  760. "24 0.581158 0.671305 0.155675 \n",
  761. "10 0.609827 0.090642 0.169303 \n",
  762. "23 0.614807 0.849846 0.154907 \n",
  763. "7 0.627938 0.103787 0.169381 \n",
  764. "17 0.631614 0.489520 0.157911 \n",
  765. "15 0.665886 0.824459 0.154979 \n",
  766. "9 0.707206 0.795643 0.154416 \n",
  767. "12 0.741972 0.626063 0.155587 \n",
  768. "16 0.812435 0.438653 0.157009 \n",
  769. "0 0.814273 0.143198 0.165081 \n",
  770. "21 0.881631 0.368024 0.157931 \n",
  771. "5 0.889638 0.494838 0.156065 \n",
  772. "3 0.938010 0.351768 0.158349 \n",
  773. "14 0.954645 0.185670 0.163750 \n",
  774. "2 0.998496 0.389903 0.157566 "
  775. ]
  776. },
  777. "execution_count": 30,
  778. "metadata": {},
  779. "output_type": "execute_result"
  780. }
  781. ],
  782. "source": [
  783. "results = pd.DataFrame()\n",
  784. "for conn, x in df.groupby(\"connection\"):\n",
  785. " model = sf.ols(\"timed_walk_trial1_time ~ r00_correlation + r00_mt_map +age_at_visit +gender\", x).fit()\n",
  786. " foo = model.pvalues\n",
  787. " foo[\"r2\"] = model.rsquared\n",
  788. " foo[\"measure\"] = conn\n",
  789. " results = results.append(foo,ignore_index=True)\n",
  790. "results.sort(\"r00_correlation\")"
  791. ]
  792. },
  793. {
  794. "cell_type": "code",
  795. "execution_count": 31,
  796. "metadata": {
  797. "collapsed": false
  798. },
  799. "outputs": [
  800. {
  801. "data": {
  802. "text/html": [
  803. "<div style=\"max-height:1000px;max-width:1500px;overflow:auto;\">\n",
  804. "<table border=\"1\" class=\"dataframe\">\n",
  805. " <thead>\n",
  806. " <tr style=\"text-align: right;\">\n",
  807. " <th></th>\n",
  808. " <th>Intercept</th>\n",
  809. " <th>age_at_visit</th>\n",
  810. " <th>gender[T.M]</th>\n",
  811. " <th>measure</th>\n",
  812. " <th>r00_correlation</th>\n",
  813. " <th>r00_mt_map</th>\n",
  814. " <th>r2</th>\n",
  815. " </tr>\n",
  816. " </thead>\n",
  817. " <tbody>\n",
  818. " <tr>\n",
  819. " <th>18</th>\n",
  820. " <td> 0.001050</td>\n",
  821. " <td> 0.027799</td>\n",
  822. " <td> 0.054800</td>\n",
  823. " <td> R_S-M1_R_PFC</td>\n",
  824. " <td> 0.042501</td>\n",
  825. " <td> 0.219691</td>\n",
  826. " <td> 0.136168</td>\n",
  827. " </tr>\n",
  828. " <tr>\n",
  829. " <th>6 </th>\n",
  830. " <td> 0.000064</td>\n",
  831. " <td> 0.023402</td>\n",
  832. " <td> 0.043774</td>\n",
  833. " <td> L_PFC_R_PFC</td>\n",
  834. " <td> 0.056876</td>\n",
  835. " <td> 0.027625</td>\n",
  836. " <td> 0.149855</td>\n",
  837. " </tr>\n",
  838. " <tr>\n",
  839. " <th>16</th>\n",
  840. " <td> 0.003051</td>\n",
  841. " <td> 0.009895</td>\n",
  842. " <td> 0.049720</td>\n",
  843. " <td> R_S-M1_L_S-M1</td>\n",
  844. " <td> 0.229823</td>\n",
  845. " <td> 0.076170</td>\n",
  846. " <td> 0.128553</td>\n",
  847. " </tr>\n",
  848. " <tr>\n",
  849. " <th>19</th>\n",
  850. " <td> 0.002014</td>\n",
  851. " <td> 0.015761</td>\n",
  852. " <td> 0.050947</td>\n",
  853. " <td> R_S-M1_R_S2</td>\n",
  854. " <td> 0.232018</td>\n",
  855. " <td> 0.163115</td>\n",
  856. " <td> 0.124931</td>\n",
  857. " </tr>\n",
  858. " <tr>\n",
  859. " <th>17</th>\n",
  860. " <td> 0.000062</td>\n",
  861. " <td> 0.014243</td>\n",
  862. " <td> 0.051321</td>\n",
  863. " <td> R_S-M1_R_AS-Sens-C</td>\n",
  864. " <td> 0.237646</td>\n",
  865. " <td> 0.067674</td>\n",
  866. " <td> 0.133613</td>\n",
  867. " </tr>\n",
  868. " <tr>\n",
  869. " <th>15</th>\n",
  870. " <td> 0.012735</td>\n",
  871. " <td> 0.002615</td>\n",
  872. " <td> 0.018871</td>\n",
  873. " <td> R_PFC_L_deepgray</td>\n",
  874. " <td> 0.376680</td>\n",
  875. " <td> 0.820765</td>\n",
  876. " <td> 0.108716</td>\n",
  877. " </tr>\n",
  878. " <tr>\n",
  879. " <th>3 </th>\n",
  880. " <td> 0.000012</td>\n",
  881. " <td> 0.008224</td>\n",
  882. " <td> 0.054157</td>\n",
  883. " <td> L_AS-Sens-C_R_S2</td>\n",
  884. " <td> 0.412252</td>\n",
  885. " <td> 0.061450</td>\n",
  886. " <td> 0.128726</td>\n",
  887. " </tr>\n",
  888. " <tr>\n",
  889. " <th>7 </th>\n",
  890. " <td> 0.000000</td>\n",
  891. " <td> 0.005829</td>\n",
  892. " <td> 0.044702</td>\n",
  893. " <td> L_S-M1_L_PFC</td>\n",
  894. " <td> 0.451307</td>\n",
  895. " <td> 0.006729</td>\n",
  896. " <td> 0.148246</td>\n",
  897. " </tr>\n",
  898. " <tr>\n",
  899. " <th>2 </th>\n",
  900. " <td> 0.000037</td>\n",
  901. " <td> 0.005917</td>\n",
  902. " <td> 0.044495</td>\n",
  903. " <td> L_AS-Sens-C_R_AS-Sens-C</td>\n",
  904. " <td> 0.462218</td>\n",
  905. " <td> 0.106427</td>\n",
  906. " <td> 0.123201</td>\n",
  907. " </tr>\n",
  908. " <tr>\n",
  909. " <th>23</th>\n",
  910. " <td> 0.000828</td>\n",
  911. " <td> 0.006420</td>\n",
  912. " <td> 0.021273</td>\n",
  913. " <td> R_deepgray_R_PFC</td>\n",
  914. " <td> 0.521677</td>\n",
  915. " <td> 0.171136</td>\n",
  916. " <td> 0.118901</td>\n",
  917. " </tr>\n",
  918. " <tr>\n",
  919. " <th>14</th>\n",
  920. " <td> 0.000064</td>\n",
  921. " <td> 0.011467</td>\n",
  922. " <td> 0.038808</td>\n",
  923. " <td> R_AS-Sens-C_R_S2</td>\n",
  924. " <td> 0.588705</td>\n",
  925. " <td> 0.041388</td>\n",
  926. " <td> 0.132721</td>\n",
  927. " </tr>\n",
  928. " <tr>\n",
  929. " <th>12</th>\n",
  930. " <td> 0.016344</td>\n",
  931. " <td> 0.002152</td>\n",
  932. " <td> 0.029329</td>\n",
  933. " <td> L_S2_R_AS-Sens-C</td>\n",
  934. " <td> 0.595380</td>\n",
  935. " <td> 0.963862</td>\n",
  936. " <td> 0.105163</td>\n",
  937. " </tr>\n",
  938. " <tr>\n",
  939. " <th>8 </th>\n",
  940. " <td> 0.000037</td>\n",
  941. " <td> 0.004482</td>\n",
  942. " <td> 0.021598</td>\n",
  943. " <td> L_S-M1_L_deepgray</td>\n",
  944. " <td> 0.603762</td>\n",
  945. " <td> 0.015219</td>\n",
  946. " <td> 0.136280</td>\n",
  947. " </tr>\n",
  948. " <tr>\n",
  949. " <th>5 </th>\n",
  950. " <td> 0.000197</td>\n",
  951. " <td> 0.012670</td>\n",
  952. " <td> 0.019629</td>\n",
  953. " <td> L_PFC_L_deepgray</td>\n",
  954. " <td> 0.657948</td>\n",
  955. " <td> 0.063127</td>\n",
  956. " <td> 0.124463</td>\n",
  957. " </tr>\n",
  958. " <tr>\n",
  959. " <th>21</th>\n",
  960. " <td> 0.000007</td>\n",
  961. " <td> 0.009961</td>\n",
  962. " <td> 0.028479</td>\n",
  963. " <td> R_deepgray_L_PFC</td>\n",
  964. " <td> 0.681705</td>\n",
  965. " <td> 0.023235</td>\n",
  966. " <td> 0.133790</td>\n",
  967. " </tr>\n",
  968. " <tr>\n",
  969. " <th>9 </th>\n",
  970. " <td> 0.002011</td>\n",
  971. " <td> 0.004132</td>\n",
  972. " <td> 0.025987</td>\n",
  973. " <td> L_S2_L_AS-Sens-C</td>\n",
  974. " <td> 0.716290</td>\n",
  975. " <td> 0.289003</td>\n",
  976. " <td> 0.112135</td>\n",
  977. " </tr>\n",
  978. " <tr>\n",
  979. " <th>13</th>\n",
  980. " <td> 0.000002</td>\n",
  981. " <td> 0.003724</td>\n",
  982. " <td> 0.038902</td>\n",
  983. " <td> R_AS-Sens-C_L_deepgray</td>\n",
  984. " <td> 0.730586</td>\n",
  985. " <td> 0.020673</td>\n",
  986. " <td> 0.134931</td>\n",
  987. " </tr>\n",
  988. " <tr>\n",
  989. " <th>22</th>\n",
  990. " <td> 0.000032</td>\n",
  991. " <td> 0.007084</td>\n",
  992. " <td> 0.017682</td>\n",
  993. " <td> R_deepgray_R_AS-Sens-C</td>\n",
  994. " <td> 0.780205</td>\n",
  995. " <td> 0.066916</td>\n",
  996. " <td> 0.123963</td>\n",
  997. " </tr>\n",
  998. " <tr>\n",
  999. " <th>11</th>\n",
  1000. " <td> 0.000081</td>\n",
  1001. " <td> 0.009850</td>\n",
  1002. " <td> 0.015995</td>\n",
  1003. " <td> L_S2_L_deepgray</td>\n",
  1004. " <td> 0.796521</td>\n",
  1005. " <td> 0.031031</td>\n",
  1006. " <td> 0.129467</td>\n",
  1007. " </tr>\n",
  1008. " <tr>\n",
  1009. " <th>1 </th>\n",
  1010. " <td> 0.001578</td>\n",
  1011. " <td> 0.004778</td>\n",
  1012. " <td> 0.017942</td>\n",
  1013. " <td> L_AS-Sens-C_L_deepgray</td>\n",
  1014. " <td> 0.813290</td>\n",
  1015. " <td> 0.215427</td>\n",
  1016. " <td> 0.112836</td>\n",
  1017. " </tr>\n",
  1018. " <tr>\n",
  1019. " <th>24</th>\n",
  1020. " <td> 0.002798</td>\n",
  1021. " <td> 0.002842</td>\n",
  1022. " <td> 0.017301</td>\n",
  1023. " <td> R_deepgray_R_S2</td>\n",
  1024. " <td> 0.825346</td>\n",
  1025. " <td> 0.203149</td>\n",
  1026. " <td> 0.112634</td>\n",
  1027. " </tr>\n",
  1028. " <tr>\n",
  1029. " <th>0 </th>\n",
  1030. " <td> 0.000000</td>\n",
  1031. " <td> 0.008066</td>\n",
  1032. " <td> 0.036150</td>\n",
  1033. " <td> L_AS-Sens-C_L_S-M1</td>\n",
  1034. " <td> 0.828228</td>\n",
  1035. " <td> 0.001947</td>\n",
  1036. " <td> 0.159483</td>\n",
  1037. " </tr>\n",
  1038. " <tr>\n",
  1039. " <th>4 </th>\n",
  1040. " <td> 0.000001</td>\n",
  1041. " <td> 0.002742</td>\n",
  1042. " <td> 0.025568</td>\n",
  1043. " <td> L_AS-Sens-C_R_deepgray</td>\n",
  1044. " <td> 0.922103</td>\n",
  1045. " <td> 0.009649</td>\n",
  1046. " <td> 0.141273</td>\n",
  1047. " </tr>\n",
  1048. " <tr>\n",
  1049. " <th>10</th>\n",
  1050. " <td> 0.000002</td>\n",
  1051. " <td> 0.003864</td>\n",
  1052. " <td> 0.021712</td>\n",
  1053. " <td> L_S2_L_S-M1</td>\n",
  1054. " <td> 0.988556</td>\n",
  1055. " <td> 0.005241</td>\n",
  1056. " <td> 0.147528</td>\n",
  1057. " </tr>\n",
  1058. " <tr>\n",
  1059. " <th>20</th>\n",
  1060. " <td> 0.000133</td>\n",
  1061. " <td> 0.004122</td>\n",
  1062. " <td> 0.033329</td>\n",
  1063. " <td> R_S-M1_R_deepgray</td>\n",
  1064. " <td> 0.998675</td>\n",
  1065. " <td> 0.084061</td>\n",
  1066. " <td> 0.120437</td>\n",
  1067. " </tr>\n",
  1068. " </tbody>\n",
  1069. "</table>\n",
  1070. "</div>"
  1071. ],
  1072. "text/plain": [
  1073. " Intercept age_at_visit gender[T.M] measure \\\n",
  1074. "18 0.001050 0.027799 0.054800 R_S-M1_R_PFC \n",
  1075. "6 0.000064 0.023402 0.043774 L_PFC_R_PFC \n",
  1076. "16 0.003051 0.009895 0.049720 R_S-M1_L_S-M1 \n",
  1077. "19 0.002014 0.015761 0.050947 R_S-M1_R_S2 \n",
  1078. "17 0.000062 0.014243 0.051321 R_S-M1_R_AS-Sens-C \n",
  1079. "15 0.012735 0.002615 0.018871 R_PFC_L_deepgray \n",
  1080. "3 0.000012 0.008224 0.054157 L_AS-Sens-C_R_S2 \n",
  1081. "7 0.000000 0.005829 0.044702 L_S-M1_L_PFC \n",
  1082. "2 0.000037 0.005917 0.044495 L_AS-Sens-C_R_AS-Sens-C \n",
  1083. "23 0.000828 0.006420 0.021273 R_deepgray_R_PFC \n",
  1084. "14 0.000064 0.011467 0.038808 R_AS-Sens-C_R_S2 \n",
  1085. "12 0.016344 0.002152 0.029329 L_S2_R_AS-Sens-C \n",
  1086. "8 0.000037 0.004482 0.021598 L_S-M1_L_deepgray \n",
  1087. "5 0.000197 0.012670 0.019629 L_PFC_L_deepgray \n",
  1088. "21 0.000007 0.009961 0.028479 R_deepgray_L_PFC \n",
  1089. "9 0.002011 0.004132 0.025987 L_S2_L_AS-Sens-C \n",
  1090. "13 0.000002 0.003724 0.038902 R_AS-Sens-C_L_deepgray \n",
  1091. "22 0.000032 0.007084 0.017682 R_deepgray_R_AS-Sens-C \n",
  1092. "11 0.000081 0.009850 0.015995 L_S2_L_deepgray \n",
  1093. "1 0.001578 0.004778 0.017942 L_AS-Sens-C_L_deepgray \n",
  1094. "24 0.002798 0.002842 0.017301 R_deepgray_R_S2 \n",
  1095. "0 0.000000 0.008066 0.036150 L_AS-Sens-C_L_S-M1 \n",
  1096. "4 0.000001 0.002742 0.025568 L_AS-Sens-C_R_deepgray \n",
  1097. "10 0.000002 0.003864 0.021712 L_S2_L_S-M1 \n",
  1098. "20 0.000133 0.004122 0.033329 R_S-M1_R_deepgray \n",
  1099. "\n",
  1100. " r00_correlation r00_mt_map r2 \n",
  1101. "18 0.042501 0.219691 0.136168 \n",
  1102. "6 0.056876 0.027625 0.149855 \n",
  1103. "16 0.229823 0.076170 0.128553 \n",
  1104. "19 0.232018 0.163115 0.124931 \n",
  1105. "17 0.237646 0.067674 0.133613 \n",
  1106. "15 0.376680 0.820765 0.108716 \n",
  1107. "3 0.412252 0.061450 0.128726 \n",
  1108. "7 0.451307 0.006729 0.148246 \n",
  1109. "2 0.462218 0.106427 0.123201 \n",
  1110. "23 0.521677 0.171136 0.118901 \n",
  1111. "14 0.588705 0.041388 0.132721 \n",
  1112. "12 0.595380 0.963862 0.105163 \n",
  1113. "8 0.603762 0.015219 0.136280 \n",
  1114. "5 0.657948 0.063127 0.124463 \n",
  1115. "21 0.681705 0.023235 0.133790 \n",
  1116. "9 0.716290 0.289003 0.112135 \n",
  1117. "13 0.730586 0.020673 0.134931 \n",
  1118. "22 0.780205 0.066916 0.123963 \n",
  1119. "11 0.796521 0.031031 0.129467 \n",
  1120. "1 0.813290 0.215427 0.112836 \n",
  1121. "24 0.825346 0.203149 0.112634 \n",
  1122. "0 0.828228 0.001947 0.159483 \n",
  1123. "4 0.922103 0.009649 0.141273 \n",
  1124. "10 0.988556 0.005241 0.147528 \n",
  1125. "20 0.998675 0.084061 0.120437 "
  1126. ]
  1127. },
  1128. "execution_count": 31,
  1129. "metadata": {},
  1130. "output_type": "execute_result"
  1131. }
  1132. ],
  1133. "source": [
  1134. "results = pd.DataFrame()\n",
  1135. "for conn, x in df.groupby(\"connection\"):\n",
  1136. " model = sf.ols(\"nhpt_dominant_time ~ r00_correlation + r00_mt_map +age_at_visit +gender\", x).fit()\n",
  1137. " foo = model.pvalues\n",
  1138. " foo[\"r2\"] = model.rsquared\n",
  1139. " foo[\"measure\"] = conn\n",
  1140. " results = results.append(foo,ignore_index=True)\n",
  1141. "results.sort(\"r00_correlation\")"
  1142. ]
  1143. },
  1144. {
  1145. "cell_type": "code",
  1146. "execution_count": null,
  1147. "metadata": {
  1148. "collapsed": true
  1149. },
  1150. "outputs": [],
  1151. "source": []
  1152. }
  1153. ],
  1154. "metadata": {
  1155. "kernelspec": {
  1156. "display_name": "Python 2",
  1157. "language": "python",
  1158. "name": "python2"
  1159. },
  1160. "language_info": {
  1161. "codemirror_mode": {
  1162. "name": "ipython",
  1163. "version": 2
  1164. },
  1165. "file_extension": ".py",
  1166. "mimetype": "text/x-python",
  1167. "name": "python",
  1168. "nbconvert_exporter": "python",
  1169. "pygments_lexer": "ipython2",
  1170. "version": "2.7.11"
  1171. }
  1172. },
  1173. "nbformat": 4,
  1174. "nbformat_minor": 0
  1175. }
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement