R: Credible set

cs {gap} | R Documentation

## Credible set

### Description

The function implements credible set as in fine-mapping.

### Usage

    cs(tbl, b = "Effect", se = "StdErr", log_p = NULL, cutoff = 0.95)

### Arguments

`tbl` | Input data.
`b` | Effect size.
`se` | Standard error.
`log_p` | if not NULL it will be used to derive z-statistic
`cutoff` | Threshold for inclusion.

### Value

Credible set.

### Examples

    ## Not run:
    \preformatted{
      zcat METAL/4E.BP1-1.tbl.gz | \
      awk 'NR==1 || ($1==4 && $2 >= 187158034 - 1e6 && $2 < 187158034 + 1e6)' > 4E.BP1.z
    }
      tbl <- within(read.delim("4E.BP1.z"),{logp <- logp(Effect/StdErr)})
      z <- cs(tbl)
      l <- cs(tbl,log_p="logp")

    ## End(Not run)

[Package gap version 1.2.3-6 Index]