Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

The smallest value of ICVF always larger than 0.1 #97

Open
Shintai516 opened this issue Apr 14, 2020 · 2 comments
Open

The smallest value of ICVF always larger than 0.1 #97

Shintai516 opened this issue Apr 14, 2020 · 2 comments

Comments

@Shintai516
Copy link

Hello,

I was trying to fit the NODDI model with AMICO in tumor patients with edema.
In the results of the model fitting, I found that most of the voxels with edema were the same value in ICVF which was 0.1.
I've tried three different datasets and came-out with the same result.
Screen Shot 2020-04-14 at 7 01 29 PM
This is the histogram of FIT_ICVF of three different datasets. The smallest value(zero excluded) was always equal to 0.1.
Screen Shot 2020-04-14 at 7 04 31 PM
This is the same data but fitting with AMICO(left) and NODDI MATLAB toolbox(right).
Is there any clue or solution to fix this problem?

Thanks!

Cheers,
Shin Tai

@daducci
Copy link
Owner

daducci commented May 5, 2020

Hi @Shintai516 ,

sorry for the delay. To create the linear operator AMICO discretizes the parameters like ICVF and OD. When we developed AMICO the aim was to reproduce the results of the original implementation and we set those parameters to given values (see

self.IC_VFs = np.linspace(0.1,0.99,12)
). You can modify those parameters using the set() function (
def set( self, dPar, dIso, IC_VFs, IC_ODs, isExvivo ):
) to suit your needs. Let us know if this solves your problem.

@Shintai516
Copy link
Author

Thanks, @daducci
This helps a lot.
I will try to modify the parameter and run the function again.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants