### Isopam.wrap

isopam.wrap<-function(dat, c.max, ...)
{
c.range<-2:c.max
res.table<-data.frame(K=2:c.max, Ind.n=NA)
for (i in 1:length(c.range)){
  iso.res<-isopam(dat=dat, c.fix=c.range[i])
  res.table$Ind.n[i]<-sum(as.numeric(as.character(iso.res$analytics[7,])))
}
print(res.table)
}


