11

次のようなデータがあります。

#d  TRUE    FALSE   Cutoff
4   28198   0   0.1
4   28198   0   0.2
4   28198   0   0.3
4   28198   13  0.4
4   28251   611 0.5
4   28251   611 0.6
4   28251   611 0.7
4   28251   611 0.8
4   28251   611 0.9
4   28251   611 1
6   19630   0   0
6   19630   0   0.1
6   19630   0   0.2
6   19630   0   0.3
6   19630   0   0.4
6   19636   56  0.5
6   19636   56  0.6
6   19636   56  0.7
6   19636   56  0.8
6   19636   56  0.9
6   19636   56  1

そこで、True(Y軸)とFalse(X軸)に基づいてプロットしたいと思います。

これが私が大まかに見せたい方法です。 これは私がそれを大まかに見せたい方法です

それを行う正しい方法は何ですか?以下の私のコードは失敗します

dat<-read.table("mydat.txt", header=F);
dis     <- c(4,6);
linecols <-c("red","blue");
plot(dat$V2 ~ dat$V3, data = dat,  xlim = c(0,611),ylim =c(0,28251), type="l")

for (i in 1:length(dis)){
datax <- subset(dat, dat$V1==dis[i], select = c(dat$V2,dat$V3))
lines(datax,lty=1,type="l",col=linecols[i]);
}
4

2 に答える 2

7

あなたのデータはすでに長いフォーマットであり、とにかくggplotグラフィックが好きなので、そのパスをお勧めします。データを読み込んだ後(とは有効な名前ではないため、Rが列名にaを追加したことに注意してくださいTRUEFALSE.次のように機能するはずです。

require(ggplot2)
ggplot(dat, aes(FALSE., TRUE., colour = as.factor(d), group = as.factor(d))) + 
  geom_line()

ggplot2のWebサイトには、優れたヒントが満載です。また、関連トピックに関する他の多くの優れたヒントについては、SOでのこの検索クエリに注意してください。

ここに画像の説明を入力してください

そして、記録のために、元のコードを変更する際の問題へのアプローチ方法は次のとおりです。

colnames(dat)[2:3] <- c("T", "F")

dis <- unique(dat$d)

plot(NA, xlim = c(0, max(dat$F)), ylim = c(0, max(dat$T)))
for (i in seq_along(dis)){
  subdat <- subset(dat, d == dis[i])
  with(subdat, lines(F,T, col = linecols[i]))
}
legend("bottomright", legend=dis, fill=linecols)

ここに画像の説明を入力してください

于 2012-05-24T02:38:06.510 に答える
6

datこの例でデータが呼び出されると仮定すると、基本的な R メソッドは次のようになります。

plot(1:max(dat$false), xlim = c(0,611),ylim =c(19000,28251), type="n")

apply(
rbind(unique(dat$d),1:2),
#the 1:2 here are your chosen colours
2,
function(x) lines(dat$false[dat$d==x[1]],dat$true[dat$d==x[1]],col=x[2])
)

結果:

ここに画像の説明を入力

編集 - 変数名に小文字の true/false を使用することは受け入れられますが、おそらくまだ最高のコーディング プラクティスではありません。

于 2012-05-24T02:59:08.043 に答える