R/DepNHPPqueueI.R

In IndTestPP: Tests of Independence and Analysis of Dependence Between Point Processes in Time

DepNHPPqueueI <-
function(lambda, d,   T=NULL,  type, 
	ent, Clambda=TRUE, ddist, argd, fixed.seed=NULL)
{

posNHs<-list(PP1=ent)
lambdaM<-cbind(lambda)
tt<-c(1:T)

for (i in c(2:d))
{
	ent<-posNHs[[i-1]]
	nEv<-length(ent)
#	ser<-do.call(paste('r',ddist,sep=''), append(list(n=nEv),append(argd[[1]][i-1],argd[-1])) )
	ser<-do.call(paste('r',ddist,sep=''), append(list(n=nEv),argd[i-1,]) )
	sal<-ent+ser
	posNHs[[i]]<-sort(sal[sal<T])

	if ((type=='NH') &(Clambda==TRUE))
	{
	dendist<-do.call(paste('d',ddist,sep=''), append(list(x=tt),argd[i-1,]) )
	convlambda<-convolve(lambdaM[,(i-1)],rev(dendist),type='o')[1:T]
	lambdaM<-cbind(lambdaM, convlambda)
	}
}
names(posNHs)<-paste('N',c(1:d), sep='')


return(list(posNHs=posNHs, lambdaM=lambdaM))
}