<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40">
<head>
<meta http-equiv=Content-Type content="text/html; charset=us-ascii">
<meta name=Generator content="Microsoft Word 12 (filtered medium)">
<style>
<!--
/* Font Definitions */
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
@font-face
{font-family:Consolas;
panose-1:2 11 6 9 2 2 4 3 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0in;
margin-bottom:.0001pt;
font-size:11.0pt;
font-family:"Calibri","sans-serif";}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:blue;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:purple;
text-decoration:underline;}
p.MsoPlainText, li.MsoPlainText, div.MsoPlainText
{mso-style-priority:99;
mso-style-link:"Plain Text Char";
margin:0in;
margin-bottom:.0001pt;
font-size:10.5pt;
font-family:Consolas;}
span.EmailStyle17
{mso-style-type:personal-compose;
font-family:"Calibri","sans-serif";
color:windowtext;}
span.PlainTextChar
{mso-style-name:"Plain Text Char";
mso-style-priority:99;
mso-style-link:"Plain Text";
font-family:Consolas;}
.MsoChpDefault
{mso-style-type:export-only;}
@page WordSection1
{size:8.5in 11.0in;
margin:1.0in 1.0in 1.0in 1.0in;}
div.WordSection1
{page:WordSection1;}
-->
</style>
<!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]-->
</head>
<body lang=EN-US link=blue vlink=purple>
<div class=WordSection1>
<p class=MsoPlainText>Dear EEGLAB experts,<o:p></o:p></p>
<p class=MsoPlainText><b><u>If you have a background in EEGLAB scripting and
can help me with this, I would greatly appreciate this</u>. <o:p></o:p></b></p>
<p class=MsoPlainText><b><o:p> </o:p></b></p>
<p class=MsoPlainText><b>I have been referring to this, but <u>not</u> able to
solve the problem:<u> </u></b><a
href="http://sccn.ucsd.edu/wiki/Chapter_02:_Writing_EEGLAB_Scripts">http://sccn.ucsd.edu/wiki/Chapter_02:_Writing_EEGLAB_Scripts</a><o:p></o:p></p>
<p class=MsoPlainText><o:p> </o:p></p>
<p class=MsoPlainText>Recently, I have had great trouble running ICA on my data
using the GUI due to memory constraints. We have not been able to fix this
otherwise. So, I have moved to the script level using a low memory script
level version of ICA. I load my dataset into the GUI and then run the
script through the GUI. Here is the script I have been using.<o:p></o:p></p>
<p class=MsoPlainText><o:p> </o:p></p>
<p class=MsoPlainText>[EEG.icaweights EEG.sphere] =
runicalowmem(reshape(EEG.data, EEG.nbchan, EEG.trials*EEG.pnts));<o:p></o:p></p>
<p class=MsoPlainText>EEG = eeg_checkset(EEG);<o:p></o:p></p>
<p class=MsoPlainText>[ALLEEG EEG CURRENTSET] = eeg_store(ALLEEG, EEG,
CURRENTSET);<o:p></o:p></p>
<p class=MsoPlainText>eeglab redraw;<o:p></o:p></p>
<p class=MsoPlainText><o:p> </o:p></p>
<p class=MsoPlainText>Here is the command window output:<o:p></o:p></p>
<p class=MsoPlainText>“Input data size [62,1048575] = 62 channels,
1048575 frames<o:p></o:p></p>
<p class=MsoPlainText>Finding 62 ICA components using logistic ICA.<o:p></o:p></p>
<p class=MsoPlainText>Decomposing 272 frames per ICA weight ((3844)^2 = 1048575
weights, Initial learning rate will be 0.000157494, block size 70.<o:p></o:p></p>
<p class=MsoPlainText>Learning rate will be multiplied by 0.9 whenever
angledelta >= 60 deg.<o:p></o:p></p>
<p class=MsoPlainText>More than 32 channels: default stopping weight change
1E-7<o:p></o:p></p>
<p class=MsoPlainText>Training will end when wchange < 1e-007 or after 512
steps.<o:p></o:p></p>
<p class=MsoPlainText>Online bias adjustment will be used.<o:p></o:p></p>
<p class=MsoPlainText>Removing mean of each channel ...<o:p></o:p></p>
<p class=MsoPlainText>Final training data range: -2048.65 to 1315.39<o:p></o:p></p>
<p class=MsoPlainText>Computing the sphering matrix...<o:p></o:p></p>
<p class=MsoPlainText>Starting weights are the identity matrix ...<o:p></o:p></p>
<p class=MsoPlainText>Sphering the data ...<o:p></o:p></p>
<p class=MsoPlainText>Beginning ICA training ...<o:p></o:p></p>
<p class=MsoPlainText>Warning: Using 'state' to set RAND's internal state
causes RAND, RANDI, and RANDN<o:p></o:p></p>
<p class=MsoPlainText>to use legacy random number generators. <o:p></o:p></p>
<p class=MsoPlainText>> In runicalowmem at 896<o:p></o:p></p>
<p class=MsoPlainText> In pop_runscript at 50<o:p></o:p></p>
<p class=MsoPlainText>step 1 - lrate 0.000157, wchange 5985964.09709237,
angledelta 0.0 deg<o:p></o:p></p>
<p class=MsoPlainText>.<o:p></o:p></p>
<p class=MsoPlainText>.<o:p></o:p></p>
<p class=MsoPlainText>.<o:p></o:p></p>
<p class=MsoPlainText>step 127 - lrate 0.000000, wchange 0.00000008, angledelta
50.9 deg<o:p></o:p></p>
<p class=MsoPlainText>Sorting components in descending order of mean projected
variance ...”<o:p></o:p></p>
<p class=MsoPlainText><o:p> </o:p></p>
<p class=MsoPlainText><o:p> </o:p></p>
<p class=MsoPlainText>It appears from the command window that ICA ran on the
data. However, the current loaded dataset that ICA was run on says there
are no ICA weights (<b>GUI attached</b>). When I go to plot components,
it says ICA needs to be run first. <b><u>So, clearly it ran but did not
save the ICA to the dataset.</u></b> <o:p></o:p></p>
<p class=MsoPlainText><o:p> </o:p></p>
<p class=MsoPlainText><b>I thought that this following command would overwrite the
variables in the workspace:</b> [ALLEEG EEG CURRENTSET] = eeg_store(ALLEEG,
EEG, CURRENTSET);<o:p></o:p></p>
<p class=MsoPlainText><o:p> </o:p></p>
<p class=MsoPlainText><b>I thought this following command would update the
display in the GUI:</b> eeglab redraw;<o:p></o:p></p>
<p class=MsoPlainText><o:p> </o:p></p>
<p class=MsoPlainText><b><u>However, this is not working.<o:p></o:p></u></b></p>
<p class=MsoNormal><o:p> </o:p></p>
<p class=MsoNormal>Please help me if you can.<o:p></o:p></p>
<p class=MsoNormal>Thanks,<br>
Ed Modestino<o:p></o:p></p>
<p class=MsoNormal><o:p> </o:p></p>
<p class=MsoNormal>Edward Justin Modestino, Ph.D.<o:p></o:p></p>
<p class=MsoNormal>Postdoctoral Research Associate<o:p></o:p></p>
<p class=MsoNormal>Ray Westphal Neuroimaging Laboratory<o:p></o:p></p>
<p class=MsoNormal>Division of Perceptual Studies<o:p></o:p></p>
<p class=MsoNormal>Department of Psychiatry and Neurobehavioral Sciences<o:p></o:p></p>
<p class=MsoNormal>University of Virginia<o:p></o:p></p>
<p class=MsoNormal>Email: ejm9f@virginia.edu<o:p></o:p></p>
<p class=MsoNormal><o:p> </o:p></p>
</div>
</body>
</html>