「インデックス」に基づいてデータ フレームを分割し、比較のために各グループを並べてプロットします。
Grp <- split(TOC, TOC$Index)
$`1`
Site Index depth_ft TOC_mg.g IC_mg.g
1 Z 1 5 12 NA
$`1`
Site Index depth_ft TOC_mg.g IC_mg.g
1 A 2 2 11 NA
...
...
1つずつ行けば、データを簡単にプロットできます。
plot(Grp$`1`$TOC_mg.g, Grp$`1`$depth_ft)
しかし、「lapply」を使用して一度にすべてのグループをプロットしたい場合、そこにデータ ポイントがない XY 軸が表示されます ( https://docs.google.com/file/d/0B6GUNg-8d30vdmZBMVhKVlR0TkE/edit?usp=共有)!誰が何が間違っているのか教えてもらえますか??
#plot
par(mfrow=c(1,5))
lapply(1:length(Grp), function(i)
plot(Grp$`i`$TOC_mg.g, Grp$`i`$depth_ft, ylim=c(0, max(TOC$depth_ft)),
xlim= c(min(TOC$TOC_mg.g, na.rm=T), max(TOC$TOC_mg.g, na.rm=T)), lwd=2, col=2 ))