Включить список IRanges в качестве столбца в data.frame - PullRequest
Новая
       19

Включить список IRanges в качестве столбца в data.frame

1 голос
/ 13 февраля 2020

У меня есть некоторые данные, структурированные примерно так: 

x01 <- c("94633X94644Y95423X96130", "124240X124494Y124571X124714", "135654X135660Y136226X136786")

, которые я использую позже как объект IRanges, через несколько шагов, которые выглядят следующим образом: 

x02 <- sapply(x01,
              function(x) do.call(rbind,
                                  strsplit(strsplit(x,
                                                    split = "Y",
                                                    fixed = TRUE)[[1]],
                                           split = "X",
                                           fixed = TRUE)),
              simplify = FALSE,
              USE.NAMES = FALSE)

x03 <- sapply(x02,
              function(x) IRanges(start = as.integer(x[, 1L]),
                                  end = as.integer(x[, 2L])),
              simplify = FALSE,
              USE.NAMES = FALSE)

> x03
[[1]]
IRanges object with 2 ranges and 0 metadata columns:
          start       end     width
      <integer> <integer> <integer>
  [1]     94633     94644        12
  [2]     95423     96130       708

[[2]]
IRanges object with 2 ranges and 0 metadata columns:
          start       end     width
      <integer> <integer> <integer>
  [1]    124240    124494       255
  [2]    124571    124714       144

[[3]]
IRanges object with 2 ranges and 0 metadata columns:
          start       end     width
      <integer> <integer> <integer>
  [1]    135654    135660         7
  [2]    136226    136786       561

Теперь я хотел бы иметь возможность хранить x03 в виде столбца в data.frame с некоторой связанной информацией с чем-то простым, например: 

> x04 <- data.frame("col1" = 1:3,
                    "col2" = x01,
                    "col3" = x03)

Это неудивительно, что у меня различное количество строк, однако я чувствую, что я видел JSON импорт в R mimi c вид структуры, которую я хочу, где рваный список размещается в столбце data.frame. Это возможная операция?

1 Ответ

1 голос
/ 13 февраля 2020

Это очень хороший вопрос, я видел его раньше с другими объектами, такими как объекты данных, но я думаю, что вышеупомянутое не работает, потому что есть as.data.frame, который можно использовать в матрице, или IRanges , он испортит размеры и не вставит его (я могу быть совершенно не прав).

Один из вариантов - использовать тиббл: 

x04 = tibble::tibble(x01=x01,x02=x02,x03=x03)
# A tibble: 3 x 3
  a                           b                 c        
  <chr>                       <list>            <list>   
1 94633X94644Y95423X96130     <chr[,2] [2 x 2]> <IRanges>
2 124240X124494Y124571X124714 <chr[,2] [2 x 2]> <IRanges>
3 135654X135660Y136226X136786 <chr[,2] [2 x 2]> <IRanges>

x04$x03
[[1]]
IRanges object with 2 ranges and 0 metadata columns:
          start       end     width
      <integer> <integer> <integer>
  [1]     94633     94644        12
  [2]     95423     96130       708

[[2]]
IRanges object with 2 ranges and 0 metadata columns:
          start       end     width
      <integer> <integer> <integer>
  [1]    124240    124494       255
  [2]    124571    124714       144

[[3]]
IRanges object with 2 ranges and 0 metadata columns:
          start       end     width
      <integer> <integer> <integer>
  [1]    135654    135660         7
  [2]    136226    136786       561

Другой вариант: 

library(S4Vectors)
DataFrame(x01=x01,x02=List(x02),x03=IRangesList(x03))
                          x01                             x02
                  <character>                          <List>
1     94633X94644Y95423X96130     94633:94644,95423:96130,...
2 124240X124494Y124571X124714 124240:124494,124571:124714,...
3 135654X135660Y136226X136786 135654:135660,136226:136786,...
                          x03
                <IRangesList>
1     94633-94644,95423-96130
2 124240-124494,124571-124714
3 135654-135660,136226-136786
...