98
edits
Changes
From Genome Analysis Wiki
Jump to navigationJump to search→A Simple Tutorial for Using the rareMETALS.single.group function
alternative = c("two.sided"), callrate.cutoff = 0,
alternative = c("two.sided"), callrate.cutoff = 0,
hwe.cutoff = 0, correctFlip = TRUE, analyzeRefAltListOnly = TRUE)
hwe.cutoff = 0, correctFlip = TRUE, analyzeRefAltListOnly = TRUE)
print(res31$res.out)
###result can be explored as below###
> names(res31)
[1] "p.value" "ref" "alt" "integratedData" "raw.data"
[6] "clean.data" "statistic" "direction.by.study" "anno" "maf"
[11] "maf.byStudy" "maf.maxdiff.vec" "ix.maf.maxdiff.vec" "maf.sd.vec" "no.sample.mat"
[16] "no.sample" "beta1.est" "beta1.sd" "QC.by.study" "hsq.est"
[21] "nearby" "cochranQ.stat" "cochranQ.df" "cochranQ.pVal" "I2"
[26] "log.mat" "pos"
> print(res31$pos)
[1] "19:11200093" "19:11200213" "19:11200235" "19:11200272" "19:11200282" "19:11200309" "19:11200412" "19:11200419"
[9] "19:11200431" "19:11200442" "19:11200475" "19:11200508" "19:11200514" "19:11200557" "19:11200579" "19:11200728"
[17] "19:11200753" "19:11200754" "19:11200806" "19:11200839" "19:11200840" "19:11200896" "19:11201124" "19:11201259"
[25] "19:11201274" "19:11201275"
> print(res31$p.value)
[1] NA NA NA NA 0.9223267 NA NA NA NA NA NA NA
[13] NA NA NA NA NA NA NA NA NA NA NA NA
[25] NA NA
== A Simple Tutorial for Using the rareMETALS.range function ==
== A Simple Tutorial for Using the rareMETALS.range function ==
