Streudiagramm in R

Question

Ich unterrichte mich selbst, indem ich einige Aufgaben in meiner Monte-Carlo-Klasse re-mache. Ich baute einen RANDU-Generator in R, und ich bin für alle Triplets in der Sequenz zuständig (u_i, u_ {i + 1}, u_ {i + 2}), wobei 0,5 < = u_ {i + 1} < = 0,51 , plot u_i gegen u_ {i + 2}.

Ich habe folgendes in MATLAB während der Klasse und ich bekomme diese:

Hier ist der Code, den ich in Latex schrieb, beachten Sie nur, dass die in der Variablen gespeicherten Daten u einfach einheitliche Zufallszahlen zwischen 0 und 1.

Nun zu meiner Frage, bekomme ich diesen Fehler in R

Error in if (0.5 <= u[i + 1] & u[i + 1] <= 0.51) { : 
    missing value where TRUE/FALSE needed 

Ich bin nicht sicher, was das Problem ist hier, hier ist meine volle R-Code:

n = 20002 
x <- 1:n 
x[1] = 1 

for(i in 2:n){ 
    x[i] = (65539*x[i-1]) %% 2^31 
} 

u <- 1:n 
u[1] = 1/(2^31) 
for(i in 2:n){ 
    u[i] = x[i]/(2^31) 
} 


for(i in 1:length(u)){ 
    if(0.5 <= u[i+1] & u[i+1] <= 0.51){ 
    plot(u[i],u[i+2]) 
    } 
} 

Irgendwelche Vorschläge sehr geschätzt werden.

Answer 1

Ich denke, das Problem ist für Ihre letzte Beobachtung in u. Wie für Ihre for Schleife, die Sie bitten, u[i+1] auszuwerten. Versuchen Sie dies:

for(i in 1:length(u)-1){ 
    if(0.5 <= u[i+1] & u[i+1] <= 0.51){ 
    plot(u[i],u[i+2]) 
    } 
}