我第一次尝试使用循环生成图表。 我必须为400种物种生成双y轴图。 由于它们必须是双y轴图,我没有使用ggplot2。 我设法在Rstudio中看到图形,但是当我尝试打开PDF(我也尝试过TIFF,JPEG,...)时,它们无法打开。 (免费翻译错误代码:由于不支持此类文件或文件已损坏,因此无法打开)
我阅读了很多关于使用dev.off()结束代码的帖子,但这没有帮助。 这是代码。使用相同代码(但没有循环)生成的单个文件打开正常...
for (i in categories){
plot_data<-subset(af, nummer_soort_naam == i)
par(mar = c(5,5,2,5))
with(plot_data,
plot(jaar,
aantal_waargenomen_individuen_per_duizend_daghokbezoeken,
type="l"))# pch=16))}
par(new = T)
with(plot_data,
plot(jaar, ruwe_aantal_ind_waargenomen,
pch=16, axes=F, xlab=NA, ylab=NA, cex=1.2,
main=paste(i,sep=" ")))
axis(side = 4)
mtext(side = 4, line = 3, 'ruwe aantal individuen waargenomen')
legend("topleft",
legend=c("gecorrigeerd aantal","ruwe aantallen"),
lty=c(1,0),
pch=c(NA, 16),
col=c("black", "black"))
dev.print(file=paste(i,".pdf",sep=""))
dev.off()
}
数据如下所示:
soort_id nummer_soort_naam jaar ruwe_aantal_ind_waargenomen aantal_waargenomen_individuen_per_duizend_daghokbezoeken
1 1_Roerdomp_Botaurus_stellaris 1990 6 3.84
1 1_Roerdomp_Botaurus_stellaris 1991 12 7.89
1 1_Roerdomp_Botaurus_stellaris 1992 10 6.9
1 1_Roerdomp_Botaurus_stellaris 1993 7 3.7
1 1_Roerdomp_Botaurus_stellaris 1994 37 20.19
2 2_Dodaars_Tachybaptus_ruficollis 1990 399 248.6
2 2_Dodaars_Tachybaptus_ruficollis 1991 569 365.9
2 2_Dodaars_Tachybaptus_ruficollis 1992 560 376.1
2 2_Dodaars_Tachybaptus_ruficollis 1993 312 163
2 2_Dodaars_Tachybaptus_ruficollis 1994 373 198.3
有人可以根据此代码帮助我并查看错误吗?
非常感谢你!
Kristijn