<html>
<head>
<style><!--
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
font-size: 10pt;
font-family:Verdana
}
--></style>
</head>
<body class='hmmessage'>
Thanks a lot, I got the trial by trial data from '*.dattimef' files. In '.ersp' files individual trials are not available, Are these values averaged across all the trial values?<BR>In '.datitc' values are visible as complex numbers in Matlab workspace, I write the data into excell file, there it write only real part, not the complex part ?<BR>I have to use the data (after conversion into another format) into LIBSVM for classification of my visual conditions. I dont know, How important is the complex part for ITC computations?<BR><BR>Saim <BR><BR>> Date: Wed, 30 Jun 2010 11:34:10 -0700<BR>> Subject: RE: [Eeglablist] Numerical Values in STUDY Structure<BR>> From: julie@sccn.ucsd.edu<BR>> To: saim_rasheed@hotmail.com<BR>> CC: arno@ucsd.edu; eeglablist@sccn.ucsd.edu<BR>> <BR>> There is an option in the precompute menu to save single-trial data, but you<BR>> will have to recompute if you did not check this the first time around<BR>> (default is NOT to save single trials because this occupies a remarkable<BR>> amount of disk space). ITC, of course, is a value that is computed across<BR>> trials.<BR>> <BR>> Julie<BR>> <BR>> -- <BR>> Julie Onton, PhD<BR>> http://sccn.ucsd.edu/~julie<BR>> <BR>> ><BR>> ><BR>> > Thankyou very much Julie and Arno,<BR>> ><BR>> > Your help really worked. I have got the ERSP and itc data into Matlab<BR>> > workspace. I have got another problem.<BR>> > I have recorded the data for three different visual stimuli conditions from 4<BR>> > electrode positions at 256 samples per second. There are certain number of<BR>> > time-locked trials, e.g. 60, 58 55 for three conditions respectively. Each<BR>> > trial consists of 3 seconds of data containing 768 data points, which were<BR>> > automatically processed into 200 time points and 100 frequency points from 3Hz<BR>> > to 128Hz. Now I can see in workspace, 100*200 variable each for all electrode<BR>> > positions as chan1_ersp, chan2_ersp, chan3_ersp and chan4_ersp. Trial<BR>> > information is bit confused for me now. Is it lost now? Is it possible to<BR>> > compute ERSP and ITC trial by trial because I need to process a single-trial<BR>> > 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<BR>> >> functions after copying them from the history.<BR>> >> For instance, after creating or loading an EEGLAB STUDY, then precomputing<BR>> >> the ERP for data channel, you may use the channel plotting interface to plot<BR>> >> the ERP for electrode CZ (for instance). Then, the EEGLAB history contains<BR>> >> 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<BR>> >> 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<BR>> >> ERP of each condition/group. If one subject only is present, each array will<BR>> >> be a vector of values for the ERP of this subject in a given condition. If<BR>> >> several subjects are present, each column of the array will contain the ERP<BR>> >> 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),';<BR>> >> ',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<BR>> >> them<BR>> >> >> by a single click and analyse them.<BR>> >> >><BR>> >> >> Where I can find the numerical values for ERSP and ITC ? Where these<BR>> >> 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<BR>> >> 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<BR>> >> 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>> > _________________________________________________________________<BR>> > Hotmail: Trusted email with powerful SPAM protection.<BR>> > https://signup.live.com/signup.aspx?id=60969<BR>> <BR><BR> <br /><hr />Your E-mail and More On-the-Go. Get Windows Live Hotmail Free. <a href='https://signup.live.com/signup.aspx?id=60969' target='_new'>Sign up now.</a></body>
</html>