Goose/Moon Sensors/Moon Temp Results

Moon Temp Statistical Analysis — 2026-07-31

2026-07-31_Moon-Temp_003

H₀ (null): The rate of change of the difference signal (ADC2 − reference) does not change at the shadow event — any observed trend differences are explained by sensor noise and environmental drift alone.

H₁ (alternate): The rate of change differs after shadow onset — ADC2 progressively diverges from the moonlit reference channels (ADC0, ADC1) in a sustained, accumulating way.

Analysis approach: Due to sensor thermal mass acting as a thermal sink, the moonlight effect is expected to manifest as a gradual slope change in the difference signal, not an instantaneous step. The analysis looks for a change in trend rate, with possible lag relative to the recorded shadow event time. Note: higher ADC counts = cooler; lower ADC counts = warmer.

Parameter	Value
Experiment	2026-07-31_Moon-Temp_003
Date	2026-07-31
Shadow event (ADC2 enters full shadow)	01:11 MDT
Sunrise (analysis cutoff −1 hr)	06:04 MDT
GoPro start	22:03 MDT

import pandas as pd
import numpy as np
import scipy.stats as stats
import statsmodels.api as sm
import plotly.graph_objects as go
import plotly.io as pio; pio.renderers.default = 'notebook'
from plotly.subplots import make_subplots
from pathlib import Path
from datetime import datetime, timedelta, date

# ── rooster.ninja colour palette ─────────────────────────────────────────────
BG       = '#0d1117'
PLOT_BG  = '#161b22'
GRID     = '#30363d'
TEXT     = '#e6edf3'
MUTED    = '#8b949e'
ACCENT   = '#58a6ff'
C0       = '#5b9cf6'   # ADC0 — blue
C1       = '#f6a542'   # ADC1 — amber
C2       = '#63d68c'   # ADC2 — green (shadow channel)
C_SHADOW = '#e05c5c'   # shadow event line — red
C_SUN    = '#ffd166'   # sunrise line — gold
FONT_FAM = "-apple-system, BlinkMacSystemFont, 'Segoe UI', sans-serif"

def rooster_layout(**kw):
    """Base Plotly layout matching rooster.ninja theme."""
    base = dict(
        paper_bgcolor=BG, plot_bgcolor=PLOT_BG,
        font=dict(color=TEXT, family=FONT_FAM),
        xaxis=dict(gridcolor=GRID, linecolor=GRID, zerolinecolor=GRID),
        yaxis=dict(gridcolor=GRID, linecolor=GRID, zerolinecolor=GRID),
        legend=dict(bgcolor='rgba(0,0,0,0)', bordercolor=GRID, borderwidth=1),
        margin=dict(l=60, r=40, t=60, b=60),
    )
    base.update(kw)
    return base

# ── experiment parameters (pre-filled from template) ─────────────────────────
EXPERIMENT_DIR  = Path('.')
experiment_id   = '2026-07-31_Moon-Temp_003'
LOG_DATE        = '2026-07-31'     # YYYY-MM-DD (morning / end date)
SHADOW_TIME_STR = '01:11'    # HH:MM:SS MDT
SUNRISE_STR     = '06:04'        # HH:MM MDT
gopro_start_str = '22:03'         # HH:MM:SS MDT (prev evening)

LOG_DATE_PREV = (date.fromisoformat(LOG_DATE) - timedelta(days=1)).isoformat()

def parse_local(date_str, time_str):
    t = time_str.replace(' MDT','').replace(' MST','').strip()
    fmt = '%Y-%m-%d %H:%M:%S' if t.count(':') == 2 else '%Y-%m-%d %H:%M'
    return datetime.strptime(f'{date_str} {t}', fmt)

gopro_start = parse_local(LOG_DATE_PREV, gopro_start_str)
shadow_dt   = parse_local(LOG_DATE_PREV, SHADOW_TIME_STR)
if shadow_dt < gopro_start:
    shadow_dt += timedelta(days=1)  # shadow event logged after midnight
sunrise_dt  = parse_local(LOG_DATE,      SUNRISE_STR)
cutoff_dt   = sunrise_dt + timedelta(hours=1)

NOISE_FLOOR = 3.0   # counts, from calibration residual std

print(f'Experiment  : {experiment_id}')
print(f'GoPro start : {gopro_start}')
print(f'Shadow event: {shadow_dt}')
print(f'Sunrise     : {sunrise_dt}')
print(f'Cutoff      : {cutoff_dt}')

