0

我正在使用一个数据框，其中每个观察值都链接到一个特定的 ID，并且我有一组变量来定义“值”，就好像我有一个因子变量一样。但是，“单元格”中的值是频率。这是一个简化版本：

ID  1  2  3
A   2  3  2
B   1  4  1

我想获得两个扩展频率的向量，以便我可以计算每个 ID 的插值中位数。也就是说，我想要以下形式的东西：

A  B
1  1
1  2
2  2
2  2
2  2
3  3
3

该psych包有一个函数，该函数interp.median可以获取每个向量并返回每个 ID 的插值中位数，我想将其作为新变量包含在原始数据帧中。我检查了vcdExtra可能通过它的expand.dft功能做到这一点的包，但我不确定它究竟是如何工作的。

任何帮助将不胜感激！

编辑：为了进一步完善，interp.median如果最终结果是一个数据框，并在最后填充 NA，则效果最好。也就是说，某种形式：

A  B
1  1
1  2
2  2
2  2
2  2
3  3
3  NA
4

3 回答 3

0

如果dat是数据集

  lst <- by(dat[,-1], dat[,1], function(x) rep(seq_along(x), x))
  lst
  #dat[, 1]: A
  #[1] 1 1 2 2 2 3 3
  #------------------------------------------------------------ 
 #dat[, 1]: B
 #[1] 1 2 2 2 2 3

 indx <- max(sapply(lst,length))
 dat2 <- do.call(data.frame,lapply(lst, function(x) c(x,rep(NA,indx-length(x)))))
 dat2
 #  A  B
 #1 1  1
 #2 1  2
 #3 2  2
 #4 2  2
 #5 2  2
 #6 3  3
 #7 3 NA

或者

  lst2 <- lapply(split(dat[,-1], dat$ID), function(x) rep(seq_along(unlist(x)), unlist(x)))

  do.call(data.frame,lapply(lst2, function(x) c(x,rep(NA,indx-length(x)))))

数据

 dat <-  structure(list(ID = c("A", "B"), `1` = c(2L, 1L), `2` = 3:4, 
`3` = c(2L, 1L)), .Names = c("ID", "1", "2", "3"), class = "data.frame", row.names = c(NA, 
 -2L))
于 2014-08-27T18:20:18.267 回答
0

这里有一种方法：

# your data
df <- data.frame(ID=c(1,2,3), A=c(2,3,2), B=c(1,4,1))

# function to repeat each ID a given number of times,
# as specified in 'colname' of df
rep_id <- function(colname) {
  unname(unlist(apply(df[, c('ID',colname)], 1, function(x) rep(x[1], x[2]))))
}

# apply this function to all columns (except the first, which is ID)
sapply(names(df)[-1], rep_id)

产量：

$A
[1] 1 1 2 2 2 3 3

$B
[1] 1 2 2 2 2 3
于 2014-08-27T18:29:25.593 回答
0

样本数据：

df <- read.table(text="
ID  1  2  3
A   2  3  2
B   1  4  1", header=TRUE, check.names=FALSE)

使用apply：

(newlist <- apply(df[2:4], 1, function(x) rep(names(x), x)))
#[[1]]
#[1] "1" "1" "2" "2" "2" "3" "3"
#
#[[2]]
#[1] "1" "2" "2" "2" "2" "3"

names(newlist) <- df$ID
#$A
#[1] "1" "1" "2" "2" "2" "3" "3"
#
#$B
#[1] "1" "2" "2" "2" "2" "3"

这会输出字符，但您可以像这样输出数字：

newlist <- apply(df[2:4], 1, function(x) rep(as.numeric(names(x)), x))
names(newlist) <- df$ID

编辑：

为了解决 OP 将向量放入 data.frame 并用 NA 填充的新请求，请在运行上述任一选项后调用它：

newlist <- sapply(newlist, function(x) x[1:max(sapply(newlist, length))])
#     A  B
#[1,] 1  1
#[2,] 1  2
#[3,] 2  2
#[4,] 2  2
#[5,] 2  2
#[6,] 3  3
#[7,] 3 NA
于 2014-08-27T18:32:47.387 回答