正規分布と一様分布があります。均一の密度に対する正規分布の密度の比率を計算したいと思います。次に、この比率の正規性をテストします。
ht <- runif(3000, 1, 18585056) # Uniform distribution
hm <- rnorm(35, 10000000, 5000000) # Normal distribution
hmd <- density(hm, from=0, to=18585056) # Kernel density of distributions over range
htd <- density(ht, from=0, to=18585056)
ratio <- hmd$y/htd$y # Ratio of kernel density values
上記の分布はhm、ht私の実験データが示す例です。実際に使用するベクトルは、R でランダムに生成されたものではありません。
QQ プロットの相関係数から、正規性の良いアイデアを得ることができることを知っています。
qqp <- qqnorm(hm)
cor(qqp$x,qqp$y)
正規分布する の場合hm、これは 1 に近い値になります。
密度ベクトルの正規性を判断する方法はありますか? 例えばhmdとratio.
(追加情報:hmおよびhtは、長さのゲノム全体でホモ接合およびヘテロ接合の SNP をモデル化しています18585056)