Выбор первого ненулевого значения из каждого столбца в кадре данных

Question

У меня есть кадр данных 10x100 с именем CoeNIST. Строки имеют порядок значимости (то есть значение в строке 1 более важно, чем значение в строке 2), и каждый столбец представляет собой другой образец. Я хотел бы извлечь только самое значительное ненулевое значение, то есть первое ненулевое значение для каждого образца.

Вот образец из первых 9 столбцов CoeNIST.

> CoeNIST[,1:9] 
     1  2  3  4 5 6 7  8 9 
1  0 352232  0  0 0 0 0 28733 0 
2 332829  0  0 380109 0 0 0 380343 0 
3  0  0  0 380111 0 0 0 380409 0 
4  0  0  0 380101 0 0 0  0 0 
5  0  0 299211 380112 0 0 0  0 0 
6  0  0  0 380103 0 0 0  0 0 
7  0  0  0 380100 0 0 0 71899 0 
8  0  0  0 24812 0 0 0  0 0 
9  0  0  0  0 0 0 0 380410 0 
10  0 332958  0  0 0 0 0 380440 0 

И вот что я хотел бы результат, чтобы выглядеть

> NIST 
[1] 332829 352232 299211 380109 NA NA NA 28733 NA 

ИЛИ ... как список ...

> NIST 
[[1]] 
[1] 332829 

[[2]] 
[1] 352232 

[[3]] 
[1] 299211 

[[4]] 
[1] 380109 

[[5]] integer(0) 

[[6]] integer(0) 

[[7]] integer(0) 

[[8]] 
[1] 28733 

[[9]] integer(0) 

Answer 1

CoeNIST <- read.table(header=TRUE,text=" 
1  2  3  4 5 6 7  8 9 
1  0 352232  0  0 0 0 0 28733 0 
2 332829  0  0 380109 0 0 0 380343 0 
3  0  0  0 380111 0 0 0 380409 0 
4  0  0  0 380101 0 0 0  0 0 
5  0  0 299211 380112 0 0 0  0 0 
6  0  0  0 380103 0 0 0  0 0 
7  0  0  0 380100 0 0 0 71899 0 
8  0  0  0 24812 0 0 0  0 0 
9  0  0  0  0 0 0 0 380410 0 
10  0 332958  0  0 0 0 0 380440 0") 

Я бы описать вашу проблему, как «выбор первого ненулевого значения в каждом столбце.» Мое решение дает NA значению, когда Есть только нули в колонке ...

apply(CoeNIST,2,function(x) (x[x>0])[1]) 
##  X1  X2  X3  X4  X5  X6  X7  X8  X9 
## 332829 352232 299211 380109  NA  NA  NA 28733  NA 

Answer 2

CoeNIST = matrix(c(0, 352232, 0, 0, 0, 0, 0, 28733, 0, 332829, 0, 0, 380109, 0, 0, 0, 380343, 0, 0, 0, 0, 380111, 0, 0, 0, 380409, 0, 0, 0, 0, 380101, 0, 0, 0, 0, 0, 0, 0, 299211, 380112, 0, 0, 0, 0, 0, 0, 0, 0, 380103, 0, 0, 0, 0, 0, 0, 0, 0, 380100, 0, 0, 0, 71899, 0, 0, 0, 0, 24812, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 380410, 0, 0, 332958, 0, 0, 0, 0, 0, 380440, 0), nrow=10, ncol=10, byrow=T) 

> CoeNIST 
     [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] 
[1,]  0 352232  0 0 0  0  0 28733  0 332829 
[2,]  0  0 380109 0 0  0 380343  0  0  0 
[3,]  0 380111  0 0 0 380409  0  0  0  0 
[4,] 380101  0  0 0 0  0  0  0 299211 380112 
[5,]  0  0  0 0 0  0  0  0 380103  0 
[6,]  0  0  0 0 0  0  0 380100  0  0 
[7,]  0 71899  0 0 0  0 24812  0  0  0 
[8,]  0  0  0 0 0  0  0  0  0 380410 
[9,]  0  0 332958 0 0  0  0  0 380440  0 
[10,]  0 352232  0 0 0  0  0 28733  0 332829 

Это для макс. каждого столбца:

apply(CoeNIST, 2, function(x){x_max = max(x); if(x_max == 0) NULL else x_max})}) 

Это первый ненулевое значение в столбце:

apply(CoeNIST, 2, function(x){x_top_non_zero = min(which(x>0)); if(x_top_non_zero == Inf) {NaN} else {x_top_non_zero}})