Rで試してハンドルエラー

Question

データで幾何平均を返す関数を作成しようとしています。 try（）でループを作り、有効なデータを渡したいのですが、実際に試したものはうまくいきませんでした。

これは

R=function(g) 
{ 
k=1 
n=length(g) 
for(i in 1 : n) 
{ 
ifelse(g[i]>0, k<-k*g[i], stop("Negative component")) 
k 
} 
t=k^(1/n) 
t 
} 

機能であると私は、私は有効なデータのための手段を取得したいが、それは0

にすべての値になり、このループで

set.seed(123) 
data <- matrix(rnorm(10000, mean=3), ncol=25, dimnames=list(NULL, paste("X", 
1:25, sep="."))) 


v=rep(0,400) 
for(i in 1 : 400) 
{ 
try("v[i]=R(data[,i])",TRUE) 
} 
v 

を、この機能を使用したいです

[1] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 

[37] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 

[73] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 

[109] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 

[145] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 

[181] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 

[217] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 

[253] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 

[289] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 

[325] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 

[361] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 

[397] 0 0 0 0 

あなたはどこに間違っていたか教えていただけますか？ ありがとう

Answer 1

数値ベクトルとしてのtの宣言が欠落していて、try()という構文が問題だと思います。次のコードは動作するはずです。

R=function(g) 
{ 
t = numeric(length(g)) 
k=1 
n=length(g) 
for(i in 1 : n) 
{ 
    ifelse(g[i]>0, k<-k*g[i], stop("Negative component")) 
} 
t=k^(1/n) 
t 
} 
set.seed(123) 
data <- matrix(rnorm(10000, mean=3), ncol=25, dimnames=list(NULL,paste("X", 
                    1:25, sep="."))) 

v=numeric(400) 
for(i in 1 : 400) 
{ 
try(v[i]<-R(data[i,]),TRUE) 
} 
v