R中如何实现这样的图

https://bbs.pinggu.org/thread-1242928-1-1.html

qoiqpwqr 发表于 2012-11-24 23:29
你是想要这样的label吗？

scatterplot3d(ma.dat,type="p", lwd=2, pch=46, box=F, scale.y=1, x.tickl ...

@qoiqpwqr
dxystata 的佣金 5 5分吧 呵呵。。。
我再做下和他发布的图 再近似些

require(scatterplot3d)

x1 <-1;x2 <-2;x3 <-3
n <- 10000
y1 <-rnorm(n,0,x1)
y2 <-rnorm(n,0,x2)
y3 <-rnorm(n,0,x3)

x <-c(rep(x1,n),rep(x2,n),rep(x3,n))
y <- c(y1,y2,y3)
z <-c(dnorm(y1,0,x1),dnorm(y2,0,x2),dnorm(y3,0,x3))
ma.dat <- data.frame(x=x,y=y,z=z)
var3d <- scatterplot3d(ma.dat,type="p",lwd=2,pch=46 ,box=F,x.ticklabs=c(1,NA,2,NA,3),scale.y=1)
#均值时的，概念密度函数值
var3d$points3d(c(1,2,3),c(0,0,0),c(1/sqrt(2*pi),1/(2*sqrt(2*pi)),1/(3*sqrt(2*pi))),col="blue",type="h",pch=" ")
#平行X，y为均值，z为无限接近0的线
var3d$points3d(c(1,2,3),c(0,0,0),c(0,0,0),col="blue",type="l",pch=" ")

dxystata,根据那个帖子，试了下 scatterplot3d 包

#查看均值相同，  标准差不同的（1,2,3） 正态分布情况
x1 <-1;x2 <-2;x3 <-3

y1 <-rnorm(10000,0,x1)
y2 <-rnorm(10000,0,x2)
y3 <-rnorm(10000,0,x3)

x <-c(rep(x1,10000),rep(x2,10000),rep(x3,10000))
y <- c(y1,y2,y3)
z <-c(dnorm(y1,0,x1),dnorm(y2,0,x2),dnorm(y3,0,x3))
ma.dat <- data.frame(x=x,y=y,z=z)
scatterplot3d(ma.dat,type="p",lwd=2,pch=46 ,box=F,xlim=c(1,2,3),scale.y=1)

@dxystata,@qoiqpwqr
其实这个图，有个问题，一直没有解决xlim的问题，我想表示为1,2,3.  而不是1.0,1.5,2.0，2.5,3.0，
想用xaxp=c(1,3,2) 也不行。看帮助，作者承认xlim上有bug

jmpamao 发表于 2012-11-24 22:34
@dxystata,@qoiqpwqr
其实这个图，有个问题，一直没有解决xlim的问题，我想表示为1,2,3.  而不是1.0,1.5,2 ...

你是想要这样的label吗？

scatterplot3d(ma.dat,type="p", lwd=2, pch=46, box=F, scale.y=1, x.ticklabs = c(1,NA,2,NA,3))

qoiqpwqr 发表于 2012-11-24 23:29
你是想要这样的label吗？

scatterplot3d(ma.dat,type="p", lwd=2, pch=46, box=F, scale.y=1, x.tickl ...

哇塞， 还可以这样。

可是 能把中间的 那个线条 也去掉吗？

jmpamao 发表于 2012-11-25 00:22
@qoiqpwqr
dxystata 的佣金 5 5分吧 呵呵。。。
我再做下和他发布的图 再近似些

图是你做的，币当然全部都要给你。

jmpamao 发表于 2012-11-25 00:22
@qoiqpwqr
dxystata 的佣金 5 5分吧 呵呵。。。
我再做下和他发布的图 再近似些

用grid = FALSE来去掉那些格子

学习了~

厉害厉害，保存学习