人口の何パーセント見つけるための私のR関数は平均値のkの標準偏差内にある、

Question

whatpercent <- function(x) { 
i<-1 
for(i in length(x)) { 
if(x[i]<= mean(x)+sd(x)) { 
vec <- c(vec,x[i]) 
} 
else{ 
next 
} 
} 
percent <- length(vec)/length(x) 
percent <- percent * 100 
return(percent) 
} 

さて、私は、私が使用したときと同じものを取得するように見える動作しません

whatpercent(mydata,1) whatpercent(mydata,2) whatpercent(mydata,3) ... 

何ができますか問題？

Answer 1

ここでは、あなたが望むものを達成するための簡単な機能があります。このようなタスクでは、ループを回避するためにベクトル化された関数を使用する必要があります。また、コードをはるかにクリーンで読みやすくします。

whatpercent <- function(x, k) { 
    mean(abs(x - mean(x)) <= sd(x) * k) 
} 

#Example 
x <- rnorm(100) 
whatpercent(x,2) 
[1] 0.95