密度プロットを作成して通常の曲線を追加しようとしています。
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
このエラーメッセージの意味がわからないので心配です..