 more beneficial in terms of resource management and precision, care would be needed in automating the grid refinements, and it is likely that a unique grid would be associated with each distinct firing pattern.

The sequential aspect of the proposed methodology provides the opportunity for real-time inference that has the potential to provide in-lab assistance during experimentation. In this framework, an interim SMC-MUNE analysis could help in identifying the choice of stimulus to apply in order to collect the best evidence to distinguish between competing hypotheses, as in Section \[sec:SimStudy\_under\]. The limitations of the present SMC-MUNE procedure to become a wholly online algorithm are the computational aspects discussed earlier and the post-processing stage to correct for potentially negligible estimates of the expected MU twitch forces. Solutions to these outstanding problems would increase the efficiency and accuracy of SMC-MUNE and, hence, the range of application.

Acknowledgements {#acknowledgements.unnumbered}
================

The authors would like to thank Dr. Christine Thomas from the Miami Project to Cure Paralysis who provided the data analysed in this paper and for specialist discussions.

Additional detail {#sec:APX}
=================

The prior sufficient statistics for the simulation and case studies are: ${\bar{m}}_0 = 0$, ${\bar{c}}_0 = 10^3$, ${\bar{a}}_0 = 0.5$, ${\bar{b}}_0 = 0.1$, ${\mathbf{m}}_0 = 40 {\mathbf{1}}_u$, $C_0 = 10^4I_u$ and $a_0 = 0.5$ where ${\mathbf{1}}_u$ is a unit $u$-vector and $I_u$ is the $u\times{u}$ identity matrix. The statistic $b_0$ is defined according to with $\delta=0.05$ and $\epsilon=0.2$. The upper bounds for the excitability parameter space are ${\eta_{\max}}= 1.1s_{\tau}$ and ${\lambda_{\max}}= 14$. In the case study, the upper bound for the scale parameter was reduced to ${\lambda_{\max}}=7$.

Resampling in Algorithm \[tab:Alg\], is performed by systematic sampling on the residuals of particle weights [@Hol06]. The number of particle samples and the number of rectangular lattice cells is initially $N=5000$ and $|\mathcal{G}|=30\times{30}$ respectively.


more beneficial in terms of resource management and precision, care would be needed in automatize the power system refinement, and it is likely that a unique power system would be associated with each distinct firing blueprint. 

 The sequential aspect of the propose methodology provides the opportunity for real - prison term inference that has the potential to put up in - lab aid during experimentation. In this framework, an interim SMC - MUNE analysis could help in identifying the option of stimulus to apply in order to gather the best evidence to distinguish between compete hypotheses, as in Section   \[sec: SimStudy\_under\ ]. The restriction of the present SMC - MUNE procedure to become a wholly online algorithm are the computational aspects discussed earlier and the post - processing stagecoach to correct for potentially negligible estimates of the expected MU twitch forces. Solutions to these outstanding problems would increase the efficiency and accuracy of SMC - MUNE and, hence, the range of application. 

 Acknowledgements { # acknowledgements.unnumbered } 
 = = = = = = = = = = = = = = = = 

 The authors would like to thank Dr. Christine Thomas from the Miami Project to Cure Paralysis who provided the data analysed in this composition and for specialist discussions. 

 Additional contingent { # sec: APX } 
 = = = = = = = = = = = = = = = = = 

 The anterior sufficient statistics for the simulation and case studies are: $ { \bar{m}}_0 = 0 $, $ { \bar{c}}_0 = 10 ^ 3 $, $ { \bar{a}}_0 = 0.5 $, $ { \bar{b}}_0 = 0.1 $, $ { \mathbf{m}}_0 = 40 { \mathbf{1}}_u$, $ C_0 = 10 ^ 4I_u$ and $ a_0 = 0.5 $ where $ { \mathbf{1}}_u$ is a unit of measurement $ u$-vector and $ I_u$ is the $ u\times{u}$ identity matrix. The statistic $ b_0 $ is defined according to with $ \delta=0.05 $ and $ \epsilon=0.2$. The upper bounds for the excitability parameter quad are $ { \eta_{\max}}= 1.1s_{\tau}$ and $ { \lambda_{\max}}= 14$. In the case study, the upper bound for the scale parameter was reduced to $ { \lambda_{\max}}=7$. 

 Resampling in Algorithm   \[tab: Alg\ ], is perform by systematic sampling on the residuals of particle weight unit [ @Hol06 ]. The number of atom sample and the number of rectangular lattice cells is initially $ N=5000 $ and $ |\mathcal{G}|=30\times{30}$ respectively.
 moge beneficial in terms on resource managgmwnt anv preciaion, card would be needed in automatmng rhe geid refinements, and it is likelj that a yniqne grid would be associabzd wifm eack vistinct firing pattern.

Tha sequential avpdcc of the proposed methodology providqs the ppoortunity for teal-tpmq inrvrtnce that has the potential to prkvide ii-lab assistance during experimentation. In thid framework, an intfrim SMC-MUNG anwoysis could felp in idtncifying the choice of stimulus to apply in ordex to collecj fhf best evideice to distinguish between competong hypotheses, as in Section \[sec:SimStudy\_unver\]. The limitations jf the prasznt SMC-MUNE procedure ti becmme d whuoly onmiie zlgorihhm are the ckmputationao aspects discussed ewgkier and the post-pwosessing stage to correct for potentiallj nefligible estimates of tye expected MU twitch forces. Sjlutions to these outstanding problems would incraase vhd eynlcievxy and accuracy of SMC-MUNE and, hence, the range jr spilication.

Acknowlebgements {#acknowlrdhekgnts.unnumbered}
================

Jhe autkkra would like to thwnk Dr. Shrisrine Thomws ftom the Miami Project to Cuee Paralysis qho provided the dcta analysed in jhis psper and for specialist disdussions.

Addltional dsgail {#sec:APX}
=================

The pfiog suxficient statistics for thq simulatmon aud case rtudoes arq: ${\bar{m}}_0 = 0$, ${\bwr{c}}_0 = 10^3$, ${\bar{a}}_0 = 0.5$, ${\bar{b}}_0 = 0.1$, ${\mathhf{m}}_0 = 40 {\mdthbf{1}}_u$, $C_0 = 10^4L_u$ and $a_0 = 0.5$ where ${\mathbf{1}}_u$ is a uimt $u$-vector anc $H_u$ ps the $u\tnmes{u}$ ldentity matrix. The statistic $b_0$ is dzfined according to witi $\delta=0.05$ and $\qpsilon=0.2$. The u[ker bounds foc the excytabulitt paramdger space are ${\rta_{\max}}= 1.1s_{\tcb}$ and ${\lamvda_{\max}}= 14$. In the casc stuag, the upper bouub dor the scale psraoetqr wes reqgced to ${\lambga_{\mab}}=7$.

Rersmplivg in Algorltho \[tab:Slg\], is performed by vystsmatic sampling on tme residuqls of pwrticle weighys [@Hol06]. The number lf pacticle sampkes and the number of rectangular lattice feljs is initiajly $K=5000$ anq $|\mathcal{G}|=30\cimes{30}$ respectively.


more beneficial in terms of resource management care be needed automating the grid that unique grid would associated with each firing pattern. The sequential aspect of proposed methodology provides the opportunity for real-time inference that has the potential to in-lab assistance during experimentation. In this framework, an interim SMC-MUNE analysis could help identifying choice stimulus apply in order to collect the best evidence to distinguish between competing hypotheses, as in Section The limitations of the present SMC-MUNE procedure to a wholly online algorithm the computational aspects discussed earlier the stage to for negligible of the expected twitch forces. Solutions to these outstanding problems would increase the efficiency and accuracy of SMC-MUNE and, hence, range of {#acknowledgements.unnumbered} ================ authors like thank Dr. Christine the Miami Project to Cure Paralysis data analysed in this paper and for specialist Additional detail ================= The prior sufficient statistics for simulation and case studies are: ${\bar{m}}_0 = 0$, = 10^3$, ${\bar{a}}_0 = 0.5$, ${\bar{b}}_0 = 0.1$, ${\mathbf{m}}_0 = 40 {\mathbf{1}}_u$, $C_0 = 10^4I_u$ = 0.5$ where ${\mathbf{1}}_u$ a unit $u$-vector $I_u$ the identity The statistic is defined according to with $\delta=0.05$ and $\epsilon=0.2$. The upper bounds the excitability parameter space are ${\eta_{\max}}= 1.1s_{\tau}$ and ${\lambda_{\max}}= 14$. case the upper bound the scale parameter was to Resampling in Algorithm \[tab:Alg\], by sampling of weights The number of particle and the number of rectangular cells is initially $N=5000$
 more beneficial in terms of reSource manaGemenT anD prEcIsioN, carE would be needed IN autOmating the grid refinemeNts, anD iT Is liKElY that A unique GRiD WOulD bE aSsoCiATeD with EacH distinCt firing paTteRn.

the sequentiaL AsPect of the pRopOsed methodolOgy ProvidEs The OPportUniTy for Real-tiME inferEnce that hAs THe poteNTial to pROViDe in-Lab assistance duriNG eXPerimentation. IN this fRaMEwORK, an IntErim SMC-MUNe aNalysIS could hELp IN IDenTIfying the choiCe of stimuluS To aPply in OrDer TO colleCt the BeST evIdence to disTingUish betweEn compETing hypOTheses, aS in SecTioN \[seC:SimsTuDy\_UndEr\]. tHe lIMiTatIOns Of the preSeNt sMC-MUnE prOCEDUre tO beCome A wholLy online algorIthM are THe cOmputAtionAl asPeCts diScusseD earlIeR and the post-procEssiNg stage to CorReCt fOr PotenTIally nEglIgiBle estiMates of THe eXpECTEd mU twitch forces. SoluTiONS tO these ouTstandINg PrOBlems wouLd IncReasE THe effIcieNCy And accurAcy of Smc-MuNe and, henCe, The ranGe Of aPplIcatiON.

AckNowledGements {#aCknowLEdgements.unnumBEred}
================

The authorS WoULD lIKe to ThaNk Dr. ChristiNe ThOMas fRom tHE MIamI projeCt to CUrE paRAlysis who provided thE dAta anaLysed In this paper anD for speciaLIST discussIons.

aDdITional detail {#seC:APX}
=================

THe prior sufFIcient stAtistIcs for thE simulatiON And case sTudIes Are: ${\Bar{M}}_0 = 0$, ${\BAr{C}}_0 = 10^3$, ${\bar{a}}_0 = 0.5$, ${\bar{b}}_0 = 0.1$, ${\mathBF{M}}_0 = 40 {\matHbF{1}}_u$, $C_0 = 10^4I_u$ anD $a_0 = 0.5$ wHere ${\matHbf{1}}_U$ is A unIt $u$-VeCtor and $I_u$ Is the $u\tiMeS{u}$ IdEnTitY matrIX. The statIsTic $B_0$ iS deFined ACcordiNg to wIth $\dElTa=0.05$ ANd $\ePsilon=0.2$. THE uPPEr boUnDs For tHe eXcItabiLity PAraMeter spAce are ${\eta_{\Max}}= 1.1S_{\Tau}$ aNd ${\LaMbda_{\max}}= 14$. in the case studY, tHe upper bouNd For The scaLE ParameteR was reduced to ${\lambda_{\max}}=7$.

RESamplinG in algorIthm \[Tab:Alg\], is pErfOrmed bY sySTematiC samplIng on ThE reSIDuals OF PaRtiClE weights [@HoL06]. tHe nUmber Of PartIcle samPles and the number of REctAngular latticE ceLls iS INiTiaLLy $n=5000$ And $|\MaTHcaL{g}|=30\Times{30}$ respectiveLy.


 more beneficial in termsof resourc e man age men tandprec ision, care wo u ld b e needed in automating  thegr i d re f in ement s, andi ti s  li ke ly  th at auniqu e g rid wou ld be asso cia te d with eachd is tinct firi ngpattern.

Th e s equent ia l a s pectofthe p ropose d  metho dology pr ov i des th e  opport u n it y fo r real-time infer e nc e  that has thepotent ia l  t o pro vid e in-lab a ss istan c e durin g  e x p e rim e ntation. In t his framewo r k,an int er imS MC-MUN E ana ly s iscould helpin i dentifyin g thec hoice o f  stimul us toapp lyin o r de rtoco l lec t  t heb est  evidenc eto  dist ingu i s h betw een  com petin g hypotheses,  as  inS ect ion \ [sec: SimS tu dy\_u nder\] . The  l imitations of t he p resent SM C-M UN E p ro cedur e  to be com e a  wholly  online alg or i t h mare the computatio na l as pects di scusse d  e ar l ier andth e p ost- p r ocess ings ta ge to co rrectf or  p otentia ll y negl ig ibl e e stima t es o f theexpected  MU t w itch forces. S o lutions to th e se o ut s tand ing  problems w ould incr ease th e e f ficie ncy a nd ac c uracy of SMC-MUNE a nd , henc e, th e range of ap plication. A cknowled geme n ts {#acknowledgem ents. unnumbered } 
======= ===== ====

Th e authors w ould lik e t o t han k D r .  C hristine Thom a s  fro mthe Mia miProject  to  Cu rePar al ysis whoprovided  t he  d at a a nalys e d in thi spap er  an d for specia listdisc us si o ns. 

Addit i on a l  det ai l{#se c:A PX }
=== ==== = === ======
The prio r s u ffic ie nt  statis tics for thesi mulation a nd  ca se stu d i es are:${\bar{m}}_0 = 0$, ${\b a r{c}}_0  =10^3$ , ${ \bar{a}}_ 0 =  0.5$,  ${ \ bar{b} }_0 =0.1$,  $ {\m a t hbf{m } } _0  =40  {\mathbf{ 1 } }_u $, $C _0  = 1 0^4I_u$  and $a_0 = 0.5$ w h ere  ${\mathbf{1} }_u $ is a  u nit $u $ -ve ct o r a n d  $I_u$ is the $ u\times{u} $i de ntity matr i x.Th e stati stic $b _0$ i s  define d accordi ng to wit h$\de l t a=0 .05$ and $ \epsilon =0.2$. Th e  uppe r  b ounds  fo r theex cit abili ty par a met er sp ace ar e${\eta _{\ma x} }= 1.1s_ {\tau}$ and ${\lambda_{ \max}} = 14$ . I n the cas e s t udy , the upp er b ound for t hesca le pa ram e ter w as r e du ced to ${ \lam b da_{\max} } =7 $.R es ampling inA l g ori thm \ [ta b :Alg\] , is  performed by sys t ematic samplin g on t heres i dual sof particle we igh ts [ @Hol06].  T he number o f partic le sampl es and  the n umber o f re c tangul ar l att ice cells  is  i n itially  $ N= 5 000$ a nd $ |\ mathca l{G}|= 3 0\ti m e s{30}$ respectiv ely.
 more_beneficial in_terms of resource management_and precision,_care_would be_needed_in automating the_grid refinements, and_it is likely that_a unique grid_would_be associated with each distinct firing pattern.

The sequential aspect of the proposed methodology provides_the_opportunity for_real-time_inference_that has the potential to_provide in-lab assistance during experimentation._In this_framework, an interim SMC-MUNE analysis could help in_identifying_the choice of_stimulus to apply in order to collect the best_evidence to distinguish between competing hypotheses,_as in Section \[sec:SimStudy\_under\]._The_limitations_of the present SMC-MUNE_procedure to become a wholly online_algorithm are the computational aspects discussed_earlier and the post-processing stage to correct_for potentially negligible estimates of the_expected MU twitch forces. Solutions_to these_outstanding problems would increase the_efficiency and accuracy_of SMC-MUNE_and, hence, the_range of application.

Acknowledgements {#acknowledgements.unnumbered}
================

The authors would_like to thank_Dr. Christine Thomas from the Miami_Project_to Cure Paralysis_who_provided_the data_analysed in this_paper_and for_specialist_discussions.

Additional detail {#sec:APX}
=================

The prior sufficient statistics_for_the simulation and case studies are: ${\bar{m}}_0_= 0$, ${\bar{c}}_0 =_10^3$,_${\bar{a}}_0 = 0.5$, ${\bar{b}}_0_= 0.1$, ${\mathbf{m}}_0 = 40_{\mathbf{1}}_u$, $C_0 = 10^4I_u$ and $a_0_= 0.5$_where ${\mathbf{1}}_u$_is a unit $u$-vector and $I_u$ is the $u\times{u}$ identity matrix._The statistic $b_0$ is defined according_to with $\delta=0.05$ and_$\epsilon=0.2$. The_upper_bounds for the_excitability_parameter space_are ${\eta_{\max}}= 1.1s_{\tau}$ and ${\lambda_{\max}}= 14$. In_the case_study, the upper bound for the_scale parameter was reduced_to_${\lambda_{\max}}=7$.

Resampling in Algorithm \[tab:Alg\], is performed by_systematic sampling on the residuals of_particle weights [@Hol06]. The number_of_particle_samples and the number of_rectangular lattice cells is initially $N=5000$_and $|\mathcal{G}|=30\times{30}$ respectively.

