--
---
You received this message because you are subscribed to the Google Groups "Deducer" group.
To unsubscribe from this group and stop receiving emails from it, send an email to deducer+u...@googlegroups.com.
For more options, visit https://groups.google.com/d/optout.
# This works
>Corbre_2 <- read.table("E:/Acoustic
Parameters/~Emballonuridae/Data by
genus/~Corbre/Corbre_2.txt",header=T,sep="\t",quote="")
Corbre_2 [c("TBC")] <- recode.variables(Corbre_2[c("TBC")] , "0 -> NA;")
# Create new BW variable then change data file name
BW<-within(Corbre_2,BW<-Fmax-Fmin)
# Create 2 new variables for fundamental H1 and third harmonic H3 #based on Fc that is generally loudest at the 2nd harmonic.
BW<-within(BW,FcH1<-Fc*.5)
BW<-within(BW,FcH3<-FcH1*3)
BatStats<-BW
#added the H1 and H3 in the output section
BatStats<-Deducer::descriptive.table (vars = d (Dur,TBC,Fmin,Fmax,BW,Fmean,Fk,FcH1,Fc,FcH3,Sc,Pmc),data= BatStats,func.names =c("Valid N","Minimum","Maximum","Mean","St. Deviation"),func.additional= list(p10=function(x) quantile(x, c(0.10), na.rm=TRUE),p25=function(x) quantile(x, c(0.25), na.rm=TRUE),p75=function(x) quantile(x, c(0.75),na.rm=TRUE),p90=function(x) quantile(x, c(0.90),na.rm=TRUE)))
## Write the results Change output name and directory for stats.
write.csv(BatStats,file="C:\\=Bat data working\\Acoustic Parameters\\Corbre_Stats.csv")