https://github.com/cran/mvtBinaryEP
Tip revision: e03023db232593f436717e7be27d1c52d9865644 authored by Kunthel By on 18 February 2009, 00:00:00 UTC
version 1.0
version 1.0
Tip revision: e03023d
xch.R
`xch` <-
function(p,rho)
{
if (p <= 0) {stop('n must be at least 1')}
else if (p==1) {r = 1}
else if (p > 1)
{
c1 = rep(rho,p-1)
r = toeplitz(c(1,c1))
}
return(r)
}