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