[Freesurfer] Fake thickness map 2

2007-08-06 Thread Inês Souta
-- Forwarded message --
From: Bruce Fischl <[EMAIL PROTECTED]>
Date: 06/08/2007 14:23
Subject: Re: [Freesurfer] Fake thickness map 2
To: Inês Souta <[EMAIL PROTECTED]>

they don't matter. For a real volume they would be the ras2vox transform
and the mr parameters. For your purposes they can be any 4x4 and you can
leave the mri_parms one out (it will detect and fill it in with zeros).

Send these to the list so other people can answer!

Bruce
On
Mon, 6 Aug 2007, Inês Souta wrote:

> Hi again,
>
> I still don't understand what all the input objects of save_mgh function
> stand for. 'vol' is the vector of thickness values, 'fname' is the name of
> the output .mgh file, and what about M and ? What do they stand
> for and how should I define them?
>
> Thanks again,
>
> Ines
>
> 2007/7/31, Bruce Fischl <[EMAIL PROTECTED]>:
>>
>> Hi Ines,
>>
>> use the mgh route. The "vol" is actually the vector of thickness values,
>> with one entry per vertex on the surface.
>>
>> cheers,
>> Bruce
>> On Tue, 31 Jul 2007, Inês Souta
>> wrote:
>>
>>> Hello,
>>>
>>> I was trying to create a cortical thickness map with constant values in
>>> order to perform some tests. At first I tried using mris_convert. I
>>> converted my thickness file into ascii:
>>>
>>> mris_convert lh.thickness lh.thickess.asc
>>>
>>> then I used matlab to change all the values in the last column to 2,0. I
>>> tried to convert the ascii file back to surface file by doing
>>>
>>> mris_convert lh.thicness.asc lh.thickness
>>>
>>> but I got "Segmentation fault". How can I convert may ascii file into a
>>> surface one?
>>>
>>> After failling to acomplish this, I tried using the save_mgh and
>> load_mgh
>>> file, nevertheless I'm not quite sure which input values are needed.
>> Where
>>> can I put my fake values in the save_mgh file? And what is "vol"?
>>>
>>> Thanks,
>>>
>>> Ines
>>>
>
___
Freesurfer mailing list
Freesurfer@nmr.mgh.harvard.edu
https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer

Re: [Freesurfer] Fake thickness map 2

2007-08-06 Thread Inês Souta
freesurfer v.3.0.5 for Linux-Centos4


2007/8/6, Bruce Fischl <[EMAIL PROTECTED]>:
>
> yes, you should be able to open it in tksurfer. What version of
> freesurfer do you have?
> On Mon, 6 Aug 2007, Inês Souta wrote:
>
> > And after creating the mgh volume, how can I read it? Is it possible to
> open
> > it in tksurfer? Should I load it as an overlay?
> > I've used one of my already pre-processed subjects, created a 'vol' with
> the
> > same number of vertices of one of the hemisferes, and saved the output
> of
> > save_mgh file as an rh.thickness.mgh file. However, when I don't seem to
> be
> > able to open it in tksurfer...
> >
> > Thanks
> >
> > Ines
> >
> > 2007/8/6, Bruce Fischl <[EMAIL PROTECTED]>:
> >>
> >> they don't matter. For a real volume they would be the ras2vox
> transform
> >> and the mr parameters. For your purposes they can be any 4x4 and you
> can
> >> leave the mri_parms one out (it will detect and fill it in with zeros).
> >>
> >> Send these to the list so other people can answer!
> >>
> >> Bruce
> >> On
> >> Mon, 6 Aug 2007, Inês Souta wrote:
> >>
> >>> Hi again,
> >>>
> >>> I still don't understand what all the input objects of save_mgh
> function
> >>> stand for. 'vol' is the vector of thickness values, 'fname' is the
> name
> >> of
> >>> the output .mgh file, and what about M and ? What do they
> >> stand
> >>> for and how should I define them?
> >>>
> >>> Thanks again,
> >>>
> >>> Ines
> >>>
> >>> 2007/7/31, Bruce Fischl <[EMAIL PROTECTED]>:
> 
>  Hi Ines,
> 
>  use the mgh route. The "vol" is actually the vector of thickness
> >> values,
>  with one entry per vertex on the surface.
> 
>  cheers,
>  Bruce
>  On Tue, 31 Jul 2007, Inês Souta
>  wrote:
> 
> > Hello,
> >
> > I was trying to create a cortical thickness map with constant values
> >> in
> > order to perform some tests. At first I tried using mris_convert. I
> > converted my thickness file into ascii:
> >
> > mris_convert lh.thickness lh.thickess.asc
> >
> > then I used matlab to change all the values in the last column to
> 2,0.
> >> I
> > tried to convert the ascii file back to surface file by doing
> >
> > mris_convert lh.thicness.asc lh.thickness
> >
> > but I got "Segmentation fault". How can I convert may ascii file
> into
> >> a
> > surface one?
> >
> > After failling to acomplish this, I tried using the save_mgh and
>  load_mgh
> > file, nevertheless I'm not quite sure which input values are needed.
>  Where
> > can I put my fake values in the save_mgh file? And what is "vol"?
> >
> > Thanks,
> >
> > Ines
> >
> >>>
> >
___
Freesurfer mailing list
Freesurfer@nmr.mgh.harvard.edu
https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer

