出口Cronbach Alpha结果 - R

时间:2017-04-15 21:20:50

标签: r psych

我按照THIS R Blogger教程来计算Cronbach alpha,它完美无缺。我想学习如何将结果导出到data.frame或文本文件中。我有什么想法可以从以下代码中导出结果:psych::alpha(d)?请注意,我查看了观星包,但是无法使用Cronbach输出仅输出回归和描述性统计数据。谢谢。

2 个答案:

答案 0 :(得分:1)

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

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

答案 1 :(得分:0)

正如所有R一样,有很多方法可以做你想做的事。第一件事是查看函数的帮助菜单(在这种情况下?alpha)。在那里你会看到alpha函数返回了许多对象。 (这是帮助文件的值部分中列出的内容。)

打印alpha输出时,只显示这些对象的子集。但是,要查看返回的整个对象列表,请使用&#34; str&#34;命令

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是描述比例结构的特别有用的统计量。您可能希望阅读有关如何使用

中的心理包找到系数omega的教程

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