Hôm nay mình mới thử vẽ biểu đồ với package
(plot3D) cho hình đẹp, sinh động với các góc nhìn đa chiều. Bạn hãy thử với
plot3D để cảm thấy sự khác biết khi vẽ biểu đồ thông thường. Mình lấy ví dụ sau
làm minh chứng.
#
install.packages(plot3D)
#
library(plot3D)
> scatter3D(x,y,z, phi=0,
bty="g",cex=0.7,xlab="D1.3,
cm",ylab="Hvn,m",zlab="M,m3",main="Fig.1",clab=c("M,
m3"))
# result Fig. 1
Hoặc
>
scatter3D(x,y,z, phi=0,
bty="g",type="h",cex=0.7,xlab="D1.3,
cm",ylab="Hvn,m",zlab="M,m3",main="Fig.2",clab=c("M,
m3"))
# result Fig. 2
# Tiếp theo
>
fit=lm(z~x+y)
>
grid.lines=16
>
grid.lines=26
>
x.pred=seq(min(x),max(x),length.out=grid.lines)
>
y.pred=seq(min(y),max(y),length.out=grid.lines)
>
xy=expand.grid(x=x.pred,y=y.pred)
>
z.pred=matrix(predict(fit,newdata=xy),nrow=grid.lines, ncol=grid.lines)
>
fitpoints=predict(fit)
>
scatter3D(x,y,z, phi=20,cex=0.7,xlab="D1.3, cm",ylab="Hvn,m",zlab="M,m3",main="Fig.2",clab=c("M,
m3"),ticktype="detailed",theta=20,surf=list(x=x.pred,
y=y.pred,z=z.pred, facets=NA, fit=fitpoints))
>
scatter3D(x,y,z, phi=20,cex=0.7,xlab="D1.3,
cm",ylab="Hvn,m",zlab="M,m3",main="Fig.3",clab=c("M,
m3"),ticktype="detailed",theta=20,surf=list(x=x.pred,
y=y.pred,z=z.pred, facets=NA, fit=fitpoints))
# result Fig. 3
# Tiếp theo
>
scatter3D(x,y,z,pch=16,ticktype="detailed",
theta=15,d=2,xlab="D1.3, cm", ylab="Hvn, m", zlab="M,
m3", main="Fig. 4", clab="M, m3")
# result Fig. 4