Experiment  : 2026-07-31_Moon-Temp_003
GoPro start : 2026-07-30 22:03:00
Shadow event: 2026-07-31 01:11:00
Sunrise     : 2026-07-31 06:04:00
Cutoff      : 2026-07-31 07:04:00

# Load raw overnight CSVs and trim to analysis window
import subprocess

def ensure_csv(date_str):
    path = EXPERIMENT_DIR / f'{date_str}.csv'
    if not path.exists():
        remote = f'rooster@10.0.10.20:~/claude-code/data/logs/{date_str}.csv'
        print(f'Fetching {remote} ...')
        subprocess.run(['scp', remote, str(path)], check=True)
    return path

df = pd.concat([
    pd.read_csv(ensure_csv(LOG_DATE_PREV), parse_dates=['timestamp']),
    pd.read_csv(ensure_csv(LOG_DATE),      parse_dates=['timestamp']),
]).sort_values('timestamp').drop_duplicates('timestamp').reset_index(drop=True)

# Trim to analysis window
df = df[(df.timestamp >= gopro_start) & (df.timestamp <= cutoff_dt)].copy()

# Temperature: use the logger's own rail-voltage-compensated t*_c columns
# (a fixed-VCC=3.3V recompute from raw counts was systematically wrong once
# v_rail started drifting — the logger now measures v_rail and corrects for it)
for ch in ['adc0', 'adc1', 'adc2']:
    df[ch.replace('adc','temp')] = df[ch.replace('adc','t') + '_c']

# Difference signal: ADC2 vs moonlit reference
df['ref']       = (df.adc0 + df.adc1) / 2
df['delta']     = df.adc2 - df.ref
df['temp_ref']  = (df.temp0 + df.temp1) / 2
df['temp_delta']= df.temp2 - df.temp_ref

# Period labels
df['period'] = np.where(df.timestamp < shadow_dt, 'pre-shadow (moonlit)', 'post-shadow')

pre  = df[df.timestamp < shadow_dt]
post = df[df.timestamp >= shadow_dt]

print(f'Loaded {len(df)} rows  |  {df.timestamp.iloc[0]}  →  {df.timestamp.iloc[-1]}')
print(f'Pre-shadow rows : {len(pre)}')
print(f'Post-shadow rows: {len(post)}')
df[['adc0','adc1','adc2','temp0','temp1','temp2','delta']].describe().round(2)

Loaded 3246 rows  |  2026-07-30 22:03:02  →  2026-07-31 07:03:52
Pre-shadow rows : 1128
Post-shadow rows: 2118

	adc0	adc1	adc2	temp0	temp1	temp2	delta
count	3246.00	3246.00	3246.00	3246.00	3246.00	3246.00	3246.00
mean	17616.73	17559.71	17447.60	10.55	10.74	11.13	-140.61
std	928.07	943.97	933.75	3.19	3.24	3.19	16.90
min	13655.00	13601.00	13476.00	6.83	7.01	7.34	-183.50
25%	16646.50	16567.25	16472.00	7.79	7.94	8.37	-151.50
50%	17808.00	17761.00	17634.00	9.88	10.04	10.48	-137.50
75%	18418.00	18374.00	18256.00	13.85	14.11	14.43	-127.62
max	18696.00	18645.00	18552.00	23.77	23.95	24.37	-107.00

# Interactive overview: all 3 ADC channels + temperature
fig = make_subplots(rows=2, cols=1, shared_xaxes=True, vertical_spacing=0.06,
                    subplot_titles=['ADC counts (offset-corrected)', 'Temperature (°C)'])

for col, color, name in [('adc0',C0,'ADC0'),('adc1',C1,'ADC1'),('adc2',C2,'ADC2 (shadow ch)')]:
    fig.add_trace(go.Scatter(x=df.timestamp, y=df[col], name=name,
                             line=dict(color=color, width=1.2), legendgroup=name), row=1, col=1)

for col, color, name in [('temp0',C0,'ADC0'),('temp1',C1,'ADC1'),('temp2',C2,'ADC2')]:
    fig.add_trace(go.Scatter(x=df.timestamp, y=df[col], name=name,
                             line=dict(color=color, width=1.2),
                             legendgroup=name, showlegend=False), row=2, col=1)

