r - 如何在 R 语言服从高斯分布的数据上绘制高斯曲线？

Question

我有一些看起来像服从高斯分布的数据。所以我用

my.glm<- glm(b1~a1,family=Gaussian)

然后使用命令

summary(my.glm).

结果是：

Call:
glm(formula = b1 ~ a1, family = gaussian)

Deviance Residuals: 
      Min         1Q     Median         3Q        Max  
-0.067556  -0.029598   0.002121   0.030980   0.044499  

Coefficients:
             Estimate Std. Error t value Pr(>|t|)    
(Intercept)  0.433697   0.018629   23.28 1.36e-12 ***
a1          -0.027146   0.001927  -14.09 1.16e-09 ***
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 

(Dispersion parameter for gaussian family taken to be 0.001262014)

Null deviance: 0.268224  on 15  degrees of freedom
Residual deviance: 0.017668  on 14  degrees of freedom
AIC: -57.531

Number of Fisher Scoring iterations: 2

我认为他们很合适。但是我怎样才能在这些数据上画出高斯曲线呢？

Answer 1

假设截距具有正态分布，您可以像这样绘制它的分布：

x <- seq(0.3,0.6,by =0.001)
plot(x, dnorm(x, 0.433697, 0.018629), type = 'l')

你可能想添加你的数据：

rug(b1)

由于您没有提供数据，我们可以弥补一些（在示例中进行一些转换以匹配统计信息）：

set.seed(0)
b <- rnorm(15)
b1 <- ((b - mean(b))/sd(b)  * 0.018629) + 0.433697
rug(b1)

您还可以覆盖数据的内核密度估计

lines(density(b1), col = 'red')

给出以下情节：

Answer 2

简单的： ?dnorm

使用 dnorm 创建所需均值和 sd 的高斯曲线，而无需将自己绑定到任何数值拟合函数。这是显示您的数据如何“拟合”理论曲线的简单且良好的方法。与绘制拟合数据并试图找出它与高斯的“接近程度”不同。