W=as.matrix(r1)
P <- sum(W*W)
S <- apply(W,2,function(x){
  x <- sum(x)
  x <- x*x
  return(x)
})
h<- S/P#计算入样概率
我有如上一段代码，现在想根据h的值大小，从W中抽出入样概率最大的前10列，请问应该怎么写这样的R代码？