# Event markers
for row in [1, 2]:
    fig.add_vline(x=shadow_dt, line_dash='dash', line_color=C_SHADOW,
                  annotation_text='Shadow event', annotation_font_color=C_SHADOW,
                  annotation_position='top right', row=row, col=1)
    fig.add_vline(x=sunrise_dt, line_dash='dot', line_color=C_SUN,
                  annotation_text='Sunrise', annotation_font_color=C_SUN,
                  annotation_position='top left', row=row, col=1)

fig.update_layout(**rooster_layout(
    title=f'{experiment_id} — Full session overview',
    height=600,
    xaxis2=dict(gridcolor=GRID, linecolor=GRID),
    yaxis2=dict(gridcolor=GRID, linecolor=GRID),
))
fig.show()

# Difference signal: ADC2 − mean(ADC0, ADC1)
# H₀ implies a constant offset. H₁ implies a slope change (progressive divergence)
# after shadow onset — not a step. Rolling slope panel shows rate of change over time.

ROLL_MEAN = 30   # rolling mean window (minutes)
ROLL_SLOPE = 45  # rolling slope window (minutes) — wider for stable slope estimate

df['delta_roll'] = df['delta'].rolling(ROLL_MEAN, center=True, min_periods=10).mean()

# Rolling OLS slope: fit delta ~ time in sliding window, report counts/min
def compute_rolling_slope(df_in, col, window_minutes):
    from datetime import timedelta
    half = timedelta(minutes=window_minutes / 2)
    slopes = []
    ts = df_in.timestamp.values
    vals = df_in[col].values
    for i, t in enumerate(df_in.timestamp):
        mask = (df_in.timestamp >= t - half) & (df_in.timestamp <= t + half)
        seg = df_in[mask]
        if len(seg) >= 5:
            x = (seg.timestamp - t).dt.total_seconds().values / 60
            b = np.polyfit(x, seg[col].values, 1)[0]
        else:
            b = np.nan
        slopes.append(b)
    return slopes

df['delta_slope'] = compute_rolling_slope(df, 'delta', ROLL_SLOPE)

fig = make_subplots(
    rows=2, cols=1, shared_xaxes=True, vertical_spacing=0.08,
    subplot_titles=[f'ADC2 − reference (Δ counts)  ·  {ROLL_MEAN}-min rolling mean',
                    f'Rolling slope of Δ  ({ROLL_SLOPE}-min window, counts/min)']
)

# Panel 1: delta signal
fig.add_trace(go.Scatter(
    x=df.timestamp, y=df.delta,
    name='ADC2 − ref (raw)', mode='lines',
    line=dict(color=C2, width=0.8), opacity=0.30,
), row=1, col=1)
fig.add_trace(go.Scatter(
    x=df.timestamp, y=df.delta_roll,
    name=f'{ROLL_MEAN}-min mean', mode='lines',
    line=dict(color=C2, width=2.5),
), row=1, col=1)

# Panel 2: rolling slope
fig.add_trace(go.Scatter(
    x=df.timestamp, y=df.delta_slope,
    name='slope (counts/min)', mode='lines',
    line=dict(color=ACCENT, width=2.0),
    showlegend=True,
), row=2, col=1)
fig.add_hline(y=0, line_dash='dot', line_color=MUTED, row=2, col=1)

# Event lines on both panels
for r in [1, 2]:
    fig.add_vline(x=shadow_dt,  line_dash='dash', line_color=C_SHADOW,
                  annotation_text='Shadow event', annotation_font_color=C_SHADOW,
                  annotation_position='top right', row=r, col=1)
    fig.add_vline(x=sunrise_dt, line_dash='dot',  line_color=C_SUN,
                  annotation_text='Sunrise', annotation_font_color=C_SUN,
                  annotation_position='top left', row=r, col=1)

fig.update_layout(**rooster_layout(
    title=f'{experiment_id} — Difference signal and rolling slope',
    height=560,
    xaxis2=dict(gridcolor=GRID, linecolor=GRID),
    yaxis=dict(title='Δ ADC counts', gridcolor=GRID, linecolor=GRID),
    yaxis2=dict(title='counts/min', gridcolor=GRID, linecolor=GRID,
                zeroline=True, zerolinecolor=MUTED),
))
fig.show()

# Pre-shadow baseline: inter-channel statistics
# Describes the calibration state and ADC noise floor during the pre-shadow window.

