パッケージgamから使用しているモデルをフィッティングしています。mgcv結果を に保存しています。modelこれまでのところ、 を使用して滑らかなコンポーネントを見てきましたplot(model)。私は最近使い始めlattice、その出力が気に入っています。を使用してこれらのグラフをプロットすることは可能latticeですか?
この私のデータセット: https://gist.github.com/plxsas/fcef4a228c18c772b4f3
m2<- gam(TotalInd ~ s(dayinyear, by=as.numeric(Site=="1"), bs="cr")
+s(dayinyear, by=as.numeric(Site=="2"), bs="cr") + s(dayinyear,
by=as.numeric(Site=="3"), bs="cr"), random=list(Replicate=~ 1), data=data)
latticeこのモデルを、私の 3 つのサイトをよりスムーズに表す 3 つのパネルを含むパッケージでプロットするにはどうすればよいですか?
dayinyearまた、適切な月の形式 (データの最初の列)の代わりに を使用していることにお気付きかもしれません。これは、一般化加法モデルがカテゴリ変数を処理しないためです。ただし、グラフの時間を月の名前で表現したいと思います(最初の列のように)。latticeプロットでそれを進める方法を知っている人はいますか?