R：使用特定于行的中断应用剪切

Question

我想将潜在分数矩阵转换为观测分数。

可以通过对原始矩阵应用 break points/thresholds 来实现，从而最终得到一个新的分类矩阵。这样做很简单，例如：

#latent variable matrix
true=matrix(c(1.45,2.45,3.45,
              0.45,1.45,2.45,
              3.45,4.45,5.45)
,ncol=3,byrow=TRUE)

#breaks for the cut function
br=c(-Inf,1,2,3,4,Inf)

#apply cut function to latent variable
observed=apply(true,c(1,2),cut,breaks=br,labels=FALSE,include.lowest=TRUE)

但是，我需要做的是对原始矩阵的每一行应用不同的中断。这些阈值存储在一个矩阵中：

#matrix of breaks for the cut function
br=matrix(c(-Inf,1,2,3,4,Inf,
             -Inf,1.5,2.5,3.5,4.5,Inf,
             -Inf,2,3,4,5,Inf)
,ncol=6,byrow=TRUE)

也就是说，br矩阵的第1行应该作为true矩阵第1行的中断并且仅针对该行，br的第2行是true的第2行的中断，等等

使用以下方法似乎不起作用：

for (i in 1:nrow(true)) {
  observed[i,]=apply(true[i,],c(1,2),cut,breaks=br[i,],labels=FALSE,include.lowest=TRUE)
}

你有什么想法吗？有什么方法可以将相应的 br 行应用到相应的 true 行并将其保存在观察到的同一行中？

非常感谢！

KH

Answer 1

一些函数式编程和 Map 可以解决问题：

splitLines = function(m) split(m, rep(1:nrow(m), ncol(m)))

do.call(rbind, Map(cut, splitLines(true), splitLines(br), labels=F, include.lowest=T))
#  [,1] [,2] [,3]
#1    2    3    4
#2    1    1    2
#3    3    4    5

Answer 2

对行数使用 sapply，（基本上只是隐藏 for 循环）给你你想要的：

values = sapply(1:nrow(true), function(i) 
     cut(true[i,], br[i,], labels=FALSE, include.lowest=TRUE)))
values = t(values)

不幸的是，我们需要一个额外的转置步骤才能以正确的方式获得矩阵。

关于你问题中的 for 循环，当你只是对一行进行子集化时，即 true[i,] 我们只得到一个向量。这会导致 apply 中断。为了避免向量，你需要一个额外的参数

true[i,, drop=FALSE]

R：使用特定于行的中断应用剪切

R: Apply cut using row-specific breaks

r

matrix

apply