« [Octave]標準正規分布に従う確率変数(乱数)を得る | トップページ | [R]ベクトルの要素を逆順に並び替える »

2019年7月 1日 (月)

[R]階乗を計算する

factorial関数を使う。以下、実行例。なお、4の階乗(4!)は24である。

> factorial(4)
[1] 24

ガンマ関数の値を計算するgamma関数を使っても計算できる。ガンマ関数の以下の性質を利用する。

Γ(n + 1) = n! (nは整数かつn≧1)

> gamma(4 + 1)
[1] 24

factorial関数とgamma関数はその実行速度が異なる。gamma関数のほうが断然早いので、gamma関数の使用を推奨する。

> system.time(for (i in 1:1000000) factorial(4))
ユーザ システム 経過
0.39 0.00 0.39
> system.time(for (i in 1:1000000) gamma(4 + 1))
ユーザ システム 経過
0.11 0.00 0.11

« [Octave]標準正規分布に従う確率変数(乱数)を得る | トップページ | [R]ベクトルの要素を逆順に並び替える »

R(数値計算)」カテゴリの記事

コメント

コメントを書く

コメントは記事投稿者が公開するまで表示されません。

(ウェブ上には掲載しません)

« [Octave]標準正規分布に従う確率変数(乱数)を得る | トップページ | [R]ベクトルの要素を逆順に並び替える »

無料ブログはココログ

■■

■■■