/packages/archive/2010/04.2010/04.29.2010/spacodi/R/r.sp.plot.R

 1r.sp.plot <-
 2function(obj,n.rep=10){
 3	sp.plot=obj
 4	zz=unlist(unname(obj))
 5	nulls=length(which(zz==0))
 6	l.data=length(zz)
 7	non.nulls=l.data-nulls
 8	yy=zz[which(zz!=0)]
 9	m.log=mean(log(yy))
10	out <- list() ### prep an output file
11	for(j in 1:n.rep){
12		foo <- as.data.frame(matrix(sample(c(round(rlnorm(n=non.nulls, meanlog=m.log)),rep(0,nulls))),nrow(sp.plot),ncol(sp.plot)))
13		rownames(foo) <- row.names(sp.plot)
14		names(foo)=names(sp.plot)
15		out[[j]]=foo
16	}
17	return(out)
18}
19