# [Eeglablist] Numerical Values in STUDY Structure

Arnaud Delorme arno at ucsd.edu
Tue Jun 29 17:52:09 PDT 2010

Another simple solution is to add additional output to the STUDY plotting functions after copying them from the history.
For instance, after creating or loading an EEGLAB STUDY, then precomputing the ERP for data channel, you may use the channel plotting interface to plot the ERP for electrode CZ (for instance). Then, the EEGLAB history contains the line

STUDY = std_erpplot(STUDY, ALLEEG, 'channels', { 'cz' });

Then, looking at the std_erpplot function help, you may add additional outputs

[STUDY erp erptimes ] = std_erpplot(STUDY, ALLEEG, 'channels', { 'cz' });

Note that the "erp" output will be a cell array containing an array for the ERP of each condition/group. If one subject only is present, each array will be a vector of values for the ERP of this subject in a given condition. If several subjects are present, each column of the array will contain the ERP for a given subject.

Arno

On Jun 29, 2010, at 10:46 AM, Julie Onton wrote:

> Good question.
>
> Below is a script to help you find, load and plot your raw data.
>
> %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
> % LOAD raw data:------------------------
> %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
> basedir = '/home/you/whereyourdatais/';
>
> subjs =
> {'S01','S02','S03','S04','S05','S06','S07','S08','S09','S10','S11','S12','S13'};
>
> subj = 3; % who do you want to plot?
> cond = 1; % what condition number?
>
> % for example, load ERSP:
> [basedir,subjs{subj},'\',STUDY.datasetinfo(1).filename(1:end-4),'.icaersp'];
>
> % OR ITC
> [basedir,subjs{subj},'\',STUDY.datasetinfo(1).filename(1:end-4),'.icaitc'];
>
> % OR spectra:
> [basedir,subjs{subj},'\',STUDY.datasetinfo(1).filename(1:end-4),'.icaspec'];
>
> % OR topomaps:
> [basedir,subjs{subj},'\',STUDY.datasetinfo(1).filename(1:end-4),'.icatopo'];
>
> % load the raw data into the variable 'RAWdata':
>
> % Mask an ERSP using calculated bootstrap limits (if you calculated them):
>
> ic = 8; % choose IC to plot
>
> oneic = ['RAWdata.comp',int2str(ic),'_ersp'];
> oneic = eval(oneic);
> % in case you want to see the removed baseline:
> onebase = ['RAWdata.comp',int2str(ic),'_erspbase'];
> onebase = eval(onebase);
> % load the bootstrap significance limits:
> oneboot = ['RAWdata.comp',int2str(ic),'_erspboot'];
> oneboot = eval(oneboot);
>
> maskERSP = oneic;
> % zero out non-sig values:
> maskERSP(find(oneic > repmat(oneboot(:,1),[1 size(oneic,2)])& oneic <
> repmat(oneboot(:,2),[1 size(oneic,2)]))) = 0;
>
> clim = 4; % set +/- color limits
> figure; imagesc(RAWdata.times,RAWdata.freqs,maskERSP,[-clim clim]);
> set(gca,'ydir','norm');
> title(['Subj ',int2str(subj),'; IC ',int2str(ic),'; ',STUDY.condition{cond}]);
> cbar;
>
>
> Hope this helps,  Julie
>
> --
> Julie Onton, PhD
> http://sccn.ucsd.edu/~julie
>
>>
>>
>>
>> Hi,
>>
>>
>>
>> I am working with STUDY structure and computed ERSP, ITC and. I can plot them
>> by a single click and analyse them.
>>
>> Where I can find the numerical values for ERSP and ITC ? Where these values
>> are stored?
>>
>>
>>
>>
>>
>>
>> Thanks.
>>
>>
>>
>> Saim
>>
>> _________________________________________________________________
>> Your E-mail and More On-the-Go. Get Windows Live Hotmail Free.
>> https://signup.live.com/signup.aspx?id=60969_______________________________________________
>> Eeglablist page: http://sccn.ucsd.edu/eeglab/eeglabmail.html
>> To unsubscribe, send an empty email to eeglablist-unsubscribe at sccn.ucsd.edu
>> For digest mode, send an email with the subject "set digest mime" to
>> eeglablist-request at sccn.ucsd.edu
>
> _______________________________________________
> Eeglablist page: http://sccn.ucsd.edu/eeglab/eeglabmail.html
> To unsubscribe, send an empty email to eeglablist-unsubscribe at sccn.ucsd.edu
> For digest mode, send an email with the subject "set digest mime" to eeglablist-request at sccn.ucsd.edu