R软件中的总和

Question

我需要计算出总和 （XI-M）^ 3

• 我取值1到25
• xi是向量x
的元素
• m是向量x
的平均值

以下代码是否正确？

s<-sum((xi-m)^3)

Answer 1

首先是一些数据：

set.seed(1)
x = rnorm(100, 20, 3)

您可以使用x - mean(x)，而不是使用scale，避免使用scale=F作为参数除以标准偏差：

sum(scale(x[1:25], scale=F)^3)

#> sum((x[1:25]-mean(x[1:25]))^3)==sum(scale(x[1:25], scale=F)^3)
#[1] TRUE

Answer 2

我认为你的代码是正确的，如果你正在寻找的东西是残差的总和（值减去平均值）由3提供，只需考虑这些例子：

y <- rnorm(25)
> sum(y - mean(y))^3
[1] -1.710569e-49
> sum(y - mean(y)^3)
[1] -4.262774
> sum((y - mean(y))^3)
[1] 9.449698