我有一个条形图,下半部分应该适合这个公式:
. 现在我想绘制整个条形图并在条形图的最后一部分显示拟合模型,因为它仅适用于该部分。但是,我找不到仅在后半部分显示折线图的方法。如果我只是做类似的事情y~axexp(-b*x^2)
submitted=c(1L, 1L, 1L, 1L, 1L, 1L, 2L, 1L, 2L, 2L, 3L, 2L, 1L, 1L, 4L,
3L, 2L, 11L, 6L, 2L, 16L, 7L, 17L, 36L, 27L, 39L, 41L, 33L, 42L,
66L, 92L, 138L, 189L, 249L, 665L, 224L, 309L, 247L, 641L, 777L,
671L, 532L, 749L, 506L, 315L, 292L, 281L, 130L, 137L, 91L, 40L,
27L, 34L, 19L, 1L)
x=seq(0:(length(submitted)-1))
y1=rs$submitted[30:(length(submitted)-1)]
x1=seq(0:(length(y1)-1))
fit1=nls(y1~a*x1*exp(-b*x1^2),start=list(a=500,b=.01),trace=TRUE)
barplot(submitted,names.arg=x, las=2, cex.axis=0.8, cex=0.8)
lines(predict(fit1))
该行已显示,但位置错误。那么如何控制画线的位置呢?