เป็นไปได้ไหมที่จะสร้าง“ ชุดขนาน” โดยใช้ R


16

ขอบคุณคำถาม Tormod (โพสต์ที่นี่ ) ฉันเจอพล็อตชุดแบบขนาน นี่คือตัวอย่างสำหรับลักษณะที่ปรากฏ: ป้อนคำอธิบายรูปภาพที่นี่ (เป็นภาพของชุดข้อมูลไททานิกการแสดงยกตัวอย่างเช่นผู้หญิงส่วนใหญ่ที่ไม่รอดชีวิตอยู่ในชั้นที่สาม ... )

ฉันชอบที่จะสามารถทำซ้ำพล็อตดังกล่าวกับอาร์ได้หรือไม่?

ขอบคุณ Tal


1
สำหรับแนวคิดเกี่ยวกับกราฟิกฉันจะตรวจสอบแกลเลอรีกราฟ R เสมอ นี่คือสิ่งจากที่นั่นที่เป็นเหมือนสิ่งที่คุณขอ: R กราฟแกลลอรี่ขนาน ฉันพบมันโดยคลิกขนานใน tag cloud แต่อาจมีตัวเลือกที่ดีกว่า
Nick Sabbe

1
ขอบคุณนิค แต่สิ่งนี้จะไม่ทำงานสำหรับข้อมูลที่เป็นหมวดหมู่หากไม่มีการปรับแต่งโค้ดให้ใหญ่ (อาจเป็นฟังก์ชั่นพื้นฐานที่ดีที่สุดในการสร้างด้วย) ฉันหวังว่าบางคนอาจจะได้ทำสิ่งที่คล้ายกันอยู่แล้ว ...
Tal Galili

คำตอบ:


25

นี่คือรุ่นที่ใช้กราฟิกพื้นฐานเท่านั้นขอบคุณความคิดเห็นของ Hadley (สำหรับรุ่นก่อนหน้าดูประวัติการแก้ไข)

ลองครั้งที่สาม

parallelset <- function(..., freq, col="gray", border=0, layer, 
                             alpha=0.5, gap.width=0.05) {
  p <- data.frame(..., freq, col, border, alpha, stringsAsFactors=FALSE)
  n <- nrow(p)
  if(missing(layer)) { layer <- 1:n }
  p$layer <- layer
  np <- ncol(p) - 5
  d <- p[ , 1:np, drop=FALSE]
  p <- p[ , -c(1:np), drop=FALSE]
  p$freq <- with(p, freq/sum(freq))
  col <- col2rgb(p$col, alpha=TRUE)
  if(!identical(alpha, FALSE)) { col["alpha", ] <- p$alpha*256 }
  p$col <- apply(col, 2, function(x) do.call(rgb, c(as.list(x), maxColorValue = 256)))
  getp <- function(i, d, f, w=gap.width) {
    a <- c(i, (1:ncol(d))[-i])
    o <- do.call(order, d[a])
    x <- c(0, cumsum(f[o])) * (1-w)
    x <- cbind(x[-length(x)], x[-1])
    gap <- cumsum( c(0L, diff(as.numeric(d[o,i])) != 0) )
    gap <- gap / max(gap) * w
    (x + gap)[order(o),]
  }
  dd <- lapply(seq_along(d), getp, d=d, f=p$freq)
  par(mar = c(0, 0, 2, 0) + 0.1, xpd=TRUE )
  plot(NULL, type="n",xlim=c(0, 1), ylim=c(np, 1),
       xaxt="n", yaxt="n", xaxs="i", yaxs="i", xlab='', ylab='', frame=FALSE)
  for(i in rev(order(p$layer)) ) {
     for(j in 1:(np-1) )
     polygon(c(dd[[j]][i,], rev(dd[[j+1]][i,])), c(j, j, j+1, j+1),
             col=p$col[i], border=p$border[i])
   }
   text(0, seq_along(dd), labels=names(d), adj=c(0,-2), font=2)
   for(j in seq_along(dd)) {
     ax <- lapply(split(dd[[j]], d[,j]), range)
     for(k in seq_along(ax)) {
       lines(ax[[k]], c(j, j))
       text(ax[[k]][1], j, labels=names(ax)[k], adj=c(0, -0.25))
     }
   }           
}

data(Titanic)
myt <- subset(as.data.frame(Titanic), Age=="Adult", 
              select=c("Survived","Sex","Class","Freq"))
myt <- within(myt, {
  Survived <- factor(Survived, levels=c("Yes","No"))
  levels(Class) <- c(paste(c("First", "Second", "Third"), "Class"), "Crew")
  color <- ifelse(Survived=="Yes","#008888","#330066")
})

with(myt, parallelset(Survived, Sex, Class, freq=Freq, col=color, alpha=0.2))

Aaron, wow, คำตอบที่ยอดเยี่ยม - ฉันหวังว่าฉันสามารถทำเครื่องหมาย V ได้สองครั้ง ขอขอบคุณ!
Tal Galili

2
ดีใจที่คุณชอบมัน มันสนุก. :) ส่วนที่ยุ่งยากเท่านั้นคือการรับสถานที่ที่บาร์ควรเริ่มต้นและสิ้นสุด (ซึ่งอยู่ในgetpฟังก์ชั่นย่อย); ส่วนที่เหลือเป็นเพียงการวาดรูปหลายเหลี่ยม
แอรอน - คืนสถานะโมนิก้า

1
แค่อีกpanel.textบรรทัด ดูการแก้ไข
แอรอน - คืนสถานะโมนิก้า

1
คุณสามารถโปร่งใสในกราฟิกพื้นฐานได้เช่นกัน
hadley

2
ใช่คุณ ฉันลืมไปอย่างสิ้นเชิงเกี่ยวกับเรื่องนี้การคุ้นเคยกับวิธีขัดแตะในการทำสิ่งต่าง ๆ สำหรับคนอื่น ๆ #FF000080ที่มีความสนใจคุณเพิ่มตัวละครคู่มากขึ้นบนสตริงสีของคุณตัวอย่างเช่น ?rgbมีรายละเอียด
Aaron - Reinstate Monica

12

ตามรหัส @Aaron ฉันพัฒนาบางสิ่งที่เรียกว่า "แผนภาพลุ่มน้ำ" ดูhttp://bc.bojanorama.pl/2014/03/alluvial-diagrams/ตัวอย่างด้านล่าง:

ป้อนคำอธิบายรูปภาพที่นี่


ฉันเห็นโพสต์ของคุณบน r-blogges แล้ว ขอบคุณมีคาล
Tal Galili

1
แพ็คเกจอย่างเป็นทางการ: cran.r-project.org/web/packages/alluvial
Orion
โดยการใช้ไซต์ของเรา หมายความว่าคุณได้อ่านและทำความเข้าใจนโยบายคุกกี้และนโยบายความเป็นส่วนตัวของเราแล้ว
Licensed under cc by-sa 3.0 with attribution required.