{"action":"created","issue":{"url":"https://api.github.com/repos/shogun-toolbox/shogun/issues/4267","repository_url":"https://api.github.com/repos/shogun-toolbox/shogun","labels_url":"https://api.github.com/repos/shogun-toolbox/shogun/issues/4267/labels{/name}","comments_url":"https://api.github.com/repos/shogun-toolbox/shogun/issues/4267/comments","events_url":"https://api.github.com/repos/shogun-toolbox/shogun/issues/4267/events","html_url":"https://github.com/shogun-toolbox/shogun/issues/4267","id":318896768,"number":4267,"title":"Problem with FastICA algorithm","user":{"login":"fyoda","id":37776179,"avatar_url":"https://avatars1.githubusercontent.com/u/37776179?v=4","gravatar_id":"","url":"https://api.github.com/users/fyoda","html_url":"https://github.com/fyoda","followers_url":"https://api.github.com/users/fyoda/followers","following_url":"https://api.github.com/users/fyoda/following{/other_user}","gists_url":"https://api.github.com/users/fyoda/gists{/gist_id}","starred_url":"https://api.github.com/users/fyoda/starred{/owner}{/repo}","subscriptions_url":"https://api.github.com/users/fyoda/subscriptions","organizations_url":"https://api.github.com/users/fyoda/orgs","repos_url":"https://api.github.com/users/fyoda/repos","events_url":"https://api.github.com/users/fyoda/events{/privacy}","received_events_url":"https://api.github.com/users/fyoda/received_events","type":"User","site_admin":false},"labels":[],"state":"open","locked":false,"assignee":null,"assignees":[],"milestone":null,"comments":24,"created_at":"2018-04-30T13:50:21Z","updated_at":"2018-05-14T09:46:49Z","closed_at":null,"author_association":"NONE","body":"Hello,\r\n\r\nI'm having problem when using the FastICA algorithm available within shogun.\r\nI'm elaborating the same data and compare shogun results with the matlab implementation available here: https://research.ics.aalto.fi/ica/fastica/code/dlcode.shtml\r\n\r\nThe test data I use are: https://www.dropbox.com/s/2j7diqad2zx7t77/test_data.txt?dl=0\r\n\r\nHere is the sample code, `X` contains 16x30000 samples:\r\n```\r\n    // 16xN sample matrix\r\n    CDenseFeatures< float64_t >* mixed_signals = new CDenseFeatures<float64_t >(X);\r\n    CFastICA* ica = new CFastICA();\r\n    ica->set_max_iter(1000);        // matlab default\r\n    ica->set_tol(0.0001);           // matlab default\r\n    CFeatures* ica_signals = ica->apply(mixed_signals);\r\n\r\n    auto FeatureMatrix = CDenseFeatures<float64_t>::obtain_from_generic(ica_signals)->get_feature_matrix();\r\n    SGMatrix< float64_t> mixing_matrix = ica->get_mixing_matrix();  // Matrix A 16x16\r\n\r\n    int nr = FeatureMatrix.num_rows;\r\n    int nc = FeatureMatrix.num_cols;\r\n```\r\n\r\nNow in `mixing_matrix` I have the A matrix and I should get the inverse to have the W matrix.\r\nInstead in `FeatureMatrix` I should get the independent components. \r\nRight?\r\n\r\nHowever, when I try to plot the 16 independent components there is no information, signals are almost zero everywhere. Am I doing something wrong?\r\n\r\nThanks a lot."},"comment":{"url":"https://api.github.com/repos/shogun-toolbox/shogun/issues/comments/388759603","html_url":"https://github.com/shogun-toolbox/shogun/issues/4267#issuecomment-388759603","issue_url":"https://api.github.com/repos/shogun-toolbox/shogun/issues/4267","id":388759603,"user":{"login":"vigsterkr","id":30353,"avatar_url":"https://avatars1.githubusercontent.com/u/30353?v=4","gravatar_id":"","url":"https://api.github.com/users/vigsterkr","html_url":"https://github.com/vigsterkr","followers_url":"https://api.github.com/users/vigsterkr/followers","following_url":"https://api.github.com/users/vigsterkr/following{/other_user}","gists_url":"https://api.github.com/users/vigsterkr/gists{/gist_id}","starred_url":"https://api.github.com/users/vigsterkr/starred{/owner}{/repo}","subscriptions_url":"https://api.github.com/users/vigsterkr/subscriptions","organizations_url":"https://api.github.com/users/vigsterkr/orgs","repos_url":"https://api.github.com/users/vigsterkr/repos","events_url":"https://api.github.com/users/vigsterkr/events{/privacy}","received_events_url":"https://api.github.com/users/vigsterkr/received_events","type":"User","site_admin":false},"created_at":"2018-05-14T09:46:49Z","updated_at":"2018-05-14T09:46:49Z","author_association":"OWNER","body":"@fyoda that's a good question. the negative coefficient is probably coming from the different eigen solver - namely that the eigenvectors are pointing to the opposite direction..... the amplitude difference is a good question though.... is the amplitude difference is consistent? i mean it looks like, meaning they are all scaled with the same factor, or?"}}
{"actor":{"display_login":"vigsterkr"}}