Question

用R语言估计参数值 请帮忙解释下面这段R语言程序每句的含义，谢谢！

> N <- 10^6
> x1 <- numeric(N)
> x2 <- numeric(N)
> n <- 0
> pai <- 0
> for(i in 1:N){
+ x1 <- runif(1,-1,1)
+ x2<- runif(1,-1,1)
+ if( x1^2+x2^2 <= 1 ){
+ n<- n+1}
+ }
> pai <- 4*n/N
> pai

Answer 1

从第一行开始，N直到pai，都是赋值语句。其中x1和x2是长度为N，类型为十进制小数的向量。

runif是生成一个随机数，取值在-1到1之间。for循环语句生成具体的两个向量，即x1和x2，其中的每个数都用runif来生成。

紧接之后的if语句用来计数n，其条件为如果对于数值x1[i]和x2[i]，如果点（x1[i]，x2[i]）位于单位圆内，n就加1。一共循环N次。

最后就得出值pai = 4*n/N。

看下来不像是一个参数值估计程序。

Answer 2

jjjjjjjjjjjjjjjjjjjjjjjjjjjjjjjjjjjjjjjjjjjj