Hello Doug,
Thanks for your help on this. Ive run it again and there doesnt seem to be a problem. Its just that when I smooth it by 15 now the residual is about 23. When I ran it the other time, the residual was 47 for some unknown reason.
Regards,
Joshua
I think when you smooth it that much strange things happen. Try just one smoothing of 15mm.
On 06/06/2013 03:36 AM, Joshua deSouza wrote:
Hi Doug,
I did not specify anything in the mris_preproc, so it would be using the unsmoothed data. I only specify smoothing during the mri_surf2surf, which is 15 and then run glmfit. So what could be the reasons why i'm getting a residual fwhm of 47?
Regards,
Joshua
On Thu, Jun 6, 2013 at 12:16 AM, Douglas N Greve <greve@nmr.mgh.harvard.edu <mailto:greve@nmr.mgh.harvard.edu>> wrote:<mailto:desouzajoshua00@gmail.com
Unless you specify otherwise, mris_preproc will use unsmoothed
data. If you then smooth it with mri_surf2surf and run run
mri_glmfit with --fwhm then you will be smoothing it twice.
doug
On 06/04/2013 09:49 PM, Joshua deSouza wrote:
Hello Doug,
Each individuals subject data has already been smoothed by 5
10 15 20 and 25 through the qache function during the initial
processing.
I'm a little confused now. I thought that the mri_preproc
would use the non-smoothed values of the subjects and create
an mgh file. After which I smooth the mgh file and then use
mri_glmfit to analyse it according to the contrasts.
If the above is not the case and the subject data has already
been smoothed, then what should I do if I wish my final mgh
file to have a smoothing value of 15 before doing a glmfit on it?
Thank you,
Joshua
On Wed, Jun 5, 2013 at 12:40 AM, Joshua deSouza
<desouzajoshua00@gmail.com <mailto:desouzajoshua00@gmail.com><mailto:greve@nmr.mgh.harvard.edu
<mailto:desouzajoshua00@gmail.com>>> wrote:
Hello Doug,
Each individuals subject data has already been smoothed by
5 10 15
20 and 25 through the qache function during the initial
processing.
I'm a little confused now. I thought that the mri_preproc
would
use the non-smoothed values of the subjects and create an mgh
file. After which I smooth the mgh file and then use
mri_glmfit to
analyse it according to the contrasts.
If the above is not the case and the subject data has
already been
smoothed, then what should I do if I wish my final mgh file to
have a smoothing value of 15 before doing a glmfit on it?
Thank you,
Joshua
On Mon, Jun 3, 2013 at 11:13 PM, Douglas Greve
<greve@nmr.mgh.harvard.edu
<mailto:greve@nmr.mgh.harvard.edu><mailto:greve@nmr.mgh.harvard.edu
<mailto:greve@nmr.mgh.harvard.edu>>> wrote:
Hi Joshua, when you add --fwhm 15 to the command line
you are
instructing mri_glmfit to smooth by 15mm. Is this what you
want? Or is the data already smoothed?
doug
On 6/3/13 12:26 AM, Joshua deSouza wrote:
Hi Doug,
The command line is mri_glmfit --y lh.t.15.mgh
--fsgd t.txt
dods --C t1.txt --surface fsaverage lh --cortex
--fwhm 15
--glmdir lh.t.15.glmdir
I am doing a whole hemisphere analysis and am
measuring
cortical thickness. I do not have my own mask but
I have
doubts about the mask that has been produced
because from
previous studies, there has been significant
thinning in the
area around the hippocampus and parts of it are
being left
out in the mask, such as the para-hippocampal region.
Regards
Joshua
On Mon, Jun 3, 2013 at 10:54 AM, Douglas Greve
<greve@nmr.mgh.harvard.edu
<mailto:greve@nmr.mgh.harvard.edu><mailto:greve@nmr.mgh.harvard.edu
<mailto:greve@nmr.mgh.harvard.edu>>> wrote:
The residual FWHM is measured from the data
itself. You
apply a certain amount of smoothing, but there
is already
some smoothness to the data, so the final fwhm
is a
combination of what is already there and what
you apply.
A fwhm of 47 is huge and may indicate that
something is
wrong. What is the input? cortical thickness?
Are you
using your own mask or is this a whole
hemisphere analysis?
doug
On 6/2/13 10:29 PM, Joshua deSouza wrote:
Hello Doug,
Thank you for your quick replies to my
queries. Really
appreciate it.
I have found the files containing the
cortical thickness
values. Thank you. I assume that the order
in which they
follow are based on the fsgd file and the
cluster
annotation numbers.
I was just slightly confused about the
residual fwhm. I
have already done a smoothing of 15 on the
file and have
recieved a residual value of 47, which I
assume is
extremely large. Since I'm finding out
values for the
entire cortex I cannot do a glmfit-sim as
there are no
fwhm cache values of 47. Could you please
tell me what
the residual values mean? What do you
recommend I do?
Is there anything I can read up on to find
out more
about this? (Please pardon my questions if
they seem a
little trivial)
Thank you
Joshua
On Sun, Jun 2, 2013 at 11:33 PM, Douglas Greve
<greve@nmr.mgh.harvard.edu
<mailto:greve@nmr.mgh.harvard.edu><mailto:greve@nmr.mgh.harvard.edu
<mailto:greve@nmr.mgh.harvard.edu>>> wrote:
This file is always created by
mri_glmfit-sim. It
might be named something different
depending on your
parameters (should always end in
y.ocn.dat)
doug
On 6/2/13 3:57 AM, Joshua deSouza wrote:
Hi Doug,
This file that your saying, will
be created after
glmfit-sim is done using the cache
function and a
threshold of 1.3 and absolute
,with y being the
fsgd file? Am I understanding this
correctly?
Thank you
Joshua
On Sat, Jun 1, 2013 at 12:14 AM,
Douglas N Greve
<greve@nmr.mgh.harvard.edu
<mailto:greve@nmr.mgh.harvard.edu><mailto:Freesurfer@nmr.mgh.harvard.edu
<mailto:greve@nmr.mgh.harvard.edu>>> wrote:
There should be a file called
something like
cache.th13.abs.y.ocn.dat.
This will have a column for
each cluster and a
row for each subject. The
value in the row will be the
input for each
subject.
doug
On 05/30/2013 11:21 PM, Joshua
deSouza wrote:
> Hello everyone,
>
> Please pardon my rather
simple question
>
> I have done a GLM fit
anlysis and have
obtained a file sig.mgh alonf
> with a few others after the
analysis.As I
know, I do not have to
> create ROIs in order to do
this. I also
understand that doing a
> multiple comparison analysis
just compares
the regions of importance
> and evaluates them giving
the most resilient
cluster depending on the
> threshold. So, this does not
give a thickness
estimate. I would like
> to know how do I get the
cortical thickness
values from the results?
> If so, can I do it after glmfit?
>
> And thanks Doug for your
help in the other
questions that I had prior
> to this!
>
> Thank you
>
>
>
_______________________________________________
> Freesurfer mailing list
>
Freesurfer@nmr.mgh.harvard.edu
<mailto:Freesurfer@nmr.mgh.harvard.edu><mailto:greve@nmr.mgh.harvard.edu <tel:617-724-2358> <tel:617-724-2358
<mailto:Freesurfer@nmr.mgh.harvard.edu>>
>
https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer
--
Douglas N. Greve, Ph.D.
MGH-NMR Center
greve@nmr.mgh.harvard.edu
<mailto:greve@nmr.mgh.harvard.edu>
<tel:617-724-2358>>
Fax: 617-726-7422
<tel:617-726-7422> <tel:617-726-7422<mailto:Freesurfer@nmr.mgh.harvard.edu
<tel:617-726-7422>>
Bugs:
surfer.nmr.mgh.harvard.edu/fswiki/BugReporting
<http://surfer.nmr.mgh.harvard.edu/fswiki/BugReporting>
<http://surfer.nmr.mgh.harvard.edu/fswiki/BugReporting>
FileDrop:
https://gate.nmr.mgh.harvard.edu/filedrop2
www.nmr.mgh.harvard.edu/facility/filedrop/index.html
<http://www.nmr.mgh.harvard.edu/facility/filedrop/index.html>
<http://www.nmr.mgh.harvard.edu/facility/filedrop/index.html>
Outgoing:
ftp://surfer.nmr.mgh.harvard.edu/transfer/outgoing/flat/greve/
_______________________________________________
Freesurfer mailing list
Freesurfer@nmr.mgh.harvard.edu
<mailto:Freesurfer@nmr.mgh.harvard.edu>greve@nmr.mgh.harvard.edu <mailto:greve@nmr.mgh.harvard.edu>
<mailto:Freesurfer@nmr.mgh.harvard.edu>>
https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer
The information in this e-mail
is intended only
for the person to whom it is
addressed. If you believe this
e-mail was sent
to you in error and the e-mail
contains patient information,
please contact
the Partners Compliance
HelpLine at
http://www.partners.org/complianceline . If the
e-mail was sent to you in error
but does not contain patient
information,
please contact the sender and
properly
dispose of the e-mail.
-- Douglas N. Greve, Ph.D.
MGH-NMR Center
Phone Number: 617-724-2358 <tel:617-724-2358>
Fax: 617-726-7422 <tel:617-726-7422>
Bugs: surfer.nmr.mgh.harvard.edu/fswiki/BugReporting
<http://surfer.nmr.mgh.harvard.edu/fswiki/BugReporting>
FileDrop: https://gate.nmr.mgh.harvard.edu/filedrop2
www.nmr.mgh.harvard.edu/facility/filedrop/index.html
<http://www.nmr.mgh.harvard.edu/facility/filedrop/index.html>
Outgoing:
ftp://surfer.nmr.mgh.harvard.edu/transfer/outgoing/flat/greve/
--
Douglas N. Greve, Ph.D.
MGH-NMR Center
greve@nmr.mgh.harvard.edu
Phone Number: 617-724-2358
Fax: 617-726-7422
Bugs: surfer.nmr.mgh.harvard.edu/fswiki/BugReporting
FileDrop: https://gate.nmr.mgh.harvard.edu/filedrop2
www.nmr.mgh.harvard.edu/facility/filedrop/index.html
Outgoing: ftp://surfer.nmr.mgh.harvard.edu/transfer/outgoing/flat/greve/