Send Freesurfer mailing list submissions to
 freesurfer@nmr.mgh.harvard.edu

To subscribe or unsubscribe via the World Wide Web, visit
     https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer
or, via email, send a message with subject or body 'help' to
   freesurfer-request@nmr.mgh.harvard.edu

You can reach the person managing the list at
 freesurfer-owner@nmr.mgh.harvard.edu

When replying, please edit your Subject line so it is more specific
than "Re: Contents of Freesurfer digest..."


Today's Topics:

1. Longitudinal LME contrast vectors (K)
2. Re: Longitudinal LME contrast vectors (Martin Reuter)
3. Question about NULL event with Optseq2 (Feng Zhou)
4. Cerebral volume without ventricles (Kaiming Yin)
5. Re: FSFAST error: ERROR: fast_selxavg3() failed\n (Douglas Greve)


----------------------------------------------------------------------

Message: 1
Date: Fri, 8 Jul 2016 15:02:58 -0700
From: K <kam.ntnu@gmail.com>
Subject: [Freesurfer] Longitudinal LME contrast vectors
To: freesurfer@nmr.mgh.harvard.edu
Message-ID:
       <CAH54DzgqKWSCwm0gu_ecbB3jrQmLc=Y6Nun5j4DujKwaxn78UA@mail.gmail.com>
Content-Type: text/plain; charset="utf-8"

Hello! We are trying to run a longitudinal morphometry analysis with a
control group and clinical group scanned twice. We'd like to look at group
x time interactions and would like help figuring out the best way to do so.
As the time variable we use centered age -- the deviation from the mean age
across all subjects and timepoints. To do so, we use the following
lme model:

intercept + centeredage.controls + centeredage.clinical + group(0,1)
+ sex(0,1)

Contrast vector: [0 1 -1 0 0]

where the centeredage column contains only the centered age for that group
and zeros for the opposite group. The thinking behind this is to
investigate which cortical regions show opposite effects of time
between the groups.

However, the FreeSurfer lme software requests a specified "time
column". Obviously, we cannot add a centered age for all subjects, nor does
it seem to make sense to use one of the two group-specific
centeredage columns as the time column.

Our alternative model for investigating interatction is:

intercept + centeredage + group(0,1) + sex(0,1) + groupxcenteredage

where centeredage is the time columns and is a common regressor for both
groups and thus, groupxcenteredage becomes 0 for control subjects.
Contrast vector: [0 0 0 0 1]
This would test whether there is an added effect of centered age in the
clinical group.

Is there a better way to do an group x time interaction analysis here, or
is there a solution for the time column? Thanks!
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mail.nmr.mgh.harvard.edu/pipermail/freesurfer/attachments/20160708/8506ba3b/attachment-0001.html

------------------------------

Message: 2
Date: Sat, 9 Jul 2016 00:39:59 +0200
From: Martin Reuter <mreuter@nmr.mgh.harvard.edu>
Subject: Re: [Freesurfer] Longitudinal LME contrast vectors
To: Freesurfer support list <freesurfer@nmr.mgh.harvard.edu>
Message-ID: <A8A075CB-740E-4032-8A0E-6BC4326636A7@nmr.mgh.harvard.edu>
Content-Type: text/plain; charset="us-ascii"

Hi K,

the models are equivalent (they are linear combinations of each other). The preferred way (at least for the Matlab LME stuff) is the second one.

Note that usually people want to control for age (or test for an age effect). To do so they would have

intercept + time_from_baseline + age_at_baseline + group + sex + groupXtime_from_baseline

Best, Martin

On Jul 9, 2016, at 12:02 AM, K <kam.ntnu@gmail.com> wrote:

Hello! We are trying to run a longitudinal morphometry analysis with a control group and clinical group scanned twice. We'd like to look at group x time interactions and would like help figuring out the best way to do so. As the time variable we use centered age -- the deviation from the mean age across all subjects and timepoints. To do so, we use the following lme model:

intercept + centeredage.controls + centeredage.clinical + group(0,1) + sex(0,1)

Contrast vector: [0 1 -1 0 0]

where the centeredage column contains only the centered age for that group and zeros for the opposite group. The thinking behind this is to investigate which cortical regions show opposite effects of time between the groups.

However, the FreeSurfer lme software requests a specified "time column". Obviously, we cannot add a centered age for all subjects, nor does it seem to make sense to use one of the two group-specific centeredage columns as the time column.

Our alternative model for investigating interatction is:

