我想及时绘制甘蔗、椰子和水稻的产量。但我得到了这个输出:
发布于 2022-04-29 17:36:36
您的数据首先应该是aggregated,因为它具有相同日期和裁剪的多个条目,所以我假设您希望sum作为聚合函数。在data参数中,我使用transform应用ISOdate (参见other answer),并为pch和color定义一个数字向量以进行绘图。
a <- aggregate(Yield ~ .,
transform(data, Year=ISOdate(Year, 1, 1, 0),
clr=as.integer(factor(Crop)) + 1), FUN=sum) 然后,plot和空基图,我建议log算法可能适合你的数据。接下来,在作物上使用lines循环by。最后,放置一个漂亮的legend,仅此而已。
这里有一个基本版本:
plot(Yield ~ Year, a, log='y', type='n')
by(a, a$Crop, \(x) with(x, lines(Yield ~ Year, type='b', col=clr, pch=clr)))
legend('right', legend=unique(data$Crop), col=unique(a$clr), pch=unique(a$clr))

发布于 2022-04-29 19:48:37
(数据,aes(x =年,y=产量,颜色=作物,组=作物)+ geom_smooth() + xlab("")
https://stackoverflow.com/questions/72059310
复制相似问题