Hi, 

I believe the spia() function code contains a bug in the line where it checks whether de list of ids is contained within all list of ids. Specifically, it checks 'de' names against 'all' values instead of names:

IDsNotP<-names(de)[!names(de)%in%all]

And it should be:

IDsNotP<-names(de)[!names(de)%in%names(all)]

Thanks!