intercept + centeredage + group(0,1) + sex(0,1) + groupxcenteredage

where centeredage is the time columns and is a common regressor for both groups and thus, groupxcenteredage becomes 0 for control subjects.

Contrast vector: [0 0 0 0 1]

This would test whether there is an added effect of centered age in the clinical group.

Is there a better way to do an group x time interaction analysis here, or is there a solution for the time column? Thanks!

_______________________________________________

Freesurfer mailing list

Freesurfer@nmr.mgh.harvard.edu

https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer

-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mail.nmr.mgh.harvard.edu/pipermail/freesurfer/attachments/20160709/bc5295b7/attachment-0001.html

------------------------------

Message: 3
Date: Sat, 9 Jul 2016 07:39:43 +0800
From: Feng Zhou <fengzhou.nj@gmail.com>
Subject: [Freesurfer] Question about NULL event with Optseq2
To: freesurfer@nmr.mgh.harvard.edu
Message-ID:
       <CAKKeuhMfNDthQ-ZLiORyVRqy7eMi4TvWxDJHScfCV+img_21xA@mail.gmail.com>
Content-Type: text/plain; charset="utf-8"

Dear experts,
I'd like to use optseq2 to optimize the order of events as well as inter
stimulus intervals. My experiment contains 3 runs and each run contains 3
events, say Cp, Cn and Pobe events
Each stimulu will be presented 4, 4 and 1 sec and the numbers of times that
each event type will be presented during the course of the run are 12 12
and 4 respectively, and interstimulus interval ranges from 8?12 sec. So
each run would last around 380 sec. I used the following code:
optseq2 --ntp 380 --tr 1 --psdwin 0 20 1 --ev Cp 4 12 --ev Cn 4 12 --ev
Prob 1 3 --evc 1 -1 0 --tnullmin 8 --tnullmax 12 --nkeep 3 --o IAPS
--nsearch 10000
However, the durations of some NULL events were more than 12 sec and the
last NULL events were more than 30 sec.
Did I use optseq2 in the right way? How could I fix it?
Thanks in advance & Best regards,
Feng Zhou
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mail.nmr.mgh.harvard.edu/pipermail/freesurfer/attachments/20160709/ca4b874f/attachment-0001.html

------------------------------

Message: 4
Date: Sat, 9 Jul 2016 12:32:42 +0100
From: Kaiming Yin <yinkaiming6987@gmail.com>
Subject: [Freesurfer] Cerebral volume without ventricles
To: FreeSurfer Mailing List <freesurfer@nmr.mgh.harvard.edu>
Message-ID:
     <CAEOOxG0uyz9qT3Wzib7++ZkYQoRYZFY_uXK-5LztPH5MyuogwA@mail.gmail.com>
Content-Type: text/plain; charset="utf-8"

Dear FS experts,

Could I ask which value in the aseg.stats file represent the whole cerebral
volume without ventricles?

Thanks,
Kaiming
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mail.nmr.mgh.harvard.edu/pipermail/freesurfer/attachments/20160709/d20cb77c/attachment-0001.html

------------------------------

Message: 5

Date: Sat, 9 Jul 2016 11:27:32 -0400
From: Douglas Greve <greve@nmr.mgh.harvard.edu>
Subject: Re: [Freesurfer] FSFAST error: ERROR: fast_selxavg3()
     failed\n
To: Erik Jahner <erik.jahner@icloud.com>,      free surfer
     <freesurfer@nmr.mgh.harvard.edu>
Message-ID: <db1bb81e-4481-03db-b261-af94a86480c4@nmr.mgh.harvard.edu>
Content-Type: text/plain; charset="windows-1252"

what do you mean nothing happens? After you hit the return key, it
should either give you the path to the command or it should say Command
not found.

doug


ps. Please remember to post to the list and not to us individually


On 7/8/16 7:16 PM, Erik Jahner wrote:

When I type "which matlab? in the command window, nothing happens. My

$PATH structure is a s follows:

echo $PATH

/Applications/freesurfer/bin:/Applications/freesurfer/fsfast/bin:/Applications/freesurfer/tktools:/usr/local/fsl/bin:/Applications/freesurfer/mni/bin:/opt/local/bin:/opt/local/sbin:/opt/local/bin:/opt/local/sbin:/usr/local/bin:/usr/bin:/bin:/usr/sbin:/sbin:/opt/X11/bin

Sorry if my response seems na?ve, I?am new to running command

lines although I?m trying learn.

Erik

Message: 6

