R语言统计
文章平均质量分 54
clancy_wu
北京中医药大学博士生,方向:针刺对脑卒中、偏头痛的机制研究。
展开
-
Mediation Analysis and causal mediation analysis
中介分析原创 2023-02-02 12:04:52 · 830 阅读 · 3 评论 -
统计学进阶之临床建立模型(一)
临床建模原创 2022-10-10 22:15:48 · 918 阅读 · 0 评论 -
Time New Roman字体
windows下,字体放在:C:\Windows\Fonts其中time new roman包含4个字体,如下:linux下,字体一般可以放在/usr/share/fonts下,也可以使用fc-list列出所有加载的字体.原创 2021-12-23 12:46:41 · 1226 阅读 · 0 评论 -
R语言相关性分析与偏相关分析
原始数据:结果用R实现:统计学里,比较重要的方法,一个是协方差分析,另一个就是偏相关分析了。协方差分析,就是在回归分析的基础上进行建模,找出去除混杂变量后,两个变量是否有显著关系,在R语言中公式如下:summary(aov(y ~ x1+x2)) # x1是混淆变量,x2是组别。相关性分析的原理也很简单,就是用y~x建立一元回归模型f(x),然后用实际的x减去模型预测的x,再根据差值分布求得p值。偏相关分析的原理:如上图,如果我们想考察X1和X2的偏相关系数,控制变量为z1, z2原创 2021-12-19 14:29:45 · 8698 阅读 · 1 评论 -
常用统计方法对比
参数方法非参数方法Welch TestWilcox test (one way or two way)Student TestMann-Whitney’ test (two way )ANOVA TestKruskal-Wallis ( U test )ANCOVA TestFriedman TestLSD Test 等Nemenyi Test原创 2021-03-11 18:18:25 · 508 阅读 · 1 评论 -
病例对照研究中—两组组间比较—的统计方法选择,基于R语言
医学中最常设计的试验就是病例对照研究,以探究某一干预措施是否有改善性。需要根据基线的情况,选择相应的方法。试验数据如下:(声明:该数据是随机自动生成的,虚拟的,该计算结果不代表任何真实的事情,该数据不适用于现实世界。)数据由试验组长病程10名,试验组短病程10名,对照组10名组成,干预手段为减肥。一、两组病例对照研究常用的三种方法(1)T检验法,将结果进行T检验。因为基线齐,说明两组干预前组间差异不明显,直接用T检验对比干预后的结果。long_after <- filter(dat原创 2021-03-11 17:11:49 · 20376 阅读 · 0 评论 -
多元回归
多元线性回归方程:y = b0 + b1X1 + b2X2 +b3X3…建立零假设:b1.b2.b3…为零,b1=b2=b3=0spss中的多元回归分析流程和一元的差不多,跳过。最后算出来p<0.05,就是b1.b2.b3…不为零;若>0.05,就是有可能为零。如果接受零假设的话,下一步就要进行逐步回归,筛选出有效的变量。在spss上操作较简单。下面是R的操作流程,使用数据为ISwR的cy...原创 2020-02-21 20:35:22 · 1938 阅读 · 0 评论 -
R中的假设检验(四)卡方检验
最近几天没碰编程,浑身难受,日子过的浑浑噩噩,还是得每天接触下编程。相较于核磁分析,我感觉统计建模更对我的胃口。今天写卡方检验(χ2检验,chi-square test)。医学中最最常见的就是四格表的卡方检验,spss教材上计算四格表资料的卡方检验结果值有:Pearson chi-square 、Continuity correction 、Likelihood ratio 、 Fisher’s...原创 2020-02-19 20:34:29 · 13286 阅读 · 0 评论 -
R中的假设检验(二)方差分析画图
再有协方差,逻辑回归,卡方,统计这一块就结束了。感觉我在统计这一块花了好多的时间啊,最近疯狂钻。等我把这弄完,R这一块我就不更新了,数据清洗也不要了。博士前都不来更新了!还有好多事情等着我去做,还有论文没写,还有fMRI没研究,还有matlab没有学,我真是个不务正业到极点了,555,希望以后会用到这些知识和手段,不然我巨亏。要赶紧抓时间学习专业知识和领域技术,能发文章的技术才是好技术。今天就写...原创 2020-02-04 23:32:42 · 8555 阅读 · 1 评论 -
R中的假设检验(二)方差分析
方差分析的原理:原本想跳过这个部分的,然而突然发现自己看懂了这些公式,就忍不住写出来。(1)Xij=X均值+(X组内均值-X均值)+(Xij-X组内均值) #Xij代表第i组第j个观察值,X均值代表所有观察值的均值,X组内均值就是每组内部的均值。可以看到打开括号就是成立的,只是巧妙的引入了组内均值和组间均值。(2)SSD总=SSD组内+SSD组间 #把X组内均值-X均值的方差计算出来,...原创 2020-02-04 00:54:31 · 7577 阅读 · 6 评论 -
R中的假设检验(三)相关与回归
学到哪儿写到哪儿,这篇后面会随时更新。一、相关分析研究变量与变量间的关系就需要分析其相关性,就需要使用相关分析,相关分析比较简单。它分为皮尔逊相关、斯皮尔曼相关、肯德尔等级相关,并且他们有各自的适用范围:(1)Pearson,皮尔逊积差相关系数,适用于两个变量为连续变量,且两个变量需要严格符合正态性分布,由于其严格的条件所以其具有较强的统计效能。(2)Spearman,斯皮尔曼相关系数,对两个...原创 2020-02-02 20:02:46 · 2666 阅读 · 0 评论 -
R-中文期刊散点图
最近开始准备间断更新R作图代码,以方便自己回顾命令,不然容易忘记。中文期刊大部分要求的是灰度图,无彩色,所以今天的图是黑灰相间散点图。如下:个人觉得画的还不错,主要是不能使用彩色,所以只能运用线条样式和颜色会读来区分不同的组别。代码如下:修改相应的代码区可以实现不同的样式,最后调整成你觉得最顺眼的。今天的分享END...原创 2019-12-25 00:22:18 · 284 阅读 · 0 评论