请教这样的统计图是怎么画出来的

相似文件 换一批

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏2 回帖

关键词：统计图 统计图

理论上用TIKZ应该是可以的

shizidushu 发表于 2017-4-26 11:19
理论上用TIKZ应该是可以的

latex吗

是的。理论上是latex的tikz包可以制作出一模一样的，tikz的官方文档有1k多页，慢慢看吧。另外如果只是要画出一个类似的图，可以有很多其他解决方案的，比如gnuplot（听说过，不了解）；r语言的ggplot2包；

1. # Invent some dummy data according to the plots
   
2. dfDmmy <-
   
3.   merge(
   
4.     data.frame(
   
5.       treatment =
   
6.         factor( x = c(1:4), levels = c(1:4),
   
7.                 labels = c("C", "T1", "T2", "T3") )
   
8.               ),
   
9.     data.frame(
   
10.       requirement =
    
11.         factor(
    
12.           x = c(1:5), levels = c(1:5),
    
13.           labels =
    
14.             paste(
    
15.               "Requires", c("Name", "Phone Number", "Email", "Address", "ID"),
    
16.               sep = "\n"
    
17.                  )
    
18.               )
    
19.               )
    
20.        )
    
21. 
    
22. set.seed(1705171044)
    
23. 
    
24. dfDmmy$mid <-
    
25.   mapply(
    
26.     FUN = rnorm,
    
27.     mean = c(0.8, 0.6, 0.5, 0.3, 0.1)[ as.integer(dfDmmy$requirement) ],
    
28.     MoreArgs = list( n = 1, sd = 0.01 )
    
29.         )
    
30. 
    
31. dfDmmy$dev <-
    
32.   c(0.05, 0.05, 0.05, 0.05, 0.025)[ as.integer(dfDmmy$requirement) ]
    
33. 
    
34. dfDmmy$lower <- with(dfDmmy, mid - dev)
    
35. 
    
36. dfDmmy$upper <- with(dfDmmy, mid + dev)
    
37. 
    
38. with(
    
39.   dfDmmy,
    
40.   {
    
41.     # Set up plot
    
42.     plot(
    
43.       NA, type = "n",
    
44.       xlim = c( 0, ( nlevels(treatment) + 1 )*nlevels(requirement) ),
    
45.       ylim = c(0, 1), axes = FALSE,
    
46.       xlab = NA_character_, ylab = "Proportion of Counties with Forums"
    
47.         )
    
48.    
    
49.     # Draw axis
    
50.     axis(
    
51.       side = 1,
    
52.       at =
    
53.         ( as.integer(requirement) - 1 )*( nlevels(treatment) + 1 ) +
    
54.         as.integer(treatment),
    
55.       labels = as.character(treatment), mgp = c(3, 0.75, 0)
    
56.         )
    
57.    
    
58.     axis(
    
59.       side = 1,
    
60.       at =
    
61.         seq.int( nlevels(requirement) )*( nlevels(treatment) + 1 ) -
    
62.         ( nlevels(treatment) + 1 )/2,
    
63.       labels = levels(requirement), mgp = c(9, 3, 0), tick = FALSE
    
64.         )
    
65.    
    
66.     axis( side = 2, at = c(0:5/5) )
    
67.    
    
68.     box()
    
69.    
    
70.     # Separate different requirements
    
71.     abline(
    
72.       v = seq.int( nlevels(requirement) - 1 )*( nlevels(treatment) + 1 )
    
73.           )
    
74.          
    
75.     # Add intervals
    
76.     mapply(
    
77.       FUN = segments,
    
78.       x0 =
    
79.         ( as.integer(requirement) - 1 )*( nlevels(treatment) + 1 ) +
    
80.         as.integer(treatment),
    
81.       y0 = lower, y1 = upper
    
82.           )
    
83.    
    
84.     # Add middle points
    
85.     mapply(
    
86.       FUN = points,
    
87.       x =
    
88.         ( as.integer(requirement) - 1 )*( nlevels(treatment) + 1 ) +
    
89.         as.integer(treatment),
    
90.       y = mid, pch = c(18, 15:17)[ as.integer(treatment) ]
    
91.           )
    
92.   }
    
93.     )

复制代码