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 轴标签包含“xi=tau”，而不是绘图代码中提交的 xi 值。

如何获得标题和用于将绘图合并到单个文件中的参数值？

4

1 回答 1

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=" ")

如果您想要希腊字母，请参阅此问题或此问题。

于 2013-08-01T22:46:52.370 回答