<html>
<head>
<style><!--
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
font-size: 10pt;
font-family:Verdana
}
--></style>
</head>
<body class='hmmessage'>
<BR>Thankyou very much Julie and Arno,<BR> <BR>Your help really worked. I have got the ERSP and itc data into Matlab workspace. I have got another problem.<BR>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.<BR>Thanks<BR> <BR>Saim<BR><BR> <BR><BR>> Subject: Re: [Eeglablist] Numerical Values in STUDY Structure<BR>> From: arno@ucsd.edu<BR>> Date: Tue, 29 Jun 2010 17:52:09 -0700<BR>> To: saim_rasheed@hotmail.com; eeglablist@sccn.ucsd.edu<BR>> <BR>> Another simple solution is to add additional output to the STUDY plotting functions after copying them from the history.<BR>> 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<BR>> <BR>> STUDY = std_erpplot(STUDY, ALLEEG, 'channels', { 'cz' }); <BR>> <BR>> Then, looking at the std_erpplot function help, you may add additional outputs<BR>> <BR>> [STUDY erp erptimes ] = std_erpplot(STUDY, ALLEEG, 'channels', { 'cz' }); <BR>> <BR>> 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.<BR>> <BR>> Arno<BR>> <BR>> On Jun 29, 2010, at 10:46 AM, Julie Onton wrote:<BR>> <BR>> > Good question.<BR>> > <BR>> > Below is a script to help you find, load and plot your raw data.<BR>> > <BR>> > %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%<BR>> > % LOAD raw data:------------------------<BR>> > %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%<BR>> > basedir = '/home/you/whereyourdatais/';<BR>> > <BR>> > subjs =<BR>> > {'S01','S02','S03','S04','S05','S06','S07','S08','S09','S10','S11','S12','S13'};<BR>> > <BR>> > subj = 3; % who do you want to plot?<BR>> > cond = 1; % what condition number?<BR>> > <BR>> > % for example, load ERSP:<BR>> > load_string =<BR>> > [basedir,subjs{subj},'\',STUDY.datasetinfo(1).filename(1:end-4),'.icaersp'];<BR>> > <BR>> > % OR ITC<BR>> > load_string =<BR>> > [basedir,subjs{subj},'\',STUDY.datasetinfo(1).filename(1:end-4),'.icaitc'];<BR>> > <BR>> > % OR spectra:<BR>> > load_string =<BR>> > [basedir,subjs{subj},'\',STUDY.datasetinfo(1).filename(1:end-4),'.icaspec'];<BR>> > <BR>> > % OR topomaps:<BR>> > load_string =<BR>> > [basedir,subjs{subj},'\',STUDY.datasetinfo(1).filename(1:end-4),'.icatopo'];<BR>> > <BR>> > % load the raw data into the variable 'RAWdata':<BR>> > RAWdata = load('-mat',load_string);<BR>> > <BR>> > % Mask an ERSP using calculated bootstrap limits (if you calculated them):<BR>> > <BR>> > ic = 8; % choose IC to plot<BR>> > <BR>> > oneic = ['RAWdata.comp',int2str(ic),'_ersp'];<BR>> > oneic = eval(oneic);<BR>> > % in case you want to see the removed baseline:<BR>> > onebase = ['RAWdata.comp',int2str(ic),'_erspbase'];<BR>> > onebase = eval(onebase);<BR>> > % load the bootstrap significance limits:<BR>> > oneboot = ['RAWdata.comp',int2str(ic),'_erspboot'];<BR>> > oneboot = eval(oneboot);<BR>> > <BR>> > maskERSP = oneic;<BR>> > % zero out non-sig values:<BR>> > maskERSP(find(oneic > repmat(oneboot(:,1),[1 size(oneic,2)])& oneic <<BR>> > repmat(oneboot(:,2),[1 size(oneic,2)]))) = 0;<BR>> > <BR>> > clim = 4; % set +/- color limits<BR>> > figure; imagesc(RAWdata.times,RAWdata.freqs,maskERSP,[-clim clim]);<BR>> > set(gca,'ydir','norm');<BR>> > title(['Subj ',int2str(subj),'; IC ',int2str(ic),'; ',STUDY.condition{cond}]);<BR>> > cbar;<BR>> > <BR>> > <BR>> > Hope this helps, Julie<BR>> > <BR>> > -- <BR>> > Julie Onton, PhD<BR>> > http://sccn.ucsd.edu/~julie<BR>> > <BR>> >> <BR>> >> <BR>> >> <BR>> >> Hi,<BR>> >> <BR>> >> <BR>> >> <BR>> >> I am working with STUDY structure and computed ERSP, ITC and. I can plot them<BR>> >> by a single click and analyse them.<BR>> >> <BR>> >> Where I can find the numerical values for ERSP and ITC ? Where these values<BR>> >> are stored?<BR>> >> <BR>> >> <BR>> >> <BR>> >> I need them for further computations. Please help.<BR>> >> <BR>> >> <BR>> >> <BR>> >> Thanks.<BR>> >> <BR>> >> <BR>> >> <BR>> >> Saim<BR>> >> <BR>> >> _________________________________________________________________<BR>> >> Your E-mail and More On-the-Go. Get Windows Live Hotmail Free.<BR>> >> https://signup.live.com/signup.aspx?id=60969_______________________________________________<BR>> >> Eeglablist page: http://sccn.ucsd.edu/eeglab/eeglabmail.html<BR>> >> To unsubscribe, send an empty email to eeglablist-unsubscribe@sccn.ucsd.edu<BR>> >> For digest mode, send an email with the subject "set digest mime" to<BR>> >> eeglablist-request@sccn.ucsd.edu<BR>> > <BR>> > _______________________________________________<BR>> > Eeglablist page: http://sccn.ucsd.edu/eeglab/eeglabmail.html<BR>> > To unsubscribe, send an empty email to eeglablist-unsubscribe@sccn.ucsd.edu<BR>> > For digest mode, send an email with the subject "set digest mime" to eeglablist-request@sccn.ucsd.edu<BR>> <BR> <br /><hr />Hotmail: Trusted email with powerful SPAM protection. <a href='https://signup.live.com/signup.aspx?id=60969' target='_new'>Sign up now.</a></body>
</html>