History of AFNI updates  

|
February 12, 2018 01:24PM
Hi, Jun-

The total number of degrees of freedom (DOFs) of your data set initially is your number of points. With every regressor in your model or time point that you censor, one degree of freedom gets used up. 3dDeconvolve is the main function where the regression model gets applied. There are some basic things that get regressed: a few polynomials (your "polort" number above is 4), motion regressors (6 rigid body parameters, possibly also the derivative of each, which is common in resting state analysis, so that would be 12 motion regressors total). 3 initial time points get removed-- so that is three more DOFs gone-- and then it looks like you have *a lot* of time points censored-- see the line of output:
Number of time points: 237 (before censor) ; 111 (after)
which means that you have 126 points censored-- over half of the data set, which should send some warning chills. And finally, doing bandpassing is *also* something that eats up degrees of freedom-- in fact, 2 DOF get regressed out for each frequency removed (and exactly how many in total get removed depends on your TR).

In total, by the time the polorts, motion regressors (+derivatives), chopped initial time points, censored time points, and bandpassed frequencies get put into the regression model, the program detects that you are trying to remove more degrees of freedom than what you have to start with. Therein lies the problem.

If you look at your volumes and time series over time, does it look problematic/noisy? That seems to be the part causing the most trouble here.

--pt
Subject Author Posted

Resting State proc.py

JW February 12, 2018 01:03PM

Re: Resting State proc.py

rick reynolds February 12, 2018 01:11PM

Re: Resting State proc.py

ptaylor February 12, 2018 01:24PM

Re: Resting State proc.py

JW February 13, 2018 09:12AM