[関東鉄道バス]
「つくばセンター」行き．
JR 牛久駅東口．
切符売場でてすぐの階段を右におりたところから乗る．
0800 駅発，
0815 森林総研着．
240 円．
時刻表
によると，
次のつくばセンター行バスは 8 時間後の 1600 発，
である．
まあ，
谷田部車庫行き
に乗ればいいんだけど．
glm()
で得られた結果を使って作図しているところです．
plot()，
で．
樹木のサイズを cex
で表現したり，
等級を pch や col
で区別できるようにする．
あと，
text(..., pos = 1)
で個体 ID を追加，
とか．
library(lattice)
してから xyplot(y ~ x | species, data = d)
というふうに樹種ごとの地図を作らせてみたり．
library(lattice)
xyplot(log(shoot) ~ log(root) | species, data = d)
てなかんじで．
「1:1 のセン」
なるものを追加しようとすると，
とたんに見とおしの悪い書きかたになり，
たとえば
xyplot(
log(shoot) ~ log(root) | species,
data = d,
panel = function(x, y) {
panel.points(x, y)
panel.abline(0, 1, col = "gray")
}
)
となってしまう．