Date: Fri, 8 Jul 2016 14:30:43 -0400

From: Douglas N Greve <greve@nmr.mgh.harvard.edu

<mailto:greve@nmr.mgh.harvard.edu>>

Subject: Re: [Freesurfer] FSFAST error: ERROR: fast_selxavg3()

failed\n

To: freesurfer@nmr.mgh.harvard.edu <mailto:freesurfer@nmr.mgh.harvard.edu>

Message-ID: <577FF153.1080602@nmr.mgh.harvard.edu

<mailto:577FF153.1080602@nmr.mgh.harvard.edu>>

Content-Type: text/plain; charset=windows-1252; format=flowed

Is matlab in your path? What happens if you type

which matlab

at the command line?

On 07/08/2016 12:34 AM, ERIK JAHNER wrote:

Hello FreeSurfer Developers,

This is Erik Jahner, I?m attempting to complete the final step in the

tutorial concerning FSFast Functional connectivity walkthrough

(http://surfer.nmr.mgh.harvard.edu/fswiki/FsFastFunctionalConnectivityWalkthrough ),

but for I get the following error when I run the following command:

selxavg3-sess -s subj439rest -a fc.lpccseed.surf.lh

Mon Jul 4 11:01:41 PDT 2016

rawfunc2surf-sess completed

----------------------------------------

Started at Mon Jul 4 11:01:40 PDT 2016

Ended at Mon Jul 4 11:01:41 PDT 2016

preproc-sess done

-------------------------------------------

/Applications/freesurfer/LagStructure/subj439rest

Mon Jul 4 11:01:41 PDT 2016

anadir =

/Applications/freesurfer/LagStructure/subj439rest/rest/fc.lpccseed.surf.lh

DoGLMFit = 1

DoContrasts = 1

UpdateNeeded = 1

------------------------------------------

------- matlab output --------------------

matlab: Command not found.

------------------------------------------

*ERROR: fast_selxavg3() failed\n*

I've searched the list and no similar errors have been reported. I

have tried to reinstall MATLAB; I have tried following the

instructions to install the runtime environment; I have tried redoing

all analysis from scratch. I have noticed that STEP 5 does not work

unless I change the -cfg from ?L_Posteriorcingulate.config?

to ?mean.L_Posteriorcingulate.config?. I don?t think this is the

problem, but I thought I would point it out.

My session directory is ?subj439rest?

Instead of titling the sub directory ?bold?, I titled it ?rest?

1) FreeSurfer version: freesurfer-Darwin-lion-stable-pub-v5.3.0

2) Platform: OSX Yousemite

3) uname -a: Darwin ERIKs-iMac.local 15.5.0 Darwin Kernel Version

15.5.0: Tue Apr 19 18:36:36 PDT 2016;

root:xnu-3248.50.21~8/RELEASE_X86_64 x86_64

4) selxavg3-sess-rest-fc.lpccseed.surf.lh-160704110140.log: see attached

___________________________________________________

Here are the command lines I put in for steps 4-8, incase this is of

some use:

STEP 4

fcseed-config-segid 1010 -fcname L_Posteriorcingulate.dat -fsd rest

-mean-cfg mean.L_Posteriorcingulate.config

STEP 5

fcseed-sess-s subj439rest -cfg mean.L_Posteriorcingulate.config

STEP 6

* fcseed-config -wm -fcname wm.dat -fsd rest -pca -cfg wm.config

* fcseed-sess -s subj439rest -cfg wm.config

* fcseed-config -vcsf -fcname vcsf.dat -fsd rest -pca -cfg vcsf.config

* fcseed-sess -s subj439rest -cfg vcsf.config

STEP7

mkanalysis-sess-analysis fc.lpccseed.surf.lh -surface fsaverage lh

-fwhm 5 -notask-taskreg L_Posteriorcingulate.dat 1 -nuisreg vcsf.dat 5

-nuisreg wm.dat 5 -mcextreg-polyfit 5 -nskip 4 -fsd rest -TR 2 -per-run

STEP 8

selxavg3-sess-s subj439rest -a fc.lpccseed.surf.lh

_______________________________________________

Freesurfer mailing list

Freesurfer@nmr.mgh.harvard.edu <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>

Phone Number: 617-724-2358

Fax: 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

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/

------------------------------

-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mail.nmr.mgh.harvard.edu/pipermail/freesurfer/attachments/20160709/a34bb36f/attachment-0001.html

------------------------------

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

End of Freesurfer Digest, Vol 149, Issue 20
*******************************************


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.