//このページが表示される方は、URLから「action=SOURCE&」を削除してみてください
[[Ｒ備忘録 - 記事一覧]]

!!!等高線を引く
*投稿者： みゅ
*カテゴリ： なし
*優先度： 普通
*状態： 完了
*日時： 2008年07月08日 09時39分35秒
//{{bugstate}}
!!内容
パッケージ【lattice】のcontourplot関数
!!!デモ
 library(lattice)
 require(stats)
 attach(environmental)
 ozo.m <- loess((ozone^(1/3)) ~ wind * temperature * radiation,
        parametric = c("radiation", "wind"), span = 1, degree = 2)
 w.marginal <- seq(min(wind), max(wind), length = 50)
 t.marginal <- seq(min(temperature), max(temperature), length = 50)
 r.marginal <- seq(min(radiation), max(radiation), length = 4)
 wtr.marginal <- list(wind = w.marginal, temperature = t.marginal,
         radiation = r.marginal)
 grid <- expand.grid(wtr.marginal)
 grid[, "fit"] <- c(predict(ozo.m, grid))
 contourplot(fit ~ wind * temperature | radiation, data = grid,
             cuts = 10, region = TRUE,
             xlab = "Wind Speed (mph)",
             ylab = "Temperature (F)",
             main = "Cube Root Ozone (cube root ppb)")
 detach()
!!!結果
{{ref_image Ｒ備忘録6_1.png}}
!!コメント
//{{comment}}