我要算出总和
(十一-米)^3
下面的代码正确吗?
s<-sum((xi-m)^3)
发布于 2015-04-29 08:52:37
一些数据首先:
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
发布于 2015-04-29 18:45:47
我认为您的代码是正确的,如果您要寻找的是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
https://stackoverflow.com/questions/29938478
复制相似问题