[情報] 數統4.48(b)r code

[rekkler (放下需要勇氣)]

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
→ rekkler:我看錯了..好像跟X=Y沒關係