R での非心 t 分布

統計

　私が使っているのは Version 2.3.1 ですが、非心 t 乱数の関数 rt() が間違ってます。

rt <-
function (n, df, ncp = 0) 
{
  if (ncp == 0) .Internal(rt(n, df))
  else rnorm(n, ncp)/(rchisq(n, df)/sqrt(df))
}

　間違ってますね。ずいぶん悩みました。
　どこかに連絡すべきでしょうか？
　それにしても、まだこのレベルのバグが R には残っているのかと思うと怖いです。確かに、「仕事」ではまだ使いにくいです。ちゃんとバリデーションしないと。


　私のほうが間違っている？？？