计算选定列上每一行的标准偏差

Question

我不习惯 R，所以为了练习，我正在尝试做我以前在 R 上的 SPSS 上做的所有事情。

在我的数据集中，每一行都是一个案例。这些列是调查问题（每个问题 1 个）。

假设我有“A1”到“A6”、“B1”到“B6”等列

我刚刚计算完 A1 到 A6 上每个人的平均值

data$meandata <- rowMeans(subset(data, select=c(A1:A6), na.rm=TRUE))

如何计算meandata 的标准差？

Answer 1

嘿，最简单的方法是使用 apply() 函数。

假设您有 25 行数据和 6 列标记为 A1 到 A6。

data <- data.frame(A1=rnorm(25,50,4),A2=rnorm(25,50,4),A3=rnorm(25,50,4),
A4=rnorm(25,50,4),A5=rnorm(25,50,4),A6=rnorm(25,50,4))

您可以使用 apply 函数通过下面的代码找到每行第 1 到第 6 列的标准差。第一个参数是您的数据对象。第二个参数是一个整数，指定 1 表示行或 2 表示列（这是函数将应用于数据框的方向）。最后一个参数是您希望应用于数据框的函数（例如本例中的均值或标准差 (sd)。请参见下面的代码。

apply(data[,1:6],1,sd)

索引可用于限制传递给应用函数的数据的行数或列数。这是通过在数据对象后面的括号内为您感兴趣的行或列输入一个数字向量来完成的。

data[c(row.vector),c(column.vector)]

假设你只想知道前 3 列的 sd。

apply(data[,1:3],1,sd)

现在让我们看看第 4 到 6 列和第 1 到 10 行的 sd

apply(data[1:10,4:6],1,sd)

为了更好的衡量，让我们找到每列的 sd

apply(data,2,sd)

请注意，sd 接近 4，这是我在为 A1 到 A6 列生成伪随机数据时指定的值。

希望对你有帮助