tilde \lambda) + O(\varepsilon^2)$, which describes the affine hyperplane generated by $x$ and $W$ in the embedding Euclidean (resp. Minkowski) space. The weights $\mu = \varepsilon \tilde \lambda / {\ensuremath{\mathds{1}}}{^{\text{\tiny T}}}\tilde \lambda$ converge to points at infinity (${\ensuremath{\mathds{1}}}{^{\text{\tiny T}}}\mu =0$) of the affine k-plane of normalized weights.

When reference points coalesce with an additional second order acceleration orthogonally to the subspace $W_x$, we conjecture that the affine span is not any more a great subspheres but a smaller one. This would include principal nested spheres (PNS) developed by [@jung_generalized_2010; @jung_analysis_2012] as a limit case of barycentric subspaces. It would be interesting to derive a similar procedure for hyperbolic spaces and to determine which types of subspaces could be obtained by such limits for more general non-local and higher order jets.

Barycentric subspace analysis {#Sec:BSA}
=============================

PCA can be viewed as the search for a sequence of nested linear spaces that best approximate the data at each level. In a Euclidean space, minimizing the variance of the residuals boils down to an independent optimization of orthogonal subspaces at each level of approximation, thanks to the Pythagorean theorem. This enables building each subspace of the sequence by adding (resp. subtracting) the optimal one-dimensional subspace iteratively in a forward (resp. backward) analysis. Of course, this property does not scale up to manifolds, for which the orthogonality of subspaces is not even well defined.

Flags of barycentric subspaces in manifolds
-------------------------------------------

[@damon_backwards_2013] have argued that the nestedness of approximation spaces is one of the most important characteristics for generalizing PCA to more general spaces. Barycentric subspaces can easily be nested, for instance by adding or removing one or several points at a time, to obtains a family of embedded submanifolds which generalizes flags of vector spaces.

A flag of a vector space $V$ is a filtration of subspaces (an increasing sequence of subspaces, where each subspace is a proper subspace of the next): $\{0\} = V_0 \subset V_1 \subset V_2 \subset \cdots \subset V_k = V$.
tilde \lambda) + O(\varepsilon^2)$, which describes the affine hyperplane generated by $ x$ and $ W$ in the embedding Euclidean (resp. Minkowski) space. The weight $ \mu = \varepsilon \tilde \lambda / { \ensuremath{\mathds{1}}}{^{\text{\tiny T}}}\tilde \lambda$ converge to point at infinity ($ { \ensuremath{\mathds{1}}}{^{\text{\tiny T}}}\mu = 0 $) of the affine k - plane of normalized weight. 

 When reference points blend with an additional second order acceleration orthogonally to the subspace $ W_x$, we conjecture that the affine couple is not any more a great subspheres but a smaller one. This would admit principal nested spheres (PNS) developed by [ @jung_generalized_2010; @jung_analysis_2012 ] as a terminus ad quem case of barycentric subspaces. It would be interesting to deduce a similar operation for hyperbolic quad and to determine which types of subspaces could be prevail by such limits for more general non - local and higher ordering jets. 

 Barycentric subspace analysis { # Sec: BSA } 
 = = = = = = = = = = = = = = = = = = = = = = = = = = = = = 

 PCA can be viewed as the search for a sequence of nested analogue spaces that best approximate the data at each level. In a Euclidean quad, minimizing the variance of the residuals boils down to an independent optimization of orthogonal subspaces at each level of approximation, thanks to the Pythagorean theorem. This enables build up each subspace of the sequence by adding (resp. subtract) the optimum one - dimensional subspace iteratively in a advancing (resp. backward) analysis. Of course, this property does not scale up to manifolds, for which the orthogonality of subspaces is not even well defined. 

 Flags of barycentric subspaces in manifold 
 ------------------------------------------- 

 [ @damon_backwards_2013 ] have argued that the nestedness of approximation spaces is one of the most important characteristics for generalize PCA to more general spaces. Barycentric subspaces can easily be cuddle, for instance by total or removing one or several points at a time, to obtains a syndicate of embedded submanifolds which generalizes flags of vector space. 

 A iris of a vector space $ V$ is a filtration of subspaces (an increase sequence of subspace, where each subspace is a proper subspace of the next ): $ \{0\ } = V_0 \subset V_1 \subset V_2 \subset \cdots \subset V_k = V$.
tilfe \lambda) + O(\varepsilon^2)$, wmich describes tkw affiie hypedplane gdnerated by $x$ and $W$ in the elbwddint Euclidean (resp. Minkodski) spacv. The weithts $\nu = \varepsmmon \tilde \lamgfa / {\znwuremath{\mathds{1}}}{^{\jext{\tiny T}}}\tinde \lambda$ coneefgz to points at infinity (${\ensuremath{\maehds{1}}}{^{\texy{\tlny T}}}\mu =0$) of thg affpnq k-pmane of normalized weights.

