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)