Yes, you are right. Your suggestion sounds a great option! But I am still not sure how to do this. It is said in the 3dDeconvolve documentation that"3dDeconvolve does LINEAR regression'', and " If you want NONLINEAR regression, see program 3dNLfim ". Would you please provide a sample about how to achieve this? Where should add the quadratic modulators in the script please? I did see an option "-legendre " but it is for baseline model.