From Genome Analysis Wiki
Jump to navigationJump to search
1,194 bytes added
, 17:47, 13 April 2015
Line 72: |
Line 72: |
| library(rareMETALS) | | library(rareMETALS) |
| res <- rareMETALS.single(score.stat.file,cov.file=NULL,range="19:11200093-11201275",alternative="two.sided",ix.gold=1,callrate.cutoff=0,hwe.cutoff=0) | | res <- rareMETALS.single(score.stat.file,cov.file=NULL,range="19:11200093-11201275",alternative="two.sided",ix.gold=1,callrate.cutoff=0,hwe.cutoff=0) |
| + | |
| + | ###result can be explored below### |
| + | > names(res10) |
| + | [1] "p.value" "ref" "alt" "integratedData" "raw.data" |
| + | [6] "clean.data" "statistic" "direction.by.study" "anno" "maf" |
| + | [11] "QC.by.study" "no.sample" "beta1.est" "beta1.sd" "hsq.est" |
| + | [16] "nearby" "pos" |
| + | > print(res10$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(res10$p.value) |
| + | [1] 0.551263675 0.056308558 0.172481571 0.734935815 0.922326732 0.053804524 0.886985353 0.903835162 0.005280228 0.266575301 |
| + | [11] 0.196122312 0.157114376 0.951477852 0.840523624 0.759482777 0.112743041 0.414147263 0.825877149 0.006090142 0.096474975 |
| + | [21] 0.096474975 0.956407850 0.038234190 0.253512486 0.550935361 0.482315038 |
| | | |
| == A Simple Tutorial for Using the rareMETALS.single.group function == | | == A Simple Tutorial for Using the rareMETALS.single.group function == |