print('Pre-shadow baseline')
print(f'  Period: {pre.timestamp.iloc[0]}  →  {pre.timestamp.iloc[-1]}')
print(f'  n = {len(pre)} rows')
print()
print(f'  {"Signal":<20}  {"mean":>8}  {"std":>6}  {"median":>8}  {"mV equiv":>10}')
print('  ' + '-'*60)
for label, series in [
    ('ADC0 − ADC1', pre.adc0 - pre.adc1),
    ('ADC0 − ADC2', pre.adc0 - pre.adc2),
    ('ADC1 − ADC2', pre.adc1 - pre.adc2),
    ('delta (ADC2−ref)', pre.delta),
]:
    print(f'  {label:<20}  {series.mean():>+8.1f}  {series.std():>6.1f}  {series.median():>+8.1f}  {series.mean()*0.125:>+10.3f} mV')
print()
print(f'  Noise floor reference (calibration residual std): ±{NOISE_FLOOR:.1f} counts')

Pre-shadow baseline
  Period: 2026-07-30 22:03:02  →  2026-07-31 01:10:52
  n = 1128 rows

  Signal                    mean     std    median    mV equiv
  ------------------------------------------------------------
  ADC0 − ADC1              +80.4     7.2     +80.0     +10.046 mV
  ADC0 − ADC2             +177.9    17.6    +171.0     +22.242 mV
  ADC1 − ADC2              +97.6    18.8     +92.0     +12.197 mV
  delta (ADC2−ref)        -137.8    17.8    -131.5     -17.219 mV

  Noise floor reference (calibration residual std): ±3.0 counts

# Distribution comparison: pre-shadow vs post-shadow delta

all_delta = pd.concat([pre.delta, post.delta])
bins = np.linspace(all_delta.min(), all_delta.max(), 50)

fig = make_subplots(rows=1, cols=2,
                    subplot_titles=['Distribution of Δ counts', 'Box plot by period'])

# Histogram
for df_part, color, name in [(pre, MUTED, 'Pre-shadow'), (post, ACCENT, 'Post-shadow')]:
    fig.add_trace(go.Histogram(
        x=df_part.delta, name=name, opacity=0.7,
        xbins=dict(start=bins[0], end=bins[-1], size=bins[1]-bins[0]),
        marker_color=color,
    ), row=1, col=1)

# Box plots
for df_part, color, name in [(pre, MUTED, 'Pre-shadow'), (post, ACCENT, 'Post-shadow')]:
    fig.add_trace(go.Box(
        y=df_part.delta, name=name, marker_color=color,
        boxmean='sd', showlegend=False,
    ), row=1, col=2)

fig.update_layout(**rooster_layout(
    title=f'{experiment_id} — Pre vs post shadow delta distributions',
    barmode='overlay',
    height=420,
    xaxis=dict(title='Δ ADC counts', gridcolor=GRID, linecolor=GRID),
    yaxis=dict(title='count', gridcolor=GRID, linecolor=GRID),
    xaxis2=dict(gridcolor=GRID, linecolor=GRID),
    yaxis2=dict(title='Δ ADC counts', gridcolor=GRID, linecolor=GRID),
))
fig.show()

# Hypothesis test: slope comparison (rate of change of delta, pre vs post shadow)
# H₀: slope of delta is the same in both periods
# H₁: slope differs after shadow onset (sustained divergence)
#
# Autocorrelation-corrected SE: HC3 residual std scaled by sqrt(n/n_eff)

def n_eff(series):
    """Effective sample size correcting for lag-1 autocorrelation."""
    rho = series.autocorr(lag=1)
    n   = len(series)
    return max(int(n * (1 - rho) / (1 + rho)), 2)

def slope_and_corrected_se(df_seg, t0):
    x = (df_seg.timestamp - t0).dt.total_seconds() / 60
    y = df_seg.delta.values
    X = sm.add_constant(x.values)
    m = sm.OLS(y, X).fit(cov_type='HC3')
    rho  = pd.Series(m.resid).autocorr(lag=1)
    ne   = max(int(len(y) * (1 - rho) / (1 + rho)), 2)
    se_c = m.bse[1] * np.sqrt(len(y) / ne)
    return m.params[1], se_c, ne, m

t0 = df.timestamp.iloc[0]
slope_pre,  se_pre,  neff_pre,  m_pre  = slope_and_corrected_se(pre.dropna(subset=['delta']),  t0)
slope_post, se_post, neff_post, m_post = slope_and_corrected_se(post.dropna(subset=['delta']), t0)

