1
我試圖用不同的參數值創建大量圖。當我提交代碼如下:如何在循環中的繪圖標籤中加入變量標題和變量值
for(i in mars){
for(j in 1:5){
negRets <- -table1.matrix[,j]
tauSigma <- gpd(negRets,threshold=i)
tau <- tauSigma$par.ests[1]
sigma <- tauSigma$par.ests[2]
#cat("For ", reitPort[j], "portfolio and MAR=", i, "the parameter estimats are:", "\n")
#cat("Tau= ", tau, "Sigma= ", sigma, "\n")
exceedence.vector <- sort(subset(negRets, negRets > i))
returns.sorted <- sort(negRets)
pdf(paste("C:\\Users\\John Broussard\\Dropbox\\evtHandbookProject\\figuresTables\\mar_",i,"_",reitPort[j],".pdf",sep=""))
plot(returns.sorted, dgpd(returns.sorted, xi = tau, beta = sigma), type ="l", col="blue", ylim=c(0,90))
title(main=cat("Tau= ", tau, " Sigma= ", sigma), ylab=" ")
dev.off()
#print(j)
#print(i)
}
}
該地塊被寫入文件,但沒有標題,和Y軸標籤中包含「十一=頭,」不喜價值在劇情代碼中提供。
我該如何獲得一個標題,以及將這些圖表的參數值合併到單個文件中?
main = paste(文件[i],其他任何未引用的變量,「StaticWordsinTitleQuoted」,sep =「」) – cianius
是否希望希臘字母代表tau和sigma?此外,dgpd似乎是未被該代碼加載的軟件包的功能。 –