0

我正在尝试创建一个密度图并将其添加为法线曲线。

histvb <- hist(color$VB,
               probability=TRUE,
               breaks=seq(4000, 16000, by=1000),
               main=NULL,xlab=c("Visible Brigthness"), 
               ylab=("Probability Densities"),cex.lab=1.2) 
vb<-color$VB
xfit<-seq(min(0),max(16000),length=16000) 
yfit<-dnorm(xfit,mean=mean(vb),sd=sd(vb))  
yfit<-yfit*diff(hist$mids[1:2])*length(vb)
lines(xfit, yfit, col="blue",lty=5, lwd=2)

当我运行这些线时，我得到了这条线的情节，一切似乎都很好，但我也收到了这条消息

yfit <- yfit*diff(hist$mids[1:2])*length(vb)
Error in hist$mids : object of type 'closure' is not subsettable

我担心这个错误消息，因为我不知道它是什么意思..

4

1 回答 1

1

你的意思是做histvb$mids[1:2]，不是hist$mids[1:2]：histvb是你在第一行创建的结果：hist是一个绘制直方图的函数。改变你的线路

yfit<-yfit*diff(hist$mids[1:2])*length(vb)

至

yfit<-yfit*diff(histvb$mids[1:2])*length(vb)
于 2013-10-09T15:05:09.747 回答