xyplot(
EW / 252 ~ RW / 100 | factor(ID),
data = d,
panel = function(x, y, ...) {
panel.points(x, y, col = "#0000ff40");
panel.abline(0, 1, col = "#ff4000")
}
)
table()
でカウントしておくか
……
やはり,
この 3 個体ね.
> table(d$ID, d$EW / 252 > d$RW / 100)
FALSE TRUE
107 35 0
108 35 0
108.2 35 0
110.2 35 0
110.3 35 0
110.5 35 0
110.6 35 0
110.7 35 0
113 35 0
113.2 35 0
220 35 0
408 35 0
429 21 14
439 30 5
685 35 0
702 35 0
703 35 0
731 22 13
753 35 0
757 35 0