用R反向统计

Question

我想做的事听起来很简单.我想用R绘制一条正常的IQ曲线,平均值为100,标准差为15.然后,我希望能够在它上面叠加数据的散点图.

有人知道怎么做吗？

Answer 1

我猜你要做的是这样的:你想用平均值100和sd = 15 来绘制模型的正常密度,并且你想要覆盖那些据称遵循模型的一组观察的经验密度正常密度,以便您可以看到模型密度与经验密度的匹配程度.下面的代码应该这样做(在此,X是实际观测的,但为了说明目的我正在与一个混合正态分布N(100,15)+ 15*N(0,1)产生它的载体,即所述声称N(100,15)分布加噪声).

require(ggplot2)
x <- round( rnorm( 1000, 100, 15 )) + rnorm(1000)*15
dens.x <- density(x)
empir.df <- data.frame( type = 'empir', x = dens.x$x, density = dens.x$y )
norm.df <-  data.frame( type = 'normal', x = 50:150, density = dnorm(50:150,100,15))
df <- rbind(empir.df, norm.df)
m <- ggplot(data = df, aes(x,density))
m + geom_line( aes(linetype = type, colour = type))