mk突变分析fortran语言
源代码
MK突变分析是一种测量常见基因突变所引起的染色体变异程度的有效方法。它使人们能够测量两个样本中间的聚类水平及其新形成的非聚类水平空间,从而更深入地了解染色体的差异性。下面是一段关于MK突变分析的fortran语言源代码:
program MKMutationAnalysis
牛板筋怎么做好吃implicit none
致自己的短句real,intent(in)::mu1, mu2, my(10)
integer, intent(in) :: n
real, intent(out) :: score1, score2
score1 = 0.0
!
Calculate mean
mu1 = sum(my(1:n))/n
!n reprents the number of elements in the array my
mu2 = sum(my(1:n))/(n-1)
!Calculate score1
阴阳怪气什么意思
do i = 1,n春节包饺子>昆虫记读书笔记
score1 = score1 + ((my(i) - mu1)**2)
end do谏太宗>狂少年
score1 = (1/n)*score1
!Calculate score 2
score2 = 0.0
美国有多少黑人
do j = 1,n
score2 = score2 + ((my(j)-mu2)**2)
end do
score2 = (1/(n-1)) * score2
end program