,j)$ entry of the matrix is the Pearson distance $s_{ij}$ color coded according to the key on the right, i.e. blue entries correspond to pairs of nodes having different time-series while yellow entries correspond to pairs with similar time-series. Thus one can coarsely distinguish between high- and low-degree nodes. This analysis on its own leads to a gross overestimate of the structural parameter, see panel (d). However, using this coarse classification and mathematical results on coupled systems (\[eq1\]) satisfying assumptions (i), (ii) and (iii) makes it possible to find an approximation of the evolution rule at each node, see panel (e). Using this additional information, one can obtain a description of the network connectivity structure, see panel (f). The network used in this figure has 1000 nodes, but only 50 are shown for clarity.[]{data-label="fig1"}](main_fig_v5_low_3_reduced.pdf){width="\linewidth"}

[**Scale-free networks of coupled [ bursting]{} neurons**]{}.

We consider coupled bursting neurons with excitatory synapses [@rulkov2001] in scale-free networks (see Methods for details on simulation of scale-free networks and formulation of the model). Our techniques work equally well for spiking dynamics.

We generate a scale-free network with $N=10^4$ nodes such that the probability of having a node of degree $k$ is proportional to $k^{-\gamma}$, where $\gamma=2.53$. For this reconstruction we only need 2000 data points for each node. Again, to every pair of time series $y_i$ and $y_j$ we assign a *Pearson distance* $s_{ij}\ge0$ and the node intensity $S_{i} = \sum_{j}s_{ij}$. The empirical distribution of the intensities $S_i$ approximates the degree distribution of the network, see Fig \[fig1\](d). In the example here, the estimated structural exponent from the distribution of $S_i$ is $\gamma_{\rm est} = 3.1$, which yields a relative error of nearly 25% with respect to the true value of $\gamma$ (see also the plots in Figure \[fig2\] a)). The functional network therefore overestimates $\gamma$, which has drastic consequences for the predicted character of the network. For example, the number of connections of a hub for a scale-free network is concentrated
, j)$ entry of the matrix is the Pearson distance $ s_{ij}$ color coded according to the keystone on the right field, i.e. blue entries represent to couple of nodes having unlike time - series while yellow entry correspond to pairs with exchangeable time - series. therefore one can coarsely distinguish between high- and low - academic degree nodes. This psychoanalysis on its own lead to a gross overestimate of the structural parameter, witness panel (d). However, practice this coarse classification and numerical results on copulate systems (\[eq1\ ]) satisfying assumptions (iodine), (ii) and (iii) makes it possible to find an approximation of the evolution rule at each node, see panel (e). Using this additional information, one can obtain a description of the network connectivity structure, see jury (f). The network used in this name hold 1000 nodes, but only 50 are indicate for clarity.[]{data - label="fig1"}](main_fig_v5_low_3_reduced.pdf){width="\linewidth " } 

 [ * * Scale - free networks of coupled [ bursting ] { } nerve cell * * ] { }. 

 We consider coupled bursting neurons with excitatory synapses [ @rulkov2001 ] in scale - free networks (see Methods for details on model of scale - complimentary networks and formulation of the model). Our techniques work evenly well for spiking dynamics. 

 We generate a scale - free network with $ N=10 ^ 4 $ nodes such that the probability of having a node of academic degree $ k$ is proportional to $ k^{-\gamma}$, where $ \gamma=2.53$. For this reconstruction we only need 2000 data points for each node. Again, to every couple of time series $ y_i$ and $ y_j$ we assign a * Pearson distance * $ s_{ij}\ge0 $ and the node volume $ S_{i } = \sum_{j}s_{ij}$. The empirical distribution of the saturation $ S_i$ approximates the degree distribution of the network, see Fig   \[fig1\](d). In the example here, the calculate structural exponent from the distribution of $ S_i$ is $ \gamma_{\rm est } = 3.1 $, which yields a relative error of nearly 25% with respect to the true value of $ \gamma$ (see besides the plots in Figure \[fig2\ ] a) ). The functional net therefore overestimates $ \gamma$, which has drastic consequence for the predicted character of the network. For example, the number of connections of a hub for a scale - free network is concentrated
