2016-01-21 103 views
0

以下是我認爲我在與使用gmp包麻煩日常使用兩種方案:R:困難操縱bigz矢量/列表

正常行爲:

### example 1 (changing a subset of a vector) 
> v1 <- rep(1,10) 
> v1[3:7] <- 2:6 
> v1 
[1] 1 1 2 3 4 5 6 1 1 1 

### example 2 (flattening a list to a vector) 
> mylist1 <- list(1:5,6:10) 
> unlist(mylist1) 
[1] 1 2 3 4 5 6 7 8 9 10 

GMP當量:

### example 1 (when changing a subset, the vector is lost) 
> v2 <- as.bigz(rep(1,10)) 
> v2[3:7] <- as.bigz(2:6) 
> v2 
Big Integer ('bigz') 10 x 1 matrix: 
    [,1] 
[1,] 1 
[2,] 1 
[3,] 2 
[4,] 3 
[5,] 4 
[6,] 5 
[7,] 6 
[8,] 1 
[9,] 1 
[10,] 1 

### example 2 (strange behavior when flattening) 
> mylist2 <- list(as.bigz(1:5),as.bigz(6:10)) 
> unlist(mylist2) 
[1] 05 00 00 00 01 .... ##### many more "raw"-like entries 

這將是如果有GMP equival不錯(即unlist.bigz),因爲我們知道在R中「增長」的對象可能過於緩慢。大多數時候,我被迫宣佈空bigz載體,不斷添加到它想:

> v3 <- as.bigz(rep(1,2)) 
> v3 <- c(v3, as.bigz(2:6)) 
> v3 <- c(v3, as.bigz(rep(1,3))) 
> v3 
Big Integer ('bigz') object of length 10: 
    [1] 1 1 2 3 4 5 6 1 1 1 
> 

我通常聲明與給定模式和長度的載體,但是當我試圖用bigz做到這一點矢量如:v4 <- as.bigz(vector(length=n)),然後嘗試填充矢量的子集,我得到像上面那樣的結果(即有一列的矩陣)。這是非常不可取的。

回答

1

可以使用c功能結合bigz對象的列表:

bigz.list <- list(
    as.bigz(1), 
    as.bigz(1), 
    as.bigz(2:6), 
    as.bigz(rep(1,3))) 

> do.call(c, bigz.list) 
Big Integer ('bigz') object of length 10: 
[1] 1 1 2 3 4 5 6 1 1 1