Функция R pnorm с разреженной матрицей

Question

Я хотел бы найти значения большой разреженной матрицы. Все элементы в этой матрице являются стандартными нормальными z-показателями. Я хочу использовать функцию pnorm, но я столкнулся с проблемой, что pnorm не поддерживает разреженную матрицу. Кроме преобразования разреженной матрицы в полную, есть ли другой более эффективный способ?

Любые предложения приветствуются!

Answer 1

Если это разреженная матрица, вы можете легко заменить 0 значений на pnorm (0 ..). Осталось рассчитать ненулевые значения, что вы можете сделать. Например, разреженная матрица:

data <- rnorm(1e5)
zero_index <- sample(1e5)[1:9e4]
data[zero_index] <- 0
mat <- matrix(data, ncol=100)
mat_sparse <- Matrix(mat, sparse=TRUE)

Создать матрицу с помощью pnorm для 0:

mat_pnorm <- matrix(pnorm(rep(0,length(mat_sparse))),ncol=ncol(mat_sparse))
nzData <- summary(mat_sparse)
mat_pnorm[as.matrix(nzData[,1:2])] <- pnorm(nzData$x)
all.equal(mat_pnorm,pnorm(mat))
[1] TRUE

Вы не указали, как вам нужны значения p, но вы можете легко получить это приведено в векторе вместо матрицы, которая использовалась выше.