Re: [Freesurfer] Fake thickness map 2

2007-08-06 Thread Inês Souta
Sorry, I was saving it with the wrong extension. I am already able to open
it in tksurfer! Thanks a lot,

Ines

2007/8/6, Inês Souta <[EMAIL PROTECTED]>:
>
> freesurfer v.3.0.5 for Linux-Centos4
>
>
> 2007/8/6, Bruce Fischl <[EMAIL PROTECTED]>:
> >
> > yes, you should be able to open it in tksurfer. What version of
> > freesurfer do you have?
> > On Mon, 6 Aug 2007, Inês Souta wrote:
> >
> > > And after creating the mgh volume, how can I read it? Is it possible
> > to open
> > > it in tksurfer? Should I load it as an overlay?
> > > I've used one of my already pre-processed subjects, created a 'vol'
> > with the
> > > same number of vertices of one of the hemisferes, and saved the output
> > of
> > > save_mgh file as an rh.thickness.mgh file. However, when I don't seem
> > to be
> > > able to open it in tksurfer...
> > >
> > > Thanks
> > >
> > > Ines
> > >
> > > 2007/8/6, Bruce Fischl < [EMAIL PROTECTED]>:
> > >>
> > >> they don't matter. For a real volume they would be the ras2vox
> > transform
> > >> and the mr parameters. For your purposes they can be any 4x4 and you
> > can
> > >> leave the mri_parms one out (it will detect and fill it in with
> > zeros).
> > >>
> > >> Send these to the list so other people can answer!
> > >>
> > >> Bruce
> > >> On
> > >> Mon, 6 Aug 2007, Inês Souta wrote:
> > >>
> > >>> Hi again,
> > >>>
> > >>> I still don't understand what all the input objects of save_mgh
> > function
> > >>> stand for. 'vol' is the vector of thickness values, 'fname' is the
> > name
> > >> of
> > >>> the output .mgh file, and what about M and ? What do they
> > >> stand
> > >>> for and how should I define them?
> > >>>
> > >>> Thanks again,
> > >>>
> > >>> Ines
> > >>>
> > >>> 2007/7/31, Bruce Fischl <[EMAIL PROTECTED]>:
> > 
> >  Hi Ines,
> > 
> >  use the mgh route. The "vol" is actually the vector of thickness
> > >> values,
> >  with one entry per vertex on the surface.
> > 
> >  cheers,
> >  Bruce
> >  On Tue, 31 Jul 2007, Inês Souta
> >  wrote:
> > 
> > > Hello,
> > >
> > > I was trying to create a cortical thickness map with constant
> > values
> > >> in
> > > order to perform some tests. At first I tried using mris_convert.
> > I
> > > converted my thickness file into ascii:
> > >
> > > mris_convert lh.thickness lh.thickess.asc
> > >
> > > then I used matlab to change all the values in the last column to
> > 2,0.
> > >> I
> > > tried to convert the ascii file back to surface file by doing
> > >
> > > mris_convert lh.thicness.asc lh.thickness
> > >
> > > but I got "Segmentation fault". How can I convert may ascii file
> > into
> > >> a
> > > surface one?
> > >
> > > After failling to acomplish this, I tried using the save_mgh and
> >  load_mgh
> > > file, nevertheless I'm not quite sure which input values are
> > needed.
> >  Where
> > > can I put my fake values in the save_mgh file? And what is "vol"?
> > >
> > > Thanks,
> > >
> > > Ines
> > >
> > >>>
> > >
>
>
>
___
Freesurfer mailing list
Freesurfer@nmr.mgh.harvard.edu
https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer

Re: [Freesurfer] Fake thickness map 2

2007-08-06 Thread Inês Souta
And after creating the mgh volume, how can I read it? Is it possible to open
it in tksurfer? Should I load it as an overlay?
I've used one of my already pre-processed subjects, created a 'vol' with the
same number of vertices of one of the hemisferes, and saved the output of
save_mgh file as an rh.thickness.mgh file. However, when I don't seem to be
able to open it in tksurfer...

Thanks

Ines

