[R] Creating vegetation distance groups from one column

Jhope jeanwaijang at gmail.com
Thu Oct 4 21:25:09 CEST 2012


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



--
View this message in context: http://r.789695.n4.nabble.com/Creating-vegetation-distance-groups-from-one-column-tp4644970p4645066.html
Sent from the R help mailing list archive at Nabble.com.




More information about the R-help mailing list