Variance of the estimated effect size for case-control data

Var.data.cc(f, N, s)

Arguments

f

Minor allele frequencies

N

Total sample size (N0+N1)

s

Proportion of cases (N1/N0+N1)

Value

Variance of estimated effect size \(\hat{\beta}\), V.

Examples

maf = runif(100, 0.05, 0.5) N0 = 5000 # number of controls N1 = 5000 # number of cases Var.data.cc(f = maf, N = N0 + N1, s = N1/(N0+N1))
#> [1] 0.0025353881 0.0008181891 0.0009185978 0.0018838648 0.0039615189 #> [6] 0.0010398200 0.0010054095 0.0013526966 0.0008490723 0.0008349987 #> [11] 0.0008103213 0.0017832687 0.0032716992 0.0012781987 0.0011257163 #> [16] 0.0016807679 0.0011238656 0.0027600597 0.0011472589 0.0008003876 #> [21] 0.0013523677 0.0008731582 0.0008481271 0.0016794491 0.0008002455 #> [26] 0.0008457706 0.0029498620 0.0009741453 0.0008648122 0.0008682137 #> [31] 0.0033360702 0.0015558156 0.0013244121 0.0008959042 0.0010254392 #> [36] 0.0010827900 0.0008600363 0.0008017172 0.0017551329 0.0015896011 #> [41] 0.0008722764 0.0010043137 0.0008928560 0.0008824952 0.0023662309 #> [46] 0.0008372615 0.0008359197 0.0008000559 0.0008005635 0.0011464753 #> [51] 0.0010459693 0.0012899149 0.0017846766 0.0009729199 0.0010097017 #> [56] 0.0008328569 0.0016427071 0.0008570213 0.0027378034 0.0012081117 #> [61] 0.0008203023 0.0013965148 0.0009408861 0.0012449967 0.0009216941 #> [66] 0.0017001414 0.0008017721 0.0009633361 0.0009615185 0.0013829869 #> [71] 0.0010637900 0.0011763816 0.0034068508 0.0010402951 0.0011450982 #> [76] 0.0036006916 0.0011668847 0.0009488545 0.0008134581 0.0011536551 #> [81] 0.0009762270 0.0009186820 0.0014334696 0.0013519527 0.0010242762 #> [86] 0.0008041682 0.0011281913 0.0016047137 0.0008764885 0.0028351501 #> [91] 0.0008000056 0.0019349560 0.0009849161 0.0008156441 0.0008549669 #> [96] 0.0014987998 0.0009594483 0.0008180840 0.0034092500 0.0010363472