R/va_2_mat.R

In dengfei2013/learnasreml: The package is to help learn the asreml

#' Trans va in asreml to blupf90
#'
#' @param va The string copyed by asreml
#' @return matri.
#' library(learnasreml)
#' vep4 = "9885.      0.6648      0.5965E-01  0.3229     -0.3258     -0.2968    0.1409E+05  0.4545E+05 -0.3570      0.2671     -0.2211     -0.2284    1.447      -18.56      0.5949E-01 -0.3278E-01  0.7477E-02  0.8339E-02 304.1       539.4     -0.7574E-01   89.72     -0.4089     -0.3238    -16.95      -24.65      0.9539E-03  -2.026      0.2736      0.4048    -11.39      -18.79      0.7850E-03  -1.184      0.8172E-01  0.1490"
#' va_2_mat(vep4)


va_2_mat = function(va){
  library(stringr)
  library(tidyverse)
  va1 = unlist(strsplit(va,"\\s+"))
  len = length(va1)
  va1 = as.numeric(va1)
  va_mat = matrix(va1,sqrt(len),byrow = T)
  mat1 = va_mat
  mat1[upper.tri(mat1)] <- t(mat1)[upper.tri(mat1)]
  return(mat1)
}