2010-01-12から1日間の記事一覧

短くSNP-HWE検定 with R

hweExact<-function(g){ n<-sum(g) nA<-2*g[1]+g[2] na<-2*g[3]+g[2] evod<-g[2]%%2 obs<-(g[2]-evod)/2+1 maxAa<-min(nA,na)-evod Aa<-seq(from=evod,to=maxAa,by=2) AA<-(nA-Aa)/2 aa<-(na-Aa)/2 pr<-rep(0,length(Aa)) pr<-exp(n*lgamma(2+1)+lgamma(nA+1…

SNPのHWE正確検定をRで