我正在尝试在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)
偏斜度(gal)错误:找不到对象'n'