Hi all, I am revisiting using reshape2 to aggregate critter (bats) occurrences by time blocks and have the final output for each row sum a percentage such that the row sums total 1 (100%). Reshape2 library is loaded on start. The old code I had stored in NotePad++ seems to no longer work as I have the following error: > input.cast <- cast(input.melt, Species ~ Date, fun.aggregate = sum) *Error: could not find function "cast"* The code saved in NotePad++ is: require(reshape2) input.melt <- melt(input) input.cast <- cast(input.melt, Species ~ Date, fun.aggregate = sum) input.percent <- input.cast[, -1] / rowSums(input.cast[, -1]) * 100 rownames(input.percent) <- input.cast[,1] The data is read in OK: >Escameca_temporal <- read.table("C:/=Bat data working/=Nica_new/Escameca_temporal.CSV",header=T,sep=",",quote="") #then renamed : >input <-Escameca_temporal The variable view shows Species (6 letter text code) as Factor;Location code as integer; Date as Factor; Time as Factor. Data looks like this (sub sample only) Species Location Date Time Cynmex 4937 5/3/2012 23:42 Cynmex 4937 5/3/2012 23:43 Cynmex 4937 5/3/2012 23:43 Eptfur 4937 5/2/2012 0:10 Eptfur 4937 5/2/2012 0:12 Eptfur 4937 5/2/2012 0:15 Eptfur 4937 5/2/2012 0:16 Eptfur 4937 5/2/2012 0:21 Eptfur 4937 5/2/2012 0:27 Eptfur 4937 5/2/2012 0:27 Eptfur 4937 5/2/2012 0:36 Eptfur 4937 5/2/2012 0:37 Eptfur 4937 5/2/2012 0:38 Eptfur 4937 5/2/2012 0:42 Eptfur 4937 5/2/2012 0:57 Eptfur 4937 5/2/2012 0:58 Eptfur 4937 5/2/2012 0:59 Eptfur 4937 5/2/2012 1:01 Eptfur 4937 5/2/2012 1:03 Eptfur 4937 5/2/2012 1:08 Eptfur 4937 5/4/2012 2:02 Eptfur 4937 5/4/2012 2:03 Eptfur 4937 5/4/2012 2:03 Eptfur 4937 5/4/2012 2:04 Eptfur 4937 5/4/2012 2:04 Eptfur 4937 5/4/2012 2:08 Eptfur 4924 1/8/2012 2:09 Eptfur 4937 5/4/2012 2:09 Eptfur 4924 1/8/2012 2:10 Eptfur 4924 1/8/2012 2:10 Eptfur 4924 1/8/2012 2:10 Eptfur 4924 1/8/2012 2:10 Eptfur 4924 1/8/2012 2:11 Eptfur 4924 1/8/2012 2:11 Eptfur 4924 1/8/2012 2:11 Eptfur 4924 1/8/2012 2:11 Eptfur 4924 1/8/2012 2:11 Eptfur 4924 1/8/2012 2:12 Eptfur 4924 1/8/2012 2:12 Eptfur 4924 1/8/2012 2:12 Eptfur 4924 1/8/2012 2:12 Eptfur 4924 1/8/2012 2:12 Eptfur 4924 1/8/2012 2:13 Eptfur 4924 1/8/2012 2:13 Eptfur 4937 5/4/2012 2:13 Eptfur 4924 1/8/2012 2:14 Eptfur 4924 1/8/2012 2:14 Eptfur 4924 1/8/2012 2:14 Eptfur 4924 1/8/2012 2:14 Eptfur 4937 5/4/2012 2:14 Eptfur 4924 1/8/2012 2:15 Eptfur 4924 1/8/2012 2:15 Eptfur 4924 1/8/2012 2:15 Eptfur 4924 1/8/2012 2:16 Eptfur 4924 1/8/2012 2:16 Eptfur 4924 1/8/2012 2:16 Eptfur 4924 1/8/2012 2:16 Eptfur 4924 1/8/2012 2:17 Any suggestions-help appreciated. Bruce [[alternative HTML version deleted]]
Hi all, I am revisiting using reshape2 to aggregate critter (bats) occurrences by time blocks and have the final output for each row sum a percentage such that the row sums total 1 (100%). Reshape2 library is loaded on start. The old code I had stored in NotePad++ seems to no longer work as I have the following error: > input.cast <- cast(input.melt, Species ~ Date, fun.aggregate = sum) *Error: could not find function "cast"* The code saved in NotePad++ is: require(reshape2) input.melt <- melt(input) input.cast <- cast(input.melt, Species ~ Date, fun.aggregate = sum) input.percent <- input.cast[, -1] / rowSums(input.cast[, -1]) * 100 rownames(input.percent) <- input.cast[,1] The data is read in OK: >Escameca_temporal <- read.table("C:/=Bat data working/=Nica_new/Escameca_temporal.CSV",header=T,sep=",",quote="") #then renamed : >input <-Escameca_temporal The variable view shows Species (6 letter text code) as Factor;Location code as integer; Date as Factor; Time as Factor. Data looks like this (sub sample only) Species Location Date Time Cynmex 4937 5/3/2012 23:42 Cynmex 4937 5/3/2012 23:43 Cynmex 4937 5/3/2012 23:43 Eptfur 4937 5/2/2012 0:10 Eptfur 4937 5/2/2012 0:12 Eptfur 4937 5/2/2012 0:15 Eptfur 4937 5/2/2012 0:16 Eptfur 4937 5/2/2012 0:21 Eptfur 4937 5/2/2012 0:27 Eptfur 4937 5/2/2012 0:27 Eptfur 4937 5/2/2012 0:36 Eptfur 4937 5/2/2012 0:37 Eptfur 4937 5/2/2012 0:38 Eptfur 4937 5/2/2012 0:42 Eptfur 4937 5/2/2012 0:57 Eptfur 4937 5/2/2012 0:58 Eptfur 4937 5/2/2012 0:59 Eptfur 4937 5/2/2012 1:01 Eptfur 4937 5/2/2012 1:03 Eptfur 4937 5/2/2012 1:08 Eptfur 4937 5/4/2012 2:02 Eptfur 4937 5/4/2012 2:03 Eptfur 4937 5/4/2012 2:03 Eptfur 4937 5/4/2012 2:04 Eptfur 4937 5/4/2012 2:04 Eptfur 4937 5/4/2012 2:08 Eptfur 4924 1/8/2012 2:09 Eptfur 4937 5/4/2012 2:09 Eptfur 4924 1/8/2012 2:10 Eptfur 4924 1/8/2012 2:10 Eptfur 4924 1/8/2012 2:10 Eptfur 4924 1/8/2012 2:10 Eptfur 4924 1/8/2012 2:11 Eptfur 4924 1/8/2012 2:11 Eptfur 4924 1/8/2012 2:11 Eptfur 4924 1/8/2012 2:11 Eptfur 4924 1/8/2012 2:11 Eptfur 4924 1/8/2012 2:12 Eptfur 4924 1/8/2012 2:12 Eptfur 4924 1/8/2012 2:12 Eptfur 4924 1/8/2012 2:12 Eptfur 4924 1/8/2012 2:12 Eptfur 4924 1/8/2012 2:13 Eptfur 4924 1/8/2012 2:13 Eptfur 4937 5/4/2012 2:13 Eptfur 4924 1/8/2012 2:14 Eptfur 4924 1/8/2012 2:14 Eptfur 4924 1/8/2012 2:14 Eptfur 4924 1/8/2012 2:14 Eptfur 4937 5/4/2012 2:14 Eptfur 4924 1/8/2012 2:15 Eptfur 4924 1/8/2012 2:15 Eptfur 4924 1/8/2012 2:15 Eptfur 4924 1/8/2012 2:16 Eptfur 4924 1/8/2012 2:16 Eptfur 4924 1/8/2012 2:16 Eptfur 4924 1/8/2012 2:16 Eptfur 4924 1/8/2012 2:17 Any suggestions-help appreciated. Bruce -- Bruce W. Miller, PhD. Neotropical bat risk assessments If we lose the bats, we may lose much of the tropical vegetation and the lungs of the planet Using acoustic sampling to map species distributions for >15 years. Providing Interactive identification keys to the vocal signatures of New World Bats For various project details see: https://sites.google.com/site/batsoundservices/ [[alternative HTML version deleted]]
library(reshape2) ?cast "Use acast or dcast depending on whether you want vector/matrix/array output or data frame output." Jean On Wed, Jun 5, 2013 at 9:35 AM, Bruce Miller <batsncats@gmail.com> wrote:> Hi all, > > I am revisiting using reshape2 to aggregate critter (bats) occurrences > by time blocks and have the final output for each row sum a percentage > such that the row sums total 1 (100%). > > Reshape2 library is loaded on start. > > The old code I had stored in NotePad++ seems to no longer work as I have > the following error: > > > input.cast <- cast(input.melt, Species ~ Date, fun.aggregate = sum) > *Error: could not find function "cast"* > > The code saved in NotePad++ is: > > require(reshape2) > input.melt <- melt(input) > input.cast <- cast(input.melt, Species ~ Date, fun.aggregate = sum) > input.percent <- input.cast[, -1] / rowSums(input.cast[, -1]) * 100 > rownames(input.percent) <- input.cast[,1] > > The data is read in OK: > >Escameca_temporal <- read.table("C:/=Bat data > working/=Nica_new/Escameca_temporal.CSV",header=T,sep=",",quote="") > #then renamed : >input <-Escameca_temporal > The variable view shows Species (6 letter text code) as Factor;Location > code as integer; Date as Factor; Time as Factor. > > Data looks like this (sub sample only) > > Species Location Date Time > Cynmex 4937 5/3/2012 23:42 > Cynmex 4937 5/3/2012 23:43 > Cynmex 4937 5/3/2012 23:43 > Eptfur 4937 5/2/2012 0:10 > Eptfur 4937 5/2/2012 0:12 > Eptfur 4937 5/2/2012 0:15 > Eptfur 4937 5/2/2012 0:16 > Eptfur 4937 5/2/2012 0:21 > Eptfur 4937 5/2/2012 0:27 > Eptfur 4937 5/2/2012 0:27 > Eptfur 4937 5/2/2012 0:36 > Eptfur 4937 5/2/2012 0:37 > Eptfur 4937 5/2/2012 0:38 > Eptfur 4937 5/2/2012 0:42 > Eptfur 4937 5/2/2012 0:57 > Eptfur 4937 5/2/2012 0:58 > Eptfur 4937 5/2/2012 0:59 > Eptfur 4937 5/2/2012 1:01 > Eptfur 4937 5/2/2012 1:03 > Eptfur 4937 5/2/2012 1:08 > Eptfur 4937 5/4/2012 2:02 > Eptfur 4937 5/4/2012 2:03 > Eptfur 4937 5/4/2012 2:03 > Eptfur 4937 5/4/2012 2:04 > Eptfur 4937 5/4/2012 2:04 > Eptfur 4937 5/4/2012 2:08 > Eptfur 4924 1/8/2012 2:09 > Eptfur 4937 5/4/2012 2:09 > Eptfur 4924 1/8/2012 2:10 > Eptfur 4924 1/8/2012 2:10 > Eptfur 4924 1/8/2012 2:10 > Eptfur 4924 1/8/2012 2:10 > Eptfur 4924 1/8/2012 2:11 > Eptfur 4924 1/8/2012 2:11 > Eptfur 4924 1/8/2012 2:11 > Eptfur 4924 1/8/2012 2:11 > Eptfur 4924 1/8/2012 2:11 > Eptfur 4924 1/8/2012 2:12 > Eptfur 4924 1/8/2012 2:12 > Eptfur 4924 1/8/2012 2:12 > Eptfur 4924 1/8/2012 2:12 > Eptfur 4924 1/8/2012 2:12 > Eptfur 4924 1/8/2012 2:13 > Eptfur 4924 1/8/2012 2:13 > Eptfur 4937 5/4/2012 2:13 > Eptfur 4924 1/8/2012 2:14 > Eptfur 4924 1/8/2012 2:14 > Eptfur 4924 1/8/2012 2:14 > Eptfur 4924 1/8/2012 2:14 > Eptfur 4937 5/4/2012 2:14 > Eptfur 4924 1/8/2012 2:15 > Eptfur 4924 1/8/2012 2:15 > Eptfur 4924 1/8/2012 2:15 > Eptfur 4924 1/8/2012 2:16 > Eptfur 4924 1/8/2012 2:16 > Eptfur 4924 1/8/2012 2:16 > Eptfur 4924 1/8/2012 2:16 > Eptfur 4924 1/8/2012 2:17 > > > Any suggestions-help appreciated. > > Bruce > > [[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. >[[alternative HTML version deleted]]