Question

我正在尝试使用ggplot2制作雷达图。

> Test1
   Country fertility gengap   LEB
1 Colombia      1.92   0.29 79.30
2     Peru      1.94   0.37 78.40
3      Min      9.23   3.83  0.00
4      Max      1.59  23.70 78.43

#normalizing the data:

Test1_norm<- data.frame(sapply(Test1[1:4, 2:4], scales::rescale))

    > Test1_norm
   fertility      gengap       LEB
1 0.04319372 0.000000000 1.0000000
2 0.04581152 0.003417343 0.9886507
3 1.00000000 0.151217428 0.0000000
4 0.00000000 1.000000000 0.9890290

我发现此代码用于绘制雷达图： http://www.cmap.polytechnique.fr/~lepennec/R/Radar/RadarAndParallelPlots.html 对于mtcars-example：

   ggplot(mtcarsmelted, aes(x = variable, y = value)) +
  geom_path(aes(group = model, color = model),  size = 2) +
  theme(strip.text.x = element_text(size = rel(0.8)),
        axis.text.x = element_text(size = rel(0.8)),
        axis.ticks.y = element_blank(),
        axis.text.y = element_blank()) +
  xlab("") + ylab("") +
  guides(color = guide_legend(ncol=2)) +
  coord_polar()

使用mtcars示例一切正常，但是当我尝试重写数据代码时，我只是得到错误。我现在正在尝试几个小时，我不确定我需要交换哪些变量以及哪些变量必须保留在示例中。

我很感激一点帮助！我还是R的新手;）

谢谢大家：）

Answer 1

首先，使用普通数据，我尝试根据您显示给我的方法绘制雷达图。 reshape2包是你需要的。这是代码：

Test1_norm <- data.frame(group=c("Colombia","Peru","Min","Max"),
              fertility=c(0.04319372,0.04581152,1,0),
              gengap=c(0,0.003417343,0.151217428,1),
              LEB=c(1,0.9886507,0,0.9890290))
Test1_norm$model <- Test1_norm$group
library(reshape2)
Test1_plot<- melt(Test1_norm)

##Define a radar coordinate system
coord_radar <- function (theta = "x", start = 0, direction = 1) 
{
    theta <- match.arg(theta, c("x", "y"))
    r <- if (theta == "x") 
        "y"
    else "x"
    ggproto("CordRadar", CoordPolar, theta = theta, r = r, start = start, 
        direction = sign(direction),
        is_linear = function(coord) TRUE)
}

#Plot
library(ggplot2)
ggplot(Test1_plot, aes(x = variable, y = value)) +
geom_polygon(aes(group = model, color = model), fill = NA, size = 2, show.legend = FALSE) +
  geom_line(aes(group = model, color = model), size = 2) +
  theme(strip.text.x = element_text(size = rel(0.8)),
        axis.text.x = element_text(size = rel(0.8)),
        axis.ticks.y = element_blank(),
        axis.text.y = element_blank()) +
  xlab("") + ylab("") +
  guides(color = guide_legend(ncol=2)) +
  coord_radar()

幸运的是，Ricardo Bion提供了一种使用ggplot2创建雷达图表的有效方法。这可能对您有所帮助：https://github.com/ricardo-bion/ggradar 您可以按照手册中的说明安装软件包并调用ggradar()功能。这是我的原始数据代码。

Test1 <- data.frame(fertility=c(1.92,1.94,9.23,1.59),
                    gengap=c(0.29,0.37,3.83,23.70),
                     LEB=c(79.30,78.40,0.00,78.43))
rownames(Test1)<-c("Colombia","Peru","Min","Max")

source("ggradar.R")
library(dplyr)
library(ggplot2)
library(scales)
Test1 %>%
      add_rownames( var = "group" ) %>%
      mutate_each(funs(rescale), -group)  -> Test1_radar
ggradar( Test1_radar)
ggsave("Test1_radar.tiff",width = 8, height = 8)

雷达图/蜘蛛图与ggplot2

1 个答案: