From Genome Analysis Wiki
Jump to navigationJump to search
501 bytes added
, 15:49, 13 April 2015
Line 74: |
Line 74: |
| | | |
| == A Simple Tutorial for Using the rareMETALS.single.group function == | | == A Simple Tutorial for Using the rareMETALS.single.group function == |
| + | [[media:Example.jpg]] |
| | | |
− | | + | res.site<-read.table("groupFile.txt",header=T) |
− | | + | refaltList <- list(pos=paste(res.site[,1],res.site[,2],sep=":"),ref=res.site$AF,alt=res.site$ALT,af=res.site$AF,af.diff.max=0.5,checkAF=T) |
− | | + | res31<-rareMETALS.single.group(score.stat.file,cov.file=NULL, range="19:11200093-11201275", refaltList, |
− | | + | alternative = c("two.sided"), callrate.cutoff = 0, |
− | | + | hwe.cutoff = 0, correctFlip = TRUE, analyzeRefAltListOnly = TRUE) |
| + | print(res31$res.out) |
| | | |
| == A Simple Tutorial for Using the rareMETALS.range function == | | == A Simple Tutorial for Using the rareMETALS.range function == |