中央値が置換され、数値データが必要

Question

グループに基づいて欠損値を代入しようとしています。 median（）関数が数値データを必要とするエラーが発生していますが、データがすべて数値であるため問題が見えません。ここには最小限再現可能な例があります。

set.seed(123) 
cluster = sample(seq(1,10),1000,replace=TRUE) 
V1 = sample(c(runif(100),NA),1000,replace=TRUE) 
V2 = sample(c(runif(100),NA),1000,replace=TRUE) 

df = as.data.frame(cbind(cluster,V1,V2)) 

df_fixed = by(df,df$cluster,function(x){replace(x,is.na(x),median(x, na.rm=TRUE))}) 

エラーが返さ：

median.default（X、na.rm = TRUE）でエラーが発生しました：必要な数値データ

は、このコードはので、しかし動作します問題は中央値関数である。

df_fixed = by(df,df$cluster,function(x){replace(x,is.na(x),1)}) 

Answer 1

df_fixed <- apply(df[,2:3], 2, function(x) { 
    md <- sapply(sort(unique(df$cluster)), function(k) median(x[df$cluster==k], na.rm=TRUE)) 
    x[is.na(x)] <- md[df$cluster][is.na(x)] 
    return(x) 
}) 

any(is.na(df_fixed)) 
# [1] FALSE