psych::describe 函数中的标准错误 - 它指的是什么?
Standard error in psych::describe function - what is it referring to?
假设我们有一个名为 "df" 的数据框,我们想知道 df 中一个名为 "x" 的变量的偏度和峰度值。
假设我们使用:
psych::describe(df$x)
并得到如下结果:
vars n mean sd median trimmed mad min max range skew kurtosis se
1 1 478 98.54 19 102.5 100.57 18.53 34 125 91 -0.94 0.47 0.87
最后一个值 se 指的是什么?偏斜或峰度的标准误差?
se
指的是"standard error of the mean"
注意:您可以在 R 终端中输入 describe
来阅读源代码。
无论如何,作为双重检查,这里是虹膜数据集
上describe
的输出
unlist(sapply(iris[1:4], describe)[13,])
#output
Sepal.Length Sepal.Width Petal.Length Petal.Width
0.06761132 0.03558833 0.14413600 0.06223645
这里是我手写的均值标准误差函数的输出
sapply(iris[1:4], sem)
#output
Sepal.Length Sepal.Width Petal.Length Petal.Width
0.06761132 0.03558833 0.14413600 0.06223645
p.s。我的 sem
函数
function(x) {
sqrt(var(x)/length(x))
}
假设我们有一个名为 "df" 的数据框,我们想知道 df 中一个名为 "x" 的变量的偏度和峰度值。 假设我们使用:
psych::describe(df$x)
并得到如下结果:
vars n mean sd median trimmed mad min max range skew kurtosis se
1 1 478 98.54 19 102.5 100.57 18.53 34 125 91 -0.94 0.47 0.87
最后一个值 se 指的是什么?偏斜或峰度的标准误差?
se
指的是"standard error of the mean"
注意:您可以在 R 终端中输入 describe
来阅读源代码。
无论如何,作为双重检查,这里是虹膜数据集
上describe
的输出
unlist(sapply(iris[1:4], describe)[13,])
#output
Sepal.Length Sepal.Width Petal.Length Petal.Width
0.06761132 0.03558833 0.14413600 0.06223645
这里是我手写的均值标准误差函数的输出
sapply(iris[1:4], sem)
#output
Sepal.Length Sepal.Width Petal.Length Petal.Width
0.06761132 0.03558833 0.14413600 0.06223645
p.s。我的 sem
函数
function(x) {
sqrt(var(x)/length(x))
}