Displaying 1 result from an estimated 1 matches for "datacod".
2009 May 20
1
Comparing spatial distributions - permutation test implementation
..."var2")])
# Compute the statistic (i.e. sum of mean squared differences)
return(sum(meanSqDiff))
}
# Get data online : http://dl.getdropbox.com/u/1047321/syrjala_data_cod.csv
system("curl http://dl.getdropbox.com/u/1047321/syrjala_data_cod.csv >
syrjala_data_cod.csv")
dataCod = read.csv(file = "syrjala_data_cod.csv", header = TRUE)
# Normalize abundances
dataCod$var1 = dataCod$var1/sum(dataCod$var1)
dataCod$var2 = dataCod$var2/sum(dataCod$var2)
# Number of permutations
nperm = 1000
# Create nperm-1 replicates of the data (one is the original
observation)...