write.relationshipMatrix {synbreed} | R Documentation |
This function can be used to write an object of class "relationshipMatrix" in the table format used by other software, i.e. WOMBAT or ASReml. The resulting table has three columns, the row, the column and the entry of the (inverse) relationshipMatrix.
write.relationshipMatrix(x, file = NULL, sorting=c("WOMBAT","ASReml"), type=c("ginv", "inv", "none"), digits = 10)
x |
Object of class "relationshipMatrix" |
file |
Path where the output should be written . If |
sorting |
Type of sorting. Use "WOMBAT" for 'row-wise' sorting of the table and "ASReml" for 'column-wise' sorting. |
type |
A character string indicating which form of |
digits |
Numeric. The result is rounded to |
Note that "WOMBAT" only uses the generalized inverse relationship matrix and expects a file with the name "ranef.gin", where 'ranef' is the name of the random effect with option 'GIN' in the 'MODEL' part of the parameter file. For ASREML, either the relationship could be saved as "*.grm" or its generalized inverse as "*.giv".
Valentin Wimmer
Meyer, K. (2006) WOMBAT - A tool for mixed model analyses in quantitative genetics by REML, J. Zhejinag Uni SCIENCE B 8: 815-821.
Gilmour, A., Cullis B., Welham S., and Thompson R. (2000) ASREML. program user manual. NSW Agriculture, Orange Agricultural Institute, Forest Road, Orange, Australia .
## Not run: # example with 9 individuals id <- 1:9 par1 <- c(0,0,0,0,1,1,1,4,7) par2 <- c(0,0,0,0,2,3,2,5,8) gener <- c(0,0,0,0,1,1,1,2,3) ped <- create.pedigree(id,par1,par2,gener) gp <- create.gpData(pedigree=ped) A <- kin(ped,ret="add") write.relationshipMatrix(A,type="ginv") ## End(Not run)