,j)$ ejtry of the matrix is tht Pearson distancg $w_{ij}$ conor cosed accofding to the key on the righv, i.e. blue entries correspond to pairs of nodes hqvinj different time-series wmnle ysplow xntries corresppnd to paiss with similas gile-series. Thus one can coarsely distynguish bftween high- anq loe-qegrsv kodes. This analysis on its own lsads to a gross overrstimate of the structural parwmeter, see panel (d). However, usung eyis coarse cuassificatpmn and matgematical results on coupled syrtems (\[eq1\]) satisfinbg wvsumptions (m), (ii) agd (iii) makes it posshble to find an approwimatmon if the evolution rule at each node, see kanel (e). Ushny this additional infirnatiot, ona cav obgaih e dsscriphioi of the nefwork connextivity structure, ste [qnel (f). The nefwork tsqd in this figure has 1000 nodes, but only 50 dre shown for clarity.[]{data-lqbel="fig1"}](main_fig_v5_low_3_redoced.pdf){wideh="\linewidth"}

[**Scale-free networks of coupled [ burstinc]{} neucovs**]{}.

Wt gjvwifer coupled bursting neurons with excitatory fgnspxes [@rulkov2001] in fcale-free nrtaoths (see Methodr for bstzils on simulation of scaje-frew networkf anc formulation of the model). Iur techniquvs wirk equally well flr spiking bynamivs.

We benerate a scale-free necwork sith $N=10^4$ noded such thzg the probabilith on hdving a node of degree $k$ if proportmonal to $k^{-\gaoma}$, ehere $\damma=2.53$. For hhis vaconstruction we ojly nged 2000 ddta points for each node. Again, to every pemr of time setiev $y_p$ and $y_j$ ce asslgn a *Pearson dystance* $s_{ij}\ge0$ cnd the uode ivtensity $S_{p} = \sum_{j}s_{io}$. The empirisal distributhln of the invensities $S_i$ qppriximater the degree dixtributiou of rhe network, see Fin \[fig1\](d). Jn the example ktre, the estimated xtrjcttrwl ev[onent from dhe aisgtibutkon of $S_i$ iw $\gamka_{\rm est} = 3.1$, which yiends z relative error og kearly 25% wuth respqct to the trie value of $\gamma$ (dee anso the llojs in Figure \[fig2\] a)). The functionzl networn tmerefore overqstinates $\gamma$, chich has drastic consequences for the pcedicted character of jhe network. For exampkt, the number pf cognections of a hub for a scalw-free network is goncentrated
,j)$ entry of the matrix is the $s_{ij}$ coded according the key on correspond pairs of nodes different time-series while entries correspond to pairs with similar Thus one can coarsely distinguish between high- and low-degree nodes. This analysis on own leads to a gross overestimate of the structural parameter, see panel (d). using coarse and results on coupled systems (\[eq1\]) satisfying assumptions (i), (ii) and (iii) makes it possible to find approximation of the evolution rule at each node, panel (e). Using this information, one can obtain a of network connectivity see (f). network used in figure has 1000 nodes, but only 50 are shown for clarity.[]{data-label="fig1"}](main_fig_v5_low_3_reduced.pdf){width="\linewidth"} [**Scale-free networks of coupled [ bursting]{} We consider neurons with synapses in networks (see Methods on simulation of scale-free networks and model). Our techniques work equally well for spiking We generate scale-free network with $N=10^4$ nodes such the probability of having a node of degree is proportional to $k^{-\gamma}$, where $\gamma=2.53$. For this reconstruction we only need 2000 data points node. Again, to every of time series and we a distance* $s_{ij}\ge0$ the node intensity $S_{i} = \sum_{j}s_{ij}$. The empirical distribution of the $S_i$ approximates the degree distribution of the network, see Fig the here, the estimated exponent from the distribution $S_i$ $\gamma_{\rm est} = 3.1$, a error with to true value of $\gamma$ also the plots in Figure a)). The functional network drastic consequences for the predicted character of the For example, the number of connections of hub for a scale-free network is concentrated
,j)$ entry of the matrix is the PeaRson distanCe $s_{ij}$ ColOr cOdEd acCordIng to the key on tHE rigHt, i.e. blue entries correspOnd to PaIRs of NOdEs havIng diffEReNT TimE-sErIes WhILe YelloW enTries coRrespond to PaiRs With similar tIMe-Series. Thus One Can coarsely dIstInguisH bEtwEEn higH- anD low-dEgree nODes. ThiS analysis On ITs own lEAds to a gROSs OverEstimate of the struCTuRAl parameter, see Panel (d). hoWEvER, UsiNg tHis coarse cLaSsifiCAtion anD MaTHEMatICal results on cOupled systeMS (\[eq1\]) SatisfYiNg aSSumptiOns (i), (iI) aND (iiI) makes it posSiblE to find an ApproxIMation oF The evolUtion rUle At eAch nODe, SeE paNeL (E). UsINg ThiS AddItional iNfOrMatioN, one CAN OBtaiN a dEscrIptioN of the network ConNectIVitY struCture, See pAnEl (f). ThE netwoRk useD iN this figure has 1000 nOdes, But only 50 arE shOwN foR cLaritY.[]{Data-laBel="Fig1"}](Main_fig_V5_low_3_redUCed.PdF){WIDtH="\linewidth"}

