I want to replace my NA values from a matrix acquired by :
read.table(…)
Those values should be the mean of the corresponding row.
x[is.na(x)] <- mean(x, na.rm=TRUE) # for vectors or for a matrix as a whole
t( apply(x, 1, function(xv) { xv[is.na(xv)] <-
mean(xv, na.rm=TRUE)
return(xv)}
) ) # for a row-oriented sol'n
a = c(NA, 1, 2, 3, 10)
a[which(is.na(a)==TRUE)] = mean(a,na.rm = T)
Here's some sample data.
m <- matrix(1:16, nrow=4)
m[c(1,4,6,11,16)] <- NA
And here's how I'd fill in missings with the row means.
k <- which(is.na(m), arr.ind=TRUE)
m[k] <- rowMeans(m, na.rm=TRUE)[k[,1]]
Your data will be in a data.frame
; you'll have to convert to a matrix first using as.matrix
. You may or may not want to leave it in that format; to convert back use as.data.frame
.