Usandoscatterplot3d em R, estou tentando desenhar linhas vermelhas das observações para o plano de regressão:

wh <- iris$Species != "setosa"
x  <- iris$Sepal.Width[wh]
y  <- iris$Sepal.Length[wh]
z  <- iris$Petal.Width[wh]
df <- data.frame(x, y, z)

LM <- lm(y ~ x + z, df)
library(scatterplot3d)
G  <- scatterplot3d(x, z, y, highlight.3d = FALSE, type = "p")
G$plane3d(LM, draw_polygon = TRUE, draw_lines = FALSE)

Para obter o equivalente em 3D da seguinte imagem:

Em 2D, eu poderia usarsegments:

pred  <- predict(model) 
segments(x, y, x, pred, col = 2)

Mas em 3D fiquei confuso com as coordenadas.