On Fri, May 20, 2011 at 4:04 PM, Chris Mcowen <chrismco...@gmail.com> wrote: > Thanks Peter, > > This is what i am after but i would like this ( the counts) for each ECO_NAME > not the total for the entire data set. > > A sample of my data is below > > data2 <- structure(list(ECO_NAME = structure(c(1L, 1L, 2L, 2L, 2L, 2L, > 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 3L, > 3L, 3L, 3L, 3L, 3L), .Label = c("Biak-Numfoor rain forests", > "Central Range montane rain forests", "Huon Peninsula montane rain forests" > ), class = "factor"), Order = structure(c(1L, 2L, 1L, 2L, 2L, > 2L, 2L, 2L, 2L, 2L, 2L, 3L, 4L, 4L, 4L, 4L, 4L, 4L, 5L, 5L, 5L, > 2L, 2L, 2L, 2L, 4L, 4L), .Label = c("Alismatales", "Asparagales", > "Commelinales", "Poales", "Zingiberales"), class = "factor"), > Family = structure(c(1L, 7L, 1L, 6L, 7L, 7L, 7L, 7L, 7L, > 7L, 7L, 3L, 2L, 4L, 4L, 5L, 8L, 8L, 9L, 9L, 9L, 7L, 7L, 7L, > 7L, 2L, 8L), .Label = c("Araceae", "Centrolepidaceae", "Commelinaceae", > "Cyperaceae", "Eriocaulaceae", "Lomandraceae", "Orchidaceae", > "Poaceae", "Zingiberaceae"), class = "factor"), Genus = structure(c(13L, > 3L, 13L, 6L, 19L, 9L, 20L, 4L, 4L, 9L, 9L, 15L, 5L, 14L, > 10L, 11L, 17L, 16L, 1L, 8L, 2L, 18L, 7L, 19L, 12L, 5L, 16L > ), .Label = c("Alpinia", "Amomum", "Bromheadia", "Bulbophyllum", > "Centrolepis", "Cordyline", "Corybas", "Curcuma", "Dendrobium", > "Eleocharis", "Eriocaulon", "Glomera", "Homalomena", "Machaerina", > "Murdannia", "Poa", "Schizostachyum", "Taeniophyllum", "Thelymitra", > "Vanda"), class = "factor")), .Names = c("ECO_NAME", "Order", > "Family", "Genus"), class = "data.frame", row.names = c(NA, -27L > ))
Chris, How about this? (Note that summed incidences are the same for all taxonomic levels.) > data.frame(unique.Genus=rowSums(xtabs(~ ECO_NAME + Genus, data2) > 0), + unique.Family=rowSums(xtabs(~ ECO_NAME + Family, data2) > 0), + unique.Order=rowSums(xtabs(~ ECO_NAME + Order, data2) > 0), + incidences.Genus=rowSums(xtabs(~ ECO_NAME + Genus, data2)), + incidences.Family=rowSums(xtabs(~ ECO_NAME + Family, data2)), + incidences.Order=rowSums(xtabs(~ ECO_NAME + Order, data2))) unique.Genus unique.Family unique.Order Biak-Numfoor rain forests 2 2 2 Central Range montane rain forests 16 9 5 Huon Peninsula montane rain forests 6 3 2 incidences.Genus incidences.Family Biak-Numfoor rain forests 2 2 Central Range montane rain forests 19 19 Huon Peninsula montane rain forests 6 6 incidences.Order Biak-Numfoor rain forests 2 Central Range montane rain forests 19 Huon Peninsula montane rain forests 6 Peter > > > On 20 May 2011, at 22:32, Peter Solymos wrote: > > Chris, > > Something like this should do the job: > > 1. > rowSums(xtabs(~ ECO_NAME + Genus, x) > 0) > 2. > rowSums(xtabs(~ ECO_NAME + Genus, x)) > > Cheers, > > Peter > > > > On Fri, May 20, 2011 at 3:19 PM, Chris Mcowen <chrismco...@gmail.com> wrote: >> Dear List, >> >> I am looking to calculate two things from my data frame and was after some >> advice. For the example below i want to know. >> >> 1. How many unique Orders/Families and Genera there are per eco-name >> >> 2. How many incidences are there for each Order/Family and Genus there are >> per eco-region >> >> I have 650 econame. >> >> I.e for Biak-Numfoor rain forests there are 2 orders, 2 families and two >> genera. Also, Alismatales are represented once, Asparagales once etc etc. >> >> Thanks for any advice. >> >> Chris >> >> >> ECO_NAME Order >> Family Genus >> Biak-Numfoor rain forests Alismatales >> Araceae Homalomena >> Biak-Numfoor rain forests Asparagales >> Orchidaceae Bromheadia >> Central Range montane rain forests Alismatales Araceae >> Homalomena >> Central Range montane rain forests Asparagales Lomandraceae >> Cordyline >> Central Range montane rain forests Asparagales Orchidaceae >> Thelymitra >> Central Range montane rain forests Asparagales Orchidaceae >> Dendrobium >> Central Range montane rain forests Asparagales Orchidaceae >> Vanda >> Central Range montane rain forests Asparagales Orchidaceae >> Bulbophyllum >> Central Range montane rain forests Asparagales Orchidaceae >> Bulbophyllum >> Central Range montane rain forests Asparagales Orchidaceae >> Dendrobium >> Central Range montane rain forests Asparagales Orchidaceae >> Dendrobium >> Central Range montane rain forests Commelinales Commelinaceae >> Murdannia >> Central Range montane rain forests Poales >> Centrolepidaceae Centrolepis >> Central Range montane rain forests Poales Cyperaceae >> Machaerina >> Central Range montane rain forests Poales Cyperaceae >> Eleocharis >> Central Range montane rain forests Poales >> Eriocaulaceae Eriocaulon >> Central Range montane rain forests Poales Poaceae >> Schizostachyum >> Central Range montane rain forests Poales Poaceae >> Poa >> Central Range montane rain forests Zingiberales >> Zingiberaceae Alpinia >> Central Range montane rain forests Zingiberales >> Zingiberaceae Curcuma >> Central Range montane rain forests Zingiberales >> Zingiberaceae Amomum >> Huon Peninsula montane rain forests Asparagales Orchidaceae >> Taeniophyllum >> Huon Peninsula montane rain forests Asparagales Orchidaceae >> Corybas >> Huon Peninsula montane rain forests Asparagales Orchidaceae >> Thelymitra >> Huon Peninsula montane rain forests Asparagales Orchidaceae >> Glomera >> Huon Peninsula montane rain forests Poales >> Centrolepidaceae Centrolepis >> Huon Peninsula montane rain forests Poales Poaceae >> Poa >> >> _______________________________________________ >> R-sig-ecology mailing list >> R-sig-ecology@r-project.org >> https://stat.ethz.ch/mailman/listinfo/r-sig-ecology >> >> > > _______________________________________________ R-sig-ecology mailing list R-sig-ecology@r-project.org https://stat.ethz.ch/mailman/listinfo/r-sig-ecology