Separação de fonte cega da mistura convexa?

18

Suponhamos que temos fontes independentes, X 1 , X 2 , . . . , X n e observo m misturas convexas: Y 1nX1,X2,...,Xnm

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

com para todos os i e a i j0 para todos os i , j .jaij=1iaij0i,j

Qual é o estado da arte na recuperação de de Y ?XY

O PCA está fora de questão porque preciso que os componentes sejam identificáveis. Examinei a ICA e a NMF - não consigo encontrar nenhuma maneira de impor a não-negatividade dos coeficientes de mistura da ICA, e a NMF não parece maximizar a independência.

anônimo
fonte
2
Eu pensaria que isso deveria ser chamado de "análise de componentes independentes não-negativos", mas parece que esse nome foi usado para a ACI com a restrição de não-negatividade nas fontes , não na matriz de mistura A ( eecs.qmul.ac.uk/ ~ markp / 2003 / Plumbley03-algoritmos-c.pdf ). Portanto, isso não se aplica ao seu caso. Pergunta interessante. XA
Ameba diz Reinstate Monica
Você não quer que as somas atropelem j em vez de i? Você pode assumir que as fontes são aproximadamente gaussianas? se eles são unimodais e têm decaimento suficientemente rápido, é possível que a instalação de um GMM seja suficiente.
Yair Daon
@YairDaon Ah, sim, obrigado. Infelizmente, as fontes são discretas e nem parecem misturas de gaussianos. Mas talvez eu pudesse aproximar-los grosseiramente como misturas gaussianas e depois refinar ainda mais. Mas seria bom ter algo mais geral / robusto
anônimo
1
Quais algoritmos ICA você já tentou? Estou um pouco enferrujado, mas acho que a suposição de não-negatividade dos coeficientes de mistura pode ser imposta em alguns algoritmos que assumem certos modelos para os sinais, como o algoritmo WASOBI (Identificação de cortinas de segunda ordem ajustada por pesos), pois pressupõe que você possa modele os sinais como processos de RA e, portanto, você pode impor condições nos coeficientes.
Néstor
1
As fontes são todas suportadas no conjunto {1,2, ..., 96}
anônimo

Respostas:

0

Isso poderia ser alcançado usando uma não linearidade exponencial em vez do tanh típico / padrão (), se X também não for negativo.

Fórmula 40 em https://www.cs.helsinki.fi/u/ahyvarin/papers/NN00new.pdf e disponível na maioria das implementações.

Por exemplo, no sklearn, basta usar fun = 'exp' https://scikit-learn.org/stable/modules/generated/sklearn.decomposition.FastICA.html

Henrique Mendonça
fonte
1
Bem-vindo ao Stats.SE. Você pode editar sua resposta e expandi-la para explicar as principais etapas dos links fornecidos? Dessa forma, as informações são pesquisáveis ​​aqui (e algumas vezes os links quebram). Você pode dar uma olhada em alguma ajuda de formatação . Enquanto você está nisso, você pode usar o LaTeX / MathJax .
Ertxiem - restabelece Monica