R/get.GH.z.R

In mercaldo/MMLB: Marginalized Models for Longitudinal Binary Data

get.GH.z <-
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");
  }
  if( q==3 ){
    z <- rep( NA, q )
    z[1] <- -1.224744871391589 * sqrt(2.0)
    z[2] <-  0.0
    z[3] <-  1.224744871391589 * sqrt(2.0)
  }
  if( q==5 ){
    z <- rep( NA, q )
    z[1] <- -2.020182870456086 * sqrt(2.0)
    z[2] <- -0.958572464613819 * sqrt(2.0)
    z[3] <-  0.000000000000000 * sqrt(2.0)
    z[4] <-  0.958572464613819 * sqrt(2.0)
    z[5] <-  2.020182870456086 * sqrt(2.0)
  }
  if( q==10 ){
    z <- rep( NA, q )
    z[1] <- -3.436159118837738 * sqrt(2.0)
    z[2] <- -2.532731674232790 * sqrt(2.0)
    z[3] <- -1.756683649299882 * sqrt(2.0)
    z[4] <- -1.036610829789514 * sqrt(2.0)
    z[5] <- -0.342901327223705 * sqrt(2.0)
    z[6] <-  0.342901327223705 * sqrt(2.0)
    z[7] <-  1.036610829789514 * sqrt(2.0)
    z[8] <-  1.756683649299882 * sqrt(2.0)
    z[9] <-  2.532731674232790 * sqrt(2.0)
    z[10] <-  3.436159118837738 * sqrt(2.0)
  }
  if( q==20 ){
    z <- rep( NA, q )
    z[1] <- -5.3874808900112 * sqrt(2.0)
    z[2] <- -4.6036824495507 * sqrt(2.0)
    z[3] <- -3.9447640401156 * sqrt(2.0)
    z[4] <- -3.3478545673832 * sqrt(2.0)
    z[5] <- -2.7888060584281 * sqrt(2.0)
    z[6] <- -2.2549740020892 * sqrt(2.0)
    z[7] <- -1.7385377121166 * sqrt(2.0)
    z[8] <- -1.2340762153953 * sqrt(2.0)
    z[9] <- -0.7374737285454 * sqrt(2.0)
    z[10] <- -0.2453407083009 * sqrt(2.0)
    z[11] <- 0.2453407083009 * sqrt(2.0)
    z[12] <- 0.7374737285454 * sqrt(2.0)
    z[13] <- 1.2340762153953 * sqrt(2.0)
    z[14] <- 1.7385377121166 * sqrt(2.0)
    z[15] <- 2.2549740020892 * sqrt(2.0)
    z[16] <- 2.7888060584281 * sqrt(2.0)
    z[17] <- 3.3478545673832 * sqrt(2.0)
    z[18] <- 3.9447640401156 * sqrt(2.0)
    z[19] <- 4.6036824495507 * sqrt(2.0)
    z[20] <- 5.3874808900112 * sqrt(2.0)
  }
  if( q==50 ){
    z <- rep( NA, q )
    z[1] <- -12.9858845
    z[2] <- -12.0530184 
    z[3] <- -11.2792333
    z[4] <- -10.5873817 
    z[5] <- -9.9480357 
    z[6] <- -9.3460396  
    z[7] <- -8.7722996  
    z[8] <- -8.2208159  
    z[9] <- -7.6873624
    z[10] <- -7.1688148  
    z[11] <- -6.6627754  
    z[12] <- -6.1673474  
    z[13] <- -5.6809923  
    z[14] <- -5.2024350  
    z[15] <- -4.7305986
    z[16] <- -4.2645578  
    z[17] <- -3.8035057  
    z[18] <- -3.3467278  
    z[19] <- -2.8935827  
    z[20] <- -2.4434875  
    z[21] <- -1.9959047
    z[22] <- -1.5503332  
    z[23] <- -1.1062993  
    z[24] <- -0.6633497  
    z[25] <- -0.2210452  
    z[26] <- 0.2210452  
    z[27] <- 0.6633497
    z[28] <- 1.1062993   
    z[29] <- 1.5503332   
    z[30] <- 1.9959047   
    z[31] <- 2.4434875   
    z[32] <- 2.8935827   
    z[33] <- 3.3467278
    z[34] <- 3.8035057   
    z[35] <- 4.2645578   
    z[36] <- 4.7305986   
    z[37] <- 5.2024350   
    z[38] <- 5.6809923   
    z[39] <- 6.1673474
    z[40] <- 6.6627754   
    z[41] <- 7.1688148   
    z[42] <- 7.6873624   
    z[43] <- 8.2208159   
    z[44] <- 8.7722996   
    z[45] <- 9.3460396
    z[46] <- 9.9480357   
    z[47] <- 10.5873817   
    z[48] <- 11.2792333  
    z[49] <- 12.0530184  
    z[50] <- 12.9858845  
  }
  #
  z
  #
}