2014-11-04 35 views
0

经过很长一段时间后,我只是跳回到R中,我对某些事情要做的简单程度感到惊讶。我已经在3个数组:统计一个数组中满足简单条件的元素数

Xs = runif(N, min=-1, max=1); 
Ys = runif(N, min=-1, max=1); 
Rs = sqrt(Xs^2 + Ys^2); 

其中,显然,X和Y(一起)限定(-1,1)内的N个点和正方形R是定义这些点的距离矢量。

如果我想要计算Rs中小于或等于1的元素数,是否有一个简单的inl-line命令来做到这一点?

+1

'表(Rs <= 1)'或'sum(Rs <= 1)' – user20650 2014-11-04 22:23:38

+0

分号在R中也不做任何事 – 2014-11-04 22:37:00

+1

应该接近于N(统计意义上的)N * 0.7853982 = N * pi/4。 'pi/4 [1] 0.7853982 > N = 1000 > sum(Rs <= 1) [1] 786' – 2014-11-04 22:43:10

回答

0
sum(Rs <= 1) 

Rs <= 1产生逻辑向量。 TRUE等于1; FALSE等于0.

相关问题