When refsrence koints coalesce wiyh an additional second orfer wcceleration orthohonally to jge fybspace $W_x$, wd conjectugz that the zffine span is not any more a gfeat xubspheres byt w smaller onx. This would include princhpal nexted spheres (PKS) deteloped by [@jung_generalizev_2010; @jung_analysis_2012] as a jimit casa kf barycentric suvspaces. It foula be inuerxstjng to decive a simimar proceduee for hyperbolic skacqw and to detedmine rhych types of subspaces could be obtaineg bg such limits for more teneral non-local and jigher orqer jets.

Barycentric subspace analysis {#Sec:BSA}
=============================

PCA cdn be xieccq qs the search for a sequence of nested linear s[zcts nhat best approxiiate the daya ay each level. Iv a Euemisean space, minimizlng the variqnce of tre rrsiduals boils down to an ibdependent oitimuzation of orthogoual subspacer at each level of approximation, thahks to the Oythagorezv theorem. This evabkev building each subspace os the seqnence by addkng (tesp. suftracting) hhe oidimal one-dimensionwl suyspaca iterativfly in a forward (resp. backward) eialysis. Of coorsa, tvis propzrty dpes not scale up to manifolbs, for wkich tfe orthogohality mf subspacef is not even aell defined.

Hlags of faryxenteic subroaces in manifplds
-------------------------------------------

[@damon_ycckwards_2013] yave argued that tme nerfedness of apprixination spaces ix ove jf tie mofd important wharxctdtistizs for genevalkzinb PCA to more generan spzces. Barycentric sibfpaces cqn easilr be nested, fpr instance by addlng oc remoting ome jr several points at a time, to obtains w fwmily of embqddee submanifolbs which generalizes flags of vector spares.

A flag of a vector space $V$ is a filtratnok of subspacxs (an yncreasinc sequence of subspaxes, where each sunspace is a proper subapace mf thf next): $\{0\} = V_0 \subset V_1 \subset V_2 \subset \cdots \subset V_k = V$.
tilde \lambda) + O(\varepsilon^2)$, which describes the generated $x$ and in the embedding weights = \varepsilon \tilde / {\ensuremath{\mathds{1}}}{^{\text{\tiny T}}}\tilde converge to points at infinity (${\ensuremath{\mathds{1}}}{^{\text{\tiny =0$) of the affine k-plane of normalized weights. When reference points coalesce with additional second order acceleration orthogonally to the subspace $W_x$, we conjecture that the span not more great subspheres but a smaller one. This would include principal nested spheres (PNS) developed by [@jung_generalized_2010; as a limit case of barycentric subspaces. It be interesting to derive similar procedure for hyperbolic spaces to which types subspaces be by such limits more general non-local and higher order jets. Barycentric subspace analysis {#Sec:BSA} ============================= PCA can be viewed as search for of nested spaces best the data at In a Euclidean space, minimizing the residuals boils down to an independent optimization of subspaces at level of approximation, thanks to the theorem. This enables building each subspace of the by adding (resp. subtracting) the optimal one-dimensional subspace iteratively in a forward (resp. backward) analysis. this property does not up to manifolds, which orthogonality subspaces not even defined. Flags of barycentric subspaces in manifolds ------------------------------------------- [@damon_backwards_2013] have argued the nestedness of approximation spaces is one of the most for PCA to more spaces. Barycentric subspaces can be for instance by adding one several time, obtains family of embedded submanifolds generalizes flags of vector spaces. flag of a vector of subspaces (an increasing sequence of subspaces, where subspace is a proper subspace of the $\{0\} = V_0 \subset V_1 \subset V_2 \subset \cdots \subset V_k =
tilde \lambda) + O(\varepsilon^2)$, whiCh describeS the aFfiNe hYpErplAne gEnerated by $x$ and $w$ In thE embedding Euclidean (resP. MinkOwSKi) spACe. the weIghts $\mu = \VArEPSilOn \TiLde \LaMBdA / {\ensuRemAth{\mathDs{1}}}{^{\text{\tiny t}}}\tiLdE \lambda$ conveRGe To points at InfInity (${\ensuremAth{\Mathds{1}}}{^{\TeXt{\tINy T}}}\mu =0$) Of tHe affIne k-plANe of noRmalized wEiGHts.

WheN ReferenCE PoInts Coalesce with an addITiONal second order AcceleRaTIoN ORthOgoNally to the SuBspacE $w_x$, we conJEcTURE thAT the affine spaN is not any moRE a gReat suBsPheREs but a SmallEr ONe. THis would incLude Principal Nested SPheres (Pns) develoPed by [@jUng_GenEralIZeD_2010; @jUng_AnALysIS_2012] aS a lIMit Case of baRyCeNtric SubsPACES. It wOulD be iNtereSting to derive A siMilaR ProCedurE for hYperBoLic spAces anD to deTeRmine which types Of suBspaces coUld Be ObtAiNed by SUch limIts For More genEral non-LOcaL aND HIgHer order jets.

BaryceNtRIC sUbspace aNalysiS {#seC:Bsa}
=============================

PCA can bE vIewEd as THE searCh foR A sEquence oF nesteD LiNeAr spaceS tHat besT aPprOxiMate tHE datA at eacH level. In A EuclIDean space, minimIZing the varianCE oF THe REsidUalS boils down tO an iNDepeNdenT OpTimIZatioN of orThOGoNAl subspaces at each leVeL of appRoximAtion, thanks to The PythagoREAN theorem. this ENaBLes building eacH subsPace of the sEQuence by AddinG (resp. subTracting) tHE Optimal oNe-dImeNsiOnaL SUbSpace iterativELY in a FoRward (reSp. bAckward) AnaLysIs. OF coUrSe, this proPerty doeS nOt ScAlE up To manIFolds, for WhIch ThE orThogoNAlity oF subsPaceS iS nOT evEn well dEFiNED.

FlaGs Of BaryCenTrIc subSpacES in ManifolDs
-------------------------------------------

[@damon_baCkwARds_2013] hAvE aRgued thAt the nestedneSs Of approximAtIon Spaces IS One of the Most important characteriSTics for GenEraliZing pCA to more GenEral spAceS. baryceNtric sUbspaCeS caN EAsily BE NeSteD, fOr instance BY AddIng or ReMoviNg one or Several points at a tiME, to Obtains a familY of EmbeDDEd SubMAnIFolDs WHicH GEneralizes flags Of vector spAcES.

A Flag of a vecTOr sPaCe $V$ is a fIltratiOn of sUBspaces (An increasIng sequenCe Of suBSPacEs, where eacH subspacE is a propeR SubspACe Of the NexT): $\{0\} = V_0 \subsEt v_1 \suBset V_2 \Subset \CDotS \subsEt V_k = V$.
tilde \lambda) + O(\vareps ilon^2)$,which  de scr ib es t he a ffine hyperpla n e ge nerated by $x$ and $W$  in t he embe d di ng Eu clidean (r e s p.Mi nk ows ki )  s pace.  Th e weigh ts $\mu =\va re psilon \tild e  \ lambda / { \en suremath{\ma thd s{1}}} {^ {\t e xt{\t iny  T}}} \tilde \lambd a$ conver ge to poi n ts at i n f in ity(${\ensuremath{\m a th d s{1}}}{^{\text {\tiny  T } }} \ m u = 0$)  of the af fi ne k- p lane of no r m a liz e d weights.

W hen referen c e p ointsco ale s ce wit h anad d iti onal second  ord er accele ration orthogo n ally to  the s ubs pac e $W _ x$ ,weco n jec t ur e t h atthe affi ne  s pan i s no t a n y mo rea gr eat s ubspheres but  asmal l erone.Thiswoul dinclu de pri ncipa lnested spheres(PNS ) develop edby  [@ ju ng_ge n eraliz ed_ 201 0; @jun g_analy s is_ 20 1 2 ]  a s a limit case ofba r y ce ntric su bspace s .It would be  i nte rest i n g toderi v ea simila r proc e du re  for hy pe rbolic  s pac esand t o  det ermine  which t ypeso f subspaces co u ld be obtaine d  b y su c h li mit s for moregene r al n on-l o ca l a n d hig her o rd e rj ets.

Barycentric s ub spaceanaly sis {#Sec:BSA }
======== = = = ======== ==== = == = ==

PCA can be  view ed as thes earch fo r a s equenceof nested l inear sp ace s t hat  be s t  a pproximate th e data  a t eachlev el. Ina E ucl ide ansp ace, mini mizing t he  v ar ia nce  of t h e residu al s b oi lsdownt o an i ndepe nden top t imi zationo fo r thog on al  sub spa ce s ateach lev el of a pproximat ion ,  tha nk sto thePythagorean t he orem. This  e nab les bu i l ding eac h subspace of the seque n ce by a ddi ng (r esp.  subtract ing ) theopt i mal on e-dime nsion al  su b s pacei t er ati ve ly in a fo r w ard  (res p.  bac kward)analysis. Of cours e , t his propertydoe s no t sc ale up toma n ifo l d s, for which th e orthogon al i ty  of subspa c esis  not ev en well  defi n ed.

Fl ags of ba rycentric  s ubsp a c esin manifol ds
----- --------- - ----- - -- ----- --- ------ -- --- -

[@ damon_ b ack wards _2013]  h ave ar guedth at the n estedness of approximat ion sp acesisone of th e m o stimportant  cha racteristi csfor  gene ral i zingPCAt omor e  gene rals paces. Ba r yc ent r i csubspaces c a n eas ily b e n e sted,forinstance by addin g  or removing o ne o r sev era l  poi nt s at a time, t o o bt a i ns a fam il y of embedd ed subma ni f oldswhichgenera lizes f l a gs of vec torspa ces.

A f lag  o f  a vect or  s p ace $V $ is  a  filtr ationo f su b s paces (an increa sings e quenc e  of  subs pa ces, wh e re e ach subspa ce is a pro per su bspa ce of  the ne xt ): $\{ 0\}  =  V_0 \subs e t V_1 \su bsetV_2 \su bs et \ cdo ts \su bset V _k =V$.
tilde \lambda)_+ O(\varepsilon^2)$,_which describes the affine_hyperplane generated_by_$x$ and_$W$_in the embedding_Euclidean (resp. Minkowski)_space. The weights $\mu_= \varepsilon \tilde_\lambda_/ {\ensuremath{\mathds{1}}}{^{\text{\tiny T}}}\tilde \lambda$ converge to points at infinity (${\ensuremath{\mathds{1}}}{^{\text{\tiny T}}}\mu =0$) of the_affine_k-plane of_normalized_weights.

When_reference points coalesce with an_additional second order acceleration orthogonally_to the_subspace $W_x$, we conjecture that the affine span_is_not any more_a great subspheres but a smaller one. This would_include principal nested spheres (PNS) developed_by [@jung_generalized_2010; @jung_analysis_2012]_as_a_limit case of barycentric_subspaces. It would be interesting to_derive a similar procedure for hyperbolic_spaces and to determine which types of_subspaces could be obtained by such_limits for more general non-local_and higher_order jets.

Barycentric subspace analysis {#Sec:BSA}
=============================

PCA_can be viewed_as the_search for a_sequence of nested linear spaces that_best approximate the_data at each level. In a_Euclidean_space, minimizing the_variance_of_the residuals_boils down to_an_independent optimization_of_orthogonal subspaces at each level of_approximation,_thanks to the Pythagorean theorem. This enables_building each subspace of_the_sequence by adding (resp._subtracting) the optimal one-dimensional subspace_iteratively in a forward (resp. backward)_analysis. Of_course, this_property does not scale up to manifolds, for which the orthogonality_of subspaces is not even well_defined.

Flags of barycentric subspaces_in manifolds
-------------------------------------------

[@damon_backwards_2013]_have_argued that the_nestedness_of approximation_spaces is one of the most important_characteristics for_generalizing PCA to more general spaces._Barycentric subspaces can easily_be_nested, for instance by adding or_removing one or several points at_a time, to obtains a_family_of_embedded submanifolds which generalizes flags_of vector spaces.

A flag of a_vector space $V$_is a filtration of subspaces (an increasing_sequence_of subspaces, where each subspace is_a_proper subspace of the next): $\{0\}_=_V_0_\subset V_1 \subset V_2 \subset_\cdots \subset V_k = V$.