2007/8/6, Bruce Fischl <[EMAIL PROTECTED]>:
>
> they don't matter. For a real volume they would be the ras2vox transform
> and the mr parameters. For your purposes they can be any 4x4 and you can
> leave the mri_parms one out (it will detect and fill it in with zeros).
>
> Send these to the list so other people can answer!
>
> Bruce
> On
> Mon, 6 Aug 2007, Inês Souta wrote:
>
> > Hi again,
> >
> > I still don't understand what all the input objects of save_mgh function
> > stand for. 'vol' is the vector of thickness values, 'fname' is the name
> of
> > the output .mgh file, and what about M and ? What do they
> stand
> > for and how should I define them?
> >
> > Thanks again,
> >
> > Ines
> >
> > 2007/7/31, Bruce Fischl <[EMAIL PROTECTED]>:
> >>
> >> Hi Ines,
> >>
> >> use the mgh route. The "vol" is actually the vector of thickness
> values,
> >> with one entry per vertex on the surface.
> >>
> >> cheers,
> >> Bruce
> >> On Tue, 31 Jul 2007, Inês Souta
> >> wrote:
> >>
> >>> Hello,
> >>>
> >>> I was trying to create a cortical thickness map with constant values
> in
> >>> order to perform some tests. At first I tried using mris_convert. I
> >>> converted my thickness file into ascii:
> >>>
> >>> mris_convert lh.thickness lh.thickess.asc
> >>>
> >>> then I used matlab to change all the values in the last column to 2,0.
> I
> >>> tried to convert the ascii file back to surface file by doing
> >>>
> >>> mris_convert lh.thicness.asc lh.thickness
> >>>
> >>> but I got "Segmentation fault". How can I convert may ascii file into
> a
> >>> surface one?
> >>>
> >>> After failling to acomplish this, I tried using the save_mgh and
> >> load_mgh
> >>> file, nevertheless I'm not quite sure which input values are needed.
> >> Where
> >>> can I put my fake values in the save_mgh file? And what is "vol"?
> >>>
> >>> Thanks,
> >>>
> >>> Ines
> >>>
> >
___
Freesurfer mailing list
Freesurfer@nmr.mgh.harvard.edu
https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer

Re: [Freesurfer] Fake thickness map 2

2007-08-06 Thread Bruce Fischl
yes, you should be able to open it in tksurfer. What version of 
freesurfer do you have?

On Mon, 6 Aug 2007, Inês Souta wrote:


And after creating the mgh volume, how can I read it? Is it possible to open
it in tksurfer? Should I load it as an overlay?
I've used one of my already pre-processed subjects, created a 'vol' with the
same number of vertices of one of the hemisferes, and saved the output of
save_mgh file as an rh.thickness.mgh file. However, when I don't seem to be
able to open it in tksurfer...

Thanks

Ines

2007/8/6, Bruce Fischl <[EMAIL PROTECTED]>:


they don't matter. For a real volume they would be the ras2vox transform
and the mr parameters. For your purposes they can be any 4x4 and you can
leave the mri_parms one out (it will detect and fill it in with zeros).

Send these to the list so other people can answer!

Bruce
On
Mon, 6 Aug 2007, Inês Souta wrote:


Hi again,

I still don't understand what all the input objects of save_mgh function
stand for. 'vol' is the vector of thickness values, 'fname' is the name

of

the output .mgh file, and what about M and ? What do they

stand

for and how should I define them?

Thanks again,

Ines

2007/7/31, Bruce Fischl <[EMAIL PROTECTED]>:


Hi Ines,

use the mgh route. The "vol" is actually the vector of thickness

values,

with one entry per vertex on the surface.

cheers,
Bruce
On Tue, 31 Jul 2007, Inês Souta
wrote:


Hello,

I was trying to create a cortical thickness map with constant values

in

order to perform some tests. At first I tried using mris_convert. I
converted my thickness file into ascii:

mris_convert lh.thickness lh.thickess.asc

then I used matlab to change all the values in the last column to 2,0.

I

tried to convert the ascii file back to surface file by doing

mris_convert lh.thicness.asc lh.thickness

but I got "Segmentation fault". How can I convert may ascii file into

a

surface one?

After failling to acomplish this, I tried using the save_mgh and

load_mgh

file, nevertheless I'm not quite sure which input values are needed.

Where

can I put my fake values in the save_mgh file? And what is "vol"?

Thanks,

Ines



___
Freesurfer mailing list
Freesurfer@nmr.mgh.harvard.edu
https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer

Re: [Freesurfer] How to do a surface-based group analysis with single-subject SPM5 contrast images

2007-08-06 Thread Doug Greve
That looks exactly right. Are you experiencing a problem? I assume that 
uS1_v01_06LOC_x001.img is a mean functional image of some sort (and not 
a contrast image). Did you check the registration?


