Separazione alla fonte cieca della miscela convessa?


18

Supponiamo di avere fonti indipendenti, e osservare miscele convesse: X 1 , X 2 , . . . , X n m Y 1nX1,X2,...,Xnm

Y1=a11X1+a12X2++a1nXn...Ym=am1X1+am2X2++amnXn

con per tutti e per tutti .jaij=1iaij0i,j

Qual è lo stato dell'arte nel recuperare da ?XY

PCA è fuori discussione perché ho bisogno che i componenti siano identificabili. Ho esaminato ICA e NMF: non riesco a trovare alcun modo per imporre la non negatività dei coefficienti di miscelazione per ICA e NMF non sembra massimizzare l'indipendenza.


2
Penserei che questo dovrebbe essere chiamato "analisi di componenti indipendenti non negativi", ma sembra che questo nome sia stato usato per ICA con il vincolo di non negatività sulle fonti , non sulla matrice di miscelazione ( eecs.qmul.ac.uk/ ~ markp / 2003 / Plumbley03-algoritmi-c.pdf ). Quindi questo non si applica al tuo caso. Domanda interessante. XUN
ameba dice di reintegrare Monica il

Non vuoi che le somme investano j invece di me? Puoi presumere che le fonti siano approssimativamente gaussiane? se sono unimodali e presentano un decadimento sufficientemente rapido, è possibile che sia sufficiente il montaggio di un GMM.
Yair Daon,

@YairDaon Ah sì grazie, ottima cattura. Sfortunatamente le fonti sono discrete e non sembrano nemmeno miscele di gaussiani. Ma forse potrei approssimarli grossolanamente come miscele gaussiane e quindi perfezionare ulteriormente. Ma sarebbe bello avere qualcosa di più generale / robusto
anonimo il

1
Quali algoritmi ICA hai provato? Sono un po 'arrugginito, ma penso che l'ipotesi di non negatività dei coefficienti di miscelazione possa essere imposta in alcuni algoritmi che assumono determinati modelli per i segnali come l'algoritmo WASOBI (Second Order-Adjusted Second Order Blind Identification, WASOBI), come si presume che tu possa modellare i segnali come processi AR e, quindi, è possibile imporre condizioni nei coefficienti.
Néstor,

1
Le fonti sono tutte supportate sul set {1,2, ..., 96}
anonimo

Risposte:


0

Si potrebbe ottenere usando una non linearità esponenziale invece della tipica / default tanh (), se anche X è non negativo.

Formula 40 in https://www.cs.helsinki.fi/u/ahyvarin/papers/NN00new.pdf e disponibile nella maggior parte delle implementazioni.

Ad esempio in sklearn usa semplicemente fun = 'exp' https://scikit-learn.org/stable/modules/generated/sklearn.decomposition.FastICA.html


1
Benvenuto in Stats.SE. Puoi modificare la tua risposta ed espanderla per spiegare i passaggi chiave dei link che fornisci? In questo modo, le informazioni sono ricercabili qui (e talvolta i collegamenti si interrompono). Potresti dare un'occhiata ad un aiuto di formattazione . Mentre ci sei, puoi usare LaTeX / MathJax .
Ertxiem - ripristina Monica
Utilizzando il nostro sito, riconosci di aver letto e compreso le nostre Informativa sui cookie e Informativa sulla privacy.
Licensed under cc by-sa 3.0 with attribution required.