作者rekkler (放下需要勇气)
看板CCU_SS97
标题[情报]数统4.48(b)r code
时间Wed Dec 10 23:02:05 2008
x=runif(100000,-4,4)
y=runif(100000,-4,4)
z=exp(-0.5*((x^2)*(y^2)+x^2+y^2-16*x-16*y))
scatterplot3d(x,y,z,type="l",angle=160)
因为normal(0,1)下,x和y在超过-3和3的probability已经接近0
所以x和y不能取太大才可以看到它的分布
从图中可以看出bimodal在x=y的左右2侧
证明可以从(a)的平均数下手
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 140.123.215.101
1F:→ rekkler:我看错了..好像跟X=Y没关系 12/10 23:19