平均値を逐次更新で求めてu[i]がどう変動していっているのか見たいです。
#54個のデータ(0〜100の点数）
x<-c(95,95,95,90,75,93,79,75,80,70,91,97,
　　　75,88,69,83,68,90,96,94,85,80,80,79,85,85,88,
　　　90,83,78,90,90,97,70,87,85,73,94,65,78,76,84,
　　　80,76,77,78,80,51,77,81,54,83,70,76)
n<-length(x)
#uは長さ54ですべて0のベクトル(すべて1にしてもできませんでした。）
u<-numeric(n)

u[1]<-x[1]
for(i in 1:n-1) {
#逐次更新式
u[i+1]<-(n/(n+1))*(u[i]+(1/n)*x[i+1])
}
上のコードで実行しているのですが、
Error in u[i + 1] <- (n/(n + 1)) * (u[i] + (1/n) * x[i + 1]) :
replacement has length zero
のようなエラーが出てしまいます、どのようにすれば修正できるでしょうか？
他の方法（u_h[i]=mean(head(x,n=i))）では求められましたが、上記のコードではどこがダメなのかわかりませんでした。