From 273b7a04ae43e525e56dfc709daf40136985a1dc Mon Sep 17 00:00:00 2001 From: Jan Benda Date: Fri, 16 May 2025 09:22:34 +0200 Subject: [PATCH] updated regimes figure --- .../2013-01-08-aa-invivo-1-contrastpeaks.csv | 201 ++++++++++++++++++ data/punitmodels.csv | 40 ++++ regimes.py | 34 ++- 3 files changed, 257 insertions(+), 18 deletions(-) create mode 100644 data/cells/2013-01-08-aa-invivo-1-contrastpeaks.csv create mode 100644 data/punitmodels.csv diff --git a/data/cells/2013-01-08-aa-invivo-1-contrastpeaks.csv b/data/cells/2013-01-08-aa-invivo-1-contrastpeaks.csv new file mode 100644 index 0000000..ac14f67 --- /dev/null +++ b/data/cells/2013-01-08-aa-invivo-1-contrastpeaks.csv @@ -0,0 +1,201 @@ +contrast/%,f1/Hz,f2/Hz,f1+f2/Hz,f2-f1/Hz +0.0,443.949,14910.2,3508.46,1874.87 +0.1,452.806,15246.7,3540.24,1940.44 +0.2,544.882,16645.1,3485.35,1912.72 +0.3,682.7,18324,3638.24,1899.45 +0.4,931.933,21777.1,3476.33,1901.9 +0.5,1319.9,26419.3,3603.34,1895.2 +0.6,1677.17,30113.2,3665.13,1903.58 +0.7,2051.77,37017.5,3581.3,1952.11 +0.8,2600.35,45354.6,3528.9,1932.91 +0.9,3179.39,51408.8,3586.04,1955.77 +1.0,3793.06,60046.4,3601.14,2021.38 +1.1,4560.06,69456.6,3659.88,1950.5 +1.2,5302.08,77393.3,3798.62,1954.91 +1.3,6057.7,89662,3838.78,2022.77 +1.4,6953.86,101945,4239.47,2093.94 +1.5,8036.16,119050,4385.71,2137.91 +1.6,9050.76,127272,4817.12,2125.21 +1.7,10056.1,144998,4976.52,2219.78 +1.8,11257.6,155870,5452.2,2266.68 +1.9,12394,175322,5411.3,2494.43 +2.0,13694.5,190274,6009.05,2716.64 +2.1,15194.4,202697,6881.44,2725.46 +2.2,16473.9,219269,7550.98,2874.82 +2.3,18131.1,240023,8429.76,3219.07 +2.4,19640,250636,9408.26,3396.35 +2.5,21065.9,274181,10328.6,3897.92 +2.6,22686.6,293049,11103.3,3901.32 +2.7,24053.4,312801,12373.8,4192.14 +2.8,25878.3,328418,13231.2,4699.77 +2.9,27746.2,350318,14787.7,5147.72 +3.0,29654.9,366179,15110.4,5370.47 +3.1,31801,390730,17722.5,5969.82 +3.2,33077.6,406406,19629.8,6332.34 +3.3,35290.4,431837,21897.7,6686.9 +3.4,37262.7,451934,22780.7,7177.93 +3.5,39426.2,468724,25578.5,7763.61 +3.6,41251.2,496171,27566.7,8252.81 +3.7,43538.5,512614,30030.8,8710.28 +3.8,45977.5,536852,32778.7,9579.56 +3.9,47802.2,551778,35418.6,10304.5 +4.0,50303.9,576449,37937,10937.7 +4.1,52556.7,591165,40125.9,10994.5 +4.2,54841.1,609369,42610.9,11711 +4.3,57061.6,630749,45146.3,12554.9 +4.4,59605,661253,49741.3,13176.4 +4.5,62300.8,675542,53359.6,13821.5 +4.6,64376.7,694640,55907.6,14282.3 +4.7,66649.1,715194,58066.9,15303.1 +4.8,69837.7,733166,62037.4,16068 +4.9,72221.9,756223,67006.4,16904.1 +5.0,75591.3,771586,69242.1,16966.2 +5.1,77824.1,791123,73400.7,18082.1 +5.2,80207.6,818047,79116.1,19107.5 +5.3,83177.5,832371,81536.3,19268 +5.4,85704.4,842162,85368.6,20081.5 +5.5,89069,861376,87832.8,20588.7 +5.6,91662,884650,94134.7,21494.2 +5.7,94680.7,894889,97977.7,21838.4 +5.8,97921.5,904816,101913,22191.6 +5.9,100948,931136,107127,23174.5 +6.0,103041,942873,110671,23834 +6.1,107384,960488,114345,24534.3 +6.2,110347,970176,118877,24613.3 +6.3,113582,981384,122701,25342.6 +6.4,116857,1.00145e+06,127355,25972.4 +6.5,120414,1.01405e+06,131264,26427.5 +6.6,123657,1.02774e+06,135772,26947.5 +6.7,126447,1.04146e+06,139290,27661.5 +6.8,130079,1.05614e+06,145632,28130 +6.9,133465,1.07134e+06,149271,27757.2 +7.0,137953,1.07971e+06,154455,27512.5 +7.1,141411,1.08797e+06,157678,28119.4 +7.2,144977,1.10694e+06,161216,28287.1 +7.3,148118,1.11439e+06,163634,28896 +7.4,152104,1.13169e+06,168685,28975.5 +7.5,156098,1.14385e+06,172825,29136.5 +7.6,159822,1.14756e+06,178335,29559.4 +7.7,164181,1.15736e+06,179959,29123.3 +7.8,168040,1.16624e+06,182076,28918.7 +7.9,172400,1.18103e+06,186428,29819.6 +8.0,176626,1.18451e+06,188667,29197.7 +8.1,181279,1.19173e+06,193500,28803.5 +8.2,184634,1.21214e+06,198357,30139.4 +8.3,190200,1.21182e+06,198934,29435.5 +8.4,194296,1.22856e+06,203986,30059.3 +8.5,198112,1.23816e+06,204816,29193.9 +8.6,203342,1.24264e+06,210840,29324 +8.7,209082,1.25161e+06,210281,29658 +8.8,212712,1.26118e+06,215732,29649.2 +8.9,216743,1.26224e+06,217129,29925.7 +9.0,221897,1.27861e+06,220978,30459.9 +9.1,227242,1.28384e+06,219459,30176.5 +9.2,231885,1.29032e+06,224311,30639 +9.3,235807,1.30028e+06,226280,31525.1 +9.4,242623,1.30407e+06,226463,31686.9 +9.5,247678,1.30729e+06,228923,32456.5 +9.6,253102,1.32134e+06,231688,32084.6 +9.7,258516,1.3264e+06,231052,32475.5 +9.8,264602,1.33153e+06,233155,33147 +9.9,269442,1.33847e+06,232817,34490 +10.1,275414,1.34956e+06,235268,34819.2 +10.2,282815,1.35012e+06,233680,35683.7 +10.3,286235,1.35827e+06,238995,36278.6 +10.4,292763,1.37063e+06,237883,37956.8 +10.5,298276,1.38364e+06,239149,40137.7 +10.6,305261,1.38599e+06,241933,40634.2 +10.7,310836,1.39169e+06,242297,41365.2 +10.8,316918,1.40115e+06,242502,43831.3 +10.9,323955,1.41268e+06,245134,44994.7 +11.0,330517,1.41557e+06,243202,45975.7 +11.1,336366,1.42822e+06,245586,49371.1 +11.2,342493,1.43178e+06,245122,49635.5 +11.3,350228,1.44189e+06,241073,52991.5 +11.4,356391,1.45024e+06,245472,54352.2 +11.5,362308,1.45751e+06,242601,57560.2 +11.6,370326,1.46238e+06,242903,59954.7 +11.7,378145,1.47136e+06,244009,62707.7 +11.8,384430,1.47565e+06,247069,65172.2 +11.9,390731,1.48948e+06,245279,68672.9 +12.0,398343,1.49927e+06,244648,72524 +12.1,404640,1.50622e+06,245258,74504.3 +12.2,412943,1.51634e+06,246579,79453.6 +12.3,420410,1.52766e+06,244759,83509.2 +12.4,426579,1.53338e+06,246532,86244.3 +12.5,433949,1.54249e+06,244684,90243 +12.6,441362,1.557e+06,245121,95798.5 +12.7,449260,1.56317e+06,244225,99349.7 +12.8,455866,1.57281e+06,244569,103224 +12.9,465165,1.58515e+06,245035,108078 +13.0,473881,1.59317e+06,245585,112656 +13.1,477951,1.60776e+06,247639,116865 +13.2,485697,1.61504e+06,247766,120404 +13.3,495065,1.62701e+06,246306,126477 +13.4,503346,1.6334e+06,244923,133331 +13.5,510810,1.6479e+06,247263,138423 +13.6,519136,1.65853e+06,246636,143831 +13.7,527044,1.66855e+06,247326,149587 +13.8,534907,1.67922e+06,247929,155337 +13.9,543905,1.68378e+06,247629,159740 +14.0,551637,1.70402e+06,249758,167789 +14.1,559545,1.71448e+06,250561,172134 +14.2,565747,1.7241e+06,251257,179693 +14.3,576568,1.7371e+06,250144,186764 +14.4,582864,1.75042e+06,254045,191566 +14.5,592815,1.76097e+06,254018,197652 +14.6,600921,1.77373e+06,251636,206065 +14.7,608614,1.78306e+06,252549,212579 +14.8,617621,1.79495e+06,252807,219745 +14.9,625915,1.80945e+06,255435,226765 +15.0,633795,1.82008e+06,256247,232830 +15.1,641274,1.83368e+06,257262,239334 +15.2,651776,1.84672e+06,257499,248376 +15.3,659979,1.85907e+06,259700,252787 +15.4,668210,1.87029e+06,259974,263609 +15.5,676716,1.88244e+06,260989,269961 +15.6,686524,1.89579e+06,264744,278111 +15.7,693030,1.9076e+06,264596,286325 +15.8,703286,1.91856e+06,266108,293875 +15.9,710356,1.9327e+06,267830,301787 +16.0,720078,1.94702e+06,269847,308048 +16.1,729912,1.95738e+06,270938,318588 +16.2,736469,1.97206e+06,272360,325156 +16.3,746812,1.98181e+06,276345,331016 +16.4,753788,1.99916e+06,277476,340580 +16.5,762482,2.0105e+06,278432,351356 +16.6,770243,2.02553e+06,281632,357980 +16.7,782091,2.03748e+06,282029,367124 +16.8,790682,2.05258e+06,284762,374996 +16.9,799104,2.06487e+06,288886,385716 +17.0,807696,2.07526e+06,291629,392316 +17.1,816360,2.08766e+06,293455,400056 +17.2,823582,2.10123e+06,294600,410871 +17.3,833400,2.11464e+06,296546,420457 +17.4,843474,2.1217e+06,298297,427367 +17.5,851916,2.13939e+06,301756,438379 +17.6,860716,2.15484e+06,304479,445897 +17.7,868406,2.16772e+06,308372,456795 +17.8,879438,2.17989e+06,309082,464561 +17.9,888231,2.19557e+06,312804,475801 +18.0,897437,2.20749e+06,315443,486268 +18.1,906100,2.21889e+06,317964,492645 +18.2,914157,2.23204e+06,320215,503433 +18.3,923485,2.24481e+06,324801,512714 +18.4,933308,2.25814e+06,328911,522824 +18.5,942521,2.27501e+06,332008,533261 +18.6,950541,2.28515e+06,334427,543097 +18.7,959812,2.29679e+06,337361,550246 +18.8,968856,2.31509e+06,341532,563251 +18.9,978769,2.32288e+06,343096,574810 +19.0,988305,2.33779e+06,347046,583050 +19.1,996233,2.35491e+06,350982,589428 +19.2,1.00677e+06,2.3664e+06,352951,603616 +19.3,1.01404e+06,2.38128e+06,357384,614489 +19.4,1.02433e+06,2.38848e+06,363376,622459 +19.5,1.03547e+06,2.40559e+06,365173,635075 +19.6,1.04181e+06,2.41607e+06,367379,645082 +19.7,1.05217e+06,2.42564e+06,372387,658039 +19.8,1.06182e+06,2.44668e+06,378595,666422 +19.9,1.07065e+06,2.462e+06,378154,679152 +20.0,1.07888e+06,2.46916e+06,383906,691408 diff --git a/data/punitmodels.csv b/data/punitmodels.csv new file mode 100644 index 0000000..1931e21 --- /dev/null +++ b/data/punitmodels.csv @@ -0,0 +1,40 @@ +cell,EODf,a_zero,delta_a,dend_tau,input_scaling,mem_tau,noise_strength,ref_period,deltat,tau_a,threshold,v_base,v_offset,v_zero +2010-11-08-al-invivo-1,744.66,9.450855200303527,0.060498440079361865,0.0007742334994649853,31.363843698084207,0.0017257848281706928,0.012409100812593247,0.0010273077926126455,5e-05,0.10223865531575653,1,0,-0.390625,0 +2011-10-25-ad-invivo-1,760.50,60.70561732896827,0.15752867579184138,0.004075295861125333,301.37869129345717,0.0029012439225658054,0.02968073683628557,0.0007434544009651972,5e-05,0.11607744742342971,1,0,-34.375,0 +2012-04-20-ad-invivo-1,811.57,61.15542998217432,0.18099114092280794,0.006549954587111268,317.4201379309943,0.0014460502988963755,0.03508490455171905,0.001007048276989675,5e-05,0.1718068006049748,1,0,-39.84375,0 +2012-04-20-ak-invivo-1,826.07,142.73039605053253,0.34408674592921096,0.0042905816022802655,373.7425313488243,0.0017358116367672107,0.016121166740672717,0.0009347144390364768,5e-05,0.2949535420035986,1,0,21.875,0 +2012-06-27-ah-invivo-1,752.07,28.875515577269294,0.20525468484541842,0.010612629869418248,554.4436702347741,0.0017107204697516367,0.0271472345939834,0.0012084459280246903,5e-05,0.11784186849482106,1,0,-148.4375,0 +2012-06-27-an-invivo-1,786.29,2.7860389486251593,0.029391448491264804,0.003985634602834923,26.682508784285027,0.001376526707110581,0.003604805695367027,0.0008998947426172306,5e-05,0.13036159256315194,1,0,-4.8828125,0 +2012-07-03-ak-invivo-1,928.45,1.1337603254658657,0.009636823781567081,0.0011835211027475872,10.551593612226277,0.0013790127193975233,0.0013081636418144473,0.00011600868359679133,5e-05,0.09604613888260315,1,0,-1.318359375,0 +2012-07-12-ag-invivo-1,744.95,3.695738241953972,0.03623619132514362,0.0018084515871059888,14.377436485984251,0.0011056964083457424,0.0041125372009629575,0.0011502269812458804,5e-05,0.11156737983241918,1,0,-0.09765625,0 +2012-07-12-ap-invivo-1,772.92,14.804318060119332,0.08421508191399865,0.0013221092024793385,44.13667585034951,0.0014555506316246737,0.016159465360721062,0.00097403156684423,5e-05,0.0714312935754097,1,0,1.171875,0 +2012-12-13-af-invivo-1,673.56,23.495808141911642,0.1324559007885552,0.010031950746925707,427.7951451118343,0.003903660950040095,0.018153856801854964,0.0009569332571649322,5e-05,0.13755336036983307,1,0,-111.71875,0 +2012-12-13-ah-invivo-1,664.70,28.779813371047524,0.16049265233555093,0.01314225632311085,365.45335373877225,0.0011890432540278138,0.014044744413304281,0.0014219606674122764,5e-05,0.1350142793970573,1,0,-87.5,0 +2012-12-13-an-invivo-1,657.91,4.458558193012667,0.030549781596917867,0.001378687228378128,35.834438055915825,0.0029913612464463734,0.005852518402871827,0.001256771411437993,5e-05,0.024870623174021396,1,0,-6.25,0 +2012-12-13-ao-invivo-1,657.82,3.2813550587302758,0.02199424769500702,0.001372234851698214,16.707319333864564,0.00221210217479549,0.007347109401659321,0.0008575662334429958,5e-05,0.05671867775854021,1,0,-1.26953125,0 +2012-12-20-aa-invivo-1,668.32,4.195725270757472,0.03494694984590694,0.004944939140795805,61.926544628777336,0.001774255774424637,0.006397082203748482,0.0013298764071622482,5e-05,0.06259061000559793,1,0,-14.84375,0 +2012-12-20-ab-invivo-1,738.71,9.280347979087807,0.02399411189063231,0.00120664709274584,46.57993094588944,0.001345591696621093,0.0076138573270267875,0.0011385589203656609,5e-05,0.048005025093914855,1,0,-4.78515625,0 +2012-12-20-ac-invivo-1,744.95,7.430387927489267,0.03459479087457316,0.0025268601414655655,43.133868817764956,0.0015547575080636315,0.007099446913853285,0.0008625278566343026,5e-05,0.06974558082044995,1,0,-5.6640625,0 +2012-12-20-ad-invivo-1,759.82,23.049443800356883,0.07505377496781868,0.004539003389588492,124.17804604983468,0.00106463192556606,0.010618483391053847,0.0010943556029644886,5e-05,0.0931679653142917,1,0,-16.2109375,0 +2012-12-20-ae-invivo-1,763.79,28.80721245217218,0.0698601314329484,0.005371563846161885,190.09285659771518,0.001645315861662125,0.011943115636543555,0.0009171798026672154,5e-05,0.06115045426599336,1,0,-31.8359375,0 +2012-12-21-ai-invivo-1,787.12,36.7842414315693,0.12763184695279117,0.0031488535186190362,291.18540792514136,0.0020959669717419043,0.03321451893799585,0.0012049087410326699,5e-05,0.12743644124243197,1,0,-54.6875,0 +2012-12-21-ak-invivo-1,796.83,2.4867707329904993,0.014127034477017693,0.002861628332432993,31.510428742268093,0.0003974311599786272,0.003398627675991102,0.0006089766381869961,5e-05,0.013593335235228779,1,0,-7.71484375,0 +2012-12-21-am-invivo-1,806.15,4.716159805342061,0.03667764979320955,0.004999856382483749,85.64267738935817,0.00241012573550433,0.011026662170574162,0.0011255575558147763,5e-05,0.0544681581478567,1,0,-21.484375,0 +2012-12-21-an-invivo-1,812.70,2.6579691222796833,0.012698001151842848,0.002929056579331126,30.47330327087832,0.0006812296678529029,0.006737763075104184,0.0010815182714474752,5e-05,0.029800954358211888,1,0,-6.54296875,0 +2013-01-08-aa-invivo-1,800.63,1.1104651443148155,0.008602584313899703,0.0011824418111601317,4.462416862986384,0.0011980061424283523,0.0016372239238665433,0.00038313415155155986,5e-05,0.03752455828256267,1,0,0.5859375,0 +2013-01-08-ab-invivo-1,800.25,54.67466209697357,0.23802565736953032,0.00314795259301606,401.47466098307893,0.003233481691970421,0.0435481413010326,0.0004075044135951958,5e-05,0.10601409774748785,1,0,-75.78125,0 +2013-02-21-ag-invivo-1,658.70,12.638812661345776,0.10167490092250338,0.003536770266850944,104.74966687228975,0.0016256371684299914,0.014138536979973691,0.0012702670089049608,5e-05,0.07815503514848808,1,0,-20.3125,0 +2013-04-10-aa-invivo-1,623.77,7.2029389446183565,0.052641907830485446,0.002224898046101902,62.65798249994464,0.002873209337269388,0.01811955430517219,0.0015625668737142603,5e-05,0.06354491016383691,1,0,-12.5,0 +2014-01-16-ak-invivo-1,803.26,2.3120229511581782,0.014165489782626235,0.007112977431183025,24.30632120240236,0.0002851096309430286,0.0028600590320415394,0.0011963117586515453,5e-05,0.06057228972437009,1,0,-4.98046875,0 +2014-12-11-aa-invivo-1,651.29,110.31713694365598,1.522266574846523,0.02274369798634087,499.44685404447773,0.0010908984221126568,0.04065807318024823,0.0010262164681389477,5e-05,0.8644462875159566,1,0,-50.0,0 +2017-07-18-ah-invivo-1,816.18,4.492478880802049,0.030680049206572385,0.0005510176207879538,11.010435284946723,0.0012660814616981633,0.01178076926186549,0.00010480712226888011,5e-05,0.05347128935359713,1,0,1.171875,0 +2017-07-18-ai-invivo-1,817.53,4.836173387076376,0.045960408902420334,0.0005713395854796994,19.082872790172893,0.0017648069889998111,0.024310257773158105,0.00037992426067294776,5e-05,0.021943818745769235,1,0,-2.5390625,0 +2017-07-18-aj-invivo-1,818.98,62.59548468438461,0.26688797342324944,0.0063976464842525225,344.9549845381489,0.0006269389203256706,0.03293495593748591,0.00048458937953982373,5e-05,0.08847603129712459,1,0,-48.4375,0 +2018-05-08-aa-invivo-1,643.65,23.90042739559335,0.1798294135622658,0.0014978457325112936,196.79744698301027,0.008998581511662336,0.1931709060535164,0.0012206757789855745,5e-05,0.08342443880014652,1,0,-38.28125,0 +2018-05-08-ab-invivo-1,646.62,21.143340239280292,0.17918314011340677,0.0021110702241677664,74.44012660456099,0.002104111838094462,0.031480734087501225,0.000841793287555979,5e-05,0.1649339650514045,1,0,-3.125,0 +2018-05-08-ac-invivo-1,655.16,35.32000871480389,0.3365254255760804,0.006773219307007808,139.23166556708537,0.0004350491142382615,0.028336890851257814,0.0013549977866199972,5e-05,0.1531910034411788,1,0,-12.5,0 +2018-05-08-ad-invivo-1,655.66,10.56270995846859,0.04501376693597617,0.0013544139474406026,32.871669695970056,0.001129016694314504,0.01343806066427332,0.0007800007176441471,5e-05,0.08300605346748607,1,0,-0.09765625,0 +2018-05-08-ae-invivo-1,649.48,23.236002441272873,0.1649467891961967,0.003929215662714291,139.62843570490134,0.0014895499625897,0.020705895621135995,0.0013078805846238773,5e-05,0.12368546391523849,1,0,-21.09375,0 +2018-05-08-af-invivo-1,649.93,50.67440031212505,0.15942545846977396,0.008870877837080987,266.86501353697287,0.0011885343507996296,0.01831310334923414,0.0011045939825865964,5e-05,0.061964902843727104,1,0,-35.15625,0 +2018-05-08-ai-invivo-1,653.62,17.97588058957365,0.12223368623652905,0.0022403560787985152,58.600496136896105,0.00113939954736546,0.02035182654993045,0.0014268613608903935,5e-05,0.07235575003635941,1,0,-1.5625,0 +2018-06-25-ad-invivo-1,840.79,62.49268603283534,0.23947474463722357,0.001883927876295424,286.0953844311519,0.002389671017958955,0.0698177647604567,0.0010856424997856035,5e-05,0.09634714314913008,1,0,-29.6875,0 diff --git a/regimes.py b/regimes.py index d0b8ad9..1a1d4c1 100644 --- a/regimes.py +++ b/regimes.py @@ -1,12 +1,16 @@ -import os import numpy as np -from scipy.stats import linregress import matplotlib.pyplot as plt +from pathlib import Path +from scipy.stats import linregress from numba import jit from thunderlab.tabledata import TableData from plotstyle import plot_style, lighter, darker +data_path = Path('data') +cell_path = data_path / 'cells' + + def load_models(file): """ Load model parameter from csv file. @@ -21,7 +25,7 @@ def load_models(file): For each cell a dictionary with model parameters. """ parameters = [] - with open(file, 'r') as file: + with file.open('r') as file: header_line = file.readline() header_parts = header_line.strip().split(",") keys = header_parts @@ -192,9 +196,12 @@ def peak_ampl(freqs, psd, f): def compute_peaks(name, cell, alpha_max, beatf1, beatf2, nfft, trials): - file_name = f'{name}-contrastpeaks.csv' - if os.path.exists(file_name): - data = TableData(file_name) + data_file = cell_path / f'{name}-contrastpeaks.csv' + data = TableData(data_file) + return data + """ + if data_file.exists(): + data = TableData(data_file) return data dt = 0.0001 tmax = nfft*dt @@ -217,9 +224,9 @@ def compute_peaks(name, cell, alpha_max, beatf1, beatf2, nfft, trials): data.append('f2', 'Hz', '%g', ampl_f2) data.append('f1+f2', 'Hz', '%g', ampl_sum) data.append('f2-f1', 'Hz', '%g', ampl_diff) - data.write(file_name) + data.write(data_file) return data - + """ def amplitude(power): power -= power[0] @@ -281,23 +288,14 @@ def plot_peaks(ax, s, data, alphas): if __name__ == '__main__': - parameters = load_models('models.csv') + parameters = load_models(data_path / 'punitmodels.csv') cell_name = '2013-01-08-aa-invivo-1' # 138Hz, CV=0.26: perfect! beatf1 = 40 beatf2 = 138 # cell_name = '2012-07-03-ak-invivo-1' # 128Hz, CV=0.24 # cell_name = '2018-05-08-ae-invivo-1' # 142Hz, CV=0.48 - - """ - parameters = load_models('models_big_fit_d_right.csv') - cell_name = '2013-01-08-aa-invivo-1' # 131Hz, CV=0.04: wrong! - beatf1 = 30 - beatf2 = 132 - """ cell = cell_parameters(parameters, cell_name) - for k in cell: - print(k, cell[k]) s = plot_style() s.lwmid = 1.0