2011-11-09 13 views
5

Chcę opracować szkielet podobny do wykresu z nieliczbowym na osi X, Y, jakkolwiek numeryczny na osi Z.R: Model szkieletowy podobny do wykresu 3D dla zmiennych jakościowych

# mydata 
set.seed(123) 
yv <- rnorm(20, 10, 3) 
gen <- rep(paste("G", 1:5, sep= ""), 4) 
env <- c(rep(c("CA","MN","SD", "WI"), each = 5)) 
mdf <- data.frame(yv, gen, env) 

Próbowałem za pomocą kraty:

require(lattice) 
wireframe(yv,gen, env, data = mdf) 

Error in UseMethod("wireframe") : 
    no applicable method for 'wireframe' applied to 
    an object of class "c('double', 'numeric')" 

Wszelkie sugestie mile widziane.

Odpowiedz

5

To wydaje się działać:

set.seed(123) 
mdf <- data.frame(yv=rnorm(20, 10, 3), 
        gen=rep(paste("G", 1:5, sep= ""), 4), 
        env=c(rep(c("CA","MN","SD", "WI"), each = 5))) 
library(lattice) 
wireframe(yv~gen*env,data=mdf,scales=list(arrows=FALSE)) 

enter image description here

+0

+1 więcej, jeśli można zrobić to wyglądać żurawia origami. –

3

Najprostszym sposobem, aby korzystać z funkcji w lattice jest użycie interfejsu formuły.

Zakładając, że yv to zmienna niezależna:

wireframe(yv ~ gen + env, data = mdf) 

enter image description here

Powiązane problemy