doug

Johannes Klackl wrote:


Dear all,

I'm currently trying to conduct a surface-based group analysis using 
1st level contrast images from SPM5.
I wanted to ask if somebody could tell me if what I'm doing is OK. I 
would relly appreciate your help!


I will now describe the procedure:

The functionals of each subject were realigned and unwarped using SPM5 
preprocessing commands and then entered into an SPM 1st level designs.
Additionally, of course, cortical models from all subjects were 
created using Freesurfer.


I used spmregister to realign the contrast images (con_0001.img) with 
the Freesurfer anatomicals.
spmregister --s v01 --mov ../f2/auS1_v01_06LOC_x001.img --reg 
v01_register.dat


Then I used mri_vol2surf (with the resulting registration file from 
spmregister) to bring the volumetric information from the contrast 
images onto each subject's individual surface.
mri_vol2surf --hemi lh --src con_0001.img --src_type analyze --srcreg 
v01_register.dat --out lh.con_0001.mgh --out_type mgh --projfrac 0.5 
--noreshape


After that, I used mri_surf2surf with ico as the target to bring all 
the functionals from the individual subject surface onto a common 
surface.
mri_surf2surf --hemi lh --srcsubject v01 --srcsurfval lh.con_0001.mgh 
--src_type mgh --trgsubject ico --trgicoorder 7 --trgsurfval 
lh.ico.con_0001.mgh --trg_type mgh


I was thinking about smoothing here, but I didn't do this yet.

Then, I used mri_concat to concatenate the functionals from all 
subjects into one file and mri_glmfit to do the statistics.


Thanks,

Johannes Klackl, University of Salzburg
___
Freesurfer mailing list
Freesurfer@nmr.mgh.harvard.edu
https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer




--
Douglas N. Greve, Ph.D.
MGH-NMR Center
[EMAIL PROTECTED]
Phone Number: 617-724-2358 
Fax: 617-726-7422


In order to help us help you, please follow the steps in:
surfer.nmr.mgh.harvard.edu/fswiki/BugReporting


___
Freesurfer mailing list
Freesurfer@nmr.mgh.harvard.edu
https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer


Re: [Freesurfer] optseq2 question about fixing null intervals

2007-08-06 Thread Doug Greve


what do you mean by "allowable null intervals"? If you mean that there 
needs to be a fixed amount of time between two stimuli, then you cannot 
use optseq for this. However, you can tell optseq that there are only 
two stimuli, which each stimulus is a pair. Then you'd have to go back 
and recode the timing file.


doug



Denis R wrote:


Hello,

We're designing a task that has two pairs of events (4 events total). 
Each pair of events has a pair of allowable null intervals associated 
with it. Is it possible to use optseq2 to create such a paradigm? Is 
it appropriate to use optseq2 in this case at all given these 
constraints?




___
Freesurfer mailing list
Freesurfer@nmr.mgh.harvard.edu
https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer



--
Douglas N. Greve, Ph.D.
MGH-NMR Center
[EMAIL PROTECTED]
Phone Number: 617-724-2358 
Fax: 617-726-7422


In order to help us help you, please follow the steps in:
surfer.nmr.mgh.harvard.edu/fswiki/BugReporting


___
Freesurfer mailing list
Freesurfer@nmr.mgh.harvard.edu
https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer

Re: [Freesurfer] questions

2007-08-06 Thread Doug Greve


For #2, you can use mri_volsynth with --pdf checker to create a checker 
board pattern. You'll need a new version, which I have put here:


ftp://surfer.nmr.mgh.harvard.edu/transfer/outgoing/flat/greve/mri_volsynth

doug


Bruce Fischl wrote:

hmmm, not sure i understand. You can visualize any scalar field in the 
volume, so yes. As for the volume, I guess you could create a 
functional volume with the value in each voxel the slice # that it is 
in and just paint that onto the surface as normal.


cheers,
Bruce


On Fri, 3 Aug 2007, Paul Aparicio wrote:


Hello All,


I have 2 questions I hope someone can help me with.  1) Is it 
possible to plot on a surface a set of arbitrary points in the 
coordinate frame of an anatomical volume?  In other words if I have a 
set of x,y, and z points (not voxel indices: e.g. sub-voxel 
coordinates) can I plot them onto a surface? 2) Is there a way to 
view the slice lines in a volume and a surface?  Can I see the 
anatomical location of the boundaries of the functional slices on the 
surface representation?


Thanks

paul
___
Freesurfer mailing list
Freesurfer@nmr.mgh.harvard.edu
https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer




___
Freesurfer mailing list
Freesurfer@nmr.mgh.harvard.edu
https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer




