5

我的剪切功能有问题。我有这种情况：

 codice
1 11GP2-0016
2 11GP2-0016
3 11GP2-0016
4  11OL2-074
5  11OL2-074

我想要一个新的变量“campione”，由变量“codice”分割，如下所示：

    codice campione
1 11GP2-0016    [1,3]
2 11GP2-0016    [1,3]
3 11GP2-0016    [1,3]
4  11OL2-074    (4,5]
5  11OL2-074    (4,5]

如何使用 cut 函数拆分“代码”创建一个变量，显示从 1 到 3 我有相同的代码，从 4 到 5 相同的代码等等？

我需要解决另一个问题。对于同样的问题，我想获得：

 codice campione
1 11GP2-0016    [11GP2-0016,11GP2-0016,11GP2-0016]
2 11GP2-0016    [11GP2-0016,11GP2-0016,11GP2-0016]
3 11GP2-0016    [11GP2-0016,11GP2-0016,11GP2-0016]
4  11OL2-074    (11OL2-074,11OL2-074]
5  11OL2-074    (11OL2-074,11OL2-074]

有什么解决方案可以做到这一点吗？

4

3 回答 3

3

使用您的数据：

d <- read.table(text = "1 11GP2-0016
2 11GP2-0016
3 11GP2-0016
4  11OL2-074
5  11OL2-074", row.names = 1, stringsAsFactors = FALSE)
names(d) <- "codice"

这是一个稍微复杂的示例，使用rle()：

drle <- with(d, rle(codice))

这给了我们运行长度codice：

> drle
Run Length Encoding
  lengths: int [1:2] 3 2
  values : chr [1:2] "11GP2-0016" "11OL2-074"

它是$lengths我用来创建两个指示的组件，即开始 ( ind1) 和结束 ( ind2) 位置：

ind1 <- with(drle, rep(seq_along(lengths), times = lengths) +
                     rep(c(0, head(lengths, -1) - 1), times = lengths))
ind2 <- ind1 + with(drle, rep(lengths- 1, times = lengths))

然后我把这些粘贴在一起：

d <- transform(d, campione = paste0("[", ind1, ",", ind2, "]"))

给予

> head(d)
      codice campione
1 11GP2-0016    [1,3]
2 11GP2-0016    [1,3]
3 11GP2-0016    [1,3]
4  11OL2-074    [4,5]
5  11OL2-074    [4,5]
于 2012-10-18T15:51:30.207 回答
3

这会做到的。如果需要，您可以添加括号/括号。

dat <- read.table(text='codice
1 11GP2-0016
2 11GP2-0016
3 11GP2-0016
4  11OL2-074
5  11OL2-074', header=TRUE)

within(dat, 
    campione <- with(rle(as.character(codice)), {
        starts <- which(! duplicated(codice))
        ends <- starts + lengths - 1
        inverse.rle(list(values=paste(starts, ends, sep=','), lengths=lengths))
    })
)

#       codice campione
# 1 11GP2-0016      1,3
# 2 11GP2-0016      1,3
# 3 11GP2-0016      1,3
# 4  11OL2-074      4,5
# 5  11OL2-074      4,5
于 2012-10-18T15:57:47.027 回答
2

另一种方法是使用rank：

left <- rank(factor(d$codice), ties.method = "min")
right <- rank(factor(d$codice), ties.method = "max")
d$campione <- paste("[", left, ",", right, "]", sep = "")
于 2012-10-18T22:46:12.683 回答