machine learning - Fast ICA using scikit learn- reconstruction error analysis -


i trying use fastica procedure in scikitlearn. validation purposes tried understand difference between pca , ica based signal reconstruction.

the original number of observed signals 6 , tried use 3 reconstruction independent components . problem both ica , pca result in same reconstruction errors no matter norm use. can 1 throw light happening here.

the code below:

 pca = pca(n_components=3)  icamodel = fastica(n_components=3,whiten=true)   data = trainingdatadict[yearspan][riskfactornames]   pcr_dict[yearspan] = pd.dataframe(pca.fit_transform(data),                                     columns=['pc1','pc2','pc3'],index=data.index)   icr_dict[yearspan] = pd.dataframe(icamodel.fit_transform(data),                                     columns=['ic1','ic2','ic3'],index=data.index)  '------------------------inverse transform of ic , pcs -----------'   pca_new_data_df = pd.dataframe(pca.inverse_transform(pcr_dict[yearspan]),                                    columns =['f1','f2','f3'],index = data.index)   ica_new_data_df = pd.dataframe(icamodel.inverse_transform(icr_dict[yearspan]),                                    columns =['f1','f2','f3'],index = data.index)

below way measure reconstruction error

'-----------reconstruction errors------------------'  print 'pca reconstruction error l2 norm:',np.sqrt((pca_new_data_df - data).apply(np.square).mean())   print 'ica reconstruction error l2 norm:',np.sqrt((ica_new_data_df - data).apply(np.square).mean())   print 'pca reconstruction error l1 norm:',(pca_new_data_df - data).apply(np.absolute).mean()   print 'ica reconstruction error l1 norm:',(ica_new_data_df - data).apply(np.absolute).mean()

below description of tails of pc , ics

pc stats :  ('2003', '2005')         kurtosis  skewness pcr_1 -0.001075 -0.101006 pcr_2  1.057140  0.316163 pcr_3  1.067471  0.047946   ic stats :  ('2003', '2005')         kurtosis  skewness icr_1 -0.221336 -0.204362 icr_2  1.499278  0.433495 icr_3  3.654237  0.072480

below results of reconstruction

pca reconstruction error l2 norm:  sptr        0.000601 sptrmdcp    0.001503 ru20intr    0.000788 lbustruu    0.002311 lf98truu    0.001811 nddueafe    0.000135 dtype: float64   ica reconstruction error l2 norm :  sptr        0.000601 sptrmdcp    0.001503 ru20intr    0.000788 lbustruu    0.002311 lf98truu    0.001811 nddueafe    0.000135

even l1 norms same. bit confused!


Comments

Post a Comment

Popular posts from this blog

Is there a better way to structure post methods in Class Based Views -

i'm using class based views in django , never specify url in ajax calls , omit action parameter in form because know these requests go through post method of class based view. notice view - profile page make many requests depending on user updating on page - post method become convoluted , if/else mess example: def post(self, request, *args, **kwargs): # handle user changing available date if request.post.get('availabledate') != none: self.updateavailabledate(request) return jsonresponse({'result':'success'}) if request.post.get('newprojectname') != none: creative_user = creativeuserprofile.objects.get(id = request.user.id) project = project.create_project(creative_user, request.post['newprojectname']) project.save() return jsonresponse({'projectid': project.id }) if request.post.get('projectimage') != none , request.files['file']: proj
Read more

performance - Why is XCHG reg, reg a 3 micro-op instruction on modern Intel architectures? -

i'm doing micro-optimization on performance critical part of code , came across sequence of instructions (in at&t syntax): add %rax, %rbx mov %rdx, %rax mov %rbx, %rdx i thought had use case xchg allow me shave instruction , write: add %rbx, %rax xchg %rax, %rdx however, dimay found agner fog's instruction tables , xchg 3 micro-op instruction 2 cycle latency on sandy bridge, ivy bridge, broadwell, haswell , skylake. 3 whole micro-ops , 2 cycles of latency! 3 micro-ops throws off 4-1-1-1 cadence , 2 cycle latency makes worse original in best case since last 2 instructions in original might execute in parallel. now... cpu might breaking instruction micro-ops equivalent to: mov %rax, %tmp mov %rdx, %rax mov %tmp, %rdx where tmp anonymous internal register , suppose last 2 micro-ops run in parallel latency 2 cycles. given register renaming occurs on these micro-architectures, though, doesn't make sense me done way. why wouldn't register renamer
Read more

jquery - Responsive Navbar with Sub Navbar -

good day, forgive me nav vocabulary not best. if question confusing apologize. i'm working on navigation bar needs centered on page @ 768px , above. items in nav open (or perhaps toggle) subnav lives directly below. these nav should appear centered on page. i have been able work through responsive portion of nav bar, having main items align vertically , sub nav menu of each item roll out below, , gave them indent separation. working well. however, issue not responsive side, behavior @ width greater 768px. happens toggling menu's sub navigation pushes other main items down page below sub nav. here's code: $(document).ready(function () { //stick in fixed 100% height behind navbar don't wrap $('#slide-nav.navbar-inverse').after($('<div class="inverse" id="navbar-height-col"></div>')); $('#slide-nav.navbar-default').after($('<div id="navbar-height-col"></d
Read more