0

2 つの列の完全なケースの数を合計するにはどうすればよいですか?

とc等しい：

      a  b
[1,] NA NA
[2,]  1  1
[3,]  1  1
[4,] NA  1

のようなものを適用する

rollapply(c, 2, function(x) sum(complete.cases(x)),fill=NA)

2この場合、単一の数値を取得したいと思います。これは、多くの列を持つ大規模なデータ セット用になるため、rollapply単純にsum(complete.cases(a,b)).

私はそれを考えすぎていますか？

ありがとう！

4

3 に答える 3

2

試しましたsum(complete.cases(x))か？

set.seed(123)
x <- matrix( sample( c(NA,1:5) , 15 , TRUE ) , 5 )
#     [,1] [,2] [,3]
#[1,]    1   NA    5
#[2,]    4    3    2
#[3,]    2    5    4
#[4,]    5    3    3
#[5,]    5    2   NA


sum(complete.cases(x))
#[1] 3

complete.cases()最初の 2 つの列のを見つけるには:

sum(complete.cases(x[,1:2]))
#[1] 4

そして、行列全体applyの a の 2 つの列に対して、matrix次のことができます。

#  Bigger data for example
set.seed(123)
x <- matrix( sample( c(NA,1:5) , 50 , TRUE ) , 5 )
#     [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10]
#[1,]    1   NA    5    5    5    4    5    2   NA    NA
#[2,]    4    3    2    1    4    3    5    4    2     1
#[3,]    2    5    4   NA    3    3    4    1    2     2
#[4,]    5    3    3    1    5    1    4    1    2     1
#[5,]    5    2   NA    5    3   NA   NA    1   NA     5

# Column indices
id <- seq( 1 , ncol(x) , by = 2 )
[1] 1 3 5 7 9
apply( cbind(id,id+1) , 1 , function(i) sum(complete.cases(x[,c(i)])) )
[1] 4 3 4 4 3

complete.cases()全体で行単位で機能するdata.frameか、データが欠落していない行をmatrix返します。TRUEマイナーなことはさておき、最も一般的に使用される関数の 1 つである"c"ため、変数名は不適切です。c()

于 2014-01-10T09:25:33.010 に答える
1

rollapply次のように使用して、隣接するマトリックス列の完全なケースの数を計算できます。

m <- matrix(c(NA,1,1,NA,1,1,1,1),ncol=4)
#     [,1] [,2] [,3] [,4]
#[1,]   NA    1    1    1
#[2,]    1   NA    1    1

library(zoo)

rowSums(rollapply(is.na(t(m)), 2, function(x) !any(x)))
#[1] 0 1 2
于 2014-01-10T09:43:15.677 に答える
0

これは、 と の両方matrixで機能するはずです。data.frame

> sum(apply(c, 1, function(x)all(!is.na(x))))

[1] 2

そして、大きな行列を単純に反復することができますM

for (i in 1:(ncol(M)-1) ){
    c <- M[,c(i,i+1]
    agreement <- sum(apply(c, 1, function(x)all(!is.na(x))))
}
于 2014-01-10T09:28:19.843 に答える