From Genome Analysis Wiki
Jump to navigationJump to search
506 bytes added
, 17:31, 12 April 2015
Line 116: |
Line 116: |
| example: | | example: |
| | | |
− | res<-conditional.rareMETALS.single(candidate.variant.vec=c("1:123","1:1234"), score.stat.file, cov.file, | + | res<-conditional.rareMETALS.single(candidate.variant.vec=c("19:11200282","19:11200309"), score.stat.file, cov.file, |
− | known.variant.vec=c("1:123","1:1234"), maf.cutoff=0.05, no.boot = 10000, | + | known.variant.vec=c("19:11200754","19:11200806","19:11200839"), maf.cutoff=0.05, no.boot =1000, |
| alternative = c("two.sided", "greater", "less"), ix.gold = 1, | | alternative = c("two.sided", "greater", "less"), ix.gold = 1, |
| out.digits = 4, callrate.cutoff = 0, hwe.cutoff = 0, | | out.digits = 4, callrate.cutoff = 0, hwe.cutoff = 0, |
| p.value.known.variant.vec = NA, anno.known.variant.vec = NA, | | p.value.known.variant.vec = NA, anno.known.variant.vec = NA, |
| anno.candidate.variant.vec = NA) | | anno.candidate.variant.vec = NA) |
− | print(res$res.out)
| + | print(res$res.out) |
| + | |
| + | POS REF ALT PVALUE AF BETA_EST BETA_SD DIRECTION_BY_STUDY ANNO POS_REF_ALT_ANNO_KNOWN |
| + | [1,] "19:11200282" "G" "A" "0.5825" "0.000599" "0.5616" "1.044" "-=" "N/A" "19:11200754/G/A/NA,19:11200806/C/T/NA,19:11200839/T/A/NA" |
| + | [2,] "19:11200309" "C" "A" "0.01484" "0.01538" "-0.3615" "0.02201" "+=" "N/A" "19:11200754/G/A/NA,19:11200806/C/T/NA,19:11200839/T/A/NA" |
| + | |
| | | |
| </pre> | | </pre> |