[Eeglablist] Numerical Values in STUDY Structure

Saim Rasheed saim_rasheed at hotmail.com
Wed Jun 30 09:26:02 PDT 2010



Thankyou very much Julie and Arno,
 
Your help really worked. I have got the ERSP and itc data into Matlab workspace. I have got another problem.
I have recorded the data for three different visual stimuli conditions from 4 electrode positions at 256 samples per second. There are certain number of time-locked trials, e.g. 60, 58 55 for three conditions respectively. Each trial consists of 3 seconds of data containing 768 data points, which were automatically processed into 200 time points and 100 frequency points from 3Hz to 128Hz. Now I can see in workspace, 100*200 variable each for all electrode positions as chan1_ersp, chan2_ersp, chan3_ersp and chan4_ersp. Trial information is bit confused for me now. Is it lost now? Is it possible to compute ERSP and ITC trial by trial because I need to process a single-trial for each visual condition.
Thanks
 
Saim

 

> Subject: Re: [Eeglablist] Numerical Values in STUDY Structure
> From: arno at ucsd.edu
> Date: Tue, 29 Jun 2010 17:52:09 -0700
> To: saim_rasheed at hotmail.com; eeglablist at sccn.ucsd.edu
> 
> 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:
> > load_string =
> > [basedir,subjs{subj},'\',STUDY.datasetinfo(1).filename(1:end-4),'.icaersp'];
> > 
> > % OR ITC
> > load_string =
> > [basedir,subjs{subj},'\',STUDY.datasetinfo(1).filename(1:end-4),'.icaitc'];
> > 
> > % OR spectra:
> > load_string =
> > [basedir,subjs{subj},'\',STUDY.datasetinfo(1).filename(1:end-4),'.icaspec'];
> > 
> > % OR topomaps:
> > load_string =
> > [basedir,subjs{subj},'\',STUDY.datasetinfo(1).filename(1:end-4),'.icatopo'];
> > 
> > % load the raw data into the variable 'RAWdata':
> > RAWdata = load('-mat',load_string);
> > 
> > % 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?
> >> 
> >> 
> >> 
> >> I need them for further computations. Please help.
> >> 
> >> 
> >> 
> >> 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
> 
 		 	   		  
_________________________________________________________________
Hotmail: Trusted email with powerful SPAM protection.
https://signup.live.com/signup.aspx?id=60969
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://sccn.ucsd.edu/pipermail/eeglablist/attachments/20100630/6f414be6/attachment-0001.html 


More information about the eeglablist mailing list