Loops 成圈散布

我正在尝试使用循环为文件夹中的每个.csv文件绘制一个ggscatter。它生成第一个绘图，但没有回归线和p值。但是，它不会为后续文件生成打印。我收到以下错误：“此外：警告消息：1:In read.table（file=file，header=header，sep=sep，quote=quote，： 第1行似乎包含嵌入的null“ 请看下面我的当前脚本。第一个循环自动从USGS网站下载streamflow数据，工作正常。问题在第二个循环中，我正在尝试绘制下载的数据

library(dataRetrieval)
library(EcoHydRology)
library(ggpubr)
library("readxl")
data <- read_excel("StationsID.xlsx")
stations=data.frame(data)
head(stations)
for (i in stations[,]) {
siteNumber <- i
discharge <- readNWISsite(siteNumber)
parameterCd <- "00060" #This code means mean daily 

#Daily data:
DailyData <- readNWISdv(siteNumber,parameterCd,
                           "1970-01-01","2020-01-01")
DailyData <- renameNWISColumns(DailyData) #renaming the default column names to more meaningful names
#names(DailyData)
#DailyData
df=data.frame(DailyData)
write.csv(x=df,file=paste(siteNumber,".csv",sep=""))
}
#Plotting the data
library(stringr)
files<-list.files(pattern = ".csv", full.names = T)
for (a in files){
  myData<-read.csv(a);
  fname <- basename(a);
  tiff(paste(fname,".tiff",sep=""), units="in", width=10, height=5, res=1000)
  q = ggscatter(myData, x = "Date", y = "Flow", shape = 19, xlab = "", ylab = "Discharge (ft³/s)", main=paste(fname),
                  add = "reg.line",                                
                  conf.int = TRUE,                                
                 add.params = list(color = "blue",
                                    fill = "lightblue")
  )+
    stat_cor(method = "pearson", label.x = min(myData$Date), label.y = max(myData$Flow*1.05))+
    theme(plot.title = element_text(hjust = 0.5))
  print(q)
  dev.off()
}

库（数据检索）
图书馆（生态水文学）
图书馆（ggpubr）
图书馆（“readxl”）
数据