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


トップ   編集 凍結 差分 バックアップ 添付 複製 名前変更 リロード   新規 一覧 単語検索 最終更新   ヘルプ   最終更新のRSS
Google
WWW を検索 OKADAJP.ORG を検索
Last-modified: 2015-03-01 (日) 01:15:59 (1718d)