那个网站做淘宝推广比较好,wordpress商业模板,seo是什么意思如何实现,网站建设基本常识基本原理方差分析(Analysis of variance, ANOVA)是用于两个或两个以上样本均数比较的方法#xff0c;还可以分析两个或多个研究因素的交互交互作用以及回归方程的线性假设检验等。其基本思想是将全部观察值间的变异——总变异按设计和需要分解成两个或多个组成部分#xff0c…基本原理方差分析(Analysis of variance, ANOVA)是用于两个或两个以上样本均数比较的方法还可以分析两个或多个研究因素的交互交互作用以及回归方程的线性假设检验等。其基本思想是将全部观察值间的变异——总变异按设计和需要分解成两个或多个组成部分再作分析。本文将针对单因素方差分析(One-Way ANOVA)开展R语言的实践。单因素方差分析是将同质的受试对象随机地分配到各处理组再观察对比实验效应的方法。在单因素方差分析中首先将总变异(SS_总)分解为组间变异(SS_组间)也叫处理变异和组内变异(SS_组内)也叫误差变异然后比较两者的平均变异(MS_组间MS_组内)比较时采用两者的比值F值(MS_组间/MS_组内)R语言实践为研究钙离子对体重的影响作用某研究者将36只肥胖模型大白鼠随机等分为三组每组12只分别给予高脂正常剂量钙(0.5%)、高脂中剂量钙(1.0%)和高脂高剂量钙(1.5%)三种不同的饲料喂养9周测其喂养前后体重的差值。结果如下高脂正常剂量钙(0.5%)332.96297.76312.57295.47284.25307.97292.12244.61261.46286.46322.49282.42高脂中剂量钙(1.0%)253.21235.87269.3258.9254.39200.87227.79237.05216.85238.03238.19243.49高脂高剂量钙(1.5%)232.55217.71261.15220.72219.46247.27280.75196.01208.24198.41240.35219.56首先导入数据my_data read.xlsx(./anova_1.xlsx, sheetIndex 1, encoding UTF-8)其次对通过箱型图和线图对数据进行可视化# 绘制箱型图
ggboxplot(my_data, x 剂量, y 体重差值,color 剂量, palette c(#00AFBB, #E7B800, #FC4E07),order c(高脂正常剂量钙, 高脂中剂量钙, 高脂高剂量钙),ylab 体重差值, xlab 剂量)
# 绘制线图
ggline(my_data, x 剂量, y 体重差值,add c(mean_se, jitter),order c(高脂正常剂量钙, 高脂中剂量钙, 高脂高剂量钙),ylab 体重差值, xlab 剂量)绘制结果如下然后进行方差分析# 计算方差分析
res.aov - aov(体重差值 ~ 剂量, data my_data)
# 输出分析结果
summary(res.aov)分析结果如下Df Sum Sq Mean Sq F value Pr(F)
剂量 2 28925 14462 27.21 1.04e-07 ***
Residuals 33 17540 532
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1上述结果表明经假设检验得出多个总体均数不全等。最后我们可以计算Tukey HSDTukey Honest Significant DifferencesR函数TukeyHSD()可以实现以在组均值之间进行多次成对比较# 确定特定组对之间的平均差异是否具有统计显著性
TukeyHSD(res.aov)输出结果如下Tukey multiple comparisons of means95% family-wise confidence levelFit: aov(formula 体重差值 ~ 剂量, data my_data)$剂量diff lwr upr p adj
高脂正常剂量钙-高脂高剂量钙 64.86333 41.76836 87.95831 0.0000002
高脂中剂量钙-高脂高剂量钙 10.98000 -12.11498 34.07498 0.4810804
高脂中剂量钙-高脂正常剂量钙 -53.88333 -76.97831 -30.78836 0.0000064上述结果表明说明中、高剂量钙与正常钙喂养9周前后体重差值不同但不能认为中、高剂量钙喂养9周前后体重差值不同。
