forked from cgmisc-team/cgmisc
-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathgwaa.to.bed.Rd
40 lines (36 loc) · 976 Bytes
/
gwaa.to.bed.Rd
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
% Generated by roxygen2 (4.1.1): do not edit by hand
% Please edit documentation in R/gwaa.to.bed.r
\name{gwaa.to.bed}
\alias{gwaa.to.bed}
\alias{gwaa2bed}
\title{Create a BED file with GWAS p-values in a region}
\usage{
gwaa.to.bed(chr, range, gwas, fname = "output.bed")
}
\arguments{
\item{chr}{chromosome to be considered}
\item{range}{a vector of genomic ranges}
\item{gwas}{a \code{\link[GenABEL]{scan.gwaa-class}} object produced by GenABEL}
\item{fname}{name of the output file. By default: output.bed}
}
\value{
a bed file with -log10 of corrected p-value ("Pc1df")
}
\description{
Returns a file containing GWAS p-values for a region in bed format. It can be further used as
UCSC Genome Browser custom track.
}
\examples{
\dontrun{
chr <- 36
range <- c(14e6,24e6)
gwas <- data.qc1
gwas.region.bed(chr, range, gwas, fname="my_output.bed")
}
}
\author{
Marcin Kierczak <\email{[email protected]}>
}
\keyword{bed,}
\keyword{gwas,}
\keyword{p-values}