Seed Arrival
Project coordinated by Joe Wright, Yu-Yun Chen and others from CTFS.
Objetive: Interspecific association of seed arrival
Functions
Last function: 02 september 2011 author: Joe Wright
Adjusted function to PEIC - Alexandre 29 September 2011
- small changes to adjust the functions to PEIC data. Only functions are left, so script are on a second file.
Results
Results file from seed rain analysis on PEIC
- {pairwise_fishertest.txt|}}
Ilha do Cardoso
- Script 29 Sep 2011
############################
## Script 29 de setembro de 2011
#################################
chuva<-read.table("/home/ale/Documentos/Ale2010/Ale_Projetos/Restinga/MariDados/sementes/chuvaAgo2011ale.txt", header=TRUE, as.is=TRUE ,sep="\t" )
str(chuva)
tapply(chuva$n.seed>0, chuva$trap, sum)
## excluindo plots fora da parcela e outras formas de vida
seed.plot<-chuva[chuva$forma_vida=="árvore" & chuva$ambiente!="RB",]
str(seed.plot)
tapply(seed.plot$n.seed>0, seed.plot$trap, sum)
##############################
### Dados de dispersão
################################################
disp.peic<-read.table("/home/ale/Documentos/Ale2010/Ale_Projetos/Restinga/MariDados/sementes/disp.peic.txt", header=TRUE, as.is=TRUE ,sep="\t" )
str(disp.peic)
disp.peic
trap.peic=data.frame(sp=seed.plot$code.sp,fecha=seed.plot$data.ale, census=seed.plot$censo, trap=seed.plot$trap, part=2, quantity=seed.plot$n.seeds, stringsAsFactors=FALSE)
str(trap.peic)
fruittoseed=data.frame(sp=unique(trap.peic$sp), f2s=1)
###################
### trap.peic e fruittoseed
########################
trap.peic=data.frame(sp=seed.plot$code.sp,fecha=seed.plot$data.ale, census=seed.plot$censo, trap=seed.plot$trap, part=2, quantity=seed.plot$n.seeds, stringsAsFactors=FALSE)
str(trap.peic)
fruittoseed=data.frame(sp=unique(trap.peic$sp), f2s=1)
#######################
#Codigo para rodar as funções do Joe Wright
###############
## codigo Joe Wright_> modificado para PEIC
source("/home/ale/Documentos/Ale2010/Ale_Projetos/SeedRainJWright/contagious_dispersal_20110902ale.r")
setdir="/home/ale/Documentos/Ale2010/Ale_Projetos/SeedRainJWright/resultaPeic"
setwd(setdir)
prepdata.ale()
randomization1()
randomization2()
PairwiseFisherExact()
PairwiseTimingAllYrs()
PairwiseTimingFisherYrs()
ZTransformTest()