Добавить легенду на график плотности в R - PullRequest
Новая
       10

Добавить легенду на график плотности в R

0 голосов
/ 03 октября 2011

Я использую следующий код в R для рисования двух кривых плотности на одном графике: 

mydata1<-read.csv(file="myfile1.csv",head=TRUE,sep=",")
mydata2<-read.csv(file="myfile2.csv",head=TRUE,sep=",")

pdf("comparison.pdf")

plot.multi.dens <- function(s)   
{
    junk.x = NULL
    junk.y = NULL
    for(i in 1:length(s)) {
        junk.x = c(junk.x, density(s[[i]])$x)
        junk.y = c(junk.y, density(s[[i]])$y)
    }
    xr <- range(junk.x)
    yr <- range(junk.y)
    plot(density(s[[1]]), xlim = xr, ylim = yr, xlab="Usage",main = "comparison")
    for(i in 1:length(s)) {
        lines(density(s[[i]]), xlim = xr, ylim = yr, col = i)
    }
}

plot.multi.dens( list(mydata2$usage,mydata1$usage))    
dev.off()

Теперь проблема в том, что производимый график показывает две линии, но график не 't включает в себя информацию о том, какая строка есть какая.Например, в выходных данных он должен показать, что красная линия - это «а», а черная - «b».Я новичок в R, поэтому у меня возникли некоторые трудности.любая помощь будет оценена!

1 Ответ

1 голос
/ 03 октября 2011

Ответ с сайта quickR 

# Compare MPG distributions for cars with 
# 4,6, or 8 cylinders
library(sm)
attach(mtcars)

# create value labels 
cyl.f <- factor(cyl, levels= c(4,6,8),
    labels = c("4 cylinder", "6 cylinder", "8 cylinder")) 

# plot densities 
sm.density.compare(mpg, cyl, xlab="Miles Per Gallon")
title(main="MPG Distribution by Car Cylinders")

# add legend via mouse click
colfill<-c(2:(2+length(levels(cyl.f)))) 
legend(locator(1), levels(cyl.f), fill=colfill)
...