Back to the main page.

Bug 11 - ICA gives complex results

Status CLOSED FIXED
Reported 2010-01-08 13:53:00 +0100
Modified 2010-03-22 13:48:03 +0100
Product: FieldTrip
Component: preproc
Version: unspecified
Hardware: PC
Operating System: Windows
Importance: P5 normal
Assigned to: Saskia Haegens
URL:
Tags:
Depends on:
Blocks:
See also:

Anne van Hoogmoed - 2010-01-08 13:53:37 +0100

If I perform ICA on my data, the result contains complex numbers. I've attached the data on which I performed the ICA and the result of the ICA. The script I used is cfg = []; ICA_64 = ft_componentanalysis (cfg,data);


Saskia Haegens - 2010-01-12 14:43:04 +0100

*** Bug 10 has been marked as a duplicate of this bug. ***


Saskia Haegens - 2010-01-12 14:46:11 +0100

*** Bug 12 has been marked as a duplicate of this bug. ***


Saskia Haegens - 2010-01-12 14:50:45 +0100

Hoi Anne, Ik ga kijken of ik je probleem kan oplossen. Zou je me daartoe wat data kunnen sturen waarmee ik het probleem kan reconstrueren? Groet, Saskia


Saskia Haegens - 2010-01-19 13:36:47 +0100

Hoi Anne, Volgens mij zit het probleem erin dat je de reference channels in je data hebt zitten. Dit maakt de data ongeschikt voor ICA, en, afhankelijk van de startcondities die variƫren per sessie, kom je in een soort loop terecht die blijkbaar (soms) leidt tot complexe output. Dit is duidelijk niet de bedoeling! Zou je dit eens willen proberen: cfg=[]; cfg.channel={'all', '-M1', '-M2'}; test = ft_componentanalysis(cfg, data) Op de dataset die je mij stuurde geeft het resultaat dat er normaal uitziet. Laat maar weten of het hiermee is opgelost, en of je nu interpreteerbare output krijgt. Groet, Saskia ----- From: Anne van Hoogmoed [mailto:Anne.vanHoogmoed@donders.ru.nl] Sent: dinsdag 12 januari 2010 16:02 To: saskia.haegens@donders.ru.nl Subject: ICA Hoi Saskia, Ik kan je de .mat-file niet opsturen, omdat Outlook het als een snelkoppeling ziet, maar ik heb het op common/temporary gezet. In het bestand ICA_mislukt.mat staat de data die ik gebruikt heb voor de ICA en het resultaat dat ik kreeg (ICA_64). Voor de zekerheid heb ik er ook nog even mijn hele file van inlezen tot ICA bij gedaan (en de data die ik daarbij gebruikt heb) voor als het probleem al eerder zou zitten. Fijn dat je er naar wilt kijken! Groeten Anne


Saskia Haegens - 2010-01-19 15:10:48 +0100

Hoi Anne, Volgens mij is dit de beste oplossing: cfg=[]; cfg.runica.pca = 61; test = ft_componentanalysis(cfg, data) Als het goed is krijg je op deze manier 61 componenten (aangezien M1=-M2, is dat ook het max aantal bij jouw 62 kanalen, dus een volledige beschrijving van je data). Op deze manier kun je het wel nog terug projecteren naar je oorspronkelijke 62 kanalen, en dus ook normaal plotten. Ik hoor graag of het is gelukt! Groet, Saskia