generate biSAD1 covariance matrix

Usage

get_biSAD1(par, n1, n2)

Arguments

par: vector with four number, first two for ck and the rest for stress
n1: scalar indicate length of time1
n2: scalar indicate length of time2

Value

biSAD1 covariance matrix

Examples

get_biSAD1(par=c(2,0.5,2,0.1),n1=4, n2 = 5)
#>       [,1] [,2]  [,3]  [,4] [,5] [,6] [,7] [,8] [,9]
#>  [1,] 0.25 0.50  1.00  2.00 0.00 0.00 0.00 0.00 0.00
#>  [2,] 0.50 1.25  2.50  5.00 0.00 0.00 0.00 0.00 0.00
#>  [3,] 1.00 2.50  5.25 10.50 0.00 0.00 0.00 0.00 0.00
#>  [4,] 2.00 5.00 10.50 21.25 0.00 0.00 0.00 0.00 0.00
#>  [5,] 0.00 0.00  0.00  0.00 0.01 0.02 0.04 0.08 0.16
#>  [6,] 0.00 0.00  0.00  0.00 0.02 0.05 0.10 0.20 0.40
#>  [7,] 0.00 0.00  0.00  0.00 0.04 0.10 0.21 0.42 0.84
#>  [8,] 0.00 0.00  0.00  0.00 0.08 0.20 0.42 0.85 1.70
#>  [9,] 0.00 0.00  0.00  0.00 0.16 0.40 0.84 1.70 3.41