R/ssq.delay.R

In TDCor: Gene Network Inference from Time-Series Transcriptomic Data

ssq.delay <-
function(delay,hr,ht,time_l,time_u,time_step,type,deriv)

{

u=seq(time_l,time_u,time_step)

if (deriv)

{

if (delay>=0)

{

k1=hr(u-delay,deriv=1)

k2=ht(u,deriv=1)

}else

{

k1=hr(u,deriv=1)

k2=ht(u+delay,deriv=1)

}

}else

{

if (delay>=0)

{

k1=hr(u-delay)

k2=ht(u)

}else

{

k1=hr(u)

k2=ht(u+delay)

}



}





if (type>0)

{

output=sum((k1-k2)^2)/length(u)

}else

{

if (deriv)

{

output=sum((-k1-k2)^2)/length(u)

}else

{

output=sum((1-k1-k2)^2)/length(u)

}

}

return(output)

}