Parcelle de plan 3D (véritable régression de la surface)

Je suis en train de simuler certaines données (x1 et x2 - mes variables explicatives), calculer y à l'aide d'une fonction spécifiée + bruit aléatoire et tracer les observations ET la véritable régression de la surface. Voici ce que j'ai à ce jour:

   set.seed(1)
   library(rgl)

   # Simulate some data 

   x1 <- runif(50)
   x2 <- runif(50)
   y <- sin(x1)*x2+x1*x2 + rnorm(50, sd=0.3)
   # 3D scatterplot of observations 
   plot3d(x1,x2,y, type="p", col="red", xlab="X1", ylab="X2", zlab="Y", site=5, lwd=15)

Maintenant, je ne suis pas sûr de savoir comment je peux ajouter de la "véritable" régression de l'avion. Je suis fondamentalement à la recherche de quelque chose comme de la courbe() où je peux brancher mon (vrai) modèle de la formule.

Merci!

Je pense que ce fil pourrait vous aider.
Peut-être ceci pourrait être d'intérêt.

OriginalL'auteur user2249626 | 2013-08-09

  1. 18

    Si vous voulais un avion, vous pouvez utiliser planes3d.

    Depuis votre modèle n'est pas linéaire, il n'est pas un avion: vous pouvez utiliser surface3d à la place.

    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, ...)
}
library(rgl)
f <- function(x1, x2)
  sin(x1) * x2 + x1 * x2
n <- 200
x1 <- 4*runif(n)
x2 <- 4*runif(n)
y <- f(x1, x2) + rnorm(n, sd=0.3)
plot3d(x1,x2,y, type="p", col="red", xlab="X1", ylab="X2", zlab="Y", site=5, lwd=15)
my_surface(f, alpha=.2 )

    Parcelle de plan 3D (véritable régression de la surface)

    Qui fonctionne merci! Est il possible que je peux sélectionner un point d'observation) et de dessiner une sphère autour de lui (à l'aide d'un rayon spécifié)? De préférence, un point qui a beaucoup d'observations autour de lui (qui finira par être dans le domaine).
    Une fois que vous connaissez autour de vous souhaitez une sphère, vous pouvez utiliser spheres3d, par exemple, i <- 1; spheres3d(x1[i], x2[i], y[i], radius=1, alpha=.5).
    savez-vous quelle est la différence entre persp3d et surface3d? regarde comme ce dernier avait add=TRUE par défaut.
    surface3d est conçu pour ajouter les parcelles existantes tout en persp3d est conçu pour être le tracé initial de la commande.

    OriginalL'auteur Vincent Zoonekynd

  2. 1

    Excuses: ( je n'ai pas lu la question très carefllly et maintenant que je me suis précipité dans l'estimation lorsque vous voulais parcelle de la Vérité.)

    Ici est une approche pour l'estimation suivie par la surface à dessiner à l'aide loess:

    mod2 <- loess(y~x1+x2)
grd<- data.frame(x1=seq(range(x1)[1],range(x1)[2],len=20), 
                 x2=seq(range(x2)[1],range(x2)[2],len=20))
grd$pred <- predict(mod2, newdata=grd)
grd <- grd[order(grd$x1,grd$x2),]
x1 <- unique(grd$x1)
x2 <- unique(grd$x2)   # shouldn't have used y
surface3d(x1, x2, z=matrix(grd$pred,length(x1),length(x2)) )

    Parcelle de plan 3D (véritable régression de la surface)

    Merci pour votre réponse, néanmoins, il est utile. Pouvez-vous me dire pourquoi vous avez commandé la grille et pourquoi vous avez utilisé unique? Cheers!

    OriginalL'auteur 42-