--
Douglas N. Greve, Ph.D.
MGH-NMR Center
[EMAIL PROTECTED]
Phone Number: 617-724-2358 
Fax: 617-726-7422


In order to help us help you, please follow the steps in:
surfer.nmr.mgh.harvard.edu/fswiki/BugReporting


___
Freesurfer mailing list
Freesurfer@nmr.mgh.harvard.edu
https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer


Re: [Freesurfer] possible bug found in mri_surfcluster

2007-08-06 Thread Doug Greve
Where is this data located? I cant find it under 
/space/ventzl/1/users/SUBJECTS_DIR



Robert Levy wrote:


Hello,

I found a possible bug in mri_surfcluster, but it may be linked 
specifically to average7 and no other srcsubjects (but I haven't 
checked that).  When I run mri_surfcluster with cluster simulation 
data, the size is reported at a seemingly arbitrary number, based on 
the fact that the #of vertices does not change.  And the mgh output in 
null (a blank surface is displayed). 

However when I run mri_surfcluster with cluster simulation data,  but  
specify srcsubject, it works correctly.  This seems like a bug, 
because I think when surfcluster runs csd files it reads the 
srcsubject from the csd files, and no other srcsubject should be 
compatible with the data, therefore it must know the srcsubject in 
both usages, but in one case it works and the other it doesn't. 


Thanks,
Rob


surflcluster with threshold of 1.3
# ClusterNo  Max   VtxMax   Size(mm^2)  TalX   TalY   TalZNVtxs
   1   -2.217  145639413.11 -6.4   38.45.0   418
   2   -1.836   29811141.08-11.9   53.2   19.0   178
   3   -1.532   16680111.99-11.1   19.1   22.6   178
   4   -1.400   83520 17.42 -3.2   11.4   17.523

surfcluster using csd files
# ClusterNo  Max   VtxMax   Size(mm^2)  TalX   TalY   TalZCWP
CWPLowCWPHi   NVtxs
   1   -2.217  145639275.40 -6.4   38.45.0  0.01190  
0.01050  0.01330   418
   2   -1.836   29811 94.06-11.9   53.2   19.0  0.50970  
0.50330  0.51610   178
   3   -1.532   16680 74.66-11.1   19.1   22.6  0.64580  
0.63970  0.65190   178
   4   -1.400   83520 11.61 -3.2   11.4   17.5  0.97920  
0.97740  0.9810023



Apparently this is the right way to do it:
with csd files, but specify srcsubj average7:

# ClusterNo  Max   VtxMax   Size(mm^2)  TalX   TalY   TalZCWP
CWPLowCWPHi   NVtxs
   1   -2.217  145639413.11 -6.4   38.45.0  0.00010  
0.0  0.00020   418
   2   -1.836   29811141.08-11.9   53.2   19.0  0.27000  
0.26430  0.27570   178
   3   -1.532   16680111.99-11.1   19.1   22.6  0.40500  
0.39870  0.41130   178
   4   -1.400   83520 17.42 -3.2   11.4   17.5  0.96970  
0.96750  0.9719023


--- freesurfer-Linux-centos4_x86_64-dev-20070802 
Setting up environment for FreeSurfer/FS-FAST (and FSL)
FREESURFER_HOME   /usr/local/freesurfer/dev
FSFAST_HOME   /usr/local/freesurfer/dev/fsfast
FSF_OUTPUT_FORMAT nii
SUBJECTS_DIR  /space/ventzl/1/users/SUBJECTS_DIR
MNI_DIR   /usr/local/freesurfer/dev/mni
FSL_DIR   /usr/local/freesurfer/dev/fsl

# PWD:  /space/ventzl/1/users/SUBJECTS_DIR

mri_surfcluster \
--src lh.ASvfix_8s_BwGroupError5.fsfast.glmdir/BwGroupError5/sig.mgh \
--csd ASevASc_8s_p05_z_lhMCSIM01-BwGroupError5.csd \
--csd ASevASc_8s_p05_z_lhMCSIM02-BwGroupError5.csd \
--csd ASevASc_8s_p05_z_lhMCSIM03-BwGroupError5.csd \
--csd ASevASc_8s_p05_z_lhMCSIM04-BwGroupError5.csd \
--csd ASevASc_8s_p05_z_lhMCSIM05-BwGroupError5.csd \
--csd ASevASc_8s_p05_z_lhMCSIM06-BwGroupError5.csd \
--csd ASevASc_8s_p05_z_lhMCSIM07-BwGroupError5.csd \
--csd ASevASc_8s_p05_z_lhMCSIM08-BwGroupError5.csd \
--csd ASevASc_8s_p05_z_lhMCSIM09-BwGroupError5.csd \
--csd ASevASc_8s_p05_z_lhMCSIM10-BwGroupError5.csd \
--sum 
lh.ASvfix_8s_BwGroupError5.fsfast.glmdir/BwGroupError5/sig.clusterMCSIM_lh.sum 
\
--ocp 
lh.ASvfix_8s_BwGroupError5.fsfast.glmdir/BwGroupError5/sig.clusterMCSIM_lh.mgh 
\


