Skip to content

Instantly share code, notes, and snippets.

@GabrielHoffman

GabrielHoffman/coloc.txt

Created January 22, 2024 18:46
Show Gist options
  • Save GabrielHoffman/26dbdab7834d72f89224014de12d24e1 to your computer and use it in GitHub Desktop.
Save GabrielHoffman/26dbdab7834d72f89224014de12d24e1 to your computer and use it in GitHub Desktop.
Download ZIP
Color from z-scores
Raw
coloc.txt
Gene expression Y
Genotype matrix X
TODO
1) Compute sd(Y[i,]) for gene i
2) Compute var(X[,j]) for variant j
3) beta = z / sqrt(v*(n+z^2))
sd = 1 / sqrt(v*(n+z^2))
varbeta = sd^2
where v = var(X[,j])
coloc on beta, varbeta, n, sdY
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment