我试图在 MASS 数据库中找到变量星系的偏度
我已经概述了偏度的公式,然后加载了 MASS 包,定义了星系数据集,但是由于偏度公式中未定义 n,因此出现错误。我不确定是否需要包含一个 n,因为我只希望我的样本量是星系中的所有结果。
我的代码是
myskewness = function(x) {
n = length(x)
xbar = sum(x) / n
myskewness = (sqrt(n) * sum((x - xbar)^3)) / sum((x -xbar)^2)^(3 / 2)
return(myskewness)
}
gal <- galaxies
myskewness(gal)
myskewness(gal) 中的错误:找不到对象“n”