[R] Sas by function in R
David L Carlson
dcarlson at tamu.edu
Thu Jan 3 17:02:23 CET 2013
You could try the command ?by to get information about the by() function in
R. To get many plots, it will be easier to get help, if you include sample
data (as you did) and provide that data in a way that makes it easy to bring
the data into R using dput():
> dput(example)
structure(list(plot = c(1L, 1L, 1L, 1L, 2L, 2L, 3L, 3L, 3L, 3L,
3L, 3L, 4L, 4L, 4L, 5L, 7L, 8L, 8L, 10L, 10L, 10L, 10L, 10L,
11L, 11L, 12L, 12L, 13L, 13L, 13L, 13L, 14L, 14L, 14L, 15L, 16L,
16L, 17L, 17L, 17L, 17L, 18L, 18L, 18L, 19L, 19L, 19L, 21L, 22L,
22L, 23L, 24L, 24L, 24L, 25L, 26L, 27L, 27L, 27L, 28L, 28L, 29L,
30L, 30L, 31L, 31L, 31L, 31L, 32L, 32L, 33L, 34L, 35L, 36L, 36L,
37L, 37L, 37L, 38L, 38L, 39L, 39L, 39L, 39L, 40L, 41L, 41L, 42L,
43L, 43L, 43L, 44L, 44L, 44L, 45L, 45L, 45L, 45L, 45L, 46L, 46L,
46L, 47L, 47L, 49L, 49L, 51L, 51L, 52L, 53L, 53L, 56L, 56L, 56L,
56L, 57L, 57L, 58L, 58L, 58L, 59L, 59L, 59L, 60L, 61L, 61L, 62L,
62L, 63L, 63L, 63L, 63L, 63L, 63L, 63L, 64L, 64L, 64L, 64L, 64L,
64L, 65L, 65L, 65L, 65L, 66L, 66L, 66L, 66L, 67L, 67L, 67L, 68L,
68L, 70L, 70L, 71L, 72L, 72L, 72L, 72L, 73L, 74L, 74L, 75L, 75L,
75L, 75L, 76L, 76L, 76L, 77L, 78L, 78L, 78L, 78L, 78L, 78L, 79L,
79L, 79L, 79L, 79L, 80L, 80L, 80L, 80L, 80L, 80L, 80L, 81L, 81L,
81L, 81L, 81L, 82L, 82L, 82L, 82L, 83L, 84L, 85L, 85L, 86L, 86L,
87L, 87L, 88L, 88L, 89L, 89L, 89L, 90L, 91L, 91L, 91L, 92L, 92L,
92L, 92L, 92L, 93L, 93L, 93L, 94L, 95L, 95L, 95L, 95L, 95L, 96L,
96L, 97L, 97L, 97L, 97L, 98L, 99L, 100L, 101L, 102L, 102L, 102L,
102L, 102L, 104L, 104L, 104L, 104L, 105L, 106L, 106L, 107L, 107L,
107L, 107L, 107L, 107L, 108L, 108L, 109L, 109L, 109L, 109L, 109L,
109L, 110L, 110L, 111L, 112L, 112L, 113L, 113L, 113L, 113L, 113L,
114L, 114L, 114L, 114L, 114L, 114L, 115L, 115L, 116L, 116L, 117L,
117L, 117L, 117L, 117L, 119L, 119L, 119L, 119L, 120L, 120L, 120L,
120L, 120L, 120L, 121L, 122L, 122L, 122L, 122L, 123L, 123L, 123L,
123L, 124L, 124L, 124L, 124L, 125L, 125L, 125L, 125L, 125L, 125L,
125L, 126L, 126L, 126L, 126L, 127L, 127L, 127L, 127L, 128L, 128L,
129L, 129L, 129L, 129L, 129L, 130L, 130L, 130L, 132L, 133L, 133L,
133L, 133L, 134L, 134L, 134L, 134L, 135L, 135L, 135L, 136L, 136L,
137L, 137L, 138L, 139L, 139L, 139L, 139L, 139L, 140L, 141L, 141L,
141L, 142L, 142L, 142L, 143L, 143L, 143L, 143L, 144L, 144L, 144L,
145L, 145L, 145L, 145L, 145L, 145L, 145L, 145L, 146L, 146L, 146L,
146L, 146L, 147L, 148L, 148L, 148L, 148L, 149L, 149L, 149L, 149L,
150L, 150L, 150L, 150L, 151L, 151L, 151L, 151L, 152L, 152L, 152L,
153L, 153L, 153L, 154L, 154L, 154L, 154L, 155L, 155L, 155L, 156L,
156L, 157L, 157L, 157L, 158L, 158L, 158L, 158L, 158L, 158L, 159L,
159L, 159L, 159L, 160L, 160L), d = c(16.3, 25, 57.8, 17, 10.8,
96.4, 76, 32, 11, 11, 106, 12.5, 19.3, 12, 15, 99.3, 11, 18.5,
77, 9.8, 101.5, 18, 6, 12.3, 80.8, 14, 12.5, 10.8, 14.5, 13.3,
12.5, 124.1, 92.3, 20.3, 9.5, 9, 104.1, 15, 62.8, 16.3, 95.2,
17.8, 14, 107.3, 22.5, 22.8, 16.5, 15.3, 47.8, 110.1, 14.3, 73.3,
14.5, 93.9, 68.5, 96.4, 44.5, 18, 7, 14, 13.8, 80, 43.5, 69.8,
12.8, 102.2, 15.8, 63, 63.5, 31.3, 6.3, 82.1, 112.7, 80, 5, 80,
11.5, 95.5, 31.3, 7, 54, 6.5, 20.5, 63, 60, 9, 44, 64, 70, 37,
57.3, 21.3, 10.5, 43.5, 33.8, 21.8, 19.3, 13.5, 77, 24, 16.3,
31.5, 8.3, 70.3, 17, 16.5, 7.8, 62.8, 64.8, 94.9, 94.5, 36.8,
19.5, 12.8, 62.5, 40, 17, 70.5, 10.5, 48, 59.8, 10.3, 16.8, 19.8,
108.5, 106, 22.8, 15.3, 17.3, 57.5, 15.8, 104.9, 7.8, 5, 5, 4.8,
6.3, 5, 9.5, 21.3, 5, 54, 17.5, 68, 6.8, 13.8, 80.5, 15.5, 28.8,
12.5, 27, 21.8, 25.3, 85.3, 18, 10.8, 14.8, 13, 90.1, 33.8, 7.8,
14.3, 72, 66, 24.8, 14.8, 93.3, 12.5, 15, 22, 90.4, 53, 21.3,
11.5, 60.3, 65.5, 69, 30.8, 70, 16.5, 15.8, 13.5, 22, 51.5, 55,
26.3, 55.3, 24.8, 55, 35, 13.3, 11.8, 16.5, 13.8, 15, 59.5, 26.5,
11, 7.8, 6.5, 5.3, 18.8, 29, 56.3, 71.3, 24, 12, 32.5, 96.1,
51.5, 55.8, 29.3, 16.8, 42.8, 85.9, 16.8, 15.8, 9, 12.3, 11.5,
11.8, 85.8, 13.5, 38.8, 9, 54.8, 8.8, 43.5, 12.3, 13, 80.5, 13.5,
75, 68.8, 22.8, 39.5, 43.3, 15.8, 14.8, 59.5, 91.4, 80.9, 79.6,
84.7, 15, 14.3, 16.8, 16.3, 13.3, 52, 14, 85.3, 14, 78.8, 10.8,
15.5, 67, 22.8, 15, 91.4, 19.3, 12.5, 16, 10.3, 57.8, 59, 12,
16.5, 58.5, 25.3, 53.3, 66, 9, 9.8, 9.3, 13.8, 13.5, 11, 50.8,
8, 6.8, 58.8, 54.3, 68, 80, 80, 20.3, 17.3, 17, 7.8, 4, 40.8,
21.5, 5.5, 14.5, 6, 33, 7, 11.3, 13, 7, 7, 61.8, 15.3, 18.5,
88.5, 21.8, 16.3, 16.5, 22, 14.5, 54.8, 15.3, 74.5, 51.5, 13.8,
14.8, 14, 49.5, 74.5, 13.3, 16, 56, 25.5, 14.8, 45.3, 21.3, 15.8,
17.8, 15.5, 43.3, 17, 31.5, 7, 91.7, 22.3, 5.8, 6.3, 7.5, 13.8,
95.2, 11.5, 11.3, 89.2, 78, 76.8, 17, 38.8, 48, 16.5, 126.7,
62.3, 83.7, 49.5, 31.3, 45, 101.2, 14.5, 26.5, 65.5, 12.3, 19.5,
90.4, 20.3, 7.3, 44.5, 21.3, 5.5, 4.8, 21.8, 14.5, 73.5, 51.3,
25.5, 4.5, 24.3, 5, 5, 20.5, 7.3, 5.3, 30, 5, 17, 15, 12.5, 74.5,
11, 66.5, 81.8, 14.5, 18.5, 17.5, 12.8, 13.3, 11.8, 15.3, 89.4,
22.3, 12.3, 80.2, 20.5, 12, 18.3, 17.3, 32.5, 69.3, 13.8, 20.3,
63, 37.8, 15, 49.5, 8.3, 13.5, 20.8, 14.8, 18.5, 35, 7, 9.8,
63.8, 21.5, 34, 13.5, 51.3, 30.8, 4.8, 4, 9.3, 33.5, 51.3, 83.8,
25, 13.8, 54.3)), .Names = c("plot", "d"), class = "data.frame", row.names =
c("1",
"2", "3", "4", "5", "13", "17", "18", "19", "20", "24", "25",
"21", "22", "26", "27", "32", "36", "38", "41", "42", "43", "44",
"47", "49", "50", "51", "53", "57", "60", "61", "62", "68", "69",
"70", "71", "75", "82", "83", "84", "86", "87", "89", "90", "91",
"94", "95", "96", "98", "101", "103", "107", "108", "109", "115",
"114", "124", "122", "125", "126", "127", "128", "138", "139",
"143", "144", "146", "147", "148", "149", "150", "153", "289",
"165", "170", "171", "172", "173", "176", "180", "181", "182",
"183", "187", "188", "190", "198", "199", "204", "205", "207",
"213", "209", "217", "218", "219", "225", "226", "227", "228",
"231", "236", "237", "238", "242", "249", "250", "256", "257",
"258", "264", "266", "273", "274", "277", "278", "279", "282",
"283", "284", "285", "286", "288", "293", "297", "299", "301",
"300", "302", "306", "309", "310", "311", "312", "313", "316",
"314", "315", "317", "318", "323", "326", "324", "325", "333",
"334", "335", "336", "343", "344", "350", "351", "352", "355",
"359", "363", "364", "370", "374", "377", "378", "379", "381",
"392", "393", "390", "391", "396", "397", "398", "399", "400",
"404", "405", "406", "407", "408", "409", "410", "422", "423",
"424", "425", "426", "427", "428", "429", "430", "431", "435",
"436", "6", "7", "8", "9", "10", "11", "12", "14", "15", "16",
"23", "28", "29", "30", "31", "33", "34", "35", "37", "39", "40",
"46", "45", "48", "52", "59", "54", "55", "56", "58", "63", "64",
"65", "66", "67", "72", "73", "74", "76", "77", "78", "79", "80",
"81", "85", "88", "92", "93", "97", "99", "100", "102", "104",
"105", "106", "110", "111", "112", "113", "116", "117", "118",
"119", "120", "121", "123", "130", "131", "129", "132", "133",
"134", "135", "136", "137", "140", "141", "142", "145", "151",
"152", "154", "155", "156", "157", "158", "159", "160", "161",
"162", "163", "164", "166", "167", "168", "169", "174", "175",
"177", "178", "179", "184", "185", "186", "191", "189", "192",
"193", "194", "195", "196", "197", "200", "201", "202", "203",
"206", "208", "210", "211", "212", "214", "215", "216", "220",
"221", "222", "223", "224", "229", "230", "232", "233", "234",
"235", "239", "240", "241", "243", "244", "245", "246", "247",
"248", "251", "252", "253", "254", "255", "259", "260", "261",
"262", "263", "265", "267", "268", "269", "270", "271", "272",
"275", "276", "280", "281", "287", "290", "291", "292", "294",
"296", "295", "298", "303", "305", "304", "307", "308", "319",
"320", "321", "322", "327", "328", "329", "330", "331", "332",
"337", "338", "339", "340", "341", "342", "345", "346", "347",
"348", "349", "353", "354", "356", "357", "358", "360", "361",
"362", "365", "366", "367", "368", "369", "371", "372", "373",
"375", "376", "380", "382", "383", "384", "385", "386", "387",
"388", "389", "394", "395", "401", "402", "403", "411", "412",
"413", "414", "415", "416", "417", "419", "418", "420", "421",
"432", "433", "434"))
It is easy to get the histograms you want using by():
> pdf("example.pdf")
> by(example$d, example$plot, hist)
> dev.off()
This produces a pdf file called example.pdf with 149 pages (one for each
histogram). The file is located in your default directory.
----------------------------------------------
David L Carlson
Associate Professor of Anthropology
Texas A&M University
College Station, TX 77843-4352
> -----Original Message-----
> From: r-help-bounces at r-project.org [mailto:r-help-bounces at r-
> project.org] On Behalf Of catalin roibu
> Sent: Thursday, January 03, 2013 1:50 AM
> To: r-help at r-project.org
> Subject: [R] Sas by function in R
>
> Hello,
> It's an alternative to use SAS by function in R?
> I want to plot d histograms by plot.from example bellow:
> Thank you!
> plot d
> 1 1 16.3
> 2 1 25.0
> 3 1 57.8
> 4 1 17.0
> 5 2 10.8
> 13 2 96.4
> 17 3 76.0
> 18 3 32.0
> 19 3 11.0
> 20 3 11.0
> 24 3 106.0
> 25 3 12.5
> 21 4 19.3
> 22 4 12.0
> 26 4 15.0
> 27 5 99.3
> 32 7 11.0
> 36 8 18.5
> 38 8 77.0
> 41 10 9.8
> 42 10 101.5
> 43 10 18.0
> 44 10 6.0
> 47 10 12.3
> 49 11 80.8
> 50 11 14.0
> 51 12 12.5
> 53 12 10.8
> 57 13 14.5
> 60 13 13.3
> 61 13 12.5
> 62 13 124.1
> 68 14 92.3
> 69 14 20.3
> 70 14 9.5
> 71 15 9.0
> 75 16 104.1
> 82 16 15.0
> 83 17 62.8
> 84 17 16.3
> 86 17 95.2
> 87 17 17.8
> 89 18 14.0
> 90 18 107.3
> 91 18 22.5
> 94 19 22.8
> 95 19 16.5
> 96 19 15.3
> 98 21 47.8
> 101 22 110.1
> 103 22 14.3
> 107 23 73.3
> 108 24 14.5
> 109 24 93.9
> 115 24 68.5
> 114 25 96.4
> 124 26 44.5
> 122 27 18.0
> 125 27 7.0
> 126 27 14.0
> 127 28 13.8
> 128 28 80.0
> 138 29 43.5
> 139 30 69.8
> 143 30 12.8
> 144 31 102.2
> 146 31 15.8
> 147 31 63.0
> 148 31 63.5
> 149 32 31.3
> 150 32 6.3
> 153 33 82.1
> 289 34 112.7
> 165 35 80.0
> 170 36 5.0
> 171 36 80.0
> 172 37 11.5
> 173 37 95.5
> 176 37 31.3
> 180 38 7.0
> 181 38 54.0
> 182 39 6.5
> 183 39 20.5
> 187 39 63.0
> 188 39 60.0
> 190 40 9.0
> 198 41 44.0
> 199 41 64.0
> 204 42 70.0
> 205 43 37.0
> 207 43 57.3
> 213 43 21.3
> 209 44 10.5
> 217 44 43.5
> 218 44 33.8
> 219 45 21.8
> 225 45 19.3
> 226 45 13.5
> 227 45 77.0
> 228 45 24.0
> 231 46 16.3
> 236 46 31.5
> 237 46 8.3
> 238 47 70.3
> 242 47 17.0
> 249 49 16.5
> 250 49 7.8
> 256 51 62.8
> 257 51 64.8
> 258 52 94.9
> 264 53 94.5
> 266 53 36.8
> 273 56 19.5
> 274 56 12.8
> 277 56 62.5
> 278 56 40.0
> 279 57 17.0
> 282 57 70.5
> 283 58 10.5
> 284 58 48.0
> 285 58 59.8
> 286 59 10.3
> 288 59 16.8
> 293 59 19.8
> 297 60 108.5
> 299 61 106.0
> 301 61 22.8
> 300 62 15.3
> 302 62 17.3
> 306 63 57.5
> 309 63 15.8
> 310 63 104.9
> 311 63 7.8
> 312 63 5.0
> 313 63 5.0
> 316 63 4.8
> 314 64 6.3
> 315 64 5.0
> 317 64 9.5
> 318 64 21.3
> 323 64 5.0
> 326 64 54.0
> 324 65 17.5
> 325 65 68.0
> 333 65 6.8
> 334 65 13.8
> 335 66 80.5
> 336 66 15.5
> 343 66 28.8
> 344 66 12.5
> 350 67 27.0
> 351 67 21.8
> 352 67 25.3
> 355 68 85.3
> 359 68 18.0
> 363 70 10.8
> 364 70 14.8
> 370 71 13.0
> 374 72 90.1
> 377 72 33.8
> 378 72 7.8
> 379 72 14.3
> 381 73 72.0
> 392 74 66.0
> 393 74 24.8
> 390 75 14.8
> 391 75 93.3
> 396 75 12.5
> 397 75 15.0
> 398 76 22.0
> 399 76 90.4
> 400 76 53.0
> 404 77 21.3
> 405 78 11.5
> 406 78 60.3
> 407 78 65.5
> 408 78 69.0
> 409 78 30.8
> 410 78 70.0
> 422 79 16.5
> 423 79 15.8
> 424 79 13.5
> 425 79 22.0
> 426 79 51.5
> 427 80 55.0
> 428 80 26.3
> 429 80 55.3
> 430 80 24.8
> 431 80 55.0
> 435 80 35.0
> 436 80 13.3
> 6 81 11.8
> 7 81 16.5
> 8 81 13.8
> 9 81 15.0
> 10 81 59.5
> 11 82 26.5
> 12 82 11.0
> 14 82 7.8
> 15 82 6.5
> 16 83 5.3
> 23 84 18.8
> 28 85 29.0
> 29 85 56.3
> 30 86 71.3
> 31 86 24.0
> 33 87 12.0
> 34 87 32.5
> 35 88 96.1
> 37 88 51.5
> 39 89 55.8
> 40 89 29.3
> 46 89 16.8
> 45 90 42.8
> 48 91 85.9
> 52 91 16.8
> 59 91 15.8
> 54 92 9.0
> 55 92 12.3
> 56 92 11.5
> 58 92 11.8
> 63 92 85.8
> 64 93 13.5
> 65 93 38.8
> 66 93 9.0
> 67 94 54.8
> 72 95 8.8
> 73 95 43.5
> 74 95 12.3
> 76 95 13.0
> 77 95 80.5
> 78 96 13.5
> 79 96 75.0
> 80 97 68.8
> 81 97 22.8
> 85 97 39.5
> 88 97 43.3
> 92 98 15.8
> 93 99 14.8
> 97 100 59.5
> 99 101 91.4
> 100 102 80.9
> 102 102 79.6
> 104 102 84.7
> 105 102 15.0
> 106 102 14.3
> 110 104 16.8
> 111 104 16.3
> 112 104 13.3
> 113 104 52.0
> 116 105 14.0
> 117 106 85.3
> 118 106 14.0
> 119 107 78.8
> 120 107 10.8
> 121 107 15.5
> 123 107 67.0
> 130 107 22.8
> 131 107 15.0
> 129 108 91.4
> 132 108 19.3
> 133 109 12.5
> 134 109 16.0
> 135 109 10.3
> 136 109 57.8
> 137 109 59.0
> 140 109 12.0
> 141 110 16.5
> 142 110 58.5
> 145 111 25.3
> 151 112 53.3
> 152 112 66.0
> 154 113 9.0
> 155 113 9.8
> 156 113 9.3
> 157 113 13.8
> 158 113 13.5
> 159 114 11.0
> 160 114 50.8
> 161 114 8.0
> 162 114 6.8
> 163 114 58.8
> 164 114 54.3
> 166 115 68.0
> 167 115 80.0
> 168 116 80.0
> 169 116 20.3
> 174 117 17.3
> 175 117 17.0
> 177 117 7.8
> 178 117 4.0
> 179 117 40.8
> 184 119 21.5
> 185 119 5.5
> 186 119 14.5
> 191 119 6.0
> 189 120 33.0
> 192 120 7.0
> 193 120 11.3
> 194 120 13.0
> 195 120 7.0
> 196 120 7.0
> 197 121 61.8
> 200 122 15.3
> 201 122 18.5
> 202 122 88.5
> 203 122 21.8
> 206 123 16.3
> 208 123 16.5
> 210 123 22.0
> 211 123 14.5
> 212 124 54.8
> 214 124 15.3
> 215 124 74.5
> 216 124 51.5
> 220 125 13.8
> 221 125 14.8
> 222 125 14.0
> 223 125 49.5
> 224 125 74.5
> 229 125 13.3
> 230 125 16.0
> 232 126 56.0
> 233 126 25.5
> 234 126 14.8
> 235 126 45.3
> 239 127 21.3
> 240 127 15.8
> 241 127 17.8
> 243 127 15.5
> 244 128 43.3
> 245 128 17.0
> 246 129 31.5
> 247 129 7.0
> 248 129 91.7
> 251 129 22.3
> 252 129 5.8
> 253 130 6.3
> 254 130 7.5
> 255 130 13.8
> 259 132 95.2
> 260 133 11.5
> 261 133 11.3
> 262 133 89.2
> 263 133 78.0
> 265 134 76.8
> 267 134 17.0
> 268 134 38.8
> 269 134 48.0
> 270 135 16.5
> 271 135 126.7
> 272 135 62.3
> 275 136 83.7
> 276 136 49.5
> 280 137 31.3
> 281 137 45.0
> 287 138 101.2
> 290 139 14.5
> 291 139 26.5
> 292 139 65.5
> 294 139 12.3
> 296 139 19.5
> 295 140 90.4
> 298 141 20.3
> 303 141 7.3
> 305 141 44.5
> 304 142 21.3
> 307 142 5.5
> 308 142 4.8
> 319 143 21.8
> 320 143 14.5
> 321 143 73.5
> 322 143 51.3
> 327 144 25.5
> 328 144 4.5
> 329 144 24.3
> 330 145 5.0
> 331 145 5.0
> 332 145 20.5
> 337 145 7.3
> 338 145 5.3
> 339 145 30.0
> 340 145 5.0
> 341 145 17.0
> 342 146 15.0
> 345 146 12.5
> 346 146 74.5
> 347 146 11.0
> 348 146 66.5
> 349 147 81.8
> 353 148 14.5
> 354 148 18.5
> 356 148 17.5
> 357 148 12.8
> 358 149 13.3
> 360 149 11.8
> 361 149 15.3
> 362 149 89.4
> 365 150 22.3
> 366 150 12.3
> 367 150 80.2
> 368 150 20.5
> 369 151 12.0
> 371 151 18.3
> 372 151 17.3
> 373 151 32.5
> 375 152 69.3
> 376 152 13.8
> 380 152 20.3
> 382 153 63.0
> 383 153 37.8
> 384 153 15.0
> 385 154 49.5
> 386 154 8.3
> 387 154 13.5
> 388 154 20.8
> 389 155 14.8
> 394 155 18.5
> 395 155 35.0
> 401 156 7.0
> 402 156 9.8
> 403 157 63.8
> 411 157 21.5
> 412 157 34.0
> 413 158 13.5
> 414 158 51.3
> 415 158 30.8
> 416 158 4.8
> 417 158 4.0
> 419 158 9.3
> 418 159 33.5
> 420 159 51.3
> 421 159 83.8
> 432 159 25.0
> 433 160 13.8
> 434 160 54.3
>
> --
> ---
> Catalin-Constantin ROIBU
> Forestry engineer, PhD
> Forestry Faculty of Suceava
> Str. Universitatii no. 13, Suceava, 720229, Romania
> office phone +4 0230 52 29 78, ext. 531
> mobile phone +4 0745 53 18 01
> +4 0766 71 76 58
> FAX: +4 0230 52 16 64
> silvic.usv.ro
>
> [[alternative HTML version deleted]]
>
> ______________________________________________
> R-help at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide http://www.R-project.org/posting-
> guide.html
> and provide commented, minimal, self-contained, reproducible code.
More information about the R-help
mailing list