slope_diff = slope_post - slope_pre
se_diff    = np.sqrt(se_pre**2 + se_post**2)
t_stat     = slope_diff / se_diff
df_t       = neff_pre + neff_post - 4
p_val      = 2 * stats.t.sf(abs(t_stat), df_t)

sig_label = lambda p: '*** p<0.001' if p<0.001 else '** p<0.01' if p<0.01 else '* p<0.05' if p<0.05 else '(ns)'

print('── Slope comparison test: H₀ = same rate of change pre vs post shadow ──')
print()
print(f'  Pre-shadow  slope: {slope_pre:+.4f} counts/min  SE={se_pre:.4f}  n={len(pre)}  n_eff={neff_pre}')
print(f'  Post-shadow slope: {slope_post:+.4f} counts/min  SE={se_post:.4f}  n={len(post)}  n_eff={neff_post}')
print()
print(f'  Slope change (post − pre): {slope_diff:+.4f} counts/min')
print(f'  t = {t_stat:.3f}  df = {df_t}  p = {p_val:.4f}  {sig_label(p_val)}')
print()
print(f'  Noise floor reference: ±{NOISE_FLOOR:.1f} counts total')
print()
if p_val < 0.05:
    direction = 'rising' if slope_diff > 0 else 'falling'
    print(f'  ✓ Reject H₀ — slope change of {slope_diff:+.4f} counts/min is significant.')
    print(f'    ADC2 delta is {direction} relative to reference at {abs(slope_diff):.4f} counts/min post-shadow.')
else:
    print(f'  ✗ Fail to reject H₀ — no significant slope change detected (p={p_val:.3f}).')
    print(f'    The rate of divergence is not distinguishable from pre-shadow trend.')

── Slope comparison test: H₀ = same rate of change pre vs post shadow ──

  Pre-shadow  slope: -0.2947 counts/min  SE=0.0382  n=1128  n_eff=16
  Post-shadow slope: +0.1158 counts/min  SE=0.0309  n=2118  n_eff=12

  Slope change (post − pre): +0.4105 counts/min
  t = 8.358  df = 24  p = 0.0000  *** p<0.001

  Noise floor reference: ±3.0 counts total

  ✓ Reject H₀ — slope change of +0.4105 counts/min is significant.
    ADC2 delta is rising relative to reference at 0.4105 counts/min post-shadow.

# Interrupted time-series regression
# Model: delta ~ time + shadow_on + shadow_on*time
# Primary result: interaction coefficient (slope change at shadow event)
# Secondary result: shadow_on level shift (may capture pre-shadow transition effect)

reg = df[['timestamp','delta']].dropna().copy()
t0r = reg.timestamp.iloc[0]
reg['time_min']    = (reg.timestamp - t0r).dt.total_seconds() / 60
reg['shadow_on']   = (reg.timestamp >= shadow_dt).astype(float)
reg['interaction'] = reg.time_min * reg.shadow_on

X     = sm.add_constant(reg[['time_min','shadow_on','interaction']])
model = sm.OLS(reg.delta, X).fit(cov_type='HC3')

i_coef = model.params['interaction']
i_ci   = model.conf_int().loc['interaction']
i_p    = model.pvalues['interaction']

l_coef = model.params['shadow_on']
l_ci   = model.conf_int().loc['shadow_on']
l_p    = model.pvalues['shadow_on']

print('── Segmented OLS regression (HC3 robust standard errors) ──')
print()
print('  PRIMARY — slope change (interaction term):')
print(f'    interaction coef : {i_coef:+.4f} counts/min')
print(f'    95% CI           : [{i_ci[0]:+.4f}, {i_ci[1]:+.4f}] counts/min')
print(f'    p-value          : {i_p:.4f}  {"*** p<0.001" if i_p<0.001 else "** p<0.01" if i_p<0.01 else "* p<0.05" if i_p<0.05 else "(ns)"}')
print()
print('  secondary — level shift at shadow_dt (shadow_on term):')
print(f'    shadow_on coef   : {l_coef:+.2f} counts  ({l_coef*0.125:+.3f} mV)')
print(f'    95% CI           : [{l_ci[0]:+.2f}, {l_ci[1]:+.2f}] counts')
print(f'    p-value          : {l_p:.4f}  {"*** p<0.001" if l_p<0.001 else "** p<0.01" if l_p<0.01 else "* p<0.05" if l_p<0.05 else "(ns)"}')
print()
print(f'  Pre-shadow trend   : {model.params["time_min"]:+.4f} counts/min')
print(f'  Post-shadow trend  : {model.params["time_min"] + i_coef:+.4f} counts/min')
print(f'  R²                 : {model.rsquared:.4f}')

# Visualization: two separate segment trend lines (makes slope difference visible)
pre_reg  = reg[reg.shadow_on == 0].copy()
post_reg = reg[reg.shadow_on == 1].copy()

pre_fit  = sm.OLS(pre_reg.delta,  sm.add_constant(pre_reg.time_min)).fit()
post_fit = sm.OLS(post_reg.delta, sm.add_constant(post_reg.time_min)).fit()

fig = go.Figure()
fig.add_trace(go.Scatter(x=reg.timestamp, y=reg.delta,
                         name='Δ (raw)', mode='lines',
                         line=dict(color=C2, width=0.8), opacity=0.35))
fig.add_trace(go.Scatter(x=pre_reg.timestamp,  y=pre_fit.fittedvalues,
                         name=f'Pre-shadow trend ({pre_fit.params['time_min']:+.4f} cnt/min)',
                         mode='lines', line=dict(color=MUTED, width=2.5, dash='solid')))
fig.add_trace(go.Scatter(x=post_reg.timestamp, y=post_fit.fittedvalues,
                         name=f'Post-shadow trend ({post_fit.params['time_min']:+.4f} cnt/min)',
                         mode='lines', line=dict(color=ACCENT, width=2.5, dash='solid')))
fig.add_vline(x=shadow_dt,  line_dash='dash', line_color=C_SHADOW,
              annotation_text='Shadow event', annotation_font_color=C_SHADOW)
fig.add_vline(x=sunrise_dt, line_dash='dot',  line_color=C_SUN,
              annotation_text='Sunrise', annotation_font_color=C_SUN)
fig.update_layout(**rooster_layout(
    title=f'Segmented trends  |  interaction (slope Δ) = {i_coef:+.4f} cnt/min  p={i_p:.4f}',
    yaxis_title='Δ ADC counts', height=420,
))
fig.show()

── Segmented OLS regression (HC3 robust standard errors) ──

  PRIMARY — slope change (interaction term):
    interaction coef : +0.4105 counts/min
    95% CI           : [+0.4005, +0.4205] counts/min
    p-value          : 0.0000  *** p<0.001

  secondary — level shift at shadow_dt (shadow_on term):
    shadow_on coef   : -74.26 counts  (-9.282 mV)
    95% CI           : [-76.25, -72.27] counts
    p-value          : 0.0000  *** p<0.001

  Pre-shadow trend   : -0.2947 counts/min
  Post-shadow trend  : +0.1158 counts/min
  R²                 : 0.6450

# Segmented OLS regression — temperature (°C) scale
# Same interrupted time-series model as the ADC-count version above,
# expressed in °C using the logger's own rail-compensated t*_c values.

reg_c = df[['timestamp','temp_delta']].dropna().copy()
t0rc  = reg_c.timestamp.iloc[0]
reg_c['time_min']    = (reg_c.timestamp - t0rc).dt.total_seconds() / 60
reg_c['shadow_on']   = (reg_c.timestamp >= shadow_dt).astype(float)
reg_c['interaction'] = reg_c.time_min * reg_c.shadow_on

X_c      = sm.add_constant(reg_c[['time_min','shadow_on','interaction']])
model_c  = sm.OLS(reg_c.temp_delta, X_c).fit(cov_type='HC3')

ic_coef = model_c.params['interaction']
ic_ci   = model_c.conf_int().loc['interaction']
ic_p    = model_c.pvalues['interaction']

lc_coef = model_c.params['shadow_on']
lc_ci   = model_c.conf_int().loc['shadow_on']
lc_p    = model_c.pvalues['shadow_on']

