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 に答える