［R］ガンマ関数の値の計算（「入門 統計解析 －医学・自然科学編」（東京図書）pp.120-121）

ガンマ関数の値を得るにはgamma関数を使う。

> alpha <- seq(0.1, 3.0, 0.1)
> ggamma <- gamma(alpha)
> dtf <- data.frame(alpha, ggamma)
> print(dtf)
   alpha    ggamma
1    0.1 9.5135077
2    0.2 4.5908437
3    0.3 2.9915690
4    0.4 2.2181595
5    0.5 1.7724539
6    0.6 1.4891922
7    0.7 1.2980553
8    0.8 1.1642297
9    0.9 1.0686287
10   1.0 1.0000000
11   1.1 0.9513508
12   1.2 0.9181687
13   1.3 0.8974707
14   1.4 0.8872638
15   1.5 0.8862269
16   1.6 0.8935153
17   1.7 0.9086387
18   1.8 0.9313838
19   1.9 0.9617658
20   2.0 1.0000000
21   2.1 1.0464858
22   2.2 1.1018025
23   2.3 1.1667119
24   2.4 1.2421693
25   2.5 1.3293404
26   2.6 1.4296246
27   2.7 1.5446858
28   2.8 1.6764908
29   2.9 1.8273551
30   3.0 2.0000000
> plot(alpha, ggamma, type = "n")
> lines(alpha, ggamma)