Try this

data.to.analyze$VegIndex <- cut(data.to.analyze$Veg, 
        breaks=seq(0, 70, 5), include.lowest=TRUE)
plot(data.to.analyze$VegIndex)



Jhope <jeanwaijang@gmail.com> wrote on 10/04/2012 02:25:09 PM:
> 
> Hi,
> 
> Allow me to recap my question. In plyr I am trying to group distances of
> nests to the vegetation into classes that are defined by (0-5m, 6-10m,
> 11-15m, 16-20m, 21-25m, 26-30m, 31-35m). Each row is a nest and all the
> distances to the vegetation is in a column. 
> 
> As well I would like to create classes from the HTL starting at (-5-0m,
> 1-5m, 6-10m, 11-15m, 16-20m, 21-25m, 26m-30m). 
> 
> I have included a sample of the data through dput() - see below. If that 
is
> too tiresome, I've also included head(data.to.analyze) below that. 
> 
> Best, Jean
> 
> 
> > dput(head(data.to.analyze, 30))
> structure(list(QuadratEvent = c("8.1A", "10.1A", "14.1A", "19.1A", 
> "23.1A", "24.1A", "25.1A", "29.1A", "30.1A", "31.1A", "32.1A", 
> "33.1A", "38.1A", "39.1A", "41.1A", "47.1A", "48.1A", "50.1A", 
> "52.1A", "53.1A", "58.1A", "63.1A", "68.1A", "73.1A", "75.1A", 
> "79.1A", "29BA", "90.1B", "90.2B", "82.1B"), ExDate = c("Nov. 26, 2010", 

> "Nov. 26, 2010", "Nov. 18, 2010", "Nov. 22, 2010", "Nov. 20, 2010", 
> "Nov. 28, 2010", "Nov. 28, 2010", "Nov. 27, 2010", "Nov. 20, 2010", 
> "Nov. 20, 2010", "Nov. 26, 2010", "Nov. 26, 2010", "Nov. 20, 2010", 
> "Nov. 26, 2010", "Nov. 27, 2010", "Nov. 27, 2010", "Nov. 26, 2010", 
> "Nov. 27, 2010", "Nov. 27, 2010", "Nov. 27, 2010", "Nov. 27, 2010", 
> "Nov. 26, 2010", "Nov. 26, 2010", "Nov. 28, 2010", "Nov. 28, 2010", 
> "Nov. 26, 2010", "Nov. 28, 2010", "Jan. 1, 2011", "Jan. 1, 2011", 
> "Dec. 23, 2010"), Sector = c("14", "13", "20", "79", "34", "34", 
> "34", "36", "67", "66", "67", "67", "67", "66", "70", "70", "68", 
> "57", "55", "56", "40", "27", "26", "27", "28", "19", "36", "70", 
> "70", "47"), Quadrat = c("8", "10", "14", "19", "23", "24", "25", 
> "29", "30", "31", "32", "33", "38", "39", "41", "47", "48", "50", 
> "52", "53", "58", "63", "68", "73", "75", "79", "29B", "90", 
> "90", "82"), Aeventexhumed = c("A", "A", "A", "A", "A", "A", 
> "A", "A", "A", "A", "A", "A", "A", "A", "A", "A", "A", "A", "A", 
> "A", "A", "A", "A", "A", "A", "A", "A", "B", "B", "B"), ST0 = c(6L, 
> 4L, 6L, 0L, 135L, 660L, 321L, 608L, 50L, 5L, 88L, 48L, 7L, 60L, 
> 0L, 69L, 195L, 68L, 162L, 47L, 595L, 144L, 111L, 478L, 636L, 
> 136L, 515L, 213L, 0L, 175L), ST1 = c(0L, 0L, 0L, 0L, 23L, 0L, 
> 0L, 0L, 4L, 0L, 0L, 0L, 0L, 0L, 1L, 2L, 0L, 39L, 2L, 1L, 4L, 
> 9L, 1L, 5L, 25L, 4L, 139L, 0L, 0L, 0L), ST2 = c(1L, 0L, 2L, 0L, 
> 1L, 2L, 2L, 0L, 153L, 0L, 0L, 0L, 117L, 93L, 3L, 2L, 0L, 0L, 
> 6L, 2L, 0L, 17L, 1L, 11L, 67L, 0L, 31L, 68L, 0L, 0L), ST3 = c(98L, 
> 109L, 26L, 0L, 6L, 0L, 0L, 0L, 1L, 3L, 11L, 88L, 3L, 342L, 9L, 
> 3L, 4L, 72L, 7L, 81L, 0L, 33L, 103L, 0L, 0L, 218L, 162L, 0L, 
> 74L, 0L), ST4 = c(1L, 1L, 11L, 0L, 0L, 0L, 0L, 0L, 27L, 1L, 0L, 
> 0L, 7L, 0L, 8L, 6L, 0L, 0L, 5L, 4L, 0L, 38L, 150L, 0L, 0L, 19L, 
> 5L, 0L, 0L, 0L), Shells = c(72L, 74L, 106L, 102L, 125L, 4L, 15L, 
> 0L, 31L, 111L, 0L, 0L, 86L, 48L, 21L, 35L, 0L, 0L, 12L, 13L, 
> 0L, 23L, 20L, 3L, 4L, 115L, 5L, 2L, 0L, 0L), Hatchlings = c(3L, 
> 4L, 1L, 1L, 48L, 0L, 0L, 0L, 45L, 99L, 0L, 0L, 93L, 9L, 0L, 0L, 
> 0L, 6L, 2L, 6L, 0L, 9L, 13L, 0L, 0L, 9L, 0L, 0L, 0L, 0L), MaxHatch = 
c(72L, 
> 74L, 106L, 102L, 125L, 4L, 15L, 0L, 45L, 111L, 0L, 0L, 93L, 48L, 
> 21L, 35L, 0L, 6L, 12L, 13L, 0L, 23L, 20L, 3L, 4L, 118L, 5L, 2L, 
> 0L, 0L), DeadHatch = c(2L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 4L, 0L, 
> 0L, 0L, 1L, 5L, 0L, 0L, 0L, 0L, 0L, 2L, 0L, 0L, 0L, 0L, 0L, 6L, 
> 0L, 0L, 0L, 0L), Oldeggs = c(0L, 21L, 165L, 0L, 0L, 36L, 0L, 
> 0L, 96L, 45L, 5L, 142L, 10L, 81L, 103L, 92L, 0L, 0L, 0L, 0L, 
> 108L, 0L, 0L, 0L, 93L, 0L, 0L, 0L, 0L, 0L), TotalEggs = c(178L, 
> 188L, 151L, 102L, 290L, 666L, 338L, 608L, 280L, 120L, 99L, 136L, 
> 227L, 543L, 42L, 117L, 199L, 185L, 194L, 148L, 599L, 264L, 386L, 
> 497L, 732L, 495L, 857L, 283L, 74L, 175L), QuadratID = structure(c(107L, 
> 4L, 22L, 38L, 43L, 44L, 46L, 54L, 57L, 58L, 60L, 61L, 67L, 69L, 
> 72L, 79L, 80L, 81L, 83L, 85L, 91L, 94L, 97L, 103L, 105L, 106L, 
> 56L, 109L, NA, 108L), .Label = c("", "1.1", "1.2", "10.1", "10.2", 
> "100.1", "100.2", "108.1", "108.2", "110.1", "110.2", "114.1", 
> "114.2", "121.1", "121.2", "127.1", "127.2", "128.1", "128.2", 
> "135.1", "135.2", "14.1", "15.1", "15.2", "151.1", "151.2", "152.1", 
> "152.2", "153.1", "153.2", "154.1", "154.2", "155.1", "155.2", 
> "16.1", "16.2", "18.1", "19.1", "21.1", "21.2", "22.1", "22.2", 
> "23.1", "24.1", "24.2", "25.1", "25.2", "26.1", "26.2", "27.1", 
> "27.2", "28.1", "28.2", "29.1", "29.2", "29B", "30.1", "31.1", 
> "31.2", "32.1", "33.1", "33.2", "35.1", "35.2", "37.1", "37.2", 
> "38.1", "38.2", "39.1", "40.1", "40.2", "41.1", "42.1", "42.2", 
> "43.1", "43.2", "46.1", "46.2", "47.1", "48.1", "50.1", "50.2", 
> "52.1", "52.2", "53.1", "53.2", "54.1", "54.2", "56.1", "56.2", 
> "58.1", "62.1", "62.2", "63.1", "65.1", "65.2", "68.1", "68.2", 
> "68.3", "68.4", "69.1", "70.1", "73.1", "74.1", "75.1", "79.1", 
> "8.1", "82.1", "90.1", "94.1", "94.2"), class = "factor"), Sector = 
c(14L, 
> 13L, 20L, 79L, 34L, 34L, 34L, 36L, 67L, 66L, 67L, 67L, 67L, 66L, 
> 70L, 70L, 68L, 57L, 55L, 56L, 40L, 27L, 26L, 27L, 28L, 19L, 36L, 
> 70L, NA, 47L), Aeventexhumed = structure(c(2L, 2L, 2L, 2L, 2L, 
> 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 
> 2L, 2L, 2L, 2L, 2L, 2L, 3L, NA, 3L), .Label = c("", "A", "B", 
> "C"), class = "factor"), Veg = c(9.75, 16.34, 3.15, 5.7, 12.2, 
> 4.82, 3.3, 1.85, 17.75, 27.13, 0.84, 8.45, 16.42, 20.4, 22.4, 
> 25, 25.4, 15.1, 13.7, 17.1, NA, 17.3, 24.3, 9.35, 5.4, 15.5, 
> 1.85, 11, NA, 3), Berm = c(NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> NA, NA, 10, NA, 11), HTL = c(4.66, -2.47, 16, NA, 2.35, 22.85, 
> 24, 15.7, 9.55, 4.2, 9.65, 19.8, 13.85, 8.7, 7.4, 2.9, 2.9, 4.26, 
> 1.2, 1.05, NA, 3.2, 3.8, 2.55, 6.55, 2, 15.7, 20, NA, 20), Estuary = 
c(NA, 
> NA, NA, 11.8, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 
> NA, NA, NA, NA, 3.5, NA, NA, NA, 8.05, NA, NA, NA, NA, NA), Rayos =
> c("Rayos.1", 
> "Rayos.1", "Rayos.2", "Rayos.2", "Rayos.3", "Rayos.3", "Rayos.3", 
> "Rayos.3", "Rayos.5", "Rayos.5", "Rayos.5", "Rayos.5", "Rayos.5", 
> "Rayos.5", "Rayos.5", "Rayos.5", "Rayos.5", "Rayos.4", "Rayos.4", 
> "Rayos.4", "Rayos.3", "Rayos.2", "Rayos.2", "Rayos.2", "Rayos.2", 
> "Rayos.2", "Rayos.3", "Rayos.5", "Rayos.NA", "Rayos.4"), NotHatched =
> c(106L, 
> 114L, 45L, 0L, 165L, 662L, 323L, 608L, 249L, 9L, 99L, 136L, 141L, 
> 495L, 21L, 82L, 199L, 185L, 182L, 135L, 599L, 241L, 366L, 494L, 
> 728L, 380L, 852L, 281L, 74L, 175L), DevelopIndex = c(3.86516853932584, 
> 4.92021276595745, 5.62251655629139, 1, 3.40689655172414, 
5.98198198198198, 
> 1.13313609467456, 2, 3.53571428571429, 1.15, 4.66666666666667, 
> 4.05882352941176, 2.26872246696035, 4.15469613259669, 5.21428571428571, 
> 2.97435897435897, 4.93969849246231, 4.1027027027027, 1.32474226804124, 
> 3.92567567567568, 3.02003338898164, 3.26136363636364, 4.03626943005181, 
> 5.95171026156942, 1.11885245901639, 4.06060606060606, 4.07117852975496, 
> 3.97879858657244, 2, 6)), .Names = c("QuadratEvent", "ExDate", 
> "Sector", "Quadrat", "Aeventexhumed", "ST0", "ST1", "ST2", "ST3", 
> "ST4", "Shells", "Hatchlings", "MaxHatch", "DeadHatch", "Oldeggs", 
> "TotalEggs", "QuadratID", "Sector", "Aeventexhumed", "Veg", "Berm", 
> "HTL", "Estuary", "Rayos", "NotHatched", "DevelopIndex"), row.names = 
c(NA, 
> 30L), class = "data.frame")
> 
> 
> ----------------------------------------
> 
> >head(data.to.analyze)
>   QuadratEvent        ExDate Sector Quadrat Aeventexhumed ST0 ST1 ST2 
ST3
> ST4 Shells
> 1         8.1A Nov. 26, 2010     14       8             A   6   0   1 98 
 
