데이터 프레임의 RGL 표면 플롯

Question

나는 scatter3d()과 Rcmdr을 사용하여 멋진 플롯을 만들었습니다. 이 플롯에는 두 개의 표면 매끄러움이 있습니다. 이제 에의 (즉, 관측치를 생성하는 함수로 정의 된 표면이 노이즈 성분을 뺀 표면)을 추가하여이 플롯에을 추가하고 싶습니다. 여기

지금까지 내 코드입니다 :

다른 스레드에 주어진 솔루션은 다음 함수를 정의하는 것입니다

library(car) 
set.seed(1) 

n <- 200 # number of observations (x,y,z) to be generated 
sd <- 0.3 # standard deviation for error term 
x <- runif(n) # generate x component 
y <- runif(n) # generate y component 

r <- sqrt(x^2+y^2) # used to compute z values 

z_t <- sin(x^2+3*y^2)/(0.1+r^2) + (x^2+5*y^2)*exp(1-r^2)/2 # calculate values of true regression function 

z <- z_t + rnorm(n, sd = sd) # overlay normally distrbuted 'noise' 
dm <- data.frame(x=x, y=y, z=z) # data frame containing (x,y,z) observations 
dm_t <- data.frame(x=x,y=y, z=z_t) # data frame containing (x,y) observations and the corresponding value of the *true* regression function 

# Create 3D scatterplot of: 
# - Observations (this includes 'noise') 
# - Surface given by Additive Model fit 
# - Surface given by bivariate smoother fit 

scatter3d(dm$x, dm$y, dm$z, fit=c("smooth","additive"), bg="white", 
      axis.scales=TRUE, grid=TRUE, ellipsoid=FALSE, xlab="x", ylab="z", zlab="y") 

: 이것은 그러나 오류를 얻을 수

my_surface <- function(f, n=10, ...) { 
    ranges <- rgl:::.getRanges() 
    x <- seq(ranges$xlim[1], ranges$xlim[2], length=n) 
    y <- seq(ranges$ylim[1], ranges$ylim[2], length=n) 
    z <- outer(x,y,f) 
    surface3d(x, y, z, ...) 
} 

f <- function(x, y) 
    sin(x^2+3*y^2)/(0.1+r^2) + (x^2+5*y^2)*exp(1-r^2)/2 

my_surface(f, alpha=0.2) 

을, 말 (이후 독일어 번역 이건 내 시스템 언어입니다. 사과하겠습니다.)

Error in outer(x, y, f) : 
    Dimension [Product 100] does not match the length of the object [200] 

그런 다음 다른 대안을 시도했습니다. ch :

x <- seq(0,1,length=20) 
y <- x 
z <- outer(x,y,f) 
surface3d(x,y,z) 

이것은 플롯에 서페이스를 추가하지만 전혀 보이지 않습니다 (예 : 관측치가 가까이 있지 않습니다.) 다음은 실제 표면이 보이는 모양입니다 (이것은 분명히 잘못되었습니다) : 고마워!

나는 문제가 사실로 확장 할 수있다 생각합니다. 여기에 평면 z = x + y에 앉아있는 몇 가지 점을 만들었습니다.

은 아마 당신의 일이 나를 도울 수 :

library(car) 

n <- 50 
x <- runif(n) 
y <- runif(n) 
z <- x+y 

scatter3d(x,y,z, surface = FALSE) 

f <- function(x,y) 
    x + y 

x_grid <- seq(0,1, length=20) 
y_grid <- x_grid 
z_grid <- outer(x_grid, y_grid, f) 
surface3d(x_grid, y_grid, z_grid) 

이 나에게 다음과 같은 플롯을 제공합니다 : 그럼 위의 내 방법을 사용하여 비행기를 플롯하려고하는 것을 계속?

Answer 1

의 기능은 플롯하기 전에 데이터를 재 측정하므로 본질적으로 모든 rgl 플로팅 기능 (surface3d 포함)과 호환되지 않습니다.

모든 rgl 기능을 사용하여 원하는 것과 같은 그림을 얻을 수 있습니다. scatter3d 대신 plot3d(x, y, z)을 사용하지만 물론 car- 스타일 축 대신 rgl- 스타일 축이 있습니다.