s*****n
发帖数: 2174 | 1 这不就是可重复组合数吗?
想想成5个圆盘和18个隔板拍成一列, 用圆盘和隔板的位置来代表选取的哪些小球.
最后结果就是 C(23, 18) = 33649
在R里面, 可以直接
library(gtools)
combinations(n=19, r=5, repeats.allowed = T)
就会输出一个 33649 x 5 的矩阵, 表示所有不同的组合情况 |
|
S******y
发帖数: 1123 | 2 > library(gtools)
> x <- rnorm(1000)
> quartiles <- quantcut( x )
> table(quartiles)
quartiles
[-3.41,-0.704] (-0.704,0.0272] (0.0272,0.732] (0.732,2.94]
250 250 250 250 |
|
i*e
发帖数: 352 | 3 GWAS的东西你应该看看PLINK和GTOOL先
再说有必要把genotype 和 sample/phenotype合并在一个单独文件吗 |
|