[**Scale-freE nETWoRks of couPled [ buRStInG]{} Neurons**]{}.

WE cOnsIder COUpled BursTInG neurons With exCItAtOry synaPsEs [@rulkOv2001] In sCalE-free NEtwoRks (see methods fOr detAIls on simulatioN Of scale-free neTWoRKS aND forMulAtion of the mOdel). oUr teChniQUeS woRK equaLly weLl FOr SPiking dynamics.

We genErAte a scAle-frEe network with $n=10^4$ nodes such THAT the probAbilITy OF having a node of DegreE $k$ is proporTIonal to $k^{-\Gamma}$, Where $\gamMa=2.53$. For this REConstrucTioN we OnlY neED 2000 DaTa points for eaCH Node. agAin, to evEry Pair of tIme SerIes $Y_i$ aNd $Y_j$ we assigN a *PearsoN dIsTaNcE* $s_{iJ}\ge0$ anD The node iNtEnsItY $S_{i} = \Sum_{j}s_{IJ}$. The emPiricAl diStRiBUtiOn of the INtENSitiEs $s_i$ ApprOxiMaTes thE degREe dIstribuTion of the NetWOrk, sEe fiG \[fig1\](d). In The example herE, tHe estimateD sTruCtural EXPonent frOm the distribution of $S_i$ is $\GAmma_{\rm eSt} = 3.1$, wHich yIeldS a relativE erRor of nEarLY 25% with rEspect To the TrUe vALUe of $\gAMMa$ (See AlSo the plots IN figUre \[fiG2\] a)). the fUnctionAl network therefore OVerEstimates $\gammA$, whIch hAS DrAstIC cONseQuENceS FOr the predicted cHaracter of ThE NeTwork. For exAMplE, tHe numbeR of connEctioNS of a hub For a scale-Free netwoRk Is coNCEntRated
,j)$ entry of the matrix i s the Pear son d ist anc e$s_{ ij}$  color coded a c cord ing to the key on theright ,i .e.b lu e ent ries co r re s p ond  t opai rs of  node s h aving d ifferent t ime -s eries whiley el low entrie s c orrespond to  pa irs wi th  si m ilartim e-ser ies. T h us one  can coar se l y dist i nguishb e tw eenhigh- and low-deg r ee nodes. This an alysis  o n  i t s  ow n l eads to agr oss o v erestim a te o f  th e  structural p arameter, s e e p anel ( d) . H o wever,  usin gt his  coarse cla ssif ication a nd mat h ematica l  result s on c oup led  sys t em s(\[ eq 1 \]) sa tis f yin g assump ti on s (i) , (i i ) a nd ( iii ) ma kes i t possible to  fi nd a n  ap proxi matio n of  t he ev olutio n rul eat each node, s ee p anel (e).  Us in g t hi s add i tional  in for mation,  one ca n  ob ta i n adescription of the  n e t wo rk conne ctivit y  s tr u cture, s ee  pa nel( f ). Th e ne t wo rk usedin thi s  f ig ure has  1 000 no de s,but  only 50 a re sho wn for c larit y .[]{data-label = "fig1"}](main _ fi g _ v5 _ low_ 3_r educed.pdf) {wid t h="\ line w id th" } 

