**450**wrote:

New to GWAS, been working on genomic control (GC,lambda). I see sometimes GC reported for the the 50th percentile, 10th percentile, 1st percentile, and 1/10th of a percentile.

I adapted code from this post: Calculating Genomic Inflation Factor

```
get_lambda <- function(in_path) {
S <- fread(in_path)
z <- qnorm(S$P.VALUE / 2)
lambda <- round(median(z^2) / 0.454, 3)
return(lambda)
}
```

This seems to replicate lambda for the 50th percentile, how would I go about getting the other percentiles?

**66k**• written 3 months ago by curious •

**450**