Skip to content

Instantly share code, notes, and snippets.

using Distributions
using GLM
using RCall
using Random
Random.seed!(2)
n = 10000
ses = collect(1:n)
pgs = rand(Normal(), n)
@ymer
ymer / make24.ipynb
Last active September 5, 2018 09:41
Make 24
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
@ymer
ymer / repository.md
Last active August 28, 2017 22:48
Bartholin biomedicine script repository

##Munge sumstats Takes an unprocessed GWAS sumstat file, performs QC and changes it to a standard format. Location: https://github.com/ymer/ldsc Details at: ? Contact: Jonatan

using Combinatorics
function find_matching_pairs(s)
l, starts = [], []
for i in 1:length(s)
if s[i] == "("
push!(starts, i)
elseif s[i] == ")"
push!(l, (pop!(starts), i))
end