Skip to content

Commit 172a68d

Browse files
authored
Add effective_dim.py with participation_ratio function (#152)
1 parent 91ff5f1 commit 172a68d

1 file changed

Lines changed: 11 additions & 0 deletions

File tree

Lines changed: 11 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,11 @@
1+
from jax import numpy as jnp
2+
3+
def participation_ratio(Z):
4+
Zc = Z - Z.mean(axis=0, keepdims=True)
5+
cov = (Zc.T @ Zc) / (Zc.shape[0] - 1)
6+
7+
tr = jnp.trace(cov)
8+
tr2_cov = tr * tr
9+
cov2_tr = jnp.trace(cov @ cov)
10+
11+
return tr2_cov / cov2_tr if cov2_tr > 0 else float("nan")

0 commit comments

Comments
 (0)