# overlay shows no clusters (tried at 0.5 as well in case it is a 
binary map)
tksurfer average7 lh inflated -overlay 
lh.ASvfix_8s_BwGroupError5.fsfast.glmdir/BwGroupError5/sig.clusterMCSIM_lh.mgh


mri_surfcluster \
--really-use-average7 \
--src lh.ASvfix_8s_BwGroupError5.fsfast.glmdir/BwGroupError5/sig.mgh \
--csd ASevASc_8s_p05_z_lhMCSIM01-BwGroupError5.csd \
--csd ASevASc_8s_p05_z_lhMCSIM02-BwGroupError5.csd \
--csd ASevASc_8s_p05_z_lhMCSIM03-BwGroupError5.csd \
--csd ASevASc_8s_p05_z_lhMCSIM04-BwGroupError5.csd \
--csd ASevASc_8s_p05_z_lhMCSIM05-BwGroupError5.csd \
--csd ASevASc_8s_p05_z_lhMCSIM06-BwGroupError5.csd \
--csd ASevASc_8s_p05_z_lhMCSIM07-BwGroupError5.csd \
--csd ASevASc_8s_p05_z_lhMCSIM08-BwGroupError5.csd \
--csd ASevASc_8s_p05_z_lhMCSIM09-BwGroupError5.csd \
--csd ASevASc_8s_p05_z_lhMCSIM10-BwGroupError5.csd \
--sum 
lh.ASvfix_8s_BwGroupError5.fsfast.glmdir/BwGroupError5/sig.clusterMCSIM_lh.sum 
\
--ocp 
lh.ASvfix_8s_BwGroupError5.fsfast.glmdir/BwGroupError5/sig.clusterMCSIM_lh.mgh 
\

--srcsubj average7

# overlay does show the clusters
tksurfer average7 lh inflated -overlay 
lh.ASvfix_8s_BwGroupError5.fsfast.glmdir/BwGroupError5/sig.clusterMCSIM_lh.mgh


--

Robert P. Levy, B.A.
Research Assistant, Manoach Lab
Massachusetts General Hospital
Charlestown Navy Yard
149 13th St., Room 2611
Charlestown, MA 02129
email: [EMAIL PROTECTED]
phone: 617-726-1908
fax: 617-726-4078

http://nmr.mgh.harvard.ed

[Freesurfer] Banks of the superior temporal sulcus question

2007-08-06 Thread Morgan Hough
I was just wondering why the banks of the superior temporal sulcus was 
made into its own label on the Desikan template. Does this correspond 
(roughly) to a particular functional area. I didn't see anything in the 
paper about why it was chosen. Is it just easy to define given 
sulcal/gyral boundaries? Thanks in advance.


Cheers,

-Morgan
___
Freesurfer mailing list
Freesurfer@nmr.mgh.harvard.edu
https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer


Re: [Freesurfer] Banks of the superior temporal sulcus question

2007-08-06 Thread Bruce Fischl
I think it's where Ron Killiany previously found large effects in 
prodromal AD.


On Mon, 6 Aug 2007, Morgan Hough wrote:

I was just wondering why the banks of the superior temporal sulcus was made 
into its own label on the Desikan template. Does this correspond (roughly) to 
a particular functional area. I didn't see anything in the paper about why it 
was chosen. Is it just easy to define given sulcal/gyral boundaries? Thanks 
in advance.


Cheers,

-Morgan
___
Freesurfer mailing list
Freesurfer@nmr.mgh.harvard.edu
https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer




___
Freesurfer mailing list
Freesurfer@nmr.mgh.harvard.edu
https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer


Re: [Freesurfer] Banks of the superior temporal sulcus question

2007-08-06 Thread Morgan Hough

Hi Bruce,

I see. It just seems to be surrounded by very large parcellation units 
that haven't been divided. Does it have good landmarks to place it along 
the sulcus?


Cheers,

-Morgan

Bruce Fischl wrote:
I think it's where Ron Killiany previously found large effects in 
prodromal AD.


On Mon, 6 Aug 2007, Morgan Hough wrote:

I was just wondering why the banks of the superior temporal sulcus 
was made into its own label on the Desikan template. Does this 
correspond (roughly) to a particular functional area. I didn't see 
anything in the paper about why it was chosen. Is it just easy to 
define given sulcal/gyral boundaries? Thanks in advance.


Cheers,

