如何在R中重新格式化条件密度图？

Question

以下可重现的示例创建了conditional density plot描述方面的更改。我在示例中包含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)

相关问题：

• 如何用X轴替换X轴上的数值 相关的文本值（例如N，SW等）？
• 是否可以控制文本值的顺序 X轴？我相信，默认情况下，X轴上的值对应 “aspect”中按字母顺序排列的文本值。
• 如何用“无变化”替换X轴上的“NA”？

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)  

e.g

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