require(vegan)
require(BiodiversityR)

com<-read.csv("com_aveOR.csv", sep=";", dec=",", header=T)
nevek<-read.csv("sepOR.csv", sep=";", dec=",", header=T)


years <- seq(1988,2012)

i <- j <- 1
c <- rainbow(dim(com[rowSums(com)>0,]))
s <- FALSE
yy <- years[rowSums(com)>0]

for(yr in years){
    if(sum(com[i,])>0){
        r <- renyiresult(com,y=nevek,factor='Year',method='s',level=paste("y",yr,sep=""))
        renyiplot(r,addit=s,legend=FALSE,col=c[j],ylim=c(0,3))
        j <- j+1
    }
    s <- TRUE
    i <- i+1
}

legend(x=1,y=2.2,legend=yy,col=c,lty=1,pch=1)