Distribution: Difference between revisions
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=" ")) }