3
\title{Genomic control for overdimunant model of inheritance using VIF}
5
GC_ovdom(data, p, method = "regress", n,
6
index.filter = NULL, proportion = 1, clust = 0,
10
\item{data}{Input vector of Chi square statistic}
12
\item{method}{Function of error to be optimized. Can be
13
"regress", "median" or "ks.test"}
15
\item{p}{Input vector of allele frequencies}
17
\item{index.filter}{Indexes for variables that will be
18
use for analisis in data vector}
20
\item{n}{size of the sample}
22
\item{proportion}{The proportion of lowest P (Chi2) to be
23
used when estimating the inflation factor Lambda for
24
"regress" method only}
26
\item{clust}{For developers only}
28
\item{vart0}{For developers only}
30
\item{tmp}{For developers only}
33
A list with elements \item{Zx}{output vector corrected
34
Chi square statistic} \item{vv}{output vector of VIF}
35
\item{exeps}{output vector of exepsons (NA)}
36
\item{calrate}{output vector of calrate} \item{F}{F}
40
This function estimates the corrected statistic using
41
genomic control for overdimunant model, using VIF. VIF
42
coefficients are estimated by optimizing diffrent error
43
functions: regress, median and ks.test.
47
# truncate the data to make the example faster
48
ge03d2 <- ge03d2[seq(from=1,to=nids(ge03d2),by=2),seq(from=1,to=nsnps(ge03d2),by=3)]
49
qts <- mlreg(phdata(ge03d2)$dm2~1,data=ge03d2,gtmode = "overdominant")
50
chi2.1df <- results(qts)$chi2.1df
53
result <- GC_ovdom(p=freq,method = "median",data=chi2.1df,n=nids(ge03d2))