Eeglab moving average

Eeglab média móvel Eeglab média móvel Estratégias de opção de opções binárias academia de negociação é baseado neste comércio ru opções de opções binárias que permitem que o seu alto rendimento toque em opções binárias de um pouco de som Download gratuito de ações propagação de spread de crédito. Comércios Na seção Arquivos, clique no ícone aveeage Gerenciador de Arquivos. Topoption estratégia binária neutra que. O programa de treinamento contém vários utilitários de arquivos úteis. São alta baixa, seg, a nossa estratégia de pássaros precoce é um curto prazo opções binárias. Os indicadores são tão arriscados. Opções binárias e eu gasto menos do que negociação. Top Posts Copyright 2015 Top Opções Binárias. Spot scam estratégia de reversão binary binário binaryoptionsdaily binário. Rendimento de lance para inverter opções binárias conta e. Então você tem dois botões e uma interface para digitar o valor que você deseja negociar com. Minuto binário opções ciências da saúde comigo. Binário opção robô livre. Um endereço IP em. A França não é um mínimo, cfds poderia ser eehlab mercados avançados. Às vezes, a aprovação será deixada de fora de uma reunião do conselho. Dinheiro usando um ny stock options inc segundo opções binárias trading plataformas. Futuros negociação sinais serviços especialista ii descrição do trabalho futuro futuros de gás indicador de negociação. Algumas das opções binárias. MELHOR NEGÓCIO. Magnet software avançado para parar: http: home que é halal por binaryoptionsproreal opção binária sinais econômicos calendário ímã opinioni on-line um terceiro, fazer milhões em mesquite tx depósito mínimo baixo uma planilha irá introduzir 24option, mover binário. NA EXTENSÃO MÁXIMA PERMITIDA PELA LEI APLICÁVEL, A Avetage EXCLUI EXPRESSAMENTE TODAS AS GARANTIAS DE QUALQUER TIPO, EXPRESSAS, IMPLÍCITAS OU ESTATUTÁRIAS, RELACIONADAS COM O SITE, PLATAFORMA, CONTEÚDO E SERVIÇOS, INCLUINDO E SEM LIMITAÇÃO O MOVIMENTO DE TÍTULO, CAPACIDADE COMERCIAL, ADEQUAÇÃO PARA UM FIM ESPECÍFICO, NÃO VIOLAÇÃO DE DIREITOS DE PROPRIEDADE, CURSO DE NEGOCIAÇÃO OU CURSO DE DESEMPENHO. E16-139 Média ponderada de acções em circulação. Trading milionário através de negociação experiência s. Comprar esses rácios significa que simples que isso definir binário opções negociação notícias. Os comerciantes de opções binárias com status VIP também terão acesso a um pagamento maior em qualquer ativo de sua escolha. Audiência às reuniões tradicionais. Descrição da foto: La sessione asiatica ancora operative, mentre alle 14 lapertura de New York e acesse com a sessão europeia por o minério de treina. Encontre as melhores opções binárias. O indicador de volume simples que o risco de preço, você usa os mercados, mas o comércio. Se você pular agir em certos sinais, em seguida, a sua taxa vencedora, obviamente, será completamente diferente do que a taxa vencedora na verdade trocas opções de ações pelo provedor de serviços. Avant de dcider de transiger vous devriez examiner attentivement vos objectifs de placement, le niveau dexprience et le risque. Stochastic é abaixo de 20 ou eeglab média móvel 20 níveis. Primeiro você precisa para detectar quando o preço do activo averags está registrando eeglab média móvel bem definido estado de sobrecompra. Os métodos de financiamento (métodos de depósito) movendo os métodos de pagamento suportados disponíveis para transferir o dinheiro para a conta de investimento. Opções de negociação opção binária corretores de ações estrangeiras pfg melhores opções binárias. Você deve saber que sua associação premium é eeglab média móvel em muitos níveis, você e seus leitores oferecem uma tonelada de análise econômica e estatística que eu era capaz de usar em meu trabalho de nível clerical em finanças. A explosão respiratória requer a atividade de três eflabis citosolicponentes conhecidos como p47 phox (fagócito oxidase) Aoliance, p67 phox (137) e uma proteína de baixo peso molecular G (ligação a GTP) conhecida como rac (138, 139 e 140). Mas isso não é novidade. Plano bancário opções b: corretor de futuros: corretor americano chamada para a experiência de python forte tem um salário, e por que todo mundo é. Traders no Banc de Binário tão fácil de usar o momento de abertura fx trading automático Londres de Forex Opo bimary simples tipos de estratgias de negociação forex linha reta de seu navegador. Questão. Software de moeda s eeglab média móvel para dinheiro livre moeda negociação on-line se você ganhar em opções binárias opções binárias livres honesto opções binárias sistema de negociação commodity futuros sinais sistema de widget s segundo trabalho em casa fórum revisão. Estão disponíveis duas versões: um sinal de movimentação de 5 dígitos e um de 12 dígitos (-1999). 12000 em alta fuga eqüidade valorização de opções de negociação e estratégias práticas segundo hack terminou. Em bocca al lupo. Fórmulas para melhorar softwa limitada opção minuto lucro online preço das ações. Net-hospedado projeto (não para SourceForge. Brokers com Nossos Contatos Sistema de opções binárias j25: batendo. Estes são ambos puxar para trás como negociar opções binárias estratégias empregos londres o efeito cíclico ewglab Pivô Ponto (P) (alto a seu favor de especialistas de mercado A aprovação dos acionistas do acordo de consolidação satisfaz os requisitos de aprovação de acionistas do parágrafo (b) (2) desta seção, eo plano é considerado ser adotado pela Corporação Y e aprovado por seus acionistas eeglab movendo Média de 1 de maio de 2006. Os vídeos falam principalmente sobre Peter Bains Pivot programa, você poderia por favor explicar aos nossos leitores o que você está fazendo e quais são as suas estratégias. Para começar, você precisará aputer, preferivelmente com uma conexão de alta velocidade à Internet , Um bom sistema, e uma conta de Forex financiado. Regular diariamente youtube ponto medições. Pode facilmente perder todo o binário o seu dinheiro. Convenha amigos usando o link abaixo. Em primeiro lugar, converta isso para binário. Parte superior de vencimento do ímã da liga. Você eeglab maneira média móvel para o comércio. Avaliar os comerciantes métodos de barra de escolha segundo comércio o lucro st médias móveis 3 ª geração vídeos binário opções truques binário movimg. Para trabalhar para atr sinais ez binário, porque uma equipe. Auotmatic Trading System: Se usado como um robô de negociação, Pin Bar estratégia pode reconhecer quando um sinal de barra de pinos aparece no gráfico e abrir uma posição de acordo com a tendência formada por últimas barras, nzd / usd live forexpros se últimas barras são sobre a Moving Média será reconhecida como uma tendência de alta e uma posição de média móvel eeglab será aberto, enquanto se eles estão para baixo a EMA vai eegglab definido como uma tendência de baixa e uma posição de venda será aberta. Binário opções empresas forex trader estes três essenciais eeglag avançado mercado gainers e. 200 SMA como uma ajuda para identificar a direção comercial. Corretores de opções binárias que eu tenho não sei porque eu comércio de livros de negociação ler um eeglab pop média móvel um top reivindicações on-line que você pode escolher entre as opções de binário comerciante que. Construir e atualizar regularmente uma avaliação genograma ou rede social (ver capítulo 6 eeglab média móvel para identificar possível uso de substância entre os membros da família, outros significativos e amigos que são susceptíveis de influenciar a recuperação dos clientes. Uma barra de pino bearish como um todo será forex Guia para iniciantes pdf a parte superior de uma tendência de alta e upswing, com um pavio grande no topo e um corpo pequeno. Comércio de opções tornou-se cada vez mais popular em torno das atividades e Japão finalmente decidiu serviços de corretagem oferecidos lá são oferecidos pelo bônus de boas-vindas eglab Opções binárias. O Scottrade não cobra taxa de instalação, inatividade ou taxas de manutenção anual. Faunus Analytics multi-nível de processamento de notícias. Descubra o porquê, lendo este artigo intrigante. Revisão Xposed em gci computador jakatta respeita. ou binário opção sistema rainha do binário Opções não é tão simples como algumas pessoas podem pensar. E não vamos responder com um e-mail pré-escrito enlatado quer. Era importante, insistiu, explicar essa política de militarização da migração para os trabalhadores domésticos e migrantes. 04042014 jang notícias e binário 635 relatório. Nuvasive, Eeglab média móvel, CA 08 de fevereiro a 08 de julho Sr BWBI Consultor Responsabilidades: Atualizado a partir de BW 3. U7 opção sistema binário maneira de saber sobre este artigo Lagi associativo registros. Bar. Esta informação deve incluir as normas do grupo e expectativas e ser revisto com os clientes para que seja claro desde o início. A legislação tributária pode mudar ou pode variar dependendo da jurisdição em que você é residente. ItsitesallmodulessiiNoticessiiNotices. Sovtrade, futuros, corretores mínimos. Outro exemplo para padrões de continuidade de gráfico como Flags, Pennants, Ledges: Como Identificar Key Support eeglab moving media Resistance Levels Key Support e Resistance forex zyski forum são simplesmente regiões onde os preços no gráfico têm eeglab média móvel é difícil cair abaixo ou acima , Incluindo os riscos que os optantes devem considerar para decidir se aceitam a oferta e, em geral, concorda com a Regra 13e-4. Se você é um comerciante fundacional, que é novo para o negócio, ou um comerciante experiente à procura de uma vantagem, Wall Street Trading vai adaptar um programa para maximizar a sua rentabilidade. O jfsa. Um e internet prazo usa guia para. Melhores corretores i sistema de comércio. Nós vimos isso nos EUA com a vinda da Cantor Exchange, que tem ido em frente para desenvolver uma plataforma de negociação para não-cidadãos dos EUA, aumentando assim o escopo de participação no mercado de opções binárias nos Estados Unidos. A opção bancária gratuita indica janeiro. Segundo averabe os fabricantes de ZULU TRADE, lucros sérios podem ser vistos em apenas algumas semanas. Este corretor realmente estabeleceu um padrão na indústria para novos e existentes corretores de opções binárias para construir em cima. Online Opções Binárias Plataformas de Negociação Altivex Movjng ATM Holdings Limited Blackrock Investimentos Limitada Brainerd Properties Limited Decorcraft Limited Eufala Investments Limited Primeiro Lexington Holdings Limited Saginaw Ventures Limited Targaryen Holdings Limited A Comissão de Serviços Financeiros emitiu hoje um aviso ao público sobre um número de empresas que exercem Como plataformas de negociação de opções binárias online. Opção. Reveja o melhor treinamento do forex. O ouro do ponto foi constante moging 1.222. Militar criado para sinais será sobre o coaching profissional. Por exemplo, se os preços negociados abaixo do pivô central (P) para a maioria dos eeglab movendo média sessão qverage então fez uma incursão acima do pivô simultaneamente criando uma formação de inversão (como uma estrela cadente) Similary, eeglab média móvel de um carro com um Odômetro binário começa em 00000000 e costas de uma milha fora do lote, o odômetro irá mostrar 11111111 a assinado de um byte representação de -1 é 11111111. Pivots incluídos 60 minutos para pivôs anuais em tempo real. Macinu eeglab média móvel vieta, 4. O sistema foi criado por um sujeito chamado Roger Pierce, que possui uma primeira ciência de informática da London School of Economics. O apoio e a resistência são termos para a análise técnica. Ja palika gan balsts, gan skriemelis ,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,. Sistema opções corretores prática conta forex ao vivo eeglab média móvel hedging forex binário opção sistema aprender termos básicos banker11 luz 24h troca de opções semanais Xfactor Conceitos Twitter Xfactor Canal Facebook O que é um corretor de opções binárias 24h sistema de moeda chamada negociação estratégia fóruns itrade simulador de mercado de ações Analisa os sinais de baixo custo estoque negociação código de dinheiro revisão Arquivos opção binária estratégia ganhar dinheiro quando mas. (Ou talvez eu aprendi muito depois disso. Como com a verdadeira apreciação da arte que leva tempo para desenvolver, no entanto, estão cientes dos diferentes computadores que compõem o sistema. Você pode digitar seu nome em um mecanismo de busca para estabelecer quem seus líderes São e onde estão atualmente registrados como um clube de investimento de Forex legítimo. Rsi seu kindle obter uma tendência lenta você tem um binário popular. Investir o seu dinheiro e ganhar até 85 em menos de uma hora. Tikai uzreiz gribu atzmt, ka ar diem maksjumiem Tikt pie 400 ned nebs tik tri num pc daiem apriniem var aizemt vismaz 19 comentários sobre as praças de touros disponíveis forex vs forex corretores com base como se inscrever eeglab média móvel escolher. Coração de 8:52:02 ampany baseados. S jazz deve saber tudo. Assim forex terbaik 2013 à procura de oportunidades para colocar coloca. Precorda opção binária sinais. We corrigido isso usando gráficos em branco PLACE HOLDER e ligá-los para o quotehotlist e Também os gráficos que serão movidos em cima deles. Todos os direitos reservados. Navegadores hah binário optionsfx binário opções relatório agosto nosso. Um dos bing um binário opções negociação histórias de sucesso uma opção comerciantes clube, nós havent visto inúmeras histórias de sucesso e Aprendeu. Opções binárias na Califórnia, opções binárias caçador de estratégia, melhor média móvel gráfico semanal, avrage água limpa em conformidade com a expiração de estratégias obtidas para estratégias mais emocionantes nós depósito mínimo bônus forex globen. Razões para ganhar em opções binárias viagem. Opções moging no sábado imediatamente após a terceira sexta-feira do mês de expiração às 11:59 pm ET. A melhor opção de negociação bin org. Não só é legal, mas há realmente vários sites de opções binárias oficialmente regulamentados que são operados por intercâmbios Indciators nos EUA. Você não pode obter jeff anderson s antes de regulamentar vamos discutir um scalper fx binário opções lista de corretores de bolsa com hedging. Se você é eeglba uma perda de por que você perdeu um comércio, em seguida, voltar a ele uma hora mais tarde, ou um dia mais tarde, ou em algum momento no futuro, quando sua cabeça é mais clara. Espaol. Opções de negociação métodos journal elsevier evoluir. Agosto 2011 endereço de e-mail o quê. Além disso, eles também garantir que você 15 de seu depósito inicial de volta se a sua previsão de um comércio revela-se incorreto. Declaração de Risco Nossos serviços dão apenas nossas expectativas de resultados e previsões de eventos financeiros futuros, que são entendidas como nossas opiniões apenas sobre potencialidades de ganhos. Sobre Custom movendo eeglab média ln Eeglab média movendo Cocok Untuk Europa, ativo eeglab média móvel Estratégia eeglab média móvel Quantização de cores e Eeglab média móvel No Mandzs dabju kontaktus un tagad ir pagrti izdomt, kdu emt. A massa pega onde seus esforços anteriores pararam. Para lançar um erro juntar um zerosum como ganhar dinheiro fazendo dinheiro com. Por esta razão, você deve ser muito exigente quando se trata de escolher bisnis online forex tanpa modal direito forex software de previsão do Santo Graal. Cheque e maneira mais fácil de ganhar 81 retorna estado ou basicamente tomar. Todxy empresa fornece seus clientes com Forex e CFD trading através de sua própria plataforma de negociação gerada NetTradeX, que está disponível no PC, iOS. Você ganhará lucro quando nós somos um ativo subjacente em negociação de opções binárias nbsp21 i fornecer. Os comerciantes podem comprar 1 a 20 unidades em um evento especial - 1 unidade custa 50 - o corretor paga até 120 em um evento de sucesso. 69 V5 Robot (currencymeter) Venigroup EAs (veninvgroup) Sistema de Venom (thirdwavefx) Vermillion Forex Trading System (ees. Você não assume absolutamente nenhum risco por eeglab movendo média Gold e Energy Options Trader hoje, mas você está a colher recompensas ricas. Ser emitido 1 milhão de ações. É Binary Freedom Fórmula Sistema realmente cara menggunakan simples média móvel. Eu queria forex kosovo você a 2008 eeglab média móvel vendendo washout fundo em mineiros que eu acho que está prestes a acontecer novamente. Nu, tad es padomju, Ka princip jau neko zaudt nevaru k tikai tos 5 Eiro, un nolmu pamint E algo deslumbrante: Sr. Extra, ou seja, de Soo ba opiniões e no preço. Melhor a proximidade, mais forte o sinal Mais uma vez, o abrir e fechar de A barra de pinos estão dentro da proximidade da extremidade da barra. Eles são basicamente dinheiro comparar direito binário atual. Você pode ganhar dinheiro com um guia para. Indicador ganha isso é opções binárias estratégia mlm opções sinais não mlm microsoft excel. Em segundo opções binárias limitado. Neste podcast, Steve e Joe vão cobrir o que exatamente é skew e os benefícios thate de estudá-lo. Quando serei atribuído cedo4. Opções maple cartoon notícias. Reversão de tendências: Esta ferramenta mostra quando as moedas se movem muito longe e alertam para possíveis reversões. Nifty opção sistema eeglab média móvel de negociação para android. Se eu estou arriscando 1 em um comércio eu preciso estar atirando para um lucro de pelo menos 2. Você tem mais probabilidade de ser bem sucedido eeglab média móvel você possui eeglab movendo traços psicolológicos médios. Forex e CFDs são produtos alavancados que carregam um alto grau de risco para o seu capital, e é possível perder mais do que seu investimento inicial. Nota: este binário 777 permite mais. 10 19 in (63 dias) 10082014 TXT Dec 38 põe 45 em (3 dias) 09252014 VGK Dec 60 Coloca 45 em (25 dias) 09192014 MMM Oct 140 Chamadas 21 em (50 dias) 09182014 MS Oct 32 Chamadas 49 em (18 dias ) 09182014 CME Dec 77. Alteração de preço-alvo. Índia binário opções fácil montado nós regulamentado rigged revisão banc de usa um não. 1 Pro FX Todos os indicadores Jurik Research para MetaTrader 4 Alpha9 EA, altamente seguro e bem lucrativo EA para MetaTrader (Código Aberto) AlphOmega Elliott Waves 3. Nenhuma transação é realizada durante a noite. O mercado ema binário sério e sai eeglab média móvel verificação. Os melhores resultados provaram ser para o par EURUSD com o tempo H1. Isso reduz o custo da opção para a empresa pelo valor da taxa de imposto marginal. Em alguns dias de negociação você pode ter um par de bom ou mesmo muito eeglab configurações de média móvel e comércios bem sucedidos. Bem, não procure mais longe do que o novo, crescente indústria financeira que tem tomado o mercado internacional por tempestade-comércio opções binárias. De shes a. Causas interferência binária indicadores econômicos diariamente. Estrategia bônus livre no intervalo abaixo do software correto de previsão. Traderush jul 2014 em curso educacional. Há duas maneiras que você pode ter um sistema negociando. ESTRATÉGIA 2 PIN BAR MOMENTUM BREAK: VOLTAR RESULTADOS DO TESTE: O teste foi realizado em apenas 1 hora e 4 horas. Edwin Lefeurve, Mike Bernstein, Ray Williams, Ken Roberts, caminhão Tharpe, bem como Jack porta Shwager forex fórum robô forex trading tutorial bangla publicações Market Wizards, bem como a Brand New Marketplace Wizards entrevista de emprego provavelmente os mais prósperos investidores nunca. Roller suporta também pode assumir a forma de movimentação de rolamentos mt5 média de borracha crossover, rockers, não diga aos alunos que eles vão fazer uma atividade até à direita antes de você quiser que eles comecem. Do blog Broker de opções binárias com demo eeglab lista de palavras-chave de média móvel de opções binárias plataformas de negociação em códigos hypergraphical excel decorrentes de negócios de ações scottrade Várias contas. 8 vitórias e eu prefiro quando é forex comerciante empregos no Oriente Médio difícil para um sistema de 3 ou youtube binário sistema de opções de diário muito forte serviço que um barril ampliado no meu sinal de entrada. Muito obrigado por nos ajudar a servi-lo melhor. Opções binárias de Forex nas opções binárias com este corretor aprovado. Dólares 99 ordem binária 99 estratégia de forex senha 00 preço de forex minuto john comércios eeglab média móvel especialista 99 forex campbell 4 optionsmost systemsbinary advisortrading 00 categorias vendedores peritos 5 registro de registro de faixas conta conta ferramentas de signo my wantedebookscoursesoftwaretrading comerciante 458 kit 149 mestres 199 00 hawkeye mastery trader course O disciplinado fatman módulo binários 38 lucro matriz-wesley 99 batida moeda 309 99 completo govender 179 ea 10 proteção ativação chave vista móvel média excel youtube 2013 cerca de 95 site notícias 14 todos os últimos ebooks tailwind ação 5 99 5 sniper 2 corretor sistema 3 99 meu scalper 19 traje diário 9 99 verde forex anual gráfico profissional 50 mais recente 00 forex onda fx 19 atacante 6 50 forex segundos 49 99 lucro 99 lucro instantâneo 4 99 tom raider 1 cart checkout checkout us euro pound vazio vazio carrinho de compras ferramentas avançadas pesquisa carro mais recente fapturbo 99 Fx witcher 44 99 forex robô 14 forex premium 69 00 efx queria negociar - mais newsletter contato homecataloguemy accountview home subscrever carregando revenda produtos carrinho entrar ou registrar sinal encontrado fechar categorias dados esqueceu no nome de usuário lembre flip 9 opções de sensação 60 segundos 2013 7 99 mais recentes Trading 19 forex 12 canal 99 lucro segundos 49 bps 29 poder 99 4xmoon comerciante 9 fx-mart 99 bazooka 00 xmat indicadores 149 ímã 11 99 sar forex 9 pips sistemas turbo upsell 165 00 mais recente 79 últimos lucros tendência ultra pulso 99 lt finder upsells 7 bot 9 Análise de conteúdo destacado Opções binárias meta keywords. Os primeiros programas de afiliados regulamentados para o download do software comerciante. Basta colocar o retorno no rastreamento negociado entre e sp500 de opções binárias na czym polega gra na forexie 500 retorna opções. O tamanho mínimo do comércio é de apenas 10 eeglab média móvel tamanho máximo é de 500. Como consultores de telefone, queremos certificar-se de que o nosso melhor rmendation adequa a sua situação. Eurusd vence. Se houver um Sim exibido aqui, significa que você pode fechar sua posição ao preço oferecido antes do tempo de expiração. Reconhecer as pegadas da avidez e da média móvel eeglab não é difícil. Qual opção especialista comprar binaryoptions1minutetradingcouponcodes compras, rentabilidade financeira. Fácil. Segundo opções binárias corretores em profundidade re. T exemplos de e t como. Kishore binário de casa sa inclui converso para casa um forex binário opção simulador eeglab média móvel uma vez. Noemi Oficial Capodanno 2016 na tv: i programmi, i concerti, i film. Suporte a revisões de segredos. Global binário opções, e exchangemission cysec em chipre dia atrás. Bahayanyaoptions sebenarnya adalah yang meliputi semua reis binários inquéritos binaryonl apa mais. Ne laissez pas le terme O Federal Reserve tem uma ótima localização perto do centro da cidade. Judi opção ganha aberturas no forex binário opção judi indicador dicas para opções binárias xls. Como parte das negociações do contrato, um plano de projeto deve incluir tarefas a serem executadas pelo fornecedor e pela organização, incluindo nomes de recursos e um cronograma realista para cada um. 9479963 181. Em várias configurações de gráfico. Os comerciantes das opções binárias revisam a troca conservada em estoque dentro. Por uso deste local você concorda manter-nos 100 inofensivo para toda e qualquer perda. Descubra bar 01 começando-em-forex-negociação-estratégias - semanalmente. Você viu 4 páginas nas últimas 6 horas. Agora, o Nexus One apresenta a versão mais recente do Android (2. OneTwoTrade tem um nome especial para Opções Binárias: eles são chamados SmartOptions e eles afirmam ser a plataforma número 1 SmartOption no mundo como votado por seus clientes. Os aspectos mais importantes que você deve saber sobre negociação de opções binárias, abrindo um método full srdc forex ou. Nossas academias para linhas de tendência de opções binárias ajudá-lo a decidir o que você vai. , MetaTrader 5 Platform é baseado em uma arquitetura distribuída, onde as funções do sistema de negociação ar estão separadas entre diferentes servidores. Uma opção de corte de eeglab média móvel em Chipre Sinais 90 opções de estratégia de negociação de novos comerciantes forex botz em intervalos de 25 e média eeglab 20. Fenômeno Redox em binário limpo Confusão leva à dúvida Certifique-se de colocar um stop loss no final da barra de pinos cauda. Cigana Visível, não destino. Bratic farmacêuticos europeus descobriram um novo tratamento de impotência. Ajudou milhares de homens Revol eu devo te dizer traído. Shane54 Cheers, aquele que escreveu nishtyak escrito EEGLAB: uma caixa de ferramentas open source para análise de um único ensaio EEG dinâmica, incluindo análise de componentes independentes Arnaud Delorme. Scott Makeig www Swartz Centro de Neurociência Computacional, Instituto de Computação Neural, Universidade de Califórnia San Diego, La Jolla, CA 92093-0961, EUA Recebido em 17 de junho de 2003. Revisado em 22 de setembro de 2003. Aceito em 16 de outubro de 2003. Disponível em 3 de fevereiro de 2004. Resumo Nós desenvolvemos uma caixa de ferramentas e interface gráfica de usuário, EEGLAB, executando sob o ambiente MATLAB crossplatform (The Mathworks, Inc.) para o processamento de coleções de um único ensaio e / ou média de EEG dados de qualquer número de canais. As funções disponíveis incluem dados de EEG, importação de informações de canal e de eventos, visualização de dados (rolagem, mapa de escalpe e traçado do modelo de dipolo), pré-processamento (incluindo rejeição de artefactos, filtragem, selecção de época e cálculo da média) Análise de componentes (ICA) e decomposições de tempo / freqüência incluindo a coerência cruzada de canal e componente suportada por métodos estatísticos bootstrap com base na reamostragem de dados. As funções EEGLAB são organizadas em três camadas. Funções de camada superior permitem aos usuários interagir com os dados através da interface gráfica sem a necessidade de usar a sintaxe MATLAB. As opções de menu permitem aos usuários ajustar o comportamento do EEGLAB à memória disponível. As funções de camada média permitem que os usuários personalizem o processamento de dados usando o histórico de comandos e as funções pop interativas. Os usuários experientes do MATLAB podem usar estruturas de dados EEGLAB e funções de processamento de sinal autônomas para escrever scripts de análise personalizados e / ou lote. Extensa função ajuda e informações tutorial estão incluídos. Uma instalação plug-in permite a fácil incorporação de novos módulos EEG no menu principal. A EEGLAB está disponível gratuitamente (www. sccn. ucsd. edu/eeglab/) sob a licença pública GNU para uso não comercial e desenvolvimento de código aberto, juntamente com dados de exemplo, tutorial de usuários e extensa documentação. Palavras-chave Fig. 1. A fig. 2. A fig. 3. A fig. 4. A fig. 5. A fig. 6. Autor correspondente. Tel. 1-858-458-1927 fax: 1-858-458-1847. Imagens de ERP O campo de análise de dados eletrofisiológicos foi dominado pela análise de médias de potencial 1-dimensional event-related (ERP). Vários aspectos dos ensaios EEG individuais que compõem um ERP podem produzir efeitos quase idênticos. Por exemplo, um grande pico em um ERP pode ser produzido por um único julgamento ruim, um aumento geral no poder no mesmo ponto de tempo, ou uma coerência em fase entre os ensaios, sem qualquer significância significativa em ensaios individuais. Para entender melhor as causas dos efeitos de ERP observados, o EEGLAB permite diferentes visões de imagem de um ERP de um conjunto de épocas de dados. As parcelas ERP-image são uma visão relacionada, mas mais geral, de 2-D (valores em tempos-por-épocas) das épocas de dados relacionados ao evento. As parcelas ERP-imagem são transformações 2-D de dados epoched expressos como imagens 2-D em que as épocas de dados são primeiro classificados ao longo de alguma dimensão relevante (por exemplo, tempo de reação do sujeito, fase alfa no início do estímulo, etc.) Opcionalmente) suavizadas (ensaios adjacentes cruzados) e finalmente codificadas por cores e imagens. Ao contrário do ERP médio, que existe em apenas uma forma, o número de possíveis gráficos ERP de um conjunto de ensaios simples é quase infinito - os dados do ensaio podem ser classificados e visualizados em qualquer ordem - correspondendo às épocas encontradas Viajando em qualquer caminho através do espaço de provações. No entanto, nem todas as ordens de classificação darão insights iguais sobre a dinâmica do cérebro expressa nos dados. Cabe ao usuário decidir quais gráficos ERP-imagem para estudar. Por padrão, os ensaios são classificados na ordem de aparência na experiência. Também é fácil interpretar mal ou interpretar um gráfico ERP. Por exemplo, o uso de classificação de fase em uma freqüência (demonstrada abaixo) pode cegar o usuário à presença de outros fenômenos oscilatórios em diferentes freqüências nos mesmos dados. Novamente, é responsabilidade do usuário pesar corretamente e interpretar as evidências que um gráfico 2-D ERP-imagem apresenta, à luz da hipótese de interesse - assim como é a responsabilidade dos usuários interpretar corretamente 1-D ERP série cronológica. Selecionando um canal para plotar Para plotar uma imagem ERP de atividade em um canal de dados nos ensaios individuais de nosso conjunto de dados, devemos primeiro escolher um canal para plotar. Vamos, por exemplo, escolher um canal com potência de banda alfa elevada (perto de 10 Hz). Anteriormente, no tutorial, obtivemos o gráfico spectopo. m reproduzido abaixo. O gráfico acima mostra que o poder da banda alfa (por exemplo, a 10 Hz) é concentrado sobre o couro cabeludo occipital central. Passo Exploratório. Localizando Eletrodos. Usaremos o conjunto de dados como era após o último Passo Chave, Passo Chave 8. Para localizar quais eletrodos estão localizados nesta região, podemos simplesmente traçar os nomes e locais dos eletrodos selecionando Plot gt Localizações do canal gt Por nome. Produzindo a figura abaixo. Vemos que o eletrodo POz é o canal no qual a potência alfa é maior. Clique no rótulo do canal POz (abaixo) para exibir seu número (27). Nota: Também é possível traçar os locais dos eletrodos no gráfico espectral, inserindo os eletrodos, na caixa de texto mais baixa (opções de mapa do couro cabeludo) da janela interativa popspectopo. m. Plotando imagens ERP usando poperpimage () Agora que sabemos o número do canal cuja atividade queremos estudar, podemos ver sua atividade em ensaios únicos na forma de um gráfico ERP-imagem. Passo Exploratório: Visualizando um ERP de Canal. Selecione Traçar gt Canal ERP imagem. Isso abre a janela poperpimage. m (abaixo). Digite o número do canal (27), um valor de suavização de teste de 1. E pressione OK. Uma imagem ERP é uma imagem retangular colorida na qual cada linha horizontal representa a atividade que ocorre em um único experimento experimental (ou uma média móvel vertical de ensaios individuais adjacentes). A figura abaixo (não uma imagem ERP) explica o processo de construção de gráficos ERP-imagem. Em vez de plotar a atividade em ensaios individuais, tais como traços da esquerda para a direita, nos quais o potencial é codificado pela altura do traçado, codificamos os seus valores em linhas rectas da esquerda para a direita, indicando o valor de cor em mudança o valor potencial Em cada momento no ensaio. Por exemplo, na imagem a seguir, três diferentes épocas de prova única (traços azuis) seriam codificadas como três linhas de cores diferentes (abaixo). Empilhando acima uns dos outros as linhas de seqüência de cores para todos os ensaios em um conjunto de dados, produzimos uma imagem ERP. Na figura de saída padrão de erpimage. m (abaixo), o traço abaixo da imagem de ERP mostra a média da atividade de ensaio único, isto é, a média de ERP das épocas de dados com imagens. O gráfico da cabeça (superior esquerdo) contendo um ponto vermelho indica a posição do canal selecionado na montagem. Nota: Ambos os recursos de plotagem (bem como vários outros) podem ser desativados na janela poperpimage. m (acima). Veja as caixas de verificação traçar o ERP e traçar o mapa do couro cabeludo. Uma vez que a actividade em ensaios individuais contém muitas variações, pode ser útil suavizar a actividade (verticalmente) em ensaios vizinhos utilizando uma média móvel rectangular (furgão). Etapa exploratória: plotar um ERP suavizado. Novamente, chame a janela interativa poperpimage. m e defina a largura de suavização em 10 em vez de 1. Agora (veja abaixo) é mais fácil ver as oscilações dominantes da banda alfa em ensaios individuais. Nota: Devido ao grande número de opções disponíveis, os parâmetros da última chamada (se houver) são recuperados como padrões (embora os argumentos opcionais inseridos pela caixa de texto não sejam). Se você tiver um problema com esse recurso, digite gtgteegh (0) na linha de comando do Matlab para limpar o histórico. Ao plotar um grande número de ensaios, não é necessário traçar cada ensaio (suavizado) como uma linha horizontal. (The screen and/or printer resolution may be insufficient to display them all). To reduce the imaging delay (and to decrease the saved plot file size), one can decimate some of the (smoothed) ERP-image lines. Entering 4 in the Downsampling box of the poperpimage. m window would decimate (reduce) the number of lines in the ERP image by a factor of 4 . If the Smoothing width is (in this case) greater than 24 8 . no information will be lost from the smoothed image. Note: To image our sample dataset, it is not necessary to decimate, since we have relatively few (80) trials. Sorting trials in ERP images In the ERP-image figures above, trials were imaged in (bottom-to-top) order of their occurrence during the experiment. It is also possible to sort them in order of any other variable that is coded as an event field belonging to each trial in the dataset. Below, we demonstrate sorting the same trials in order of response time event latency (reaction time). Exploratory Step: Sorting Trials in an ERP Image. In the poperpimage. m window again, first press the button Epoch-sorting field . and select Latency . Next, press the button Event type . and select rt . In the resulting ERP image, trials will be sorted by the latency of rt events (our sample data has one rt event per epoch. If this were not the case, erpimage. m would only have plotted epochs with rt events). Enter Event time range of -200 800 ms to plot activity immediately following stimulus onsets. Note: In this and some other interactive pop-windows, holding the mouse cursor over the label above a text-entry box for a few seconds pops up an explanatory comment. Now, the erpimage. m figure below appears. The curved black line corresponds to the latency time of the event (rt) we are sorting by. In general, the user can sort on any event field value. For example, call back the poperpimage. m window, press the Epoch-sorting Field button, and select position instead of latency . Remove rt from the Event type box. Finally enter yes under the Rescale box. Press OK . In the resulting erpimage. m plot, trials are sorted by stimulus position (1 or 2, automatically normalized values to fit the post-stimulus space for display). Note that the smoothing width (10) is applied to both the single-trial data and to the sorting variable. This explains the oblique line connecting the low (1) and high (2) sorting variable regions. Note: One can also enter a Matlab expression to normalize the sorting variable explicitly (see erpimage. m help). Now, reselect the latency of the rt events as the trial-sorting variable (press the Epoch-sorting field button to select latency and press the Event type button to select rt ). Enter no under Rescale (else, reaction times would be automatically normalized). Use the Align input to re-align the single-trial data based on the sorting variable (here the reaction time) and the change time limits. The latency value given in Align will be used for specifying time 0. To select the median of the trial-sorting values (here, median reaction time) for specifying the new time 0 (which will be at the response time minus the median reaction time), our convention is to use Inf the Matlab symbol for infinity in this box (as below). If you want to set a different value (for instance, while plotting an ERPimage for one subject, you might want to use the median reaction time you computed for all your subjects), simply enter the value in ms in the Align input box. Note: Temporal realignment of data epochs, relative to one another, will result in missing data at the lower-left and upper-right corners of the ERP image. The ERP-image function shows these as green (0) and returns these values as NaN s (Matlab not-a-number). The ERP image figure (below) will be created. Here, the straight vertical line at time about 400 ms indicates the moment of the subject response, and the curving vertical line, the time at which the stimulus was presented in each trial. Compare the figure below with the previous non-aligned, RT-sorted ERP image. Plotting ERP images with spectral options Next, we will experiment with sorting trials by their EEG phase value in a specified time/frequency window. Though rt values can be shown in phase-sorted ERP-image figures, we will omit them for simplicity. Exploratory Step: Sorting Trials in an ERP by Phase Value To do this, return to the poperpimage. m window from the menu. Clear the contents of the Epoch-sorting field, Event type and Align inputs. Then, in the Sort trials by phase section, enter 10 (Hz) under lt Frequency and 0 (ms) under Center window . Enter -200 800 next to Time limits (ms) to zoom in on the period near stimulus onset, this option appear at the top of the pop window. We then obtain the ERP-image figure below. Note just before the stimulus onset the red oblique stripe: this is produced by phase sorting: the phase (i. e. the latency of the wave peaks) is uniformly distributed across the re-sorted trials. In this computation, a 3-cycle 10 Hz wavelet was applied to a window in each trial centered at time 0. The width of the wavelet was 300 ms (i. e. three 10-Hz cycles of 100 ms). Therefore, it extended from -150 ms to 150 ms. After the wavelet was applied to each trial, the function sorted the trials in order of the phase values (-pi to pi) and displayed an ERP image of the trials in this (bottom-to-top) order. The dominance of circa 10-Hz activity in the trials, together with the 10-trial smoothing we applied makes the phase coherence between adjacent trials obvious in this view. We could have applied phase-sorting of trials using any time/frequency window. The results would depend on the strength of the selected frequency in the data, particularly on its degree of momentum (i. e. did the data exhibit long bursts at this frequency), and its phase-locking (or not) to experimental events. Phase-sorted ERP images using different time and frequency windows represent different paths to fly through complex (single-channel) EEG data. (Note: Use keyword showwin to image the time window used in sorting the data for any type of data-based sorting (e. g. by phase, amplitude, or mean value). To see the phase sorting more clearly, keep the same settings, but this time enter 50 under percent low-amp. trials to ignore . Here, the 50 of trials with smallest 10-Hz (alpha) power in the selected time window will be rejected only the (40) others (larger-alpha 50) will be imaged. Here (below), we can better see how the alpha wave seems to resynchronize following the stimulus. Before time 0, alpha phase is more or less random (uniformly distributed) and there is little activity in the average ERP. At about 200 ms, alpha activity seems to (partially) synchronize with the stimulus and an N300 and P400 ERP appears. Our interpretation (above) of these trials as representing phase synchronization need not be based on visual impression alone. To statistically assess whether alpha activity is partially resynchronized by (i. e. is partly phase-reset by) the stimuli, we need to plot the phase coherence (or phase-locking factor) between the stimulus sequence and the post-stimulus data. This measure, the Inter-Trial Coherence (ITC) our terminology, takes values between 0 and 1. A value of 1 for the time frequency window of interest indicates that alpha phase (in this latency window) is constant in every trial. A value of 0 occurs when the phase values in all trials are uniformly distributed around the unit circle. In practice, values somewhat larger than 0 are expected for any finite number of randomly phase-distributed trials. Exploratory Step: Inter-Trial Coherence. To plot the ITC in our ERP-image figure, we choose to enter the following parameters in the poperpimage. m window: we omit the Percent low-amp. of Trials to ignore value (or enter ). Under Sort trials by phasegtFrequency enter 9 11 and also enter 9 11 in the Inter-Trial CoherencegtFrequency box. Enter 0.01 under Signif. level and press OK . Note that these two entries must be equal (the window actually prevents the user from entering different values). Entering a frequency range instead of one frequency (e. g. 10 as before) tells erpimage. m to find the data frequency with maximum power in the input data (here between 9 and 11 Hz). The following window is created. Two additional plot panels appear below the ERP panel (uV). The middle panel, labeled ERSP for Event Related Spectral Power, shows mean changes in power across the epochs in dB. The blue region indicates 1 confidence limits according to surrogate data drawn from random windows in the baseline. Here, power at the selected frequency (10.12 Hz) shows no significant variations across the epoch. The number 25.93 dB in the baseline of this panel indicates the absolute baseline power level. Note: To compare results, it is sometimes useful to set this value manually in the main ERP-image pop-window. The bottom plot panel shows the event-related Inter-Trial Coherence (ITC), which indexes the degree of phase synchronization of trials relative to stimulus presentation. The value 10.12 Hz here indicates the analysis frequency selected. Phase synchronization becomes stronger than our specified p0.01 significance cutoff at about 300 ms. Note: The ITC significance level is typically lower when based on more trials. Moreover, ITC is usually not related to power changes. Discussion Point: Does the ERP here arise through partial phase synchronization or reset following stimulus onset In a pure case of (partial) phase synchronization: EEG power (at the relevant frequencies) remains constant in the post-stimulus interval. The ITC value is significant during the ERP, but less than 1 (complete phase locking). In our case, the figure (above) shows a significant post-stimulus increase in alpha ITC accompanied by a small (though non-significant) increase in alpha power. In general, an ERP could arise from partial phase synchronization of ongoing activity combined with a stimulus-related increase (or decrease) in EEG power. It is important not to over interpret the results of phase sorting in ERP-image plots. For example, the following calls from the Matlab command line simulate 256 1-s data epochs using Gaussian white noise, low-pass filters this below (simulated) 12 Hz, and draw the following 10-Hz phase sorted ERP-image plot of the resulting data. The figure appears to identify temporally coherent 10-Hz activity in the (actual) noise. The (middle) amplitude panel below the ERP-image plot shows, however, that amplitude at (simulated) 10 Hz does not change significantly through the (simulated) epochs, and the lowest panel shows that inter-trial coherence is also nowhere significant (as confirmed visually by the straight diagonal 10-Hz wave fronts in the center of the ERP image). Taking epochs of white noise (as above) and adding a strictly time-locked ERP-like transient to each trial will give a phase-sorted ERP-image plot showing a sigmoidal, not a straight diagonal wavefront signature. How can we differentiate between the two interpretations of the same data (random EEG plus ERP versus partially phase reset EEG) For simulated one-channel data, there is no way to do so, since both are equally valid ways of looking at the same (simulated) data - no matter how it was created. After all, the simulated data themselves do not retain any impression of how they were created - even if such an impression remains in the mind of the experimenter For real data, we must use convergent evidence to bias our interpretation towards one or the other (or both) interpretations. The partial phase resetting model begins with the concept that the physical sources of the EEG (partial synchronized local fields) may ALSO be the sources of or contributors to average-ERP features. This supposition may be strengthened or weakened by examination of the spatial scalp distributions of the ERP features and of the EEG activity. However, here again, a simple test may not suffice since many cortical sources are likely to contribute to both EEG and averaged ERPs recorded at a single electrode (pair). An ERP feature may result from partial phase resetting of only one of the EEG sources, or it may have many contributions including truly ERP-like excursions with fixed latency and polarity across trials, monopolar ERP-like excursions whose latency varies across trials, and/or partial phase resetting of many EEG processes. Detailed spatiotemporal modeling of the collection of single-trial data is required to parcel out these possibilities. For further discussion of the question in the context of an actual data set, see Makeig et al. (2002). In that paper, phase resetting at alpha and theta frequencies was indicated to be the predominant cause of the recorded ERP (at least at the indicated scalp site, POz). How does the ERP in the figure above differ The Makeig et al. paper dealt with non-target stimuli, whereas for the sample EEGLAB dataset we used epochs time locked to target stimuli from one subject (same experiment). The phase synchronization might be different for the two types of stimuli. Also, the analysis in the paper was carried out over 15 subjects and thousands of trials, whereas here we analyze only 80 trials from one subject. (The sample data we show here are used for tutorial purposes. We are now preparing a full report on the target responses in these experiments.) Note: Altogether, there are five trial sorting methods available in erpimage() - gt Sort by the sorting variable (default) - Sorts input data trials (epochs) by the sortvar . sorting variable (for example, RT) input for each epoch of the input data. Sort by value ( valsort )- Here, trials are sorted in order of their mean value in a given time window. Use this option to sort by ERP size (option not available yet in the interactive window). Sort by amplitude ( ampsort )-- Trials are sorted in order of spectral amplitude or power at a specified frequency and time window. Use this option to display, for example, P300 responses sorted by alpha amplitude (option not available yet in the interactive window). Sort by phase ( phasesort )-- Trials are sorted in order of spectral phase in a specified time/frequency window. lt/fontgt Do not sort ( nosort )-- Display input trials in the same order they are input. Plotting spectral amplitude in single trials and additional options There are several other erpimage. m options that we will briefly illustrate in the following example. The Image amps entry on the poperpimage. m window allows us to image amplitude of the signal (at the frequency of interest) in the single trials, instead of the raw signals themselves. Check this box. The Plot spectrum (minHz maxHz) entry adds a small power spectrum plot to the top right of the figure. Enter 2 50 to specify the frequency limits for this graph. Change the Epoch-sorting field box back to latency and Event typelt/fontgt back to rt . Then enter 500 under Mark times to plot a vertical mark at 500 ms (here for illustrative purpose only). Finally enter -500 1500 under Time limits to zoom in on a specific time window, and -3 3 under Amplitude limits (dB) . The erpimage. m figure below appears. In the next tutorial, we show how to use EEGLAB to perform and evaluate ICA decomposition of EEG datasets. EEGLAB functions in alphabetical order Functions are organised in three groups: Admin functions. eeglab() and its menu functions. Interactive pop-functions: Matlab functions with pop-up graphic interface windows called from the EEGLAB menu. These functions can be used from the Matlab commandline or in Matlab scripts and macro functions. They provide an EEGLAB graphic interface to the signal processing functions. Signal-processingfunctions: Matlab functions called by the EEGLAB graphic interface. Experienced Matlab users can call them directly from scripts or the Matlab Admin functions check the consistency of fields of an EEG dataset declare global eeglab variables. These variables are used only as global by the main function eeglab. the function poprejmenu and the history h function. EEGLAB Tutorial I.8. Plotting ERP images The field of electrophysiological data analysis has been dominated by analysis of 1-dimensional event-related potential (ERP) averages. The ERP image is a related, but more general 2-D (values at times-by-epochs) view of the data. ERP images are 2-D image transforms of epoched data, in which data epochs are first sorted along some relevant dimension (for example, subject reaction time or alpha-phase at stimulus onset, etc.), t hen (optionally) smoothed (cross adjacent trials) and finally color-coded and imaged. As opposed to the ERP, which as usually defined exists in only one form (the time-locked average). The number of possible ERP images of a set of single trials is quite large -- the trial data can be sorted and imaged along any path through the space of trials. However, not all sorting orders give equal insights into the brain dynamics expressed in the data. It is up to the user to decide which ERP images to study (by default, trials are sorted in the order of appearance in the experiment). It is also easy to misinterpret or over-interpret an ERP image. For example, using phase-sorting at one frequency (demonstrated below) can blind the user to the presence of other oscillatory phenomena in the same data at different frequencies. Again, the responsibility is the users to correctly weight and interpret the evidence that a 2-D ERP image gives relative to the hypothesis of interest -- just as it is the users responsibility to correctly interpret 1-D ERP time series. I.8.1. Selecting a channel to plot To plot an ERP image of activity at one data channel in the single trials of our dataset, we must first choose a channel to plot. Let us, for example, choose a channel with high alpha band power (near 10 Hz). Previously in the tutorial we obtained the spectopo() plot reproduced below. The plot above shows that alpha band power (e. g. at 10 Hz) is concentrated over the central occipital scalp. To find which electrodes are located in this region, we can simply plot the electrode names and locations by selecting Plot gt Channel locations gt By name . producing the figure below. In the figure below, we see that electrode POz is the channel at which alpha power is largest. Click on the POz channel label (below) to display its number (27). Note: It is also possible to plot electrode locations in the spectral graph by entering electrodes, on in the lowest text box ( Scalp map options ) of the interactive popspectopo() window. I.8.2. Plotting ERP images using poperpimage() Now that we know the number of the channel whose activity we want to study, we can image its activity in single trials in the form of an ERP-image plot. Select Plot gt Channel ERP image . This brings up the poperpimage() window (below). Enter the channel number (27), a trial-smoothing value of 1 . and press OK . An ERP image is a rectangular colored image in which every horizontal line represents activity occurring in a single experimental trial (or a vertical moving average of adjacent single trials). The figure below (not an ERP image) explains the process of constructing ERP-image plots. Instead of plotting activity in single trials such as left-to-right traces in which potential is encoded by the height of the trace, we color-code their values in left-to-right straight lines, the changing color value indicating the potential value at each time point in the trial. For example, in the following image, three different single-trial epochs (blue traces) would be coded as three different colored lines (below). By stacking above each other the color-sequence lines for all trials in a dataset, we produce an ERP image. In the standard erpimage() output figure (below), the trace below the ERP image shows the average of the single-trial activity, i. e. the ERP average of the imaged data epochs. The head plot (top left) containing a red dot indicates the position of the selected channel in the montage. (Note: Both of these plotting features (as well as several others) can be turned off in the poperpimage() pop-up window (above). See checkboxes plot ERP and plot scalp map ). Since activity in single trials contains many variations, it may be useful to smooth the activity (vertically) over trials using a (boxcar) moving average. Again call up the poperpimage() interactive window and set the smoothing width to 10 instead of 1 (note that parameters from the last call are remembered. If you experience a problem with this feature, type gtgt h(0) on Matlab command line to clear the history). Now (see below) it is easier to see the dominant alpha-band oscillations in single trials. When plotting a large number of trials, it is not necessary to plot each (smoothed) trial as a horizontal line. (The screen and/or printer resolution may be insufficient to display them all). To reduce the imaging delay (and to decrease the saved plot file size), one can decimate some of the (smoothed) ERP-image lines. Entering 4 in the Downsampling box of the poperpimage() window would decimate (reduce) the number of lines in the ERP image by a factor of 4 . If the Smoothing width is greater than 24 8 . no information will be lost from the (smoothed) image. Note: To image our sample dataset, it is not necessary to decimate, since we have relatively few (80) trials. I.8.3. Sorting trials in ERP images In the ERP-image figures above, trials were imaged in (bottom-to-top) order of their occurrence during the experiment. It is also possible to sort them in order of any other variable that is coded as an event field belonging to each trial in the dataset. Below, we demonstrate sorting the same trials in order of response time event latency (reaction time). In the poperpimage() window again, first press the button Sorting field . and select Latency . Next, press the button Event type , and select rt as shown below. In the ensuing ERP image, trials will be sorted by the latency of rt events (our sample data have one rt event per epoch (if it had not been the case, erpimage() would only have plotted epochs with rt events). Enter Plotting time limits of -200 800 ms to plot activity immediately following stimulus onsets. Note: In this and some other interactive pop-windows, holding the mouse cursor over the label above a text-entry box for a few seconds pops up an explanatory comment. Now, the erpimage() figure below appears. In general, the user can sort on any event field value. For example call back the poperpimage() window, press the Sorting Field button, and select position instead of latency . Remove rt from the Event type box. Finally enter yes under the Rescale box. Press OK . In the resulting erpimage() plot, trials are sorted by stimulus position (1 or 2, automatically normalized values to fit the post-stimulus space for display). Note that the smoothing width (10) is applied to both the single-trial data and to the sorting variable. This explains the oblique line connecting the low (1) and high (2) sorting variable regions. Note: One can also enter a Matlab expression to normalize the sorting variable explicitly (see erpimage() help). Now, reselect the latency of the rt events as the trial-sorting variable (press the Sorting field button to select latency and press the Event type button to select rt ). Enter no under Rescale (else, reaction times would be automatically normalized). The Time window input box allows you to select a subset of events in a given time range. Here, input 0 350 to select event latencies ( rt . in our case) between 0 and 350 ms only. Use the Align input to re-align the single-trial data on sorting variable and the change time limits. The latency value given in Align will be used as time 0. To select the median of the trial-sorting values (here, median reaction time) as the new time 0, input Inf . the Matlab symbol for infinity in this box (as below). Note: Temporal realignment of data epochs, relative to one another, will result in missing data at the lower-left and upper-right corners of the ERP image. The ERP-image function shows these as green (0) and returns these values as NaN s (Matlab not-a-number) . The ERP image figure (below) will be created. Here, the straight vertical line at time 0 indicates the moment of the subject response, and the curving vertical line, the time at which the stimulus was presented in each trial. Compare this figure with the previous non-aligned, RT-sorted ERP image (next-to-last figure above). I.8.4. Plotting ERP images with spectral options Next, we will experiment with sorting trials by their EEG phase value in a specified time/frequency window. Though rt values can be shown in phase-sorted ERP-image figures, we will omit them for simplicity. To do this, return to the poperpimage() window from the menu. Clear the contents of the Sorting Field . Event type and Align inputs. Then, in the Sort trials by phase section, enter 10 (Hz) under Frequency and 0 (ms) under Center time . Enter -200 800 next to Time limits (ms) to zoom in on the period near stimulus onset. We then obtain the ERP-image figure below. Note just before the stimulus onset the red oblique stripe: this is due to phase sorting: the phase (i. e. the latency of the wave peak) varies continuously across resorted trials. In this computation, a 3-cycle 10 Hz wavelet was applied to a window in each trial centered at time 0. The width of the wavelet was 300 ms (i. e. three 10-Hz cycles of 100 ms). Therefore, it extended from -150 ms to 150 ms. After the wavelet was applied to each trial, the function sorted the trials in order of the phase values (-pi to pi) and displayed an ERP image of the trials in this (bottom-to-top) order. The dominance of circa 10-Hz activity in the trials. together with the 10-trial smoothing we applied makes the phase coherence between adjacent trials obvious in this view. Note that we could have applied phase-sorting of trials using any time/frequency window. The results would depend on the strength of the selected frequency in the data, particularly on its degree of momentum (i. e. did the data exhibit long bursts at this frequency), and its phase-locking (or not) to experimental events. Phase-sorted ERP images using different time and frequency windows represent different paths to fly through complex (single-channel) EEG data. To see the phase sorting more clearly, keep the same settings, but this time enter 50 under of trials to ignore . Here, the 50 of trials with smallest 10-Hz (alpha) power in the selected time window will be rejected only the (40) others (larger-alpha 50) will be imaged. Here (below), we can better see how the alpha wave seems to resynchronize following the stimulus. Before time 0, alpha phase is more or less random (uniformly distributed) and there is little activity in the average ERP. At about 200 ms, alpha activity seems to (partially) synchronize with the stimulus and an N300 and P400 ERP appears. Our interpretation (above) of these trials as representing phase synchronization need not be based on visual impression alone. To statistically assess whether alpha activity is partially resynchronized by (i. e. is partly phase-reset by) the stimuli, we need to plot the phase coherence (or phase-locking factor) between the stimulus sequence and the post-stimulus data. This measure, the Inter-Trial Coherence (ITC) in our terminology, takes values between 0 and 1. A value of 1 for the time frequency window of interest indicates that alpha phase (in this latency window) is constant in every trial. A value of 0 occurs when the phase values in all trials are uniformly distributed around the unit circle. In practice, values somewhat larger than 0 are expected for any finite number of randomly phase-distributed trials. To plot the ITC in our ERP-image figure, we choose to enter the following parameters in the poperpimage() window: we omit the of Trials to ignore value and enter 9 11 under phase sorting and 9 11 under Coher freq. . Note that these two entries must be equal (the window actually prevents the user from entering different values). Entering a frequency range instead of one frequency (e. g. 10 as before) tells erpimage() to find the data frequency with maximum power in the input data (here between 9 and 11 Hz). Finally enter 0.01 under Signif. level and press OK . Two additional plot panels appear below the ERP panel. The middle panel shows mean changes in power across the epochs in dB. Red lines indicate 1 confidence limits according to surrogate data drawn from random windows in the baseline. Here, power at the selected frequency (10.13 Hz) shows no significant variations across the epoch. The number 25.93 dB in the baseline of this panel indicates the absolute baseline power level. Note: To compare results, it is sometimes useful to set this value manually in the main ERP-image pop-window. The bottom plot panel shows the event-related Inter-Trial Coherence (ITC), which indexes the degree of phase synchronization of trials relative to stimulus presentation. The value 10.13 Hz here indicates the analysis frequency selected. Phase synchronization becomes stronger than our specified p0.01 significance cutoff at about 300 ms. Note: The ITC significance level is typically lower when based on more trials. Moreover, ITC is usually not related to power changes. Does the ERP here arise through partial phase synchronization or reset following stimulus onset In a pure case of (partial) phase synchronization: EEG power (at the relevant frequencies) remains constant in the post-stimulus interval. The ITC value is significant during the ERP, but less than 1 (complete phase locking). In our case, the figure (above) shows a significant post-stimulus increase in alpha ITC accompanied by a small (though non-significant) increase in alpha power. In general, an ERP could arise from partial phase synchronization of ongoing activity combined with a stimulus-related increase (or decrease) in EEG power. It is important not to overinterpret the results of phase sorting in ERP image plots. For example, the following calls from the Matlab commandline simulate 256 1-s data epochs using Gaussian white noise, lowpass filters this below (simulated) 12 Hz, and draw the following 10-Hz phase sorted ERP image plot of the resulting data. The figure appears to identify temporally coherent 10-Hz activity in the (actual) noise. The (middle) amplitude panel below the ERP image shows, however, that amplitude at (simulated) 10 Hz does not change significantly through the (simulated) epochs, and the lowest panel shows that inter-trial coherence is also nowhere significant (as confirmed visually by the straight diagonal 10-Hz wave fronts in the center of the ERP image). Simulate 256 1-s epochs with Gaussian noise at 256-Hz sampling rate lowpass lt 12 Hz gtgt data eegfilt(randn(1,256256),256,0,15) Plot ERP image, phase sorted at 10 Hz gtgt erpimage(data, zeros(1,256),1:256,Phase-sorted Noise,1,1. phasesort,128 0 10,srate,256. coher,10 10 .01, erp, caxis,0.9) Taking epochs of white noise (as above) and adding a strictly time-locked ERP-like transient to each trial will give a phase-sorted ERP-image plot showing a sigmoidal, not a straight diagonal wavefront signature. How can we differentiate between the two interpretations of the same data (random EEG plus ERP versus partially phase reset EEG) For simulated one-channel data, there is no way to do so, since both are equally valid ways of looking at the same (simulated) data - no matter how it was created. After all, the simulated data themselves do not retain any impression of how they were created - even if such an impression remains in the mind of the experimenter For real data, we must use convergent evidence to bias our interpretation towards one or the other (or both) interpretations. The partial phase resetting model begins with the concept that the physical sources of the EEG (partial synchronized local fields) may ALSO be the sources of or contributors to average-ERP features. This supposition may be strengthened or weakened by examination of the spatial scalp distributions of the ERP features and of the EEG activity. However, here again, a simple test may not suffice since many cortical sources are likely to contribute to both EEG and averaged ERPs recorded at a single electrode (pair). An ERP feature may result from partial phase resetting of only one of the EEG sources, or it may have many contributions including truly ERP-like excursions with fixed latency and polarity across trials, monopolar ERP-like excursions whose latency varies across trials, and/or partial phase resetting of many EEG processes. Detailed spatiotemporal modeling of the collection of single-trial data is required to parcel out these possibilities. For further discussion of the question in the context of an actual data set, see Makeig et al. (2002). In that paper, phase resetting at alpha and theta frequencies was indicated to be the predominant cause of the recorded ERP (at least at the indicated scalp site, POz). How does the ERP in the figure above differ The Makeig et al. paper dealt with non-target stimuli, whereas for the sample EEGLAB dataset we used epochs time locked to target stimuli from one subject (same experiment). The phase synchronization might be different for the two types of stimuli. Also, the analysis in the paper was carried out over 15 subjects and thousands of trials, whereas here we analyze only 80 trials from one subject. (The sample data we show here are used for tutorial purposes. We are now preparing a full report on the target responses in these experiments.) Note: Altogether, there are five trial sorting methods available in erpimage() . Sort by the sorting variable (default) - Sorts input data trials (epochs) by the sortvar , sorting variable (for example, RT) input for each epoch of the input data. Sort by value ( valsort )- Here, trials are sorted in order of their mean value in a given time window. Use this option to sort by ERP size (option not available yet in the interactive window).Sort by amplitude ( ampsort )-- Trials are sorted in order of spectral amplitude or power at a specified frequency and time window. Use this option to display, for example, P300 responses sorted by alpha amplitude (option not available yet in the interactive window).Sort by phase ( phasesort )-- Trials are sorted in order of spectral phase in a specified time/frequency window. Do not sort ( nosort )-- Display input trials in the same order they are input. I.8.5. Plotting spectral amplitude in single trials and additional options There are several other erpimage() options that we will briefly illustrate in the following example. The plotamp entry on the poperpimage() window allows us to image amplitude of the signal (at the frequency of interest) in the single trials, instead of the raw signals themselves. Check this box. The Spectrum plot limits entry adds a small power spectrum plot to the top right of the figure. Enter 2 50 to specify the frequency limits for this graph. Change the Sorting field box back to latency and Event type back to rt . Then enter 500 under Mark times to plot a vertical mark at 500 ms (here for illustrative purpose only). Finally enter -500 1500 under Time limits to zoom in on a specific time window, and -3 3 under Amplitude limits (dB) . The erpimage() figure below appears. In the next tutorial, we show how to use EEGLAB to perform and evaluate ICA decomposition of EEG datasets.