2012-05-24 15 views
11

mam dane, które wygląda tak:Jak wykreślić wiele linii w R

#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 

więc chcę wykreślić je w oparciu o prawdziwe (oś Y) i false (oś X).

W ten sposób chcę, aby wyglądało mniej więcej. This is the way I want it to appear roughly

Co to jest właściwy sposób? poniżej Mój kod nie

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]); 
} 

Odpowiedz

7

Ponieważ dane są już w długim formacie i lubię grafiki ggplot każdym razie, polecam tę ścieżkę. Po odczytaniu danych w (zauważ, że TRUE i FALSE nie są poprawne nazwy, więc R dołączona . do nazw kolumn) dodaje powinno działać:

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

ggplot2 strona jest pełna dobrych porad. Zauważ także this search query on SO, aby uzyskać wiele innych dobrych wskazówek na podobne tematy.

enter image description here

A dla przypomnienia, oto jak bym podejść problemu modyfikowania oryginalnego kodu:

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) 

enter image description here

+0

dzięki. Po odczytaniu danych. Dał mi ten błąd: Błąd w eval (expr, envir, enclos): obiekt "FALSE". nie znaleziono. – neversaint

+1

@neversaint - Domyślam się, że R ma problemy, ponieważ TRUE i FALSE są zarezerwowanymi słowami w R. Zmień nazwy kolumn na coś takiego jak T i F i wszystko powinno być w porządku. Nie powinienem był nawet próbować podawać tego przykładu razem z tymi nazwami ... po prostu sprawi ci to kłopot. 'colnames()' jest funkcją, której szukasz. – Chase

6

Oto metoda baza R, zakładając, że dane jest nazywany dat w tym przykładzie:

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]) 
) 

Wynik:

enter image description here

edycja - podczas korzystania małe true/false dla nazw zmiennych zostanie zaakceptowana, to chyba jeszcze nie jest najlepsza praktyka kodowania.