手机欧美视频网站模板下载 迅雷下载地址,visualstudio网页设计作业,页面效果图,wordpress设置ssl网站打不开看过许多统计教程#xff0c;这篇是我最推荐的介 绍数值变量如果服从正态分布#xff0c;采用均数标准差进行统计描述#xff0c;采用方差分析进行组间比较#xff0c;如果组间差异有统计学意义#xff0c;进一步采用LSD法(也可以是其它方法)进行两两比较。如果不服从正… 看过许多统计教程这篇是我最推荐的介 绍数值变量如果服从正态分布采用均数±标准差进行统计描述采用方差分析进行组间比较如果组间差异有统计学意义进一步采用LSD法(也可以是其它方法)进行两两比较。如果不服从正态分布采用中位数(四分位数间距)进行统计描述组间比较采用非参数检验(Kruskal-Wallis秩和检验)当组间总的有统计学差异进一步采用Dunn法(也可以是其它方法)进行多重比较。对于分类变量采用例数(百分比)的形式描述组间比较采用χ2检验或精确概率法当组间总的有统计学差异进一步采用χ2检验或精确概率法进行两两比较两两比较的P值需采用Boferroni法(也可以选择其它方法)进行校正。对于只会SPSS软件的同学来说做两两比较的卡方检验或精确概率法时需要一次次选择个案一次次点菜单太浪费时间和精力。用R软件的优势在于做两两比较时不用再一次一次点菜单了直接用几行命令就可以了。我们想看看不同工作岗位的男女比例是否有差异。软件操作1、导入数据mydata clipboard)#查看数据str(mydata)head(mydata)2、变量因子化mydata$sex$sex,levelsc(0,1),labelsc(mydata$work$work,levelsc(1,2,3,4),labelsc(#查看数据str(mydata)head(mydata)3、对于分析涉及到的变量形成表格数据我们想研究不同工作岗位的男女比例是否有差异需要把涉及到的两个变量work和sex生成表格输入tab$work,mydatatab4、卡方检验和精确概率法#卡方检验chisq.test(tab)#精确概率法fisher.test(tab)结果解读使用chisq.test命令进行卡方检验结果显示统计量χ220.911p0.0001099在论文中需写成χ220.91p0.001说明不同岗位的男女比例有差异但具体是哪两个岗位有差异还需进一步两两比较。(注论文中卡方值一般保留两位小数p值一般保留三位小数如果p0.001可以直接写成p0.001)。如果不符合卡方检验条件需改用精确概率法其实就算符合卡方检验条件也可以直接用精确概率法。使用fisher.test命令进行精确概率法结果显示p8.359*10-5论文中需要写成p0.001说明不同岗位的男女比例有差异但具体是哪两个岗位有差异还需进一步两两比较。5、两两比较#install.packages(rcompanion)library(rcompanion)#两两比较卡方检验pairwiseNominalIndependence(tab,chisq TRUE,fisher FALSE,gtest FALSE,method bonferroni)上面红色的warning messages提示卡方检验得到的结果可能不准确为了保险起见我们采用Fisher精确概率法进行两两比较。#两两比较精确概率法pairwiseNominalIndependence(tab,chisq FALSE,fisher TRUE,gtest FALSE,method bonferroni)结果解读事后两两比较由于检验次数的增多会增加一类错误的概率所以通常需要校正p值或校正检验水准两者本质是一样的。我们通过R软件进行两两比较的好处是可以直接设定两两比较p值校正方法如本例中我们直接设置为bonferroni法p.adj.Fisher列即为校正后的p值可以直接与0.05相比如果0.05说明两者之间差异有统计学意义否则说明无统计学意义。本例中A岗和B岗的男女比例差异有统计学意义p6.42*10-5需要写成p0.001。需要说明的是有的p值经校正之后1这是很正常的。知识拓展本例中p值校正的方法我们选择的是bonferroni法只要在论文中说明所用方法即可。大家也可以尝试别的方法更多信息可以在Rstudio中输入?p.adjust获取。希望今天的文章对你有用随手点个在看吧。本文由刘老师医学统计(IDourdata)原创转载请注明来源。转载、投稿、商务合作、统计学习/咨询请在公众号后台添加作者微信。你点的每个赞我都认真当成了喜欢