Plot ROIs from PA2 StructureFunction notebook

Seed ROIs

plotting.plot_roi('/fmriDataRaw/fmri_data_raw/DSI_test/AAL626_regions/AAL626_final_195.nii')

png

plotting.plot_roi('/fmriDataRaw/fmri_data_raw/DSI_test/AAL626_regions/AAL626_final_198.nii')

<nilearn.plotting.displays.OrthoSlicer at 0x7ff4e9aeb9b0>

png

plotting.plot_roi('/fmriDataRaw/fmri_data_raw/DSI_test/AAL626_regions/AAL626_final_200.nii')

<nilearn.plotting.displays.OrthoSlicer at 0x7ff4e98a7c88>

png

plotting.plot_roi('/fmriDataRaw/fmri_data_raw/DSI_test/AAL626_regions/AAL626_final_202.nii')

<nilearn.plotting.displays.OrthoSlicer at 0x7ff4e9665f60>

png

Predicting sedentary change

ROI 195, sedentary

stepAIC reduced model does not contain PPI terms

ROI 198, sedentary

Reduced model with PPI terms only is significantly predictive of behavior change

roilist=[199, 237, 286, 74, 76, 79]
roipathlist=[roidir+'AAL626_final_'+str(x)+'.nii' for x in roilist]
roi_imgs = image.concat_imgs(roipathlist)
plotting.plot_prob_atlas(roi_imgs, display_mode='z',colorbar=True)

/usr/local/anaconda3/lib/python3.6/site-packages/numpy/ma/core.py:2766: UserWarning: Warning: converting a masked element to nan.
  order=order, subok=True, ndmin=ndmin)





<nilearn.plotting.displays.ZSlicer at 0x7ff4f1a39be0>

png

plotting.plot_prob_atlas(roi_imgs, display_mode='y')

/usr/local/anaconda3/lib/python3.6/site-packages/numpy/ma/core.py:2766: UserWarning: Warning: converting a masked element to nan.
  order=order, subok=True, ndmin=ndmin)





<nilearn.plotting.displays.YSlicer at 0x7ff4f20192e8>

png

ROI 200, sedentary

Reduced model with PPI terms only is significantly predictive of behavior change

roilist=[282, 77, 78]
roipathlist=[roidir+'AAL626_final_'+str(x)+'.nii' for x in roilist]
roi_imgs = image.concat_imgs(roipathlist)
plotting.plot_prob_atlas(roi_imgs, display_mode='z',colorbar=True)

/usr/local/anaconda3/lib/python3.6/site-packages/numpy/ma/core.py:2766: UserWarning: Warning: converting a masked element to nan.
  order=order, subok=True, ndmin=ndmin)





<nilearn.plotting.displays.ZSlicer at 0x7ff4f1467c88>

png

ROI 202, sedentary

Reduced model with PPI terms only is significantly predictive of behavior change

roilist=[282, 314,81,84]
roipathlist=[roidir+'AAL626_final_'+str(x)+'.nii' for x in roilist]
roi_imgs = image.concat_imgs(roipathlist)
plotting.plot_prob_atlas(roi_imgs, display_mode='z',colorbar=True)

/usr/local/anaconda3/lib/python3.6/site-packages/numpy/ma/core.py:2766: UserWarning: Warning: converting a masked element to nan.
  order=order, subok=True, ndmin=ndmin)





<nilearn.plotting.displays.ZSlicer at 0x7ff4f0f3c7f0>

png

282 is the only ROI to come up in more than one list

plotting.plot_roi('/fmriDataRaw/fmri_data_raw/DSI_test/AAL626_regions/AAL626_final_282.nii')

<nilearn.plotting.displays.OrthoSlicer at 0x7ff4f84657b8>

png

Predicting mod-vig change

ROI 195, mod-vig

PPI terms alone are not a significantly predictive model; reduced model with covariates is

roilist=[14,22,38,39,40,120,500]
roipathlist=[roidir+'AAL626_final_'+str(x)+'.nii' for x in roilist]
roi_imgs = image.concat_imgs(roipathlist)
plotting.plot_prob_atlas(roi_imgs, display_mode='z',colorbar=True)

/usr/local/anaconda3/lib/python3.6/site-packages/numpy/ma/core.py:2766: UserWarning: Warning: converting a masked element to nan.
  order=order, subok=True, ndmin=ndmin)





