Как сохранить равные относительные размеры точек на графике?

Question

Я хочу сделать график, где размеры кружков указывают размер образца. если я использую сюжет в p1(), он работает нормально.

но если я попытаюсь раскрасить другой тип точек, то относительный размер будет неправильным.

Как бы получить красные и зеленые круги одинакового размера?

p1<-function() {
    plot(t$x,t$y,cex=100*t$size,xlim=c(0,1),ylim=c(0.,1.))    
}
p2<-function() {
    plot(t$x[t$r=="0"],t$y[t$r=="0"],xlim=c(0,1),ylim=c(0.,1.),cex=100*t$size,col="red")
    points(t$x[t$r=="1"],t$y[t$r=="1"],xlim=c(0,1),ylim=c(0.,1.),cex=100*t$size,col="green")
}
l<-20
x<-seq(0,1,1/l)
y<-sqrt(x)
r=round(runif(n=length(x),min=0,max=.8))
n<-1:length(x)
size=n/sum(n)
t<-data.frame(x,y,r,n,size)
t$r<-factor(r)
str(t)
p1()

Answer 1

Вы должны изменить функцию p2 немного. Ты используешь t$sizeвсе это, когда вы должны быть подмножеством фактора t$r, так как вы делаете это при построении точек.

Если вы заговор t$x[t$r == "0"] против t$y[t$r == "0"] тогда вы должны использовать размеры, соответствующие тем точкам, которые t$size[t$r == "0"], В качестве альтернативы, вы можете установить подкадр данных t сначала, а затем используйте эти два результирующих кадра данных для построения точек. Смотрите функцию p2_alt в конце.

p2 <- function() {
  plot(t$x[t$r == "0"], t$y[t$r == "0"],
       xlim = c(0, 1), ylim = c(0., 1.),
       cex = 100*t$size[t$r == "0"],
       col = "red",
       xlab = "x", ylab = "y")
  points(t$x[t$r == "1"],
         t$y[t$r == "1"],
         xlim = c(0, 1), ylim = c(0., 1.),
         cex = 100*t$size[t$r == "1"],
         col = "green")
}

set.seed(651)    # make the results reproducible

l <- 20
x <- seq(0, 1, 1/l)
y <- sqrt(x)
r <- round(runif(n = length(x), min = 0, max = 0.8))
n <- 1:length(x)
size <- n/sum(n)
t <- data.frame(x, y, r, n, size)
t$r <- factor(r)
#str(t)
#p1()
p2()

p2_alt <- function() {
  df1 <- subset(t, r == "0")
  df2 <- subset(t, r == "1")
  plot(df1$x, df1$y,
       xlim = c(0, 1), ylim = c(0., 1.),
       cex = 100*df1$size,
       col = "red",
       xlab = "x", ylab = "y")
  points(df2$x,
         df2$y,
         xlim = c(0, 1), ylim = c(0., 1.),
         cex = 100*df2$size,
         col = "green")
}

p2_alt()

График точно такой же, но, возможно, код более читабелен.

Наконец, обратите внимание, что я добавил аргументы xlab а также ylab как для p2() а также p2_alt(),