(R) 3次元散布図・等高線マップ作成メモ

 3次元散布図（固定）

library(lattice)
cloud(Porosity~Depth*Inv_Vp,data=Dataset,groups=HCS)

 3次元散布図（マウスで回転）

library(rgl)
scatter3d(Dataset$Depth,Dataset$Inv_Vp,Dataset$Porosity,xlab="Depth",ylab="Inv_Vp",zlab="Porosity")

 観測点の分布を密度関数として等高線を作成

> require(MASS)
要求されたパッケージ MASS をロード中です
> x<-Dataset$Depth
> y<-Dataset$Inv_Vp
> d<-kde2d(x,y,n=50)
> image(d)
> contour(d,add=T)
> points(x,y,pch=20)

 不等間隔のデータを補完して等高線図を作成

> library(akima)
> dd<-interp(Dataset$Depth, Dataset$Inv_Vp, Dataset$Porosity, duplicate='mean')
> image(dd$x, dd$y, dd$z, xlab="Depth", ylab="Inv_Vp")
> contour(dd$x,dd$y, dd$z, add=TRUE)
> points(Dataset$Depth, Dataset$Inv_Vp, pch=20)