无法获得与绘制的直方图分布不同的密度函数

Question

我想测试此代码定义的统计体验的正常性。

# a

y0<-50
a<-10
b<-0.4
theta<-4
T<-30
h<-0.01
t<-seq(0,T,h)
n<-length(t)
M<-1000
est<-rep(0,M)

 for( j in 1:M)
 {
   W<-c(0,cumsum(rnorm(n-1,0,1)*sqrt(h)))
   y<-rep(y0,n)


for(i in 1:(n-1))
{
  y[i+1]<-y[i]+a*t[i]*h-b*y[i]*h+theta*(W[i+1]-W[i])
}

delta<-1/n
X<-y
M<-(a/b)+(X-(a/b))*exp(-b*delta)
est[j]<-sqrt((2*b)/(1-exp(-2*b*delta))*(1/n)*sum(X[2:n]-M[1:(n-1)]))
 }

hist(sqrt(n)*(est-theta),nclass = 40,freq = FALSE)

请注意，我根本不需要获得直方图本身的折线图。我只需要添加其他正态分布来与直方图进行比较。这个新的正态分布由下式定义：

 x<-seq(-1,1,length=1000)
    y<-dnorm(x,mean = 0,sd=sqrt(theta^2/2))
    lines(x,y)

但是我没有在上面定义的密度函数图中添加任何曲线图。

 只是我无法获得密度分布

我确实明白这是什么问题？谢谢你的帮助！