pointwise()

pointwise()

旧modregパッケージ(現在はstatsに統合)で、loess()はあるんですが、pointwise()がなぜかなくなってしまったため、S-PLUSプログラムの移植で困ることがあるかもしれません。 簡単なものですが旧modreg(R-1.8.1に収録)のソースから引用します。

関数 pointwise() の解説

省略(please read source code)

pointwise <- function(results, coverage)
{
   fit <- results$fit
   lim <- qt((1 - coverage)/2, results$df, lower = FALSE) * results$se.fit
   list(fit = fit, lower = fit - lim, upper = fit + lim)
}

関数 pointwise() の引数

省略(please read source code)

関数 pointwise() の使用例

data("cars")
cars.lo <- loess(dist ~ speed, cars)
cars.pred <- predict(cars.lo, data.frame(speed = seq(5, 30, 1)), se = TRUE)
pointwise(cars.pred, coverage=0.90)

• エクスポートしてないのが意図的なのかは存じませんが,stats:::pointwise(コロン3つでstats環境を指定してpointwiseとします) で呼び出せます. -- なかま 2005-06-07 (火) 17:59:56
• おぉ、確かに呼び出せますね！ありがとうございます。計算式自体はおかしくないように思えるのですが、なぜかもともとinternal functionという扱いになっています... -- 投稿者 2005-06-07 (火) 18:17:27

最新の20件