[[projetos:pp_peic:restrito:peic|Voltar]]
====== 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

  * {{contagious_dispersal_20110902.r|}}

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.
  * {{contagious_dispersal_20110902ale.r|}}

==== Results ====

Results file from seed rain analysis on PEIC

  *{{randomization1.txt|}} 
  *{{randomization2.txt|}}
  *{pairwise_fishertest.txt|}}
  *{{pairwise_timing_allyears.txt|}}
  *{{pairwise_timing_fisheryrs.txt|}}

===== Ilha do Cardoso =====

  *  Script 29 Sep 2011
<code>
############################
## 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()

</code>