> A=matrix(floor(rnorm(20,0,4)),4,5)
> A[2,5]=NA
> A
[,1] [,2] [,3] [,4] [,5]
[1,] 6 1 -3 0 -3
[2,] -7 -1 1 -6 NA
[3,] -5 3 -10 4 -3
[4,] -5 1 -8 -2 0
> mean(A)
[1] NA
[1,] 6 1 -3 0 -3.000000
[2,] -7 -1 1 -6 -1.947368
[3,] -5 3 -10 4 -3.000000
[4,] -5 1 -8 -2 0.000000
> A[2,5]=NA
> A
[,1] [,2] [,3] [,4] [,5]
[1,] 6 1 -3 0 -3
[2,] -7 -1 1 -6 NA
[3,] -5 3 -10 4 -3
[4,] -5 1 -8 -2 0
> mean(A)
[1] NA
> A[is.na(A)]=mean(A,na.rm=T) # A [ is.na( A ) ] 选中所有NA,A [ ! is.na( A ) ] 选中所有非NA
> A
[,1] [,2] [,3] [,4] [,5][1,] 6 1 -3 0 -3.000000
[2,] -7 -1 1 -6 -1.947368
[3,] -5 3 -10 4 -3.000000
[4,] -5 1 -8 -2 0.000000