統計ソフトRの初心者です。sumの扱いに困ってます

質問

Rを使ったコードの書き方がわからず苦労しています。

(x[j]-mu[i])^2
をjが1からｎまで足し合わせたいのですが、
これをRでどう書いていいのかわかりません。

どなたか教えてください。よろしくお願いします。

回答No.2

backs

xもmuもベクトルなのでしょう？　xとmuの大きさ（データ数）が同じなら：

> x <- rnorm(10)
> mu <- rnorm(10)
> sum((x - mu)^2)
[1] 12.50165

とすればよいのでは？と思いますが・・・

ちなみに、やってみれば分かることですがNo.1が提示してくださった方法だと：

> (my.sum <- (x - mu) * t(x - mu))

ということになりますね。私は代入演算子は=ではなく、<-を使った方がよいと思いますが(^_^;)

回答No.1

sinisorsa

ｘとmuを行ベクトルとして与えます。
sum=(x-mu)*t(x-mu)
で計算できます。
t(a)はaの転置をとる関数です。
上の式は、ベクトル(x-mu)同士の内積を与えます。