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 <jeanwaij...@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]]
______________________________________________
R-help@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.
