>summary(result2)
Family: gaussian
Link function: identity
Formula:
recruit ~ s(spring, k = 4)
Parametric coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 0.36939 0.0200718.41<2e-16 ***
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Approximate significance of smooth terms:
edf Ref.df F p-value
s(spring) 1.315 1.554 1.567 0.156
R-sq. (adj) = 0.0713Deviance explained = 10.9%
GCV = 0.014294 Scale est. = 0.013292 n = 33
====================
>结果2<-gam(新成员(弹簧,k=4),族系=高斯,数据=a)
>总结(结果2)
族系:高斯
链接功能:标识
公式:
新成员(弹簧,k=4)recruit~s(spring, k=4)
参数系数:
估计标准误差t值Pr(>| t |)
(截距)0.36939 0.02007 18.41<2e-16***
显著性代码:0'**'0.001'**'0.01'*'0.05'.'0.1''1
平滑项的近似意义:
edf参考df p值
S(弹性)1.315 1.554 1.567 0.156
R平方(ADJ)=0.0713 解释偏差=10.9%
GCV=0.014294 标度估计=0.013292N=33
GCV
abbr.
总热值(Gross Caloric Valve)
最终出图如下
这里自动做统计检验的函数是 stat_compare_means()
读入数据
作图
这个函数来自于ggpubr这个包,只需要指定根据那一列来分组就可以了
默认的是Wilcoxon Rank Sum and Signed Rank Tests,如果要用t检验指定method参数
如果想把P值改成星号,直接加label=“p.signif”参数
这里如果不显著会在图上显示ns,如果不想要ns,可以加 hide.ns = TRUE 参数
星号的位置可以手动指定,用 label.y = c(26,31) 参数
使用到的是 ggsignif 这个包
小明的数据分析笔记本