как свернуть данные по интервалам? - PullRequest
0 голосов
/ 15 октября 2019

Я хотел бы взять набор данных, подобный этому:

dat <- data.frame(pos = 1:120, state = c(rep("state1", 30), rep("state2",30), rep("state3",30), rep("state1", 30)))

И свернуть его в эту форму:

dat2 <- data.frame(start = seq(1,120,30), end = seq(30,120,30), state = c("state1","state2","state3","state1"))

В заключение, я хочу знать начало и конецкаждой категории, как они появляются в данных.

Ответы [ 2 ]

0 голосов
/ 15 октября 2019

используя базу R, вы можете использовать rle:

with( rle(as.character(dat$state)),
      data.frame(state=values,end = cumsum(head(lengths))->end,start = c(1,head(end,-1)+1)))
   state end start
1 state1  30     1
2 state2  60    31
3 state3  90    61
4 state1 120    91
0 голосов
/ 15 октября 2019

Вы можете использовать group_by из dplyr для этой цели. Код ниже:

require(dplyr)

dat1 = dat %>% 
  dplyr::arrange(state,pos) %>% 
  dplyr::mutate(occurence=0)

occurence = 0

for(i in 1:nrow(dat1)){
  if((i != 1) && ((dat1$pos[i] - dat1$pos[i-1])>1)){
    occurence = occurence + 1
  }
  dat1$occurence[i] = occurence
}

dat2 = dat1 %>% 
  dplyr::group_by(state,occurence) %>% 
  dplyr::summarise(start = min(pos,na.rm=T),
                   end = max(pos,na.rm=T)) %>% 
  dplyr::arrange(start)

Дайте мне знать, если он работает.

Вывод

# A tibble: 4 x 4
# Groups:   state [3]
  state  occurence start   end
  <fct>      <dbl> <int> <int>
1 state1         0     1    30
2 state2         1    31    60
3 state3         1    61    90
4 state1         1    91   120

При необходимости вы можете удалить столбец 'occurence'. Использование:

dat2 = dat2 %>% dplyr::select(-occurence)
Добро пожаловать на сайт PullRequest, где вы можете задавать вопросы и получать ответы от других членов сообщества.
...