> 1     72
> 2        10.1A Nov. 26, 2010     13      10             A   4   0   0 
109 
> 1     74
> 3        14.1A Nov. 18, 2010     20      14             A   6   0   2 26 

> 11    106
> 4        19.1A Nov. 22, 2010     79      19             A   0   0   0 0 
> 0    102
> 5        23.1A Nov. 20, 2010     34      23             A 135  23   1 6 
> 0    125
> 6        24.1A Nov. 28, 2010     34      24             A 660   0   2 0 
> 0      4
>   Hatchlings MaxHatch DeadHatch Oldeggs TotalEggs QuadratID Sector
> Aeventexhumed   Veg Berm
> 1          3       72         2       0       178       8.1     14  
> A  9.75   NA
> 2          4       74         0      21       188      10.1     13  
> A 16.34   NA
> 3          1      106         1     165       151      14.1     20  
> A  3.15   NA
> 4          1      102         0       0       102      19.1     79  
> A  5.70   NA
> 5         48      125         0       0       290      23.1     34  
> A 12.20   NA
> 6          0        4         0      36       666      24.1     34  
> A  4.82   NA
>     HTL Estuary   Rayos NotHatched DevelopIndex
> 1  4.66      NA Rayos.1        106     3.865169
> 2 -2.47      NA Rayos.1        114     4.920213
> 3 16.00      NA Rayos.2         45     5.622517
> 4    NA    11.8 Rayos.2          0     1.000000
> 5  2.35      NA Rayos.3        165     3.406897
> 6 22.85      NA Rayos.3        662     5.981982

	[[alternative HTML version deleted]]


