0

以下可重现的示例创建了一个条件密度图，描述了各个方面的变化。我在示例中包含了 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 轴上的数值替换为相关的文本值（例如 N、SW 等）？
  • 是否可以控制 X 轴上文本值的顺序？我相信，默认情况下，X 轴上的值对应于“方面”中按字母顺序排列的文本值。
  • 如何将 X 轴上的“NA”替换为“No Change”？
4

1 回答 1

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
于 2013-02-02T18:25:05.223 回答