print('── Segmented OLS regression — °C scale (HC3 robust standard errors) ──')
print()
print('  PRIMARY — slope change (interaction term):')
print(f'    interaction coef : {ic_coef:+.6f} °C/min')
print(f'    95% CI           : [{ic_ci[0]:+.6f}, {ic_ci[1]:+.6f}] °C/min')
print(f'    p-value          : {ic_p:.4f}  {"*** p<0.001" if ic_p<0.001 else "** p<0.01" if ic_p<0.01 else "* p<0.05" if ic_p<0.05 else "(ns)"}')
print()
print('  secondary — level shift at shadow_dt (shadow_on term):')
print(f'    shadow_on coef   : {lc_coef:+.4f} °C')
print(f'    95% CI           : [{lc_ci[0]:+.4f}, {lc_ci[1]:+.4f}] °C')
print(f'    p-value          : {lc_p:.4f}  {"*** p<0.001" if lc_p<0.001 else "** p<0.01" if lc_p<0.01 else "* p<0.05" if lc_p<0.05 else "(ns)"}')
print()
print(f'  Pre-shadow trend   : {model_c.params["time_min"]:+.6f} °C/min')
print(f'  Post-shadow trend  : {model_c.params["time_min"] + ic_coef:+.6f} °C/min')
print(f'  R²                 : {model_c.rsquared:.4f}')

pre_rc  = reg_c[reg_c.shadow_on == 0].copy()
post_rc = reg_c[reg_c.shadow_on == 1].copy()

pre_fitc  = sm.OLS(pre_rc.temp_delta,  sm.add_constant(pre_rc.time_min)).fit()
post_fitc = sm.OLS(post_rc.temp_delta, sm.add_constant(post_rc.time_min)).fit()

fig = go.Figure()
fig.add_trace(go.Scatter(x=reg_c.timestamp, y=reg_c.temp_delta,
                         name='Δ °C (raw)', mode='lines',
                         line=dict(color=C2, width=0.8), opacity=0.35))
fig.add_trace(go.Scatter(x=pre_rc.timestamp,  y=pre_fitc.fittedvalues,
                         name=f'Pre-shadow trend ({pre_fitc.params["time_min"]:+.6f} °C/min)',
                         mode='lines', line=dict(color=MUTED, width=2.5, dash='solid')))
fig.add_trace(go.Scatter(x=post_rc.timestamp, y=post_fitc.fittedvalues,
                         name=f'Post-shadow trend ({post_fitc.params["time_min"]:+.6f} °C/min)',
                         mode='lines', line=dict(color=ACCENT, width=2.5, dash='solid')))
fig.add_vline(x=shadow_dt,  line_dash='dash', line_color=C_SHADOW,
              annotation_text='Shadow event', annotation_font_color=C_SHADOW)
fig.add_vline(x=sunrise_dt, line_dash='dot',  line_color=C_SUN,
              annotation_text='Sunrise', annotation_font_color=C_SUN)
fig.update_layout(**rooster_layout(
    title=f'Segmented trends (°C)  |  interaction (slope Δ) = {ic_coef:+.6f} °C/min  p={ic_p:.4f}',
    yaxis_title='Δ °C', height=420,
))
fig.show()

── Segmented OLS regression — °C scale (HC3 robust standard errors) ──

  PRIMARY — slope change (interaction term):
    interaction coef : -0.001414 °C/min
    95% CI           : [-0.001450, -0.001377] °C/min
    p-value          : 0.0000  *** p<0.001

  secondary — level shift at shadow_dt (shadow_on term):
    shadow_on coef   : +0.2619 °C
    95% CI           : [+0.2550, +0.2688] °C
    p-value          : 0.0000  *** p<0.001

  Pre-shadow trend   : +0.001066 °C/min
  Post-shadow trend  : -0.000348 °C/min
  R²                 : 0.6419

# Difference signal in °C scale
# Same slope-change analysis as ADC counts, expressed in °C-equivalent units.
# Note: higher ADC = cooler, lower ADC = warmer — °C conversion reflects this.

df['tdelta_roll'] = df.temp_delta.rolling(ROLL_MEAN, center=True, min_periods=10).mean()
df['tdelta_slope'] = compute_rolling_slope(df, 'temp_delta', ROLL_SLOPE)

fig = make_subplots(
    rows=2, cols=1, shared_xaxes=True, vertical_spacing=0.08,
    subplot_titles=['ADC2 − reference in °C  ·  rolling mean',
                    'Rolling slope (°C/min)']
)

fig.add_trace(go.Scatter(x=df.timestamp, y=df.temp_delta,
                         name='Δ °C (raw)', mode='lines',
                         line=dict(color=C2, width=0.8), opacity=0.30), row=1, col=1)
fig.add_trace(go.Scatter(x=df.timestamp, y=df.tdelta_roll,
                         name=f'{ROLL_MEAN}-min mean', mode='lines',
                         line=dict(color=C2, width=2.5)), row=1, col=1)
