R_2変量正規分布 - RとLinuxと...

RとLinuxと...


R_2変量正規分布

_ Baayen より

Analyzing Linguistic Data

pp.105-

library(MASS)
x <- mvrnorm(n = 1000, mu = c(0,0), 
Sigma = cbind(c(1, 0.8),  c(0.8,1)))

# 単変量の場合の density 関数にあたるのが kde2d

persp(kde2d(x[,1], x[,2], n = 50),
     phi = 30, theta = 20, # angle の指定
     d = 10,
     col = "lightblue",
     shade = 0.75, ltheta = -100,
     border = NA,
     expand = 0.5,
     xlab = "X",ylab ="Y", zlab = "density")
mtext("bivariate standard normal", 3, 1)

kde2d 関数の引数 n が分かりにくいが,x,y 座標の数

 lmreg による画像
 
添付ファイル: filebaayen106.png 648件 [詳細]
 
Link: Rの備忘録(1827d)
Last-modified: 2007-09-22 (土) 17:03:38 (4106d)