Dear Vladimir,
I tried using "extract" in GUI and I entered -56 -8 -24;-20 -96 -8 to extract the signals of these two sources. When I clicked on "mip" the attached figure was displayed. As you see, the sources are different from the ones I have defined. It seems that it has picked the coordinates in a zigzag form. I also changed my definition type to [-56 -8 -24];[-20 -96 -8] and [-56 -8 -24;-20 -96 -8] but the result was the same. So, I checked the D.inv{1}.source.XYZ and found out that the sources are saved correctly here. Now, I am not sure if this is just a display error or the sources are really extracted from wrong coordinates. I am trying to make some small changes to the "extract" code to be able to call it from another mfile and not just GUI.
Also, when the extract window pops up and asks for the name of sources, I try different kinds of source definitions like 's1' 's2' or 's1';'s2' or 's1','s2' but it saves the whole phrase as the name of the first source and it defines the name of the second source itself as 'Source 2'.
Best,
Pegah
--------------------------------------------------------
Pegah Tayaranian Hosseini
PhD Student
Room 4077, Tizard building (13)
Institute of Sound and Vibration Research
University of Southampton, SO17 1BJ, UK
Tel: 023 8059 2850
email: [log in to unmask]
________________________________________
From: Vladimir Litvak [[log in to unmask]]
Sent: 28 June 2011 10:49
To: Tayaranian Hosseini P.
Cc: [log in to unmask]
Subject: Re: [SPM] Inversion Matrix
Dear Pegah,
You can export your reconstruction to an image and use the attached
code to find local peaks in that image.
Best,
Vladimir
On Tue, Jun 28, 2011 at 10:36 AM, Tayaranian Hosseini P.
<[log in to unmask]> wrote:
> Dear Vladimir,
>
> I know the function but the problem is when my source reconstruction is complete, I do not know where the active sources are situated in millimeters. I need to select four or five of the most active sources and extract their signals but as I do not know their positions in millimeters, I thought maybe if I have a matrix of all the sources, I can select the rows which have higher intensities.
> Would you please tell me how I can have a rough estimation of source positions in millimeters?
>
> Best,
> Pegah
>
> --------------------------------------------------------
> Pegah Tayaranian Hosseini
> PhD Student
> Room 4077, Tizard building (13)
> Institute of Sound and Vibration Research
> University of Southampton, SO17 1BJ, UK
>
> Tel: 023 8059 2850
> email: [log in to unmask]
> ________________________________________
> From: Vladimir Litvak [[log in to unmask]]
> Sent: 23 June 2011 18:24
> To: Tayaranian Hosseini P.
> Cc: [log in to unmask]
> Subject: Re: [SPM] Inversion Matrix
>
> Dear Pegah,
>
> The functionality you are trying to develop is already available in
> spm_eeg_inv_extract in the latest SPM update. This function is
> accessible from GUI and batch.
>
> Best,
>
> Vladimir
>
>
>
> On 23 Jun 2011, at 17:49, "Tayaranian Hosseini P." <[log in to unmask]> wrote:
>
>> Dear SPMers,
>>
>> I am going to try doing some simple manipulations in extracting source signals but I have some problems.
>>
>> * I could find out that the J matrix is a conditional expectation in i conditions. What are these i conditions?
>>
>> * I can see that the T matrix has the same number of conditions. Would you please explain what this T matrix does?
>>
>> * After the inversion is complete, does M*U*data(in the desirable time course) give me the signal of the sources (in that specific time course) or there exists any other step that should be done?
>>
>> Thanks for your help.
>>
>> Best,
>> Pegah
>
|