导出 Cronbach Alpha 结果 - R

Question

我正在按照 THIS R Blogger 教程计算 Cronbach alpha，效果非常好。我想了解如何将结果导出到 data.frame 或文本文件中。关于如何从以下代码导出结果的任何想法：psych::alpha(d)？请注意，我查看了 stargazer 包，但无法让它与 Cronbach 一起工作，仅输出回归和描述性统计数据。谢谢。

Answer 1

这样输出可以保存为txt文件。您还可以使用 $ 运算符对使用 alpha 函数创建的对象进行子集化，以仅获取您感兴趣的信息。

setwd("~/Desktop")
out <- psych::alpha(d)
capture.output(out,file = "alpha.txt")

Answer 2

就像 R 的所有事物一样，有很多方法可以做您想做的事。第一件事是查看功能的帮助菜单（在本例中为 ?alpha）。在那里您会看到许多对象从 alpha 函数返回。（这是帮助文件的值部分中列出的内容。）

当您打印 alpha 的输出时，您只会看到这些对象的一个子集。但是，要查看返回的整个对象列表，请使用 "str" 命令

my.results <- alpha(my.data)
str(my.results)     #or just list the names of the objects
 names(my.alpha)
 [1] "total"         "alpha.drop"    "item.stats"    "response.freq" "keys"          "scores"        "nvar"          "boot.ci"      
 [9] "boot"          "Unidim"        "Fit"           "call"          "title"

然后您可以选择捕获这些对象中的任何一个供您自己使用。因此

my.alpha <- alpha(ability) #use the ability data set in the psych package
my.alpha #will give the normal (and nicely formatted output)
totals <- my.alpha$total   #just get one object from my.alpha
totals                     #show that object

将生成一行（没有花哨的输出）：

raw_alpha std.alpha   G6(smc) average_r      S/N         ase      mean        sd
 0.8292414 0.8307712 0.8355999 0.2347851 4.909159 0.006384736 0.5125148 0.2497765

您可以对返回的任何对象执行此操作。我们大多数编写包的人打印我们认为是函数输出的基本元素，但包括其他有用的信息。我们还允许其他功能（例如摘要）打印出其他信息。

因此，使用上面的示例，

summary(my.alpha)  #prints the rounded to 2 decimals my.alpha$total object
Reliability analysis   
 raw_alpha std.alpha G6(smc) average_r S/N    ase mean   sd
      0.83      0.83    0.84      0.23 4.9 0.0064 0.51 0.25

最后一句警告。我们中的许多人并不认为 alpha 是描述量表结构的特别有用的统计数据。您可能想阅读有关如何使用 psych 包查找系数 omega 的教程

http://personality-project.org/r/psych/HowTo/R_for_omega.pdf

导出 Cronbach Alpha 结果 - R

Export Cronbach Alpha Results - R

r

psych