<nilearn.plotting.displays.ZSlicer at 0x7ff4f0a9a908>

png

plotting.plot_prob_atlas(roi_imgs, display_mode='x',colorbar=True)

/usr/local/anaconda3/lib/python3.6/site-packages/numpy/ma/core.py:2766: UserWarning: Warning: converting a masked element to nan.
  order=order, subok=True, ndmin=ndmin)





<nilearn.plotting.displays.XSlicer at 0x7ff4ef875be0>

png

ROI 198, mod-vig

Reduced models are not significantly predictive

roilist=[14,75, 76, 79, 80, 85, 222, 286]
roipathlist=[roidir+'AAL626_final_'+str(x)+'.nii' for x in roilist]
roi_imgs = image.concat_imgs(roipathlist)
plotting.plot_prob_atlas(roi_imgs, display_mode='z',colorbar=True)

/usr/local/anaconda3/lib/python3.6/site-packages/numpy/ma/core.py:2766: UserWarning: Warning: converting a masked element to nan.
  order=order, subok=True, ndmin=ndmin)





<nilearn.plotting.displays.ZSlicer at 0x7ff4ef3b3828>

png

plotting.plot_prob_atlas(roi_imgs, display_mode='x',colorbar=True)

/usr/local/anaconda3/lib/python3.6/site-packages/numpy/ma/core.py:2766: UserWarning: Warning: converting a masked element to nan.
  order=order, subok=True, ndmin=ndmin)





<nilearn.plotting.displays.XSlicer at 0x7ff4ee4cb6d8>

png

ROI 14 shows up twice

plotting.plot_roi('/fmriDataRaw/fmri_data_raw/DSI_test/AAL626_regions/AAL626_final_14.nii')

---------------------------------------------------------------------------

NameError                                 Traceback (most recent call last)

<ipython-input-2-d6e0440484b8> in <module>()
----> 1 plotting.plot_roi('/fmriDataRaw/fmri_data_raw/DSI_test/AAL626_regions/AAL626_final_14.nii')


NameError: name 'plotting' is not defined

ROI 200, mod-vig

Reduced model (PPI terms only) is significantly predictive

roilist=[275,302,68,77]
roipathlist=[roidir+'AAL626_final_'+str(x)+'.nii' for x in roilist]
roi_imgs = image.concat_imgs(roipathlist)
plotting.plot_prob_atlas(roi_imgs, display_mode='z',colorbar=True)

/usr/local/anaconda3/lib/python3.6/site-packages/numpy/ma/core.py:2766: UserWarning: Warning: converting a masked element to nan.
  order=order, subok=True, ndmin=ndmin)





<nilearn.plotting.displays.ZSlicer at 0x7ff4ee007ef0>

png

plotting.plot_prob_atlas(roi_imgs, display_mode='x',colorbar=True)

/usr/local/anaconda3/lib/python3.6/site-packages/numpy/ma/core.py:2766: UserWarning: Warning: converting a masked element to nan.
  order=order, subok=True, ndmin=ndmin)





<nilearn.plotting.displays.XSlicer at 0x7ff4edb60588>

png

ROI 202, mod-vig

Reduced model (PPI terms only) is significantly predictive

roilist=[240,245,253,286,289,297,205,307,49,81,93,95,97]
roipathlist=[roidir+'AAL626_final_'+str(x)+'.nii' for x in roilist]
roi_imgs = image.concat_imgs(roipathlist)
plotting.plot_prob_atlas(roi_imgs, display_mode='z',colorbar=True)

/usr/local/anaconda3/lib/python3.6/site-packages/numpy/ma/core.py:2766: UserWarning: Warning: converting a masked element to nan.
  order=order, subok=True, ndmin=ndmin)
/usr/local/anaconda3/lib/python3.6/site-packages/matplotlib/contour.py:920: UserWarning: linewidths is ignored by contourf
  warnings.warn('linewidths is ignored by contourf')





<nilearn.plotting.displays.ZSlicer at 0x7ff4ed681978>

png

plotting.plot_prob_atlas(roi_imgs, display_mode='x',colorbar=True)

---------------------------------------------------------------------------

NameError                                 Traceback (most recent call last)

<ipython-input-1-9e95a262559a> in <module>()
----> 1 plotting.plot_prob_atlas(roi_imgs, display_mode='x',colorbar=True)


NameError: name 'plotting' is not defined