analysis(matlab): rate + count outputs, variation batch 3

Co-Authored-By: Claude Opus 4.8 <noreply@anthropic.com>
2026-07-24 02:32:52 -04:00
parent 6ba21a1a35
commit 0e16690c8f
72 changed files with 1875 additions and 897 deletions
@@ -0,0 +1,68 @@
+==============================================================================
+VARIATION: prev_f_full   [metric: success RATE (success/attempts)]
+==============================================================================
+model: behavior ~ stim + day + stim:day + (1|rat)   (behavior = success RATE (success/attempts))
+day = training day within window (0 = first analyzed day)
+treatment (stim=1): b2_f
+control   (stim=0): a2_f
+N = 24 rats, 231 sessions   raw day coverage: treat 0..9, control 0..9
+(equal day coverage over this window)
+
+==============================================================================
+FULL MODEL SUMMARY -- fitlme
+==============================================================================
+
+Linear mixed-effects model fit by ML
+
+Model information:
+    Number of observations             231
+    Fixed effects coefficients           4
+    Random effects coefficients         24
+    Covariance parameters                2
+
+Formula:
+    behavior ~ 1 + day*stim + (1 | rat)
+
+Model fit statistics:
+    AIC        BIC        LogLikelihood    Deviance
+    -441.62    -420.96    226.81           -453.62 
+
+Fixed effects coefficients (95% CIs):
+    Name                   Estimate    SE           tStat     DF     pValue    
+    {'(Intercept)'}         0.25306     0.024155    10.477    227    3.3213e-21
+    {'day'        }        0.017111    0.0027396     6.246    227    2.0536e-09
+    {'stim'       }        0.040889     0.034129    1.1981    227       0.23214
+    {'day:stim'   }        0.006292    0.0038198    1.6472    227        0.1009
+
+
+    Lower         Upper   
+       0.20547     0.30066
+      0.011713     0.02251
+     -0.026361     0.10814
+    -0.0012348    0.013819
+
+Random effects covariance parameters (95% CIs):
+Group: rat (24 Levels)
+    Name1                  Name2                  Type           Estimate
+    {'(Intercept)'}        {'(Intercept)'}        {'std'}        0.06824 
+
+
+    Lower       Upper   
+    0.049364    0.094332
+
+Group: Error
+    Name               Estimate    Lower      Upper  
+    {'Res Std'}        0.081522    0.07404    0.08976
+
+
+effect                     t(df) / F(df1)     p (resid)    Satterthwaite: p (df)
+----------------------------------------------------------------------------
+stim x day (interaction)   t(227)=  1.65 F(1)=  2.713  p=0.1009    p=0.101 (df=209)
+day (learning)             t(227)=  6.25 F(1)= 39.012  p=2.054e-09    p=2.316e-09 (df=209)
+stim (main, window start)  t(227)=  1.20 F(1)=  1.435  p=0.2321    p=0.2378 (df=41)
+interaction 95% CI: [-0.001235, +0.01382]
+HONEST LME (per-animal random slope, day|rat): interaction F(1,22.9)=1.78, p=0.1956
+  (Satterthwaite DF ~= residual on this random-intercept model; the random-slope
+   model above is the honest learning-rate test -- DF collapses toward the animal count.)
+INTERPRETATION: stim x day interaction n.s. -- slopes parallel (no differential learning rate) (p=0.1009, slope diff=+0.006292)
+Paper (N=24, count): interaction t(227)=2.68, F(1)=7.12, p=0.008.