Data: download

#Read input data
b.data <- read.table(“3dscatterogram.csv”, header=T,fill=T)
# call scatterplot3d package
library(scatterplot3d)
attach(b.data)
# simple 3d scatterplot
scatterplot3d(SOI_O,Rainfall,SLP_O, main=”3D Scatterplot”)

# 3d scatterplot with Fit planes
s3d <-scatterplot3d(SOI_O,SLP_O,Rainfall, pch=16, highlight.3d=TRUE, type=”h”, main=”3D Scatterplot”)
fit <- lm(Rainfall ~ SOI_O+SLP_O)
s3d$plane3d(fit)