2013-08-01 22 views
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軸標籤中包含「十一=頭,」不喜價值在劇情代碼中提供。

我該如何獲得一個標題,以及將這些圖表的參數值合併到單個文件中?

+0

main = paste(文件[i],其他任何未引用的變量,「StaticWordsinTitleQuoted」,sep =「」) – cianius

+0

是否希望希臘字母代表tau和sigma?此外,dgpd似乎是未被該代碼加載的軟件包的功能。 –

回答

1

您正在尋找paste()函數,而不是cat()函數。 cat()函數輸出一個對象,而paste()在轉換爲字符後連接一個向量。

plot(returns.sorted, dgpd(returns.sorted, xi = tau, beta = sigma), 
    type ="l", col="blue", ylim=c(0,90), 
    main = paste("Tau= ", tau, " Sigma= ", sigma), ylab=" ") 

如果你想希臘字母看到this問題或this問題。