ryamada22 https://blog.hatena.ne.jp/ryamada22/ ryamadaの遺伝学・遺伝統計学メモ https://ryamada22.hatenablog.jp/ R 多項分布 サンプリング 1細胞 NGS 複数の細胞からのデータがある 細胞数が分かっている リード数を細胞数に分けて、疑似１細胞発現データを作りたい sample()関数を使って疑似作成できることを示し # 検体数 N # ある遺伝子のNGSリード数が Mだった # m1 + ... + mN = M; mi >= 0となるような{mi}をランダムに生成したい N <- 100 M <- 50 # M本のリードを1:Nの検体にランダムに(等確率で)割り付ける s <- sample(1:N,M,replace=TRUE) # 検体IDごとに集計する t <- tabulate(s) # 関数にする my.divide.reads <… 190 <iframe src="https://hatenablog-parts.com/embed?url=https%3A%2F%2Fryamada22.hatenablog.jp%2Fentry%2F20190903%2F1567469235" title="疑似1細胞発現データを作る - ryamadaの遺伝学・遺伝統計学メモ" class="embed-card embed-blogcard" scrolling="no" frameborder="0" style="display: block; width: 100%; height: 190px; max-width: 500px; margin: 10px 0px;"></iframe> Hatena Blog https://hatena.blog 2019-09-03 09:07:15 rich https://ryamada22.hatenablog.jp/entry/20190903/1567469235 1.0 100%