I would instead implement a loop inside of GaussianModel.fit to iterate
over voxels, fit the mode for each one of them and then return an object
that already has arrays of parameters for each voxel in the mask.


On Mon, Feb 24, 2014 at 2:25 PM, Kevin DeSimone <notifications@github.com>wrote:

> In popeye/prf.py:
>
> > +    """
> > +
> > +    def __init__(self, stimulus):
> > +
> > +        # this is a weird notation
> > +        PopulationModel.__init__(self, stimulus)
> > +
> > +    def fit(self, data, bounds, error_function, norm_func=utils.zscore):
> > +        """ Fit method of the DTI model class
> > +
> > +        Parameters
> > +        ----------
> > +        data : array
> > +            The measured signal from one voxel.
> > +
> > +        mask : array
>
> Yes, I took the mask out b/c as is the prf classes only handle a single
> voxel at a time. The mask would be sort of useless.
>
> BUT, this does speak to the larger issue of how to use the prf module with
> multiple voxels of data that need to be fit. I assume each voxel needs to
> have its own instantiation of GaussianFit(), which would render a mask
> useless within that class, no?
>
> —
> Reply to this email directly or view it on GitHub<https://github.com/kdesimone/popeye/pull/8/files#r10014242>
> .
>