如何在泰勒圖中標註點？

plotrix程序包有一個叫做taylor.diagram的函數，它繪製兩個向量 - 一個表示數據和另一個模型輸出。如何在泰勒圖中標註點？

下面是一個例子：

require(plotrix) 
set.seed(10) 
data <- sort(runif(100, 8,12)) 
model <- sort(rnorm(100, 10, 4)) 
taylor.diagram(data, model)

而在這個例子中，我想提高模型後更新的情節：

model2 <- sort(rnorm(100, 10,2)) 
taylor.diagram(data, model2, add = TRUE)

要產生這樣的：

enter image description here

如何添加「模型1」和「模型2」等標籤以識別t這些點？（更新：從而不做，事後的模型值確定標籤位置）

來源

2013-02-28 Abe

您的標籤在基地的圖形都用同樣的方法，用text：

text(1.5,0.5,labels = "Model2") 
text(3.5,1,labels = "Model1")

enter image description here

來源

2013-02-28 23:27:49 joran

我希望能夠繪製他們沒有他們的位置 – Abe 2013-03-01 00:07:40

@Abe的文檔的預先了解'taylor.diagram'表示它不返回點的位置（或者根本沒有任何繪圖信息），所以我不認爲你在這裏有很多選擇。 – joran 2013-03-01 00:57:08

這裏有兩種方法

example(taylor.diagram)顯示了一個體面的方式來放置一個傳說右上角（在1.5*sd(data), 1.5*sd(data)），但這需要兩點不同的顏色。
另一種選擇是基於從原來的Taylor 2001 reference的公式來計算位置 - 或從源代碼到taylor.diagram功能複製它們，近
```
dy <- 1.1 # text offset coefficient 
sd.f <- sd(model) 
R <- cor(data, model, use = 'pairwise') 
x <- sd.f * R 
y <- sd.f * sin(acos(R)) + dy * sd.f 
text(x, y, "Model") 
```
您需要計算這些各型號，但只有模型輸入和標籤會改變。你可能想保持偏移量一樣。

來源

2013-03-01 16:44:21

第三個解決方案是創建函數taylor.diagram包括文本標籤的修改後的版本。在這種情況下，所有要做的就是添加一個參數，例如text，並且在原始函數（右花括號之前的兩行）中調用points之後，行text(sd.f * R, sd.f * sin(acos(R)), labels=text, pos=3)。

taylor.diagram.modified <- function (ref, model, add = FALSE, col = "red", 
            pch = 19, pos.cor = TRUE, xlab = "", ylab = "", 
            main = "Taylor Diagram", show.gamma = TRUE, 
            ngamma = 3, gamma.col = 8, sd.arcs = 0, ref.sd = FALSE, 
            grad.corr.lines = c(0.2, 0.4, 0.6, 0.8, 0.9), pcex = 1, 
            cex.axis = 1, normalize = FALSE, mar = c(5, 4, 6, 6), 
            text, ...) #the added parameter 
{ 
    grad.corr.full <- c(0, 0.2, 0.4, 0.6, 0.8, 0.9, 0.95, 0.99,1) 
    R <- cor(ref, model, use = "pairwise") 
    sd.r <- sd(ref) 
    sd.f <- sd(model) 
    if (normalize) { 

    ... #I didn't copy here the full function because it's quite long: to obtain it 
    ... #simply call `taylor.diagram` in the console or `edit(taylor.diagram)`. 

      } 
      S <- (2 * (1 + R))/(sd.f + (1/sd.f))^2 
     } 
    } 
    points(sd.f * R, sd.f * sin(acos(R)), pch = pch, col = col, 
      cex = pcex) 
    text(sd.f * R, sd.f * sin(acos(R)), #the line to add 
     labels=text, cex = pcex, pos=3) #You can change the pos argument to your liking 
    invisible(oldpar) 
}

然後，只需在text論點提供標籤名稱：

require(plotrix) 
set.seed(10) 
data <- sort(runif(100, 8,12)) 
model <- sort(rnorm(100, 10, 4)) 
taylor.diagram.modified(data, model, text="Model 1") 
model2 <- sort(rnorm(100, 10,2)) 
taylor.diagram.modified(data, model2, add = TRUE, text="Model 2")

enter image description here

來源

2013-03-11 09:25:03 plannapus

如何在泰勒圖中標註點？

回答

相關問題