Displaying 1 result from an estimated 1 matches for "kacmat".
2013 Apr 25
2
Vectorized code for generating the Kac (Clement) matrix
Hi,
I am generating large Kac matrices (also known as Clement matrix). This a tridiagonal matrix. I was wondering whether there is a vectorized solution that avoids the `for' loops to the following code:
n <- 1000
Kacmat <- matrix(0, n+1, n+1)
for (i in 1:n) Kacmat[i, i+1] <- n - i + 1
for (i in 2:(n+1)) Kacmat[i, i-1] <- i-1
The above code is fast, but I am curious about vectorized ways to do this.
Thanks in advance.
Best,
Ravi
Ravi Varadhan, Ph.D.
Assistant Professor
The Center on Aging and Health
D...