ピークが2つ以上ある分布をRで作成してみた。関数の合成です。そしてその自作分布に従う乱数を発生させました。 まず下のような正規分布の関数fxとgxを設定します。平均と分散はそれぞれ異なる値とします。 fx <- function(x) dnorm(x, mean=1, sd = 1)gx <-…

Rで遊ぶのがすっかり楽しくなってきました。 その１に続いて、ファットヘッドミノーの生命表を使います（Miller and Ankley, 2004）。今度は初期個体数による絶滅確率の違いを検討してみます。その１ではw(9,2,0)の初期個体数11匹で計算していましたが、その…