3
\title{Genomic control for various model of inheritance using VIF}
5
GC(data, p, x, method = "regress", n,
6
index.filter = NULL, proportion = 1, clust = 0,
7
vart0 = 0, tmp = 0, CA = FALSE, p.table = 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{x}{Model of inheritance (0 for recessive,0.5 for
18
additive, 1 for dominant)}
20
\item{index.filter}{Indexes for variables that will be
21
use for analisis in data vector}
23
\item{n}{size of the sample}
25
\item{proportion}{The proportion of lowest P (Chi2) to be
26
used when estimating the inflation factor Lambda for
27
"regress" method only}
29
\item{clust}{For developers only}
31
\item{vart0}{For developers only}
33
\item{tmp}{For developers only}
35
\item{CA}{For developers only}
37
\item{p.table}{For developers only}
40
A list with elements \item{Zx}{output vector corrected
41
Chi square statistic} \item{vv}{output vector of VIF}
42
\item{exeps}{output vector of exepsons (NA)}
43
\item{calrate}{output vector of calrate} \item{F}{F}
47
This function estimates corrected statistic using genomic
48
control for diffrent models (recessive,dominant,additive
49
etc.), using VIF. VIF coefficients are estimated by
50
optimizing diffrent error functions: regress, median and
55
# truncate the data to make the example faster
56
ge03d2 <- ge03d2[seq(from=1,to=nids(ge03d2),by=2),seq(from=1,to=nsnps(ge03d2),by=3)]
57
qts <- mlreg(dm2~sex,data=ge03d2,gtmode = "dominant")
58
chi2.1df <- results(qts)$chi2.1df
61
result <- GC(p=freq,x=1,method = "median",CA=FALSE,data=chi2.1df,n=nids(ge03d2))