你真的懂了R中的stem函数是如何绘制茎叶图的么？

本文原创，转载请注明出处，本人Q1273314690（交流学习）

 

哭晕

你真的学会了stem()函数了吗？

stem()函数的使用方法是：

stem(x, scale=1,width=80, atom=le-08)

其中x是数据向量．

scale控制绘出茎叶图的长度．

width绘图的宽度．

atom是容差，如果选择scale=2，即将10个个位数俞成两段，0～4为一段，5～9为另一段。

然而事实上，我经过反复的试验，发现width最好取较大的数，他既不表示数据的取值范围也不表示最长的那片叶子的长度，也不表示所有的数据的个数

 

With the default scale, you see that numbers left of the bar go up by two - hence anything after the 4 is forty-something（四十几） or fifty-something（五十几）:

 

1. d<-c(60,85,72,59,37,75,93,7,98,63,41,90,5,17,97)
2. > stem(d,scale=1)
4.   The decimal point is 1 digit(s) to the right of the |
5.  
6.   0|577
7.   2|7
8.   4|19
9.   6|0325
10.   8|50378

Using scale=2, you'll see numbers left of the bar go up by one, so now you can get exact reconstruction（再造） of your input, since you input integers:

 

1. > stem(d,scale=2)
2.  
3.   The decimal point is 1 digit(s) to the right of the |
4.  
5.   0|57
6.   1|7
7.   2|
8.   3|7
9.   4|1
10.   5|9
11.   6|03
12.   7|25
13.   8|5
14.   9|0378

 

Going further you can even split it by first and second five within each decade:

 

1. > stem(d,scale=4)
2.  
3.   The decimal point is 1 digit(s) to the right of the |
4.  
5.   0|57
6.   1|      
7.   1|7
8.   2|
9.   2|
10.   3|
11.   3|7
12.   4|1
13.   4|
14.   5|
15.   5|9
16.   6|03
17.   6|
18.   7|2
19.   7|5
20.   8|
21.   8|5
22.   9|03
23.   9|78

Stem plots do not always guarantee you can reproduce the data by reversing the process, and that's not what they're for.

 

This is probably more of a Stack Overflow question than a CV question, because it focuses on how R works and why, rather than the statistical aspects of stem and leaf plots. Nonetheless...

The way the function is coded is designed to shorten the length of the output, so that it better fits in the console. Few people, I believe, find that terribly helpful, or at least I don't. Just always remember to start with scale=2, and you may have to play with it further or adjust the width argument. Also, know that there is a fancier偏好者，（发烧友的比较级） version stem.leaf() in Rcmdr. 

 

 

总结：

R的stem函数其实是一个比较糟糕的设计，由于其设计的初衷是让其在控制台上能尽量简短的显示（否则控制台宽度不够），所以，当数据之间的差距较大的时候，就会出问题，他会跳着提升枝干，所以一般要设置sacle，而sacle设置的越大，分茎越多，精度越高，如果你的scale较小，他甚至会自动帮你的数据做四舍五入（这样会降低精度）

> test<-c(57,122,1000)

> stem(test,scale = 2)

 

  The decimal point is 2 digit(s) to the right of the |

 

   0 | 6

   1 | 2

   2 |

   3 |

   4 |

   5 |

   6 |

   7 |

   8 |

   9 |

  10 | 0

 

> stem(test,scale = 10)

 

  The decimal point is 1 digit(s) to the right of the |

 

    4 | 7

    6 |

    8 |

   10 |

   12 | 2

   14 |

   16 |

   18 |

   20 |

   22 |

   24 |

   26 |

   28 |

   30 |

   32 |

   34 |

   36 |

   38 |

   40 |

   42 |

   44 |

   46 |

   48 |

   50 |

   52 |

   54 |

   56 |

   58 |

   60 |

   62 |

   64 |

   66 |

   68 |

   70 |

   72 |

   74 |

   76 |

   78 |

   80 |

   82 |

   84 |

   86 |

   88 |

   90 |

   92 |

   94 |

   96 |

   98 |

  100 | 0

 

> stem(test,scale = 20)

 

  The decimal point is 1 digit(s) to the right of the |

 

    5 | 7

    6 |

    7 |

    8 |

    9 |

   10 |

   11 |

   12 | 2

   13 |

   14 |

   15 |

   16 |

   17 |

   18 |

   19 |

   20 |

   21 |

   22 |

   23 |

   24 |

   25 |

   26 |

   27 |

   28 |

   29 |

   30 |

   31 |

   32 |

   33 |

   34 |

   35 |

   36 |

   37 |

   38 |

   39 |

   40 |

   41 |

   42 |

   43 |

   44 |

   45 |

   46 |

   47 |

   48 |

   49 |

   50 |

   51 |

   52 |

   53 |

   54 |

   55 |

   56 |

   57 |

   58 |

   59 |

   60 |

   61 |

   62 |

   63 |

   64 |

   65 |

   66 |

   67 |

   68 |

   69 |

   70 |

   71 |

   72 |

   73 |

   74 |

   75 |

   76 |

   77 |

   78 |

   79 |

   80 |

   81 |

   82 |

   83 |

   84 |

   85 |

   86 |

   87 |

   88 |

   89 |

   90 |

   91 |

   92 |

   93 |

   94 |

   95 |

   96 |

   97 |

   98 |

   99 |

  100 | 0

来自为知笔记(Wiz)