如何在R中重新格式化條件密度圖？

Question

以下可重現的示例創建了一個條件密度圖，該圖描述了方面的變化。 我在示例中包括了NoData，因為我的原始數據集也具有NoData。

aspect = c("N", "N", "E", "S", "SW", "SW", "E", "W", "N", "N")
change = c(-1, NA, 1, NA, 1, 1, -1, 1, -1, NA)
mydata = data.frame(aspect,change)

x = factor(mydata$change)

cdplot(x~mydata$aspect)

相關問題：

• 如何用相關的文本值（例如N，SW等）替換X軸上的數值？
• 是否可以控制X軸上文本值的順序？ 我相信默認情況下，X軸上的值對應於“方面”中按字母順序排列的文本值。
• 如何將X軸上的“ NA”替換為“ No Change”？

Answer 1

只是為了顯示@Roland在評論中所說的內容，您以錯誤的方式使用了該函數。

aspect = c("N", "N", "E", "S", "SW", "SW", "E", "W", "N", "N")
change = sample(rep(c(-1,1,NA), each = 100),100,replace=T)
mydata = data.frame(aspect,change,stringsAsFactors=T)
str(mydata)
'data.frame':   100 obs. of  2 variables:
 $aspect: Factor
$change: num 
h <- cdplot(aspect~change, data = mydata)

現在，h包含了方面方面的條件密度函數。

str(h)
List of 4
 $ E :function (v)  
 $ N :function (v)  
 $ S :function (v)  
 $ SW:function (v)  

例如

h$E(-Inf)
[1] 0.21875
h$E(Inf)
[1] 0.25