-Morgan
___
Freesurfer mailing list
Freesurfer@nmr.mgh.harvard.edu
https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer




___
Freesurfer mailing list
Freesurfer@nmr.mgh.harvard.edu
https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer


Re: [Freesurfer] Banks of the superior temporal sulcus question

2007-08-06 Thread Bruce Fischl

Hi Morgan,

I'm not really sure. You could check with Ron and Rahul (ccd).

cheers,
Bruce
On Mon, 6 Aug 
2007, Morgan Hough wrote:



Hi Bruce,

I see. It just seems to be surrounded by very large parcellation units that 
haven't been divided. Does it have good landmarks to place it along the 
sulcus?


Cheers,

-Morgan

Bruce Fischl wrote:
I think it's where Ron Killiany previously found large effects in prodromal 
AD.


On Mon, 6 Aug 2007, Morgan Hough wrote:

I was just wondering why the banks of the superior temporal sulcus was 
made into its own label on the Desikan template. Does this correspond 
(roughly) to a particular functional area. I didn't see anything in the 
paper about why it was chosen. Is it just easy to define given 
sulcal/gyral boundaries? Thanks in advance.


Cheers,

-Morgan
___
Freesurfer mailing list
Freesurfer@nmr.mgh.harvard.edu
https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer








___
Freesurfer mailing list
Freesurfer@nmr.mgh.harvard.edu
https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer


[Freesurfer] mris_anatomical_stats question

2007-08-06 Thread Morgan Hough
I am reporting thickness values from subjects native space constrained 
by mask labels using mris_anatomical_stats. I am currently getting the 
average of the unsmoothed estimates. I would like to be able to get 
smoothed values but the usage reports the smoothing option as -nsmooth 
<#> where the # is the # of smoothing iterations. If we are using 10mm 
FWHM smoothing for the group results can I estimate the number of 
smoothing iterations to obtain something equivalent? Would it be easier 
to use mri_surf2surf and create a  smoothed thickness file for each 
subject? Thanks again.


Cheers,

-Morgan
___
Freesurfer mailing list
Freesurfer@nmr.mgh.harvard.edu
https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer


[Freesurfer] editing labels

2007-08-06 Thread Sarah J Greene
Hi all,
  I just wanted to ask about editing labels in tksurfer.  I have
converted my annot files to labels using mri_annotation2label, and now
have all of the labels in the label directory.  However when I edit
the labels, and then save them using save selected label (highlighting
the label and choosing the correct label in the correct hemisphere),
nothing is saved when I pull up the image again.  I have tried to
convert the labels to annot files using mris_label2annot, but this
isn't working either.  I'm not sure what I am doing incorrectly.
Thanks,
sarah
___
Freesurfer mailing list
Freesurfer@nmr.mgh.harvard.edu
https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer


Re: [Freesurfer] editing labels

2007-08-06 Thread Kevin Teich
Hi Sarah,

When you save the label, is there any output on the shell?

Can you load and edit a single label, save it out, exit tksurfer,
restart it, and load the label in again, with the changes you made
present in the label?



On Mon, 2007-08-06 at 14:09 -0400, Sarah J Greene wrote:
> Hi all,
>   I just wanted to ask about editing labels in tksurfer.  I have
> converted my annot files to labels using mri_annotation2label, and now
> have all of the labels in the label directory.  However when I edit
> the labels, and then save them using save selected label (highlighting
> the label and choosing the correct label in the correct hemisphere),
> nothing is saved when I pull up the image again.  I have tried to
> convert the labels to annot files using mris_label2annot, but this
> isn't working either.  I'm not sure what I am doing incorrectly.
> Thanks,
> sarah


-- 
Kevin Teich
___
Freesurfer mailing list
Freesurfer@nmr.mgh.harvard.edu
https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer


Re: [Freesurfer] Thickness standard deviation map

2007-08-06 Thread Doug Greve
The easiest thing is #2. It will produce rstd.mgh as an output, and 
that's what you want.


Burmicz, Ryszarda wrote:


Hi Freesurfers,

I have tried looking this up on the list, but I don't understand 
exactly what to do in this scenario: I have a group of subjects (all 
controls).  I want to create a surface map of the standard deviations 
of the thickness values (i.e. to show how the thicknesses vary across 
subjects).  I am not sure what is the easiest way to go about this - 
can I:


1) simply convert the ?h.thickness files to .asc and calculate the 
deviations at each vertex (I am under the impression that this is not 
correct as they need to be registered to a common surface first), then 
convert back to .thickness and load into tksurfer


2) use mris_preproc then mris_glmfit?  I don't necessarily want to 
calculate any stats, but I wouldn't mind a pretty picture :)


If I have 4 groups of 10 subjects each, so if I churn these all 
through glm_fit will I be able to access separate thickness files for 
each group? And how can I map standard deviations for each group in 
this scenario?


