0

我正在尝试用 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 示例:

   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 还是很陌生;)

谢谢大家 :)

4

1 回答 1

0

首先,使用正常数据,我尝试根据您向我展示的方法绘制雷达图。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)
于 2017-01-20T06:54:22.737 回答