For-loop vs while döngü içinde R

Question

R. çalışırken ilginç bir şey fark ettim For döngüsü ve while-döngüsü kullanılarak uygulanan 1'den N'ye kadar kareleri hesaplayan basit bir programım olduğunda, davranış aynı. (Bu durumda vektörleştirmeyi umursamıyorum veya işlevleri uygula).

fn2 <- function (N) 
{ 
    i=1 
    while(i <= N) { 
     y <- i*i 
     i <- i + 1 
    } 
} 

fn1 <- function (N) 
{ 
    for(i in 1:N) { 
     y <- i*i 
    } 
} 

VE

sonuçlar şöyle:

system.time(fn1(60000)) 
    user system elapsed 
    2.500 0.012 2.493 
There were 50 or more warnings (use warnings() to see the first 50) 
Warning messages: 
1: In i * i : NAs produced by integer overflow 
. 
. 
. 

system.time(fn2(60000)) 
    user system elapsed 
    0.138 0.000 0.137 

Şimdi döngü için-hızlıdır, benim tahminim çünkü orada ön tahsisi ve optimizasyonlar sonucunda olduğunu biliyoruz. Ama neden taşıyor?

GÜNCELLEME:

fn3 <- function (N) 
{ 
    i <- 1:N 
    y <- i*i 
} 
system.time(fn3(60000)) 
    user system elapsed 
    0.008 0.000 0.009 
Warning message: 
In i * i : NAs produced by integer overflow 

Yani belki onun bir korkak bellek sorunu: Şimdi vektörler ile başka bir yol çalışıyor? OS X'de 4Gb bellek ve R'deki tüm varsayılan ayarlar ile çalışıyorum. Bu 32 ve 64-bit sürümlerde gerçekleşir (bu süreler daha hızlıdır).

Alex

Answer 1

için 1 (yani, olan bir kayan nokta sayısı) sayısal değil, tam sayı olduğu ve 1:6000 sayısal ve tamsayıdır. Ancak

> as.integer(60000)*as.integer(60000) 
[1] NA 
Warning message: 
In as.integer(60000) * as.integer(60000) : NAs produced by integer overflow 

3,6 milyar kayan nokta kolayca sunulabilen: kare

> print(class(1)) 
[1] "numeric" 
> print(class(1:60000)) 
[1] "integer" 

60000 imzalanmış 32 bitlik tamsayı, dolayısıyla bir taşma hatası alıyorum gösterilemeyecek 3,6 milyar olduğu

:

> as.single(60000)*as.single(60000) 
[1] 3.6e+09 

bir kayan nokta gösterimi için dönüştürmek sizin for kodunu çözmek için

function (N) 
{ 
    for(i in as.single(1:N)) { 
     y <- i*i 
    } 
} 

Answer 2

değişken döngü için bir tamsayı dizisi, ve böylece en sonunda, bu do:

süre döngü bir kayan nokta sayısı oluşturma oysa

> y=as.integer(60000)*as.integer(60000) 
Warning message: 
In as.integer(60000) * as.integer(60000) : NAs produced by integer overflow 

.

Onun da nedeni bu durum farklı:

> seq(0,2,1) 
[1] 0 1 2 
> seq(0,2) 
[1] 0 1 2 

Bana inanmıyor musunuz?

> identical(seq(0,2),seq(0,2,1)) 
[1] FALSE 

için:

> is.integer(seq(0,2)) 
[1] TRUE 
> is.integer(seq(0,2,1)) 
[1] FALSE 

Answer 3

Ve zamanlaması hakkında:

fn1 <- function (N) { 
    for(i in as.numeric(1:N)) { y <- i*i } 
} 
fn2 <- function (N) { 
    i=1 
    while (i <= N) { 
     y <- i*i 
     i <- i + 1 
    } 
} 

system.time(fn1(60000)) 
# user system elapsed 
# 0.06 0.00 0.07 
system.time(fn2(60000)) 
# user system elapsed 
# 0.12 0.00 0.13 

Ve şimdi For döngüsü while-döngü daha hızlı olduğunu biliyoruz. Zamanlama sırasında uyarıları göz ardı edemezsiniz.