Distribution: Difference between revisions

From ZhangLabWiki
Jump to navigation Jump to search
>Shicheng
(Created page with " table.binomial.test<-function(region,event.time){ region1<-table(region) region2<-unique(region) p=length(region)/(length(region2)*event.time) pvalue<-c() for(i in ...")
 
>Shicheng
No edit summary
Line 4: Line 4:
  region1<-table(region)
  region1<-table(region)
  region2<-unique(region)
  region2<-unique(region)
  p=length(region)/(length(region2)*event.time)
  p=length(region)/(length(region2)*event.time)
  pvalue<-c()
  pvalue<-c()
Line 11: Line 10:
   pvalue<-c(pvalue,ptmp)
   pvalue<-c(pvalue,ptmp)
  }
  }
  return(pvalue)
  return(pvalue)
  print(paste("there are",sum(pvalue<0.05/length(pvalue),"regions were over-preferred in the samples"),sep=" "))
  print(paste("there are",sum(pvalue<0.05/length(pvalue),"regions were over-preferred in the samples"),sep=" "))
  }
  }

Revision as of 18:16, 7 July 2015

table.binomial.test<-function(region,event.time){
 
region1<-table(region)
region2<-unique(region)
p=length(region)/(length(region2)*event.time)
pvalue<-c()
for(i in 1:length(table(region))){
 ptmp<-binom.test(table(region)[i],event.time,p,alternative="two.sided")$p.value
 pvalue<-c(pvalue,ptmp)
}
return(pvalue)
print(paste("there are",sum(pvalue<0.05/length(pvalue),"regions were over-preferred in the samples"),sep=" "))
}