[** Scale -f r ee networks of coupled  [  burst ing]{ } neurons**]{ }.

We con s i d er coupl ed b u rs t ing neurons wi th ex citatory s y napses [ @rulk ov2001]in scale- f r ee netwo rks  (s eeMet h o ds  for detailso n  sim ul ation o f s cale-fr eenet wor ksan d formula tion ofth emo de l).  Ourt echnique swor kequ allyw ell fo r spi king  d yn a mic s.

Weg en e r ateasc ale- fre enetwo rk w i th$N=10^4 $ nodes s uch that  t he  probab ility of havi ng  a node of  d egr ee $k$ i s propor tional to $k^{-\gamma}$ ,  where$\g amma= 2.53 $. For th isrecons tru c tion w e only  need  2 000 d ata p o i nt s f or  each node . Aga in, t oever y pairof time series $y_ i $ a nd $y_j$ we a ssi gn a * Pe ars o nd ist an c e*$ s _{ij}\ge0$ andthe node i nt e ns ity $S_{i} = \ su m_{j}s_ {ij}$.The e m pirical  distribu tion of t he  int e n sit ies $S_i$approxim ates thed egree di strib uti on ofth e n etwor k, see Fig  \[fi g1\](d ).  In th e exa mp le here,  the estimated structur al exp onent  fr om the di str i but ion of $S _i$is $\gamma _{\ rmest}= 3 . 1$, w hich yi eld s  a re lati v e error o f  n ear l y  2 5% with res p e c t t o the  tr u e valu e of  $\gamma$ (see al s o the plots in  Fig u r e \ [fi g 2\]a) ). The functio nal  n e t work the re fore overes timates$\ g amma$ , whic h hasdrastic c on s equenc es f orthe predi cte dc haracte rof the ne twor k.  For e xample ,  the n umber of connect ionso f  a hu b  fo r a s ca le-free netw ork is con centrated
,j)$ entry_of the_matrix is the Pearson_distance $s_{ij}$_color_coded according_to_the key on_the right, i.e._blue entries correspond to_pairs of nodes_having_different time-series while yellow entries correspond to pairs with similar time-series. Thus one can_coarsely_distinguish between_high-_and_low-degree nodes. This analysis on_its own leads to a_gross overestimate_of the structural parameter, see panel (d). However,_using_this coarse classification_and mathematical results on coupled systems (\[eq1\]) satisfying assumptions_(i), (ii) and (iii) makes it_possible to find_an_approximation_of the evolution rule_at each node, see panel (e)._Using this additional information, one can_obtain a description of the network connectivity_structure, see panel (f). The network_used in this figure has_1000 nodes,_but only 50 are shown_for clarity.[]{data-label="fig1"}](main_fig_v5_low_3_reduced.pdf){width="\linewidth"}

[**Scale-free networks_of coupled_[ bursting]{} neurons**]{}.

We_consider coupled bursting neurons with excitatory_synapses [@rulkov2001] in_scale-free networks (see Methods for details_on_simulation of scale-free_networks_and_formulation of_the model). Our_techniques_work equally_well_for spiking dynamics.

We generate a scale-free_network_with $N=10^4$ nodes such that the probability_of having a node_of_degree $k$ is proportional_to $k^{-\gamma}$, where $\gamma=2.53$. For_this reconstruction we only need 2000_data points_for each_node. Again, to every pair of time series $y_i$ and $y_j$_we assign a *Pearson distance* $s_{ij}\ge0$_and the node intensity_$S_{i} =_\sum_{j}s_{ij}$._The empirical distribution_of_the intensities_$S_i$ approximates the degree distribution of the_network, see_Fig \[fig1\](d). In the example here, the_estimated structural exponent from_the_distribution of $S_i$ is $\gamma_{\rm est}_= 3.1$, which yields a relative_error of nearly 25% with_respect_to_the true value of $\gamma$_(see also the plots in Figure_\[fig2\] a)). The_functional network therefore overestimates $\gamma$, which has_drastic_consequences for the predicted character of_the_network. For example, the number of_connections_of_a hub for a scale-free_network is concentrated