R: Genomic control for over-dominant model of inheritance using...
VIFGC_ovdom
R Documentation
Genomic control for over-dominant model of inheritance using VIF
Description
This function estimates the corrected statistic using
genomic control for the over-dominant model, using VIF.
VIF coefficients are estimated by optimizing different
error functions: regress, median and ks.test.
Function of error to be optimized. Can be
"regress", "median" or "ks.test"
p
Input vector of allele frequencies
index.filter
Indexes for variables that will be
use for analysis in data vector
n
size of the sample
proportion
The proportion of lowest P (Chi2) to be
used when estimating the inflation factor Lambda for
"regress" method only
plot
If TRUE, plot of lambda will be produced
lmax
The threshold for lambda for plotting
(optional)
color
The color of the plot
clust
For developers only
vart0
For developers only
tmp
For developers only
Value
A list with elements
Zx
output vector corrected
Chi square statistic
vv
output vector of VIF
exeps
output vector of exepsons (NA)
calrate
output vector of calrate
F
F
K
K
Author(s)
Yakov Tsepilov
Examples
require(GenABEL.data)
data(ge03d2)
# truncate the data to make the example faster
ge03d2 <- ge03d2[seq(from=1,to=nids(ge03d2),by=2),seq(from=1,to=nsnps(ge03d2),by=3)]
qts <- mlreg(phdata(ge03d2)$dm2~1,data=ge03d2,gtmode = "overdominant")
chi2.1df <- results(qts)$chi2.1df
s <- summary(ge03d2)
freq <- s$Q.2
result <- VIFGC_ovdom(p=freq,method = "median",data=chi2.1df,n=nids(ge03d2))
Results
R version 3.3.1 (2016-06-21) -- "Bug in Your Hair"
Copyright (C) 2016 The R Foundation for Statistical Computing
Platform: x86_64-pc-linux-gnu (64-bit)
R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.
R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.
Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.
> library(GenABEL)
Loading required package: MASS
Loading required package: GenABEL.data
> png(filename="/home/ddbj/snapshot/RGM3/R_CC/result/GenABEL/VIFGC_ovdom.Rd_%03d_medium.png", width=480, height=480)
> ### Name: VIFGC_ovdom
> ### Title: Genomic control for over-dominant model of inheritance using VIF
> ### Aliases: VIFGC_ovdom
> ### Keywords: htest
>
> ### ** Examples
>
> require(GenABEL.data)
> data(ge03d2)
> # truncate the data to make the example faster
> ge03d2 <- ge03d2[seq(from=1,to=nids(ge03d2),by=2),seq(from=1,to=nsnps(ge03d2),by=3)]
> qts <- mlreg(phdata(ge03d2)$dm2~1,data=ge03d2,gtmode = "overdominant")
Warning message:
In mlreg(phdata(ge03d2)$dm2 ~ 1, data = ge03d2, gtmode = "overdominant") :
You analysed X chromosome without adjustment for sex
> chi2.1df <- results(qts)$chi2.1df
> s <- summary(ge03d2)
> freq <- s$Q.2
> result <- VIFGC_ovdom(p=freq,method = "median",data=chi2.1df,n=nids(ge03d2))
>
>
>
>
>
> dev.off()
null device
1
>
providing 
.