R/get.GH.w.R

In mercaldo/MMLB: Marginalized Models for Longitudinal Binary Data

get.GH.w <-
function( q ){
  #
  if( (q!=3)*(q!=5)*(q!=10)*(q!=20)*(q!=50) ){
    stop("get_GH_z():  Error.  q should be 3, 5, 10, 20, or 50.\n");
  }
  #
  PI <- 3.1415927
  #
  if( q==3 ){
    w <- rep( NA, q )
    w[1] <- ( 2.954089751509e-1 ) /sqrt(PI)
    w[2] <- ( 1.181635900604    ) /sqrt(PI)
    w[3] <- ( 2.954089751509e-1 ) /sqrt(PI)
  }
  if( q==5 ){
    w <- rep( NA, q )
    w[1] <- ( 1.995324205905e-2 ) / sqrt(PI)
    w[2] <- ( 3.936193231522e-1 ) / sqrt(PI)
    w[3] <- ( 9.453087204829e-1 ) / sqrt(PI)
    w[4] <- ( 3.936193231522e-1 ) / sqrt(PI)
    w[5] <- ( 1.995324205905e-2 ) / sqrt(PI)
  }
  if( q==10 ){
    w <- rep( NA, q )
    w[1] <- ( 7.640432855233e-6 ) / sqrt(PI)
    w[2] <- ( 1.343645746781e-3 ) / sqrt(PI)
    w[3] <- ( 3.387439445548e-2 ) / sqrt(PI)
    w[4] <- ( 2.401386110823e-1 ) / sqrt(PI)
    w[5] <- ( 6.108626337353e-1 ) / sqrt(PI)
    w[6] <- ( 6.108626337353e-1 ) / sqrt(PI)
    w[7] <- ( 2.401386110823e-1 ) / sqrt(PI)
    w[8] <- ( 3.387439445548e-2 ) / sqrt(PI)
    w[9] <- ( 1.343645746781e-3 ) / sqrt(PI)
    w[10] <- ( 7.640432855233e-6 ) / sqrt(PI)
  }
  if( q==20 ){
    w <- rep( NA, q )
    w[1] <- ( 2.229393645534e-13) / sqrt(PI)
    w[2] <- ( 4.399340992273e-10) / sqrt(PI)
    w[3] <- ( 1.086069370769e-7 ) / sqrt(PI)
    w[4] <- ( 7.802556478532e-6 ) / sqrt(PI)
    w[5] <- ( 2.283386360163e-4 ) / sqrt(PI)
    w[6] <- ( 3.243773342238e-3 ) / sqrt(PI)
    w[7] <- ( 2.481052088746e-2 ) / sqrt(PI)
    w[8] <- ( 1.090172060200e-1 ) / sqrt(PI)
    w[9] <- ( 2.866755053628e-1 ) / sqrt(PI)
    w[10] <- ( 4.622436696006e-1 ) / sqrt(PI)
    w[11] <-( 4.622436696006e-1 ) / sqrt(PI)
    w[12] <-( 2.866755053628e-1 ) / sqrt(PI)
    w[13] <-( 1.090172060200e-1 ) / sqrt(PI)
    w[14] <-( 2.481052088746e-2 ) / sqrt(PI)
    w[15] <-( 3.243773342238e-3 ) / sqrt(PI)
    w[16] <-( 2.283386360163e-4 ) / sqrt(PI)
    w[17] <-( 7.802556478532e-6 ) / sqrt(PI)
    w[18] <-( 1.086069370769e-7 ) / sqrt(PI)
    w[19] <-( 4.399340992273e-10) / sqrt(PI)
    w[20] <-( 2.229393645534e-13) / sqrt(PI)
  }
  if( q==50 ){
    w <- rep( NA, q )
    w[1] <- 1.034608e-37
    w[2] <- 9.443415e-33 
    w[3] <- 6.856281e-29 
    w[4] <- 1.206045e-25 
    w[5] <- 7.995094e-23 
    w[6] <- 2.522483e-20
    w[7] <- 4.368172e-18 
    w[8] <- 4.566698e-16 
    w[9] <- 3.083829e-14 
    w[10] <- 1.414229e-12 
    w[11] <- 4.576637e-11
    w[12] <- 1.077061e-09 
    w[13] <- 1.888226e-08 
    w[14] <- 2.514610e-07 
    w[15] <- 2.584938e-06 
    w[16] <- 2.078485e-05
    w[17] <- 1.321726e-04 
    w[18] <- 6.708281e-04 
    w[19] <- 2.738161e-03 
    w[20] <- 9.045054e-03 
    w[21] <- 2.430481e-02
    w[22] <- 5.334352e-02 
    w[23] <- 9.593054e-02 
    w[24] <- 1.416854e-01 
    w[25] <- 1.721259e-01 
    w[26] <- 1.721259e-01
    w[27] <- 1.416854e-01 
    w[28] <- 9.593054e-02 
    w[29] <- 5.334352e-02 
    w[30] <- 2.430481e-02 
    w[31] <- 9.045054e-03
    w[32] <- 2.738161e-03 
    w[33] <- 6.708281e-04 
    w[34] <- 1.321726e-04 
    w[35] <- 2.078485e-05 
    w[36] <- 2.584938e-06
    w[37] <- 2.514610e-07 
    w[38] <- 1.888226e-08 
    w[39] <- 1.077061e-09 
    w[40] <- 4.576637e-11 
    w[41] <- 1.414229e-12
    w[42] <- 3.083829e-14 
    w[43] <- 4.566698e-16 
    w[44] <- 4.368172e-18 
    w[45] <- 2.522483e-20 
    w[46] <- 7.995094e-23
    w[47] <- 1.206045e-25 
    w[48] <- 6.856281e-29 
    w[49] <- 9.443415e-33 
    w[50] <- 1.034608e-37 
  }
  #
  w
  #
}