def test_coefficient_of_determination(): """ Test the computation of coefficient of determination """ cod = ozm.coeff_of_determination(TM1, TM2) # Test for regressions: npt.assert_almost_equal(cod, np.array([[[ -3.44661890e+00, -9.56966221e-01, -2.95984527e-01, -2.14809778e-01], [ -1.13823857e+00, -1.47624354e+00, -2.49836988e-01, 4.45058522e-01]], [[ -2.12348903e+00, -1.03695127e+00, 3.43880861e-03, -5.08955429e-01], [ -2.70970026e+00, -8.62731412e-01, 3.21255708e-01, 1.52058544e-01]], [[ -1.54499435e+00, -1.12129147e+00, -1.34573166e+00, -5.70547139e-01], [ -1.31661328e+00, -1.31546355e-02, -9.43582307e-03, -4.12026495e-01]]]))
help='Bvals file (FSL format)') parser.add_argument('out_file', action='store', metavar='File', help='Output file name (.nii.gz)') parser.add_argument('--mask_file', action='store', metavar='File', help='Mask file (only the voxels within the binary mask will be analyzed (.nii.gz; default: analyze all) ', default=None) params = parser.parse_args() if __name__ == "__main__": Model1 = dti.TensorModel(params.dwi_file1, params.bvecs_file1, params.bvals_file1, mask=params.mask_file, params_file='temp') Model2 = dti.TensorModel(params.dwi_file2, params.bvecs_file2, params.bvals_file2, mask=params.mask_file, params_file='temp') # Do it and save: nib.Nifti1Image(ana.coeff_of_determination(Model1, Model2), Model1.affine).to_filename(params.out_file)
metavar='File', help='Output file name (.nii.gz)') parser.add_argument( '--mask_file', action='store', metavar='File', help= 'Mask file (only the voxels within the binary mask will be analyzed (.nii.gz; default: analyze all) ', default=None) params = parser.parse_args() if __name__ == "__main__": Model1 = dti.TensorModel(params.dwi_file1, params.bvecs_file1, params.bvals_file1, mask=params.mask_file, params_file='temp') Model2 = dti.TensorModel(params.dwi_file2, params.bvecs_file2, params.bvals_file2, mask=params.mask_file, params_file='temp') # Do it and save: nib.Nifti1Image(ana.coeff_of_determination(Model1, Model2), Model1.affine).to_filename(params.out_file)