绘制calc.poly的多项式输出

Question

我首先道歉，因为我认为这是一个非常简单的问题，但是我无法从帮助文件中收集包裹＆＃39; polynom＆＃39;如何解决这个问题呢。对于几年中的一年，我有两个向量x（一年中的d）和y（e表示产蛋指数）数据：

d=c(169,176,183,190,197,204,211,218,225,232,239,246)
e=c(0,0,0.006839425,0.027323127,0.024666883,0.005603878,0.016599262,0.002810977,0.00560387  8,0,0.002810977,0.002810977)

我希望每年使用poly.calc函数创建一个多项式函数，我可以用它来插入最大产蛋量的时间。我希望然后将该函数叠加在数据图上。首先，我对poly.calc函数没有任何问题：

egg1996<-poly.calc(d,e)
egg1996
3216904000 - 173356400*x + 4239900*x^2 - 62124.17*x^3 + 605.9178*x^4 - 4.13053*x^5 +  
0.02008226*x^6 - 6.963636e-05*x^7 + 1.687736e-07*x^8 

我可以简单地

plot(d,e)

但是当我尝试使用线条函数在图上叠加函数时，我感到困惑。帮助文件声明poly.calc的输出是类多项式的对象，因此我假设＆＃34; egg1996＆＃34;将是＆＃34; x＆＃34;在：

lines(x, len = 100, xlim = NULL, ylim = NULL, ...)

但我似乎无法根据列出的例子：

lines (poly.calc( 2:4), lty = 2)

或基于论点：

1. x类的对象＆＃34;多项式＆＃34;。
2. 要进行评估的向量的len大小。
3. xlim，ylim具有合理默认值的x和y值的范围

提出一个成功绘制多项式的命令＆＃34; egg1996＆＃34;到原始数据。 我知道这个问题在你们的下方，但我会非常感谢你们的帮助。非常感谢。

Answer 1

我不使用polynom包，但结果数据集与第一个plot()调用完全不同（X和Y轴）。如果您不介意在两个单独的面板中使用它，则会提供两个图表进行比较：

library(polynom)

d <- c(169,176,183,190,197,204,211,218,225,232,239,246)
e <- c(0,0,0.006839425,0.027323127,0.024666883,0.005603878,
       0.016599262,0.002810977,0.005603878,0,0.002810977,0.002810977)

egg1996 <- poly.calc(d,e)

par(mfrow=c(1,2))
plot(d, e)
plot(egg1996)