fig.add_trace(go.Scatter(x=df.timestamp, y=df.tdelta_slope,
                         name='slope (°C/min)', mode='lines',
                         line=dict(color=ACCENT, width=2.0)), row=2, col=1)
fig.add_hline(y=0, line_dash='dot', line_color=MUTED, row=2, col=1)

for r in [1, 2]:
    fig.add_vline(x=shadow_dt,  line_dash='dash', line_color=C_SHADOW,
                  annotation_text='Shadow event', annotation_font_color=C_SHADOW,
                  annotation_position='top right', row=r, col=1)
    fig.add_vline(x=sunrise_dt, line_dash='dot',  line_color=C_SUN,
                  annotation_text='Sunrise', annotation_font_color=C_SUN,
                  annotation_position='top left', row=r, col=1)

fig.update_layout(**rooster_layout(
    title=f'{experiment_id} — Difference signal in °C scale',
    height=560,
    xaxis2=dict(gridcolor=GRID, linecolor=GRID),
    yaxis=dict(title='Δ °C', gridcolor=GRID, linecolor=GRID),
    yaxis2=dict(title='°C/min', gridcolor=GRID, linecolor=GRID,
                zeroline=True, zerolinecolor=MUTED),
))
fig.show()

from IPython.display import Markdown, display

decision = (
    f'Reject H\u2080 \u2014 slope change of {slope_diff:+.4f} counts/min is significant ({sig_label(p_val)}). '
    f'ADC2 delta is {"rising" if slope_diff > 0 else "falling"} relative to reference post-shadow.'
    if p_val < 0.05
    else f'Fail to reject H\u2080 \u2014 no significant slope change detected (p={p_val:.3f}).'
)

md = f"""
## Conclusions

### Hypothesis test result

| Test | Result | Interpretation |
|---|---|---|
| Pre-shadow slope | {slope_pre:+.4f} counts/min (n_eff={neff_pre}) | baseline rate of change |
| Post-shadow slope | {slope_post:+.4f} counts/min (n_eff={neff_post}) | rate after shadow onset |
| Slope change (post − pre) | {slope_diff:+.4f} counts/min &nbsp; p={p_val:.4f} &nbsp; {sig_label(p_val)} | **primary test of H\u2081** |
| OLS interaction coef | {i_coef:+.4f} counts/min &nbsp; p={i_p:.4f} &nbsp; {sig_label(i_p)} | confirms slope change |
| OLS level shift (shadow_on) | {l_coef:+.2f} counts &nbsp; p={l_p:.4f} &nbsp; {sig_label(l_p)} | transition-period offset |

### Rolling slope chart

Onset timing relative to shadow_dt: *leads / coincides / lags by ~___ min*

Shape: *gradual ramp / abrupt shift / no visible change*

### Interpretation

**Decision:** {decision}

Notes: *(partial shadow transition, lag, confounders, etc.)*

### Companion video

- GoPro timelapse source: `{experiment_id}`
- DaVinci overlay: `davinci/adc_bar_chart_overlay.mov`
- Key timestamp: shadow event at **{SHADOW_TIME_STR} MDT**
"""

display(Markdown(md))

Conclusions

Hypothesis test result

Test	Result	Interpretation
Pre-shadow slope	-0.2947 counts/min (n_eff=16)	baseline rate of change
Post-shadow slope	+0.1158 counts/min (n_eff=12)	rate after shadow onset
Slope change (post − pre)	+0.4105 counts/min   p=0.0000   *** p<0.001	primary test of H₁
OLS interaction coef	+0.4105 counts/min   p=0.0000   *** p<0.001	confirms slope change
OLS level shift (shadow_on)	-74.26 counts   p=0.0000   *** p<0.001	transition-period offset

Rolling slope chart

Onset timing relative to shadow_dt: leads / coincides / lags by ~___ min

Shape: gradual ramp / abrupt shift / no visible change

Interpretation

Decision: Reject H₀ — slope change of +0.4105 counts/min is significant (*** p<0.001). ADC2 delta is rising relative to reference post-shadow.

Notes: (partial shadow transition, lag, confounders, etc.)

Companion video

• GoPro timelapse source: 2026-07-31_Moon-Temp_003
• DaVinci overlay: davinci/adc_bar_chart_overlay.mov
• Key timestamp: shadow event at 01:11 MDT