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)


トップ   編集 凍結 差分 バックアップ 添付 複製 名前変更 リロード   新規 一覧 検索 最終更新   ヘルプ   最終更新のRSS
Last-modified: 2023-03-25 (土) 11:19:16