Mathew,
Sorry to keep hassling you like this, but I found the solution to the
error i described to you this morning, assuming the relevant SPM.mat is
already loaded into the current workspace the instruction sequence is:
boxdef = struct('centre', [x y z], 'widths', [x y z]);
roi = maroi_box(boxdef);
y = getdata(roi, SPM.xY.VY);
Which as you suggested loads the data into a [number of images in
design] x [number of
voxels in box] matrix. This avoids the problems associated with trying
to get getdata() to take the mardo_5 object.
Thanks again
Lee.
Matthew Brett wrote:
> Hi,
>
>> I currently have a problem where I am trying to determine the mean
>> percent signal change of the voxels within a box shaped volume. These
>> voxels may or may not be above threshold in the statistical parametric
>> map. I can't find any way to do this in SPM5. Generally speaking I can
>> only get SPM5 to give information on one voxel at a time and never about
>> a non-significant voxel.
>
> You can get the raw data, and the mean of the data, using marsbar.
> See: http://marsbar.sourceforge.net/
>
> boxdef = struct('centre', [0 0 0], 'widths', [10 10 10]);
> roi = maroi_box(boxdef);
> D = mardo('SPM.mat');
> y = getdata(roi, D);
>
> That gets the data into a [number of images in design] x [number of
> voxels in box] matrix.
>
> Best,
>
> Matthew
>
>
--
Lee Walsh
Prince of Wales Medical Research Institute
Sydney, Australia
Ph: +61 2 9399 1014
Fax: +61 2 9399 1027
Email: [log in to unmask]
|