Many thanks,

Rysia




___
Freesurfer mailing list
Freesurfer@nmr.mgh.harvard.edu
https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer



--
Douglas N. Greve, Ph.D.
MGH-NMR Center
[EMAIL PROTECTED]
Phone Number: 617-724-2358 
Fax: 617-726-7422


In order to help us help you, please follow the steps in:
surfer.nmr.mgh.harvard.edu/fswiki/BugReporting


___
Freesurfer mailing list
Freesurfer@nmr.mgh.harvard.edu
https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer

Re: [Freesurfer] determining sizes of clusters at given CWP cutoff

2007-08-06 Thread Doug Greve

I don't have an easy way to do this, sorry.

Robert Levy wrote:


Hi,

I was trying to figure out if there is any way to determine the exact 
size at which clusters becomes significant given a specified 
significance, and it appeared from the mc-corrected overlay that in 
this overlay the clusters do not continously change in size but appear 
discretely at a given threshold and are the same size at higher 
thresholds in an all-or-nother sort of way.  I take this to mean that 
because the monte carlo simulation was performed at a given threshold, 
the clusters at that threshold size are corrected for mutiple 
comparisons based on the probability distribution of the chance data 
(though I don't know the details).  Because it is done in this way, it 
seems I would either have to run several monte carlo simulations, at 
successively higher minimum thresholds, in order to arrive at the 
critical cutoff sizes for significance.  This is impractical, so I 
guess my question is, does the way in which the monte carlo 
simulations are done preclude the possibility of getting the cluster 
sizes at which they will have a certain CWP value, or is there some 
way of doing this?


Thanks,
Rob

--

Robert P. Levy, B.A.
Research Assistant, Manoach Lab
Massachusetts General Hospital
Charlestown Navy Yard
149 13th St., Room 2611
Charlestown, MA 02129
email: [EMAIL PROTECTED]
phone: 617-726-1908
fax: 617-726-4078

http://nmr.mgh.harvard.edu/manoachlab
 




___
Freesurfer mailing list
Freesurfer@nmr.mgh.harvard.edu
https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer



--
Douglas N. Greve, Ph.D.
MGH-NMR Center
[EMAIL PROTECTED]
Phone Number: 617-724-2358 
Fax: 617-726-7422


In order to help us help you, please follow the steps in:
surfer.nmr.mgh.harvard.edu/fswiki/BugReporting


___
Freesurfer mailing list
Freesurfer@nmr.mgh.harvard.edu
https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer

[Freesurfer] dev version of mri_glmfit for intel macs

2007-08-06 Thread Catherine Hartley

Hi,

I'm getting a bus error when I run mri_glmfit on my Intel Mac.  When  
I had this problem before, I got a development version of mri_glmfit  
which solved the problem.  Unfortunately I've reinstalled a newer  
version of Freesurfer since then and overwrote that binary file and  
it's no longer on  the public distribution server.


If someone could repost the latest version of mri_glmfit for Intel  
Macs (10.4.9), that would be great.


The same problem is described here:
http://www.mail-archive.com/freesurfer@nmr.mgh.harvard.edu/msg04632.html

Thanks,
Cate
___
Freesurfer mailing list
Freesurfer@nmr.mgh.harvard.edu
https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer


Re: [Freesurfer] dev version of mri_glmfit for intel macs

2007-08-06 Thread Nick Schmansky
Cate,

There is a current dev distribution of freesurfer for the Intel Mac
here:

ftp://surfer.nmr.mgh.harvard.edu/transfer/outgoing/fsdev/nicks/dev/macintel

Its the .dmg file, which you can drag and drop to your desktop, then
double-click to start the install GUI.

A new mri_glmfit is also in that dir if you just want that.

Nick


On Mon, 2007-08-06 at 17:56 -0400, Catherine Hartley wrote:
> Hi,
> 
> I'm getting a bus error when I run mri_glmfit on my Intel Mac.  When  
> I had this problem before, I got a development version of mri_glmfit  
> which solved the problem.  Unfortunately I've reinstalled a newer  
> version of Freesurfer since then and overwrote that binary file and  
> it's no longer on  the public distribution server.
> 
> If someone could repost the latest version of mri_glmfit for Intel  
> Macs (10.4.9), that would be great.
> 
> The same problem is described here:
> http://www.mail-archive.com/freesurfer@nmr.mgh.harvard.edu/msg04632.html
> 
> Thanks,
> Cate
> ___
> Freesurfer mailing list
> Freesurfer@nmr.mgh.harvard.edu
> https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer
> 
> 

___
Freesurfer mailing list
Freesurfer@nmr.mgh.harvard.edu
https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer