[R] More

arun smartpink111 at yahoo.com
Mon Apr 22 20:03:35 CEST 2013


Hi Bruce,

This could be also done by:
dat1<- read.csv("Sample_all_3_locations.csv",header=TRUE,stringsAsFactors=FALSE,sep=",")
library(plyr)
res1<-ddply(unique(dat1[,c(1,3:4)]),.(Location_name),summarize,Survey_Time=sum(Survey_Time)) #changed here
res1
#      Location_name Survey_Time
#1 1_f3p1_La_Campana        46.4
#2 2_f4p1_La_Campana        69.5
#3   79_f2p1_Acetuna        36.4

#The rest as mentioned before.
res2<-ddply(dat1,.(Location_name,SPEC_CODE),summarize, AI=sum(AI))
 Final<- mutate(join(res1,res2),RA=10*(AI/Survey_Time))


A.K.




----- Original Message -----
From: Bruce Miller <batsncats at gmail.com>
To: arun <smartpink111 at yahoo.com>
Cc: 
Sent: Monday, April 22, 2013 12:55 PM
Subject: Re: More

Hi Arun,

My last message overlapped yours.
YES...
that seems correct.
Now if I can get this to work with teh entire data set all is be great.

Thanks again for your valuable time and help.

Bruce


On 4/22/2013 12:16 PM, arun wrote:
> Hi Bruce,
>
>  From your excel sheets, it looks like for the step1, you were just taking only the "Survey Time" of the first rows of unique Start_Day for each Location_name.
>
> Step2 results looks the same as "res2"
> If this is the case:
>
> dat1<- read.csv("Sample_all_3_locations.csv",header=TRUE,stringsAsFactors=FALSE,sep=",")
> library(plyr)
>
> res1<-ddply(aggregate(Survey_Time~Location_name+Start_Day,data=dat1,function(x) head(x,1)),.(Location_name),summarize,Survey_Time=sum(Survey_Time))
> res1
> #      Location_name Survey_Time
> #1 1_f3p1_La_Campana        46.4
> #2 2_f4p1_La_Campana        69.5
> #3   79_f2p1_Acetuna        36.4
>
> res2<-ddply(dat1,.(Location_name,SPEC_CODE),summarize, AI=sum(AI))
>
> Final<- mutate(join(res1,res2),RA=10*(AI/Survey_Time))
>   head(Final)
> #      Location_name Survey_Time SPEC_CODE  AI         RA
> #1 1_f3p1_La_Campana        46.4      Buzz  20  4.3103448
> #2 1_f3p1_La_Campana        46.4    Cencen  20  4.3103448
> #3 1_f3p1_La_Campana        46.4    Eptfur 215 46.3362069
> #4 1_f3p1_La_Campana        46.4      Frag   3  0.6465517
> #5 1_f3p1_La_Campana        46.4    Molspp  22  4.7413793
> #6 1_f3p1_La_Campana        46.4    Myoele   1  0.2155172
> dim(Final)
> #[1] 44  5
> A.K.
>
>
> ----- Original Message -----
> From: Bruce Miller <batsncats at gmail.com>
> To: arun <smartpink111 at yahoo.com>
> Cc:
> Sent: Monday, April 22, 2013 11:32 AM
> Subject: Re: More
>
> Arun,
>
> Attached are two Excel sheets that I manually did  what I am hoping I can automate in R for the larger data sets.
>
> These use the same 3 locations and data I sent earlier.
>
> Step 1 sums the total survey time for each location then step 2 sums the AI values for he SPEC_CODE by location.
>
> Then the calculation of Summed AI value/Summed Survey time can be completed for each location and SPEC_Code.
> Hopefully this clarifies it more.
>
> Bruce
>
> -- Bruce W. Miller, Ph.D.
> Conservation Ecologist
> Neotropical Bat Projects
>
>
> office details
> 11384 Alpine Road
> Stanwood, Mi. 49346
> Phone (231) 679-6059
>


-- 
Bruce W. Miller, Ph.D.
Conservation Ecologist
Neotropical Bat Projects


office details
11384 Alpine Road
Stanwood, Mi. 49346
Phone (231) 679-6059



More information about the R-help mailing list