> set.seed(1)
> n <- 50
> x <- runif(n, 0, 100)
> y <- runif(n, 0, 100)
>
> D <- as.matrix(dist(cbind(x, y)))
>
> phi <- 3/50
> sigmasq <- 50
> tausq <- 20
> mu <- 150
>
> s <- (sigmasq*exp(-phi*D))
> w <-  mvrnorm(1, rep(0, n), s)
错误: 没有"mvrnorm"这个函数
> Y <- mvrnorm(1, rep(mu, n) + w, tausq*diag(n))
错误: 没有"mvrnorm"这个函数
> X <- as.matrix(rep(1, length(Y)))
错误于as.matrix(rep(1, length(Y))) : 找不到对象'Y'
>