표시기 변수에서 값을 얻는 방법

Question

먼저 20m X 30m 직사각형에서 두 점의 좌표를 모델화해야했습니다. 물론 이러한 점은 균일 분포 그래서 내 코드의 첫 번째 부분이다 따르

X1 <- runif(1,0,20) 
    X2 <- runif(1,0,30) 
    point1 <- c(X1,X2) 
    point1 

는 I 번째 포인트 ('POINT2')에 대해 동일한 코드를 사용하지만, 각각 Y1 및 Y2와 X1과 X2를 대체했다.

내가 그 두 점 사이의 거리를 찾아야했다 : 나는 점은 서로 10m에 5m 내에있는 이벤트 같이 정의하면 이제

distance <- sqrt(((X1-Y1)^2)+((X2-Y2)^2)) 

을, 나는 표시기 변수를 찾을 필요 이 사건의 이것은 내가있어 무엇을,하지만 난 그게 맞아 잘 모르겠어요 :

x=0 
    if (distance < 10 & distance > 5) 
    { 
    x=1 
    } 
    Z <- c(distance,x) 

나는 각각의 시뮬레이션에 값을 저장 얼마나, 이러한 명령을 1,000 번 반복하고, 최소 및 최대를 찾을 수 있다면 1000의 reps의 분리 가치?

Answer 1

그들을 호출 할 수 있습니다.

points <- cbind(X1, X2) 

상황에 당신이 "병렬"배열로 X1 및 X2 값을 유지하고 달성하는 비 치수 객체가 아닌 매트릭스를 정의하고 싶다고 제안 : 당신이 사용하고자하는 확신 해요.

points1 <- matrix(runif(2000), ncol=2) 
points1 <- matrix(runif(2000), ncol=2) 
dists <- rowSums((points1-points2)^2) 
Z <- dists <10 & dists >5 

Answer 2

닫기. ifelse 또는 단순히 벡터화 된 양식을 사용하여 이벤트를 정의하십시오. 1000 샘플의 경우, runif을 사용하여 1000 개의 샘플을 생성하면됩니다. 또한 c(X1 , X2) 할 필요가 없습니다, 당신은 당신이 X1과 X2를 사용하고 나중에 정의되지 않은 Y1을 참조 한순간에 ... 당신의 거리 계산에서 벡터로

#Just make 1000 calls to runif 
X1 <- runif(1000,0,20) 
X2 <- runif(1000,0,30) 
Y1 <- runif(1000,0,20) 
Y2 <- runif(1000,0,30) 
distance <- sqrt(((X1-Y1)^2)+((X2-Y2)^2)) 

head(distance) 
#[1] 9.050522 19.512849 10.413407 7.736564 2.742174 13.729397 

# gives 1 if points are within 5 to 10m of each other 
event <- ifelse (distance >= 5 & distance <= 10 , 1 , 0) 

#Or even better, from @GavinSimpson's comment just use a vectorised form (we use as.nuemric to turn TRUE/FALSE into 1 and 0, but you could leave that out if you wish) 
event <- as.numeric(distance >= 5 & distance <= 10) 

head(event) 
#[1] 1 0 0 1 0 0 

# And the minimum and maximum sepration distance of those events 
min(distance[distance >= 5 & distance <= 10]) 
#[1] 5.017296 
max(distance[distance >= 5 & distance <= 10]) 
#[1] 9.989868