Убрать крайние точки в mvrnorm и plot3d

Question

Я хотел бы использовать plot3d, чтобы визуализировать следующее и сделать сферы более плотными, но при этом удаляются только крайние точки сферы1, но не последующие сферы. Фактически, сфера 2 и сфера 3 исчезли.

      library('rgl') 
library('MASS') 

mu1 <- c(0,0,0) 
mu2 <- c(3,3,3) 
mu3 <- c(6,6,6) 

covmat <- diag(c(1,1,1))

sphere1 <- mvrnorm(n = 1500, mu = mu1, Sigma = covmat)
sphere1 <- as.data.frame(sphere1)
sphere1 <- subset(sphere2, X1^2+X2^2+X3^2<6)

sphere2 <- mvrnorm(n = 1500, mu = mu2, Sigma = covmat)
sphere2 <- as.data.frame(sphere2)
sphere2 <- subset(sphere2, X1^2+X2^2+X3^2<6)

sphere3 <- mvrnorm(n = 1500, mu = mu2, Sigma = covmat)
sphere3 <- as.data.frame(sphere3)
sphere3 <- subset(sphere3, X1^2+X2^2+X3^2<6)

plot3d(sphere1, type = 's', add = TRUE, radius = .3, col='red')
plot3d(sphere2, type = 's', add = TRUE, radius = .3, col='red')
plot3d(sphere3, type = 's', add = TRUE, radius = .3, col='red')