"use strict"; const {jsx: _jsx, jsxs: _jsxs} = arguments[0]; const {useMDXComponents: _provideComponents} = arguments[0]; function _createMdxContent(props) { const _components = { code: "code", pre: "pre", span: "span", ..._provideComponents(), ...props.components }; return _jsx(_components.pre, { className: "language-r", children: _jsxs(_components.code, { className: "language-R code-highlight", children: [_jsxs(_components.span, { className: "code-line", children: ["ABLUP ", _jsx(_components.span, { className: "token operator", children: "<-" }), " ", _jsx(_components.span, { className: "token keyword", children: "function" }), _jsx(_components.span, { className: "token punctuation", children: "(" }), "pedfile", _jsx(_components.span, { className: "token punctuation", children: "," }), " datafile", _jsx(_components.span, { className: "token punctuation", children: "," }), " vars ", _jsx(_components.span, { className: "token operator", children: "=" }), " ", _jsx(_components.span, { className: "token number", children: "5" }), _jsx(_components.span, { className: "token punctuation", children: "," }), " vare ", _jsx(_components.span, { className: "token operator", children: "=" }), " ", _jsx(_components.span, { className: "token number", children: "55" }), _jsx(_components.span, { className: "token punctuation", children: ")" }), " ", _jsx(_components.span, { className: "token punctuation", children: "{" }), "\n"] }), _jsxs(_components.span, { className: "code-line", children: [" ", _jsx(_components.span, { className: "token comment", children: "# Read pedigree and data files" }), "\n"] }), _jsxs(_components.span, { className: "code-line", children: [" peddat ", _jsx(_components.span, { className: "token operator", children: "<-" }), " read.csv", _jsx(_components.span, { className: "token punctuation", children: "(" }), "file ", _jsx(_components.span, { className: "token operator", children: "=" }), " pedfile", _jsx(_components.span, { className: "token punctuation", children: "," }), " header ", _jsx(_components.span, { className: "token operator", children: "=" }), " ", _jsx(_components.span, { className: "token boolean", children: "TRUE" }), _jsx(_components.span, { className: "token punctuation", children: ")" }), "\n"] }), _jsxs(_components.span, { className: "code-line", children: [" wkdat ", _jsx(_components.span, { className: "token operator", children: "<-" }), " read.csv", _jsx(_components.span, { className: "token punctuation", children: "(" }), "file ", _jsx(_components.span, { className: "token operator", children: "=" }), " datafile", _jsx(_components.span, { className: "token punctuation", children: "," }), " header ", _jsx(_components.span, { className: "token operator", children: "=" }), " ", _jsx(_components.span, { className: "token boolean", children: "TRUE" }), _jsx(_components.span, { className: "token punctuation", children: ")" }), "\n"] }), _jsx(_components.span, { className: "code-line", children: " \n" }), _jsxs(_components.span, { className: "code-line", children: [" ", _jsx(_components.span, { className: "token comment", children: "# Compute the A matrix" }), "\n"] }), _jsxs(_components.span, { className: "code-line", children: [" A ", _jsx(_components.span, { className: "token operator", children: "<-" }), " Amatrix", _jsx(_components.span, { className: "token punctuation", children: "(" }), "peddat", _jsx(_components.span, { className: "token punctuation", children: "," }), " ploidy ", _jsx(_components.span, { className: "token operator", children: "=" }), " ", _jsx(_components.span, { className: "token number", children: "2" }), _jsx(_components.span, { className: "token punctuation", children: ")" }), "\n"] }), _jsx(_components.span, { className: "code-line", children: " \n" }), _jsxs(_components.span, { className: "code-line", children: [" ", _jsx(_components.span, { className: "token comment", children: "# Compute the inverse of the additive genetic matrix" }), "\n"] }), _jsxs(_components.span, { className: "code-line", children: [" IG ", _jsx(_components.span, { className: "token operator", children: "<-" }), " ", _jsx(_components.span, { className: "token punctuation", children: "(" }), "vare ", _jsx(_components.span, { className: "token operator", children: "/" }), " vars", _jsx(_components.span, { className: "token punctuation", children: ")" }), " ", _jsx(_components.span, { className: "token operator", children: "*" }), " solve", _jsx(_components.span, { className: "token punctuation", children: "(" }), "A", _jsx(_components.span, { className: "token punctuation", children: ")" }), "\n"] }), _jsx(_components.span, { className: "code-line", children: " \n" }), _jsxs(_components.span, { className: "code-line", children: [" ", _jsx(_components.span, { className: "token comment", children: "# Create model matrices" }), "\n"] }), _jsxs(_components.span, { className: "code-line", children: [" sex ", _jsx(_components.span, { className: "token operator", children: "<-" }), " as.factor", _jsx(_components.span, { className: "token punctuation", children: "(" }), "wkdat", _jsx(_components.span, { className: "token operator", children: "$" }), "sex", _jsx(_components.span, { className: "token punctuation", children: ")" }), "\n"] }), _jsxs(_components.span, { className: "code-line", children: [" X ", _jsx(_components.span, { className: "token operator", children: "<-" }), " model.matrix", _jsx(_components.span, { className: "token punctuation", children: "(" }), _jsx(_components.span, { className: "token operator", children: "~" }), " sex ", _jsx(_components.span, { className: "token operator", children: "-" }), " ", _jsx(_components.span, { className: "token number", children: "1" }), _jsx(_components.span, { className: "token punctuation", children: ")" }), "\n"] }), _jsx(_components.span, { className: "code-line", children: " \n" }), _jsxs(_components.span, { className: "code-line", children: [" sire ", _jsx(_components.span, { className: "token operator", children: "<-" }), " as.factor", _jsx(_components.span, { className: "token punctuation", children: "(" }), "wkdat", _jsx(_components.span, { className: "token operator", children: "$" }), "sire", _jsx(_components.span, { className: "token punctuation", children: ")" }), "\n"] }), _jsxs(_components.span, { className: "code-line", children: [" Z ", _jsx(_components.span, { className: "token operator", children: "<-" }), " model.matrix", _jsx(_components.span, { className: "token punctuation", children: "(" }), _jsx(_components.span, { className: "token operator", children: "~" }), " sire ", _jsx(_components.span, { className: "token operator", children: "-" }), " ", _jsx(_components.span, { className: "token number", children: "1" }), _jsx(_components.span, { className: "token punctuation", children: ")" }), "\n"] }), _jsx(_components.span, { className: "code-line", children: " \n" }), _jsxs(_components.span, { className: "code-line", children: [" ", _jsx(_components.span, { className: "token comment", children: "# Phenotype data" }), "\n"] }), _jsxs(_components.span, { className: "code-line", children: [" y ", _jsx(_components.span, { className: "token operator", children: "<-" }), " wkdat", _jsx(_components.span, { className: "token operator", children: "$" }), "yobs\n"] }), _jsx(_components.span, { className: "code-line", children: " \n" }), _jsxs(_components.span, { className: "code-line", children: [" ", _jsx(_components.span, { className: "token comment", children: "# Set up MME components" }), "\n"] }), _jsxs(_components.span, { className: "code-line", children: [" XX ", _jsx(_components.span, { className: "token operator", children: "<-" }), " t", _jsx(_components.span, { className: "token punctuation", children: "(" }), "X", _jsx(_components.span, { className: "token punctuation", children: ")" }), " ", _jsx(_components.span, { className: "token percent-operator operator", children: "%*%" }), " X\n"] }), _jsx(_components.span, { className: "code-line", children: " \n" }), _jsxs(_components.span, { className: "code-line", children: [" XZ ", _jsx(_components.span, { className: "token operator", children: "<-" }), " t", _jsx(_components.span, { className: "token punctuation", children: "(" }), "X", _jsx(_components.span, { className: "token punctuation", children: ")" }), " ", _jsx(_components.span, { className: "token percent-operator operator", children: "%*%" }), " Z\n"] }), _jsx(_components.span, { className: "code-line", children: " \n" }), _jsxs(_components.span, { className: "code-line", children: [" ZZ ", _jsx(_components.span, { className: "token operator", children: "<-" }), " t", _jsx(_components.span, { className: "token punctuation", children: "(" }), "Z", _jsx(_components.span, { className: "token punctuation", children: ")" }), " ", _jsx(_components.span, { className: "token percent-operator operator", children: "%*%" }), " Z\n"] }), _jsx(_components.span, { className: "code-line", children: " \n" }), _jsxs(_components.span, { className: "code-line", children: [" LHS1 ", _jsx(_components.span, { className: "token operator", children: "<-" }), " cbind", _jsx(_components.span, { className: "token punctuation", children: "(" }), "XX", _jsx(_components.span, { className: "token punctuation", children: "," }), " XZ", _jsx(_components.span, { className: "token punctuation", children: ")" }), "\n"] }), _jsxs(_components.span, { className: "code-line", children: [" LHS2 ", _jsx(_components.span, { className: "token operator", children: "<-" }), " cbind", _jsx(_components.span, { className: "token punctuation", children: "(" }), "t", _jsx(_components.span, { className: "token punctuation", children: "(" }), "XZ", _jsx(_components.span, { className: "token punctuation", children: ")" }), _jsx(_components.span, { className: "token punctuation", children: "," }), " ", _jsx(_components.span, { className: "token punctuation", children: "(" }), "ZZ ", _jsx(_components.span, { className: "token operator", children: "+" }), " IG", _jsx(_components.span, { className: "token punctuation", children: ")" }), _jsx(_components.span, { className: "token punctuation", children: ")" }), "\n"] }), _jsxs(_components.span, { className: "code-line", children: [" LHS ", _jsx(_components.span, { className: "token operator", children: "<-" }), " rbind", _jsx(_components.span, { className: "token punctuation", children: "(" }), "LHS1", _jsx(_components.span, { className: "token punctuation", children: "," }), " LHS2", _jsx(_components.span, { className: "token punctuation", children: ")" }), "\n"] }), _jsx(_components.span, { className: "code-line", children: " \n" }), _jsxs(_components.span, { className: "code-line", children: [" RHS1 ", _jsx(_components.span, { className: "token operator", children: "<-" }), " t", _jsx(_components.span, { className: "token punctuation", children: "(" }), "X", _jsx(_components.span, { className: "token punctuation", children: ")" }), " ", _jsx(_components.span, { className: "token percent-operator operator", children: "%*%" }), " y\n"] }), _jsxs(_components.span, { className: "code-line", children: [" RHS2 ", _jsx(_components.span, { className: "token operator", children: "<-" }), " t", _jsx(_components.span, { className: "token punctuation", children: "(" }), "Z", _jsx(_components.span, { className: "token punctuation", children: ")" }), " ", _jsx(_components.span, { className: "token percent-operator operator", children: "%*%" }), " y\n"] }), _jsxs(_components.span, { className: "code-line", children: [" RHS ", _jsx(_components.span, { className: "token operator", children: "<-" }), " rbind", _jsx(_components.span, { className: "token punctuation", children: "(" }), "RHS1", _jsx(_components.span, { className: "token punctuation", children: "," }), " RHS2", _jsx(_components.span, { className: "token punctuation", children: ")" }), "\n"] }), _jsx(_components.span, { className: "code-line", children: " \n" }), _jsxs(_components.span, { className: "code-line", children: [" ", _jsx(_components.span, { className: "token comment", children: "# Solve the mixed model equations" }), "\n"] }), _jsxs(_components.span, { className: "code-line", children: [" sol ", _jsx(_components.span, { className: "token operator", children: "<-" }), " solve", _jsx(_components.span, { className: "token punctuation", children: "(" }), "LHS", _jsx(_components.span, { className: "token punctuation", children: ")" }), " ", _jsx(_components.span, { className: "token percent-operator operator", children: "%*%" }), " RHS\n"] }), _jsxs(_components.span, { className: "code-line", children: [" return", _jsx(_components.span, { className: "token punctuation", children: "(" }), "sol", _jsx(_components.span, { className: "token punctuation", children: ")" }), "\n"] }), _jsxs(_components.span, { className: "code-line", children: [_jsx(_components.span, { className: "token punctuation", children: "}" }), "\n"] }), _jsx(_components.span, { className: "code-line", children: "\n" }), _jsx(_components.span, { className: "code-line", children: "\n" }), _jsxs(_components.span, { className: "code-line", children: ["RRBLUP ", _jsx(_components.span, { className: "token operator", children: "<-" }), " ", _jsx(_components.span, { className: "token keyword", children: "function" }), _jsx(_components.span, { className: "token punctuation", children: "(" }), "data", _jsx(_components.span, { className: "token punctuation", children: "," }), " trtNam", _jsx(_components.span, { className: "token punctuation", children: "," }), " trnIdx", _jsx(_components.span, { className: "token punctuation", children: "," }), " tstIdx", _jsx(_components.span, { className: "token punctuation", children: "," }), " snp.start", _jsx(_components.span, { className: "token punctuation", children: "," }), " snp.end", _jsx(_components.span, { className: "token punctuation", children: "," }), "\n"] }), _jsxs(_components.span, { className: "code-line", children: [" mean2pq", _jsx(_components.span, { className: "token operator", children: "=" }), _jsx(_components.span, { className: "token number", children: "3.54" }), _jsx(_components.span, { className: "token punctuation", children: "," }), "vare ", _jsx(_components.span, { className: "token operator", children: "=" }), " ", _jsx(_components.span, { className: "token number", children: "245.0" }), _jsx(_components.span, { className: "token punctuation", children: "," }), "vara ", _jsx(_components.span, { className: "token operator", children: "=" }), " ", _jsx(_components.span, { className: "token number", children: "35.24" }), _jsx(_components.span, { className: "token punctuation", children: ")" }), " ", _jsx(_components.span, { className: "token punctuation", children: "{" }), "\n"] }), _jsxs(_components.span, { className: "code-line", children: [" ", _jsx(_components.span, { className: "token comment", children: "# hyper-parameters and additional index variables" }), "\n"] }), _jsxs(_components.span, { className: "code-line", children: [" lambda ", _jsx(_components.span, { className: "token operator", children: "<-" }), " mean2pq ", _jsx(_components.span, { className: "token operator", children: "*" }), " vare ", _jsx(_components.span, { className: "token operator", children: "/" }), " vara\n"] }), _jsx(_components.span, { className: "code-line", children: " \n" }), _jsxs(_components.span, { className: "code-line", children: [" ", _jsx(_components.span, { className: "token comment", children: "# Read and preprocess data" }), "\n"] }), _jsxs(_components.span, { className: "code-line", children: [" wkdat ", _jsx(_components.span, { className: "token operator", children: "<-" }), " read.csv", _jsx(_components.span, { className: "token punctuation", children: "(" }), "file ", _jsx(_components.span, { className: "token operator", children: "=" }), " data", _jsx(_components.span, { className: "token punctuation", children: "," }), " header ", _jsx(_components.span, { className: "token operator", children: "=" }), " ", _jsx(_components.span, { className: "token boolean", children: "TRUE" }), _jsx(_components.span, { className: "token punctuation", children: ")" }), "\n"] }), _jsxs(_components.span, { className: "code-line", children: [" yidx ", _jsx(_components.span, { className: "token operator", children: "<-" }), " grep", _jsx(_components.span, { className: "token punctuation", children: "(" }), "trtNam", _jsx(_components.span, { className: "token punctuation", children: "," }), " colnames", _jsx(_components.span, { className: "token punctuation", children: "(" }), "wkdat", _jsx(_components.span, { className: "token punctuation", children: ")" }), _jsx(_components.span, { className: "token punctuation", children: ")" }), "\n"] }), _jsxs(_components.span, { className: "code-line", children: [" nInd ", _jsx(_components.span, { className: "token operator", children: "<-" }), " nrow", _jsx(_components.span, { className: "token punctuation", children: "(" }), "wkdat", _jsx(_components.span, { className: "token punctuation", children: ")" }), "\n"] }), _jsx(_components.span, { className: "code-line", children: " \n" }), _jsxs(_components.span, { className: "code-line", children: [" ", _jsx(_components.span, { className: "token comment", children: "# Incidence vector for the overall mean as the only fixed effect" }), "\n"] }), _jsxs(_components.span, { className: "code-line", children: [" X ", _jsx(_components.span, { className: "token operator", children: "<-" }), " matrix", _jsx(_components.span, { className: "token punctuation", children: "(" }), _jsx(_components.span, { className: "token number", children: "1" }), _jsx(_components.span, { className: "token punctuation", children: "," }), " nrow ", _jsx(_components.span, { className: "token operator", children: "=" }), " nInd", _jsx(_components.span, { className: "token punctuation", children: "," }), " ncol ", _jsx(_components.span, { className: "token operator", children: "=" }), " ", _jsx(_components.span, { className: "token number", children: "1" }), _jsx(_components.span, { className: "token punctuation", children: ")" }), "\n"] }), _jsx(_components.span, { className: "code-line", children: " \n" }), _jsxs(_components.span, { className: "code-line", children: [" ", _jsx(_components.span, { className: "token comment", children: "# SNP genotypes" }), "\n"] }), _jsxs(_components.span, { className: "code-line", children: [" Z ", _jsx(_components.span, { className: "token operator", children: "<-" }), " as.matrix", _jsx(_components.span, { className: "token punctuation", children: "(" }), "wkdat", _jsx(_components.span, { className: "token punctuation", children: "[" }), _jsx(_components.span, { className: "token punctuation", children: "," }), " snp.start", _jsx(_components.span, { className: "token operator", children: ":" }), "snp.end", _jsx(_components.span, { className: "token punctuation", children: "]" }), _jsx(_components.span, { className: "token punctuation", children: ")" }), "\n"] }), _jsxs(_components.span, { className: "code-line", children: [" p ", _jsx(_components.span, { className: "token operator", children: "<-" }), " apply", _jsx(_components.span, { className: "token punctuation", children: "(" }), "Z", _jsx(_components.span, { className: "token punctuation", children: "," }), " ", _jsx(_components.span, { className: "token number", children: "2" }), _jsx(_components.span, { className: "token punctuation", children: "," }), " mean", _jsx(_components.span, { className: "token punctuation", children: ")" }), " ", _jsx(_components.span, { className: "token operator", children: "/" }), " ", _jsx(_components.span, { className: "token number", children: "2" }), "\n"] }), _jsxs(_components.span, { className: "code-line", children: [" nSNP ", _jsx(_components.span, { className: "token operator", children: "<-" }), " length", _jsx(_components.span, { className: "token punctuation", children: "(" }), "p", _jsx(_components.span, { className: "token punctuation", children: ")" }), "\n"] }), _jsx(_components.span, { className: "code-line", children: " \n" }), _jsxs(_components.span, { className: "code-line", children: [" ", _jsx(_components.span, { className: "token keyword", children: "for" }), " ", _jsx(_components.span, { className: "token punctuation", children: "(" }), "j ", _jsx(_components.span, { className: "token keyword", children: "in" }), " ", _jsx(_components.span, { className: "token number", children: "1" }), _jsx(_components.span, { className: "token operator", children: ":" }), "ncol", _jsx(_components.span, { className: "token punctuation", children: "(" }), "Z", _jsx(_components.span, { className: "token punctuation", children: ")" }), _jsx(_components.span, { className: "token punctuation", children: ")" }), " ", _jsx(_components.span, { className: "token punctuation", children: "{" }), "\n"] }), _jsxs(_components.span, { className: "code-line", children: [" Z", _jsx(_components.span, { className: "token punctuation", children: "[" }), _jsx(_components.span, { className: "token punctuation", children: "," }), "j", _jsx(_components.span, { className: "token punctuation", children: "]" }), " ", _jsx(_components.span, { className: "token operator", children: "<-" }), " Z", _jsx(_components.span, { className: "token punctuation", children: "[" }), _jsx(_components.span, { className: "token punctuation", children: "," }), "j", _jsx(_components.span, { className: "token punctuation", children: "]" }), " ", _jsx(_components.span, { className: "token operator", children: "-" }), " ", _jsx(_components.span, { className: "token number", children: "2" }), " ", _jsx(_components.span, { className: "token operator", children: "*" }), " p", _jsx(_components.span, { className: "token punctuation", children: "[" }), "j", _jsx(_components.span, { className: "token punctuation", children: "]" }), "\n"] }), _jsxs(_components.span, { className: "code-line", children: [" ", _jsx(_components.span, { className: "token punctuation", children: "}" }), "\n"] }), _jsx(_components.span, { className: "code-line", children: " \n" }), _jsxs(_components.span, { className: "code-line", children: [" ", _jsx(_components.span, { className: "token comment", children: "# adjusted phenotypes or breeding varlues" }), "\n"] }), _jsxs(_components.span, { className: "code-line", children: [" y ", _jsx(_components.span, { className: "token operator", children: "<-" }), " wkdat", _jsx(_components.span, { className: "token punctuation", children: "[" }), _jsx(_components.span, { className: "token punctuation", children: "," }), " yidx", _jsx(_components.span, { className: "token punctuation", children: "]" }), "\n"] }), _jsx(_components.span, { className: "code-line", children: " \n" }), _jsxs(_components.span, { className: "code-line", children: [" ", _jsx(_components.span, { className: "token comment", children: "# Set up MME" }), "\n"] }), _jsxs(_components.span, { className: "code-line", children: [" XX ", _jsx(_components.span, { className: "token operator", children: "<-" }), " t", _jsx(_components.span, { className: "token punctuation", children: "(" }), "X", _jsx(_components.span, { className: "token punctuation", children: "[" }), "trnIdx", _jsx(_components.span, { className: "token punctuation", children: "," }), _jsx(_components.span, { className: "token punctuation", children: "]" }), _jsx(_components.span, { className: "token punctuation", children: ")" }), " ", _jsx(_components.span, { className: "token percent-operator operator", children: "%*%" }), " X", _jsx(_components.span, { className: "token punctuation", children: "[" }), "trnIdx", _jsx(_components.span, { className: "token punctuation", children: "," }), _jsx(_components.span, { className: "token punctuation", children: "]" }), "\n"] }), _jsxs(_components.span, { className: "code-line", children: [" XZ ", _jsx(_components.span, { className: "token operator", children: "<-" }), " t", _jsx(_components.span, { className: "token punctuation", children: "(" }), "X", _jsx(_components.span, { className: "token punctuation", children: "[" }), "trnIdx", _jsx(_components.span, { className: "token punctuation", children: "," }), _jsx(_components.span, { className: "token punctuation", children: "]" }), _jsx(_components.span, { className: "token punctuation", children: ")" }), " ", _jsx(_components.span, { className: "token percent-operator operator", children: "%*%" }), " Z", _jsx(_components.span, { className: "token punctuation", children: "[" }), "trnIdx", _jsx(_components.span, { className: "token punctuation", children: "," }), _jsx(_components.span, { className: "token punctuation", children: "]" }), "\n"] }), _jsxs(_components.span, { className: "code-line", children: [" ZZ ", _jsx(_components.span, { className: "token operator", children: "<-" }), " t", _jsx(_components.span, { className: "token punctuation", children: "(" }), "Z", _jsx(_components.span, { className: "token punctuation", children: "[" }), "trnIdx", _jsx(_components.span, { className: "token punctuation", children: "," }), _jsx(_components.span, { className: "token punctuation", children: "]" }), _jsx(_components.span, { className: "token punctuation", children: ")" }), " ", _jsx(_components.span, { className: "token percent-operator operator", children: "%*%" }), " Z", _jsx(_components.span, { className: "token punctuation", children: "[" }), "trnIdx", _jsx(_components.span, { className: "token punctuation", children: "," }), _jsx(_components.span, { className: "token punctuation", children: "]" }), "\n"] }), _jsxs(_components.span, { className: "code-line", children: [" Xy ", _jsx(_components.span, { className: "token operator", children: "<-" }), " t", _jsx(_components.span, { className: "token punctuation", children: "(" }), "X", _jsx(_components.span, { className: "token punctuation", children: "[" }), "trnIdx", _jsx(_components.span, { className: "token punctuation", children: "]" }), _jsx(_components.span, { className: "token punctuation", children: ")" }), " ", _jsx(_components.span, { className: "token percent-operator operator", children: "%*%" }), " y", _jsx(_components.span, { className: "token punctuation", children: "[" }), "trnIdx", _jsx(_components.span, { className: "token punctuation", children: "]" }), "\n"] }), _jsxs(_components.span, { className: "code-line", children: [" Zy ", _jsx(_components.span, { className: "token operator", children: "<-" }), " t", _jsx(_components.span, { className: "token punctuation", children: "(" }), "Z", _jsx(_components.span, { className: "token punctuation", children: "[" }), "trnIdx", _jsx(_components.span, { className: "token punctuation", children: "," }), _jsx(_components.span, { className: "token punctuation", children: "]" }), _jsx(_components.span, { className: "token punctuation", children: ")" }), " ", _jsx(_components.span, { className: "token percent-operator operator", children: "%*%" }), " y", _jsx(_components.span, { className: "token punctuation", children: "[" }), "trnIdx", _jsx(_components.span, { className: "token punctuation", children: "]" }), "\n"] }), _jsx(_components.span, { className: "code-line", children: " \n" }), _jsxs(_components.span, { className: "code-line", children: [" LHS1 ", _jsx(_components.span, { className: "token operator", children: "<-" }), " cbind", _jsx(_components.span, { className: "token punctuation", children: "(" }), "XX", _jsx(_components.span, { className: "token punctuation", children: "," }), " XZ", _jsx(_components.span, { className: "token punctuation", children: ")" }), "\n"] }), _jsxs(_components.span, { className: "code-line", children: [" LHS2 ", _jsx(_components.span, { className: "token operator", children: "<-" }), " cbind", _jsx(_components.span, { className: "token punctuation", children: "(" }), "t", _jsx(_components.span, { className: "token punctuation", children: "(" }), "XZ", _jsx(_components.span, { className: "token punctuation", children: ")" }), _jsx(_components.span, { className: "token punctuation", children: "," }), " ", _jsx(_components.span, { className: "token punctuation", children: "(" }), "ZZ ", _jsx(_components.span, { className: "token operator", children: "+" }), " diag", _jsx(_components.span, { className: "token punctuation", children: "(" }), "lambda", _jsx(_components.span, { className: "token punctuation", children: "," }), " nrow ", _jsx(_components.span, { className: "token operator", children: "=" }), " nSNP", _jsx(_components.span, { className: "token punctuation", children: "," }), " ncol ", _jsx(_components.span, { className: "token operator", children: "=" }), " nSNP", _jsx(_components.span, { className: "token punctuation", children: ")" }), _jsx(_components.span, { className: "token punctuation", children: ")" }), _jsx(_components.span, { className: "token punctuation", children: ")" }), "\n"] }), _jsxs(_components.span, { className: "code-line", children: [" LHS ", _jsx(_components.span, { className: "token operator", children: "<-" }), " rbind", _jsx(_components.span, { className: "token punctuation", children: "(" }), "LHS1", _jsx(_components.span, { className: "token punctuation", children: "," }), " LHS2", _jsx(_components.span, { className: "token punctuation", children: ")" }), "\n"] }), _jsx(_components.span, { className: "code-line", children: " \n" }), _jsxs(_components.span, { className: "code-line", children: [" RHS ", _jsx(_components.span, { className: "token operator", children: "<-" }), " rbind", _jsx(_components.span, { className: "token punctuation", children: "(" }), "Xy", _jsx(_components.span, { className: "token punctuation", children: "," }), " Zy", _jsx(_components.span, { className: "token punctuation", children: ")" }), "\n"] }), _jsx(_components.span, { className: "code-line", children: " \n" }), _jsxs(_components.span, { className: "code-line", children: [" ", _jsx(_components.span, { className: "token comment", children: "# model solutions (i.e., overal mean and SNP effects)" }), "\n"] }), _jsxs(_components.span, { className: "code-line", children: [" sol ", _jsx(_components.span, { className: "token operator", children: "<-" }), " solve", _jsx(_components.span, { className: "token punctuation", children: "(" }), "LHS", _jsx(_components.span, { className: "token punctuation", children: ")" }), " ", _jsx(_components.span, { className: "token percent-operator operator", children: "%*%" }), " RHS\n"] }), _jsx(_components.span, { className: "code-line", children: " \n" }), _jsxs(_components.span, { className: "code-line", children: [" ", _jsx(_components.span, { className: "token comment", children: "# predictions" }), "\n"] }), _jsxs(_components.span, { className: "code-line", children: [" predictions ", _jsx(_components.span, { className: "token operator", children: "<-" }), " Z", _jsx(_components.span, { className: "token punctuation", children: "[" }), "tstIdx", _jsx(_components.span, { className: "token punctuation", children: "," }), _jsx(_components.span, { className: "token punctuation", children: "]" }), " ", _jsx(_components.span, { className: "token percent-operator operator", children: "%*%" }), " sol", _jsx(_components.span, { className: "token punctuation", children: "[" }), _jsx(_components.span, { className: "token number", children: "2" }), _jsx(_components.span, { className: "token operator", children: ":" }), "length", _jsx(_components.span, { className: "token punctuation", children: "(" }), "sol", _jsx(_components.span, { className: "token punctuation", children: ")" }), _jsx(_components.span, { className: "token punctuation", children: "]" }), "\n"] }), _jsxs(_components.span, { className: "code-line", children: [" row.names", _jsx(_components.span, { className: "token punctuation", children: "(" }), "predictions", _jsx(_components.span, { className: "token punctuation", children: ")" }), _jsx(_components.span, { className: "token operator", children: "<-" }), "wkdat", _jsx(_components.span, { className: "token operator", children: "$" }), "animal", _jsx(_components.span, { className: "token punctuation", children: "[" }), "tstIdx", _jsx(_components.span, { className: "token punctuation", children: "]" }), "\n"] }), _jsx(_components.span, { className: "code-line", children: " \n" }), _jsxs(_components.span, { className: "code-line", children: [" return", _jsx(_components.span, { className: "token punctuation", children: "(" }), "list", _jsx(_components.span, { className: "token punctuation", children: "(" }), "solution ", _jsx(_components.span, { className: "token operator", children: "=" }), " round", _jsx(_components.span, { className: "token punctuation", children: "(" }), "sol", _jsx(_components.span, { className: "token punctuation", children: "," }), " ", _jsx(_components.span, { className: "token number", children: "3" }), _jsx(_components.span, { className: "token punctuation", children: ")" }), _jsx(_components.span, { className: "token punctuation", children: "," }), " predictions ", _jsx(_components.span, { className: "token operator", children: "=" }), " round", _jsx(_components.span, { className: "token punctuation", children: "(" }), "predictions", _jsx(_components.span, { className: "token punctuation", children: "," }), " ", _jsx(_components.span, { className: "token number", children: "3" }), _jsx(_components.span, { className: "token punctuation", children: ")" }), _jsx(_components.span, { className: "token punctuation", children: ")" }), _jsx(_components.span, { className: "token punctuation", children: ")" }), "\n"] }), _jsxs(_components.span, { className: "code-line", children: [_jsx(_components.span, { className: "token punctuation", children: "}" }), "\n"] }), _jsx(_components.span, { className: "code-line", children: "\n" }), _jsxs(_components.span, { className: "code-line", children: ["GBLUP ", _jsx(_components.span, { className: "token operator", children: "<-" }), " ", _jsx(_components.span, { className: "token keyword", children: "function" }), _jsx(_components.span, { className: "token punctuation", children: "(" }), "file", _jsx(_components.span, { className: "token punctuation", children: "," }), " snp.start", _jsx(_components.span, { className: "token punctuation", children: "," }), " snp.end", _jsx(_components.span, { className: "token punctuation", children: "," }), " trtNam", _jsx(_components.span, { className: "token punctuation", children: "," }), " \n"] }), _jsxs(_components.span, { className: "code-line", children: [" vare ", _jsx(_components.span, { className: "token operator", children: "=" }), " ", _jsx(_components.span, { className: "token number", children: "245.0" }), _jsx(_components.span, { className: "token punctuation", children: "," }), "vara ", _jsx(_components.span, { className: "token operator", children: "=" }), " ", _jsx(_components.span, { className: "token number", children: "35.24" }), _jsx(_components.span, { className: "token punctuation", children: ")" }), " ", _jsx(_components.span, { className: "token punctuation", children: "{" }), "\n"] }), _jsxs(_components.span, { className: "code-line", children: [" ", _jsx(_components.span, { className: "token comment", children: "# read-in the data" }), "\n"] }), _jsxs(_components.span, { className: "code-line", children: [" wkdat ", _jsx(_components.span, { className: "token operator", children: "<-" }), " read.csv", _jsx(_components.span, { className: "token punctuation", children: "(" }), "file", _jsx(_components.span, { className: "token operator", children: "=" }), "file", _jsx(_components.span, { className: "token punctuation", children: "," }), " header ", _jsx(_components.span, { className: "token operator", children: "=" }), " ", _jsx(_components.span, { className: "token boolean", children: "TRUE" }), _jsx(_components.span, { className: "token punctuation", children: ")" }), "\n"] }), _jsxs(_components.span, { className: "code-line", children: [" nInd ", _jsx(_components.span, { className: "token operator", children: "<-" }), " nrow", _jsx(_components.span, { className: "token punctuation", children: "(" }), "wkdat", _jsx(_components.span, { className: "token punctuation", children: ")" }), "\n"] }), _jsx(_components.span, { className: "code-line", children: " \n" }), _jsxs(_components.span, { className: "code-line", children: [" ", _jsx(_components.span, { className: "token comment", children: "# Daughter yield deviation for milk fat production" }), "\n"] }), _jsxs(_components.span, { className: "code-line", children: [" yidx", _jsx(_components.span, { className: "token operator", children: "<-" }), "grep", _jsx(_components.span, { className: "token punctuation", children: "(" }), "trtNam", _jsx(_components.span, { className: "token punctuation", children: "," }), "colnames", _jsx(_components.span, { className: "token punctuation", children: "(" }), "wkdat", _jsx(_components.span, { className: "token punctuation", children: ")" }), _jsx(_components.span, { className: "token punctuation", children: "," }), "ignore.case", _jsx(_components.span, { className: "token operator", children: "=" }), _jsx(_components.span, { className: "token boolean", children: "TRUE" }), _jsx(_components.span, { className: "token punctuation", children: ")" }), "\n"] }), _jsxs(_components.span, { className: "code-line", children: [" y ", _jsx(_components.span, { className: "token operator", children: "<-" }), " wkdat", _jsx(_components.span, { className: "token punctuation", children: "[" }), _jsx(_components.span, { className: "token punctuation", children: "," }), "yidx", _jsx(_components.span, { className: "token punctuation", children: "]" }), "\n"] }), _jsx(_components.span, { className: "code-line", children: " \n" }), _jsxs(_components.span, { className: "code-line", children: [" ", _jsx(_components.span, { className: "token comment", children: "# Genotypes and allelic frequencies " }), "\n"] }), _jsxs(_components.span, { className: "code-line", children: [" Z ", _jsx(_components.span, { className: "token operator", children: "<-" }), " as.matrix", _jsx(_components.span, { className: "token punctuation", children: "(" }), "wkdat", _jsx(_components.span, { className: "token punctuation", children: "[" }), _jsx(_components.span, { className: "token punctuation", children: "," }), " snp.start", _jsx(_components.span, { className: "token operator", children: ":" }), "snp.end", _jsx(_components.span, { className: "token punctuation", children: "]" }), _jsx(_components.span, { className: "token punctuation", children: ")" }), "\n"] }), _jsxs(_components.span, { className: "code-line", children: [" p ", _jsx(_components.span, { className: "token operator", children: "<-" }), " apply", _jsx(_components.span, { className: "token punctuation", children: "(" }), "wkdat", _jsx(_components.span, { className: "token punctuation", children: "[" }), _jsx(_components.span, { className: "token punctuation", children: "," }), " snp.start", _jsx(_components.span, { className: "token operator", children: ":" }), "snp.end", _jsx(_components.span, { className: "token punctuation", children: "]" }), _jsx(_components.span, { className: "token punctuation", children: "," }), " ", _jsx(_components.span, { className: "token number", children: "2" }), _jsx(_components.span, { className: "token punctuation", children: "," }), " mean", _jsx(_components.span, { className: "token punctuation", children: ")" }), " ", _jsx(_components.span, { className: "token operator", children: "/" }), " ", _jsx(_components.span, { className: "token number", children: "2" }), "\n"] }), _jsxs(_components.span, { className: "code-line", children: [" sum2pq ", _jsx(_components.span, { className: "token operator", children: "<-" }), " ", _jsx(_components.span, { className: "token number", children: "2" }), " ", _jsx(_components.span, { className: "token operator", children: "*" }), " sum", _jsx(_components.span, { className: "token punctuation", children: "(" }), "p ", _jsx(_components.span, { className: "token operator", children: "*" }), " ", _jsx(_components.span, { className: "token punctuation", children: "(" }), _jsx(_components.span, { className: "token number", children: "1" }), " ", _jsx(_components.span, { className: "token operator", children: "-" }), " p", _jsx(_components.span, { className: "token punctuation", children: ")" }), _jsx(_components.span, { className: "token punctuation", children: ")" }), "\n"] }), _jsxs(_components.span, { className: "code-line", children: [" nSNP ", _jsx(_components.span, { className: "token operator", children: "<-" }), " length", _jsx(_components.span, { className: "token punctuation", children: "(" }), "p", _jsx(_components.span, { className: "token punctuation", children: ")" }), "\n"] }), _jsx(_components.span, { className: "code-line", children: " \n" }), _jsxs(_components.span, { className: "code-line", children: [" lambda ", _jsx(_components.span, { className: "token operator", children: "<-" }), " vare ", _jsx(_components.span, { className: "token operator", children: "/" }), " vara\n"] }), _jsx(_components.span, { className: "code-line", children: " \n" }), _jsxs(_components.span, { className: "code-line", children: [" ", _jsx(_components.span, { className: "token comment", children: "# The overall mean as the only fixed effect" }), "\n"] }), _jsxs(_components.span, { className: "code-line", children: [" X ", _jsx(_components.span, { className: "token operator", children: "<-" }), " matrix", _jsx(_components.span, { className: "token punctuation", children: "(" }), _jsx(_components.span, { className: "token number", children: "1" }), _jsx(_components.span, { className: "token punctuation", children: "," }), " nrow ", _jsx(_components.span, { className: "token operator", children: "=" }), " nInd", _jsx(_components.span, { className: "token punctuation", children: "," }), " ncol ", _jsx(_components.span, { className: "token operator", children: "=" }), " ", _jsx(_components.span, { className: "token number", children: "1" }), _jsx(_components.span, { className: "token punctuation", children: ")" }), "\n"] }), _jsx(_components.span, { className: "code-line", children: " \n" }), _jsxs(_components.span, { className: "code-line", children: [" ", _jsx(_components.span, { className: "token comment", children: "# Center the SNP genotypes" }), "\n"] }), _jsxs(_components.span, { className: "code-line", children: [" ", _jsx(_components.span, { className: "token keyword", children: "for" }), " ", _jsx(_components.span, { className: "token punctuation", children: "(" }), "j ", _jsx(_components.span, { className: "token keyword", children: "in" }), " ", _jsx(_components.span, { className: "token number", children: "1" }), _jsx(_components.span, { className: "token operator", children: ":" }), "ncol", _jsx(_components.span, { className: "token punctuation", children: "(" }), "Z", _jsx(_components.span, { className: "token punctuation", children: ")" }), _jsx(_components.span, { className: "token punctuation", children: ")" }), " ", _jsx(_components.span, { className: "token punctuation", children: "{" }), "\n"] }), _jsxs(_components.span, { className: "code-line", children: [" Z", _jsx(_components.span, { className: "token punctuation", children: "[" }), _jsx(_components.span, { className: "token punctuation", children: "," }), " j", _jsx(_components.span, { className: "token punctuation", children: "]" }), " ", _jsx(_components.span, { className: "token operator", children: "<-" }), " Z", _jsx(_components.span, { className: "token punctuation", children: "[" }), _jsx(_components.span, { className: "token punctuation", children: "," }), " j", _jsx(_components.span, { className: "token punctuation", children: "]" }), " ", _jsx(_components.span, { className: "token operator", children: "-" }), " ", _jsx(_components.span, { className: "token number", children: "2" }), " ", _jsx(_components.span, { className: "token operator", children: "*" }), " p", _jsx(_components.span, { className: "token punctuation", children: "[" }), "j", _jsx(_components.span, { className: "token punctuation", children: "]" }), "\n"] }), _jsxs(_components.span, { className: "code-line", children: [" ", _jsx(_components.span, { className: "token punctuation", children: "}" }), "\n"] }), _jsx(_components.span, { className: "code-line", children: " \n" }), _jsxs(_components.span, { className: "code-line", children: [" W ", _jsx(_components.span, { className: "token operator", children: "<-" }), " diag", _jsx(_components.span, { className: "token punctuation", children: "(" }), _jsx(_components.span, { className: "token number", children: "1" }), _jsx(_components.span, { className: "token punctuation", children: "," }), " nrow ", _jsx(_components.span, { className: "token operator", children: "=" }), " nInd", _jsx(_components.span, { className: "token punctuation", children: "," }), " ncol ", _jsx(_components.span, { className: "token operator", children: "=" }), " nInd", _jsx(_components.span, { className: "token punctuation", children: ")" }), "\n"] }), _jsx(_components.span, { className: "code-line", children: " \n" }), _jsxs(_components.span, { className: "code-line", children: [" ", _jsx(_components.span, { className: "token comment", children: "# Set up MME" }), "\n"] }), _jsxs(_components.span, { className: "code-line", children: [" XX ", _jsx(_components.span, { className: "token operator", children: "<-" }), " t", _jsx(_components.span, { className: "token punctuation", children: "(" }), "X", _jsx(_components.span, { className: "token punctuation", children: ")" }), " ", _jsx(_components.span, { className: "token percent-operator operator", children: "%*%" }), " X\n"] }), _jsxs(_components.span, { className: "code-line", children: [" XW ", _jsx(_components.span, { className: "token operator", children: "<-" }), " t", _jsx(_components.span, { className: "token punctuation", children: "(" }), "X", _jsx(_components.span, { className: "token punctuation", children: ")" }), " ", _jsx(_components.span, { className: "token percent-operator operator", children: "%*%" }), " W\n"] }), _jsxs(_components.span, { className: "code-line", children: [" Xy ", _jsx(_components.span, { className: "token operator", children: "<-" }), " t", _jsx(_components.span, { className: "token punctuation", children: "(" }), "X", _jsx(_components.span, { className: "token punctuation", children: ")" }), " ", _jsx(_components.span, { className: "token percent-operator operator", children: "%*%" }), " y\n"] }), _jsxs(_components.span, { className: "code-line", children: [" Wy ", _jsx(_components.span, { className: "token operator", children: "<-" }), " t", _jsx(_components.span, { className: "token punctuation", children: "(" }), "W", _jsx(_components.span, { className: "token punctuation", children: ")" }), " ", _jsx(_components.span, { className: "token percent-operator operator", children: "%*%" }), " y\n"] }), _jsx(_components.span, { className: "code-line", children: " \n" }), _jsxs(_components.span, { className: "code-line", children: [" ", _jsx(_components.span, { className: "token comment", children: "# Add 0.01 to the diagonal elements " }), "\n"] }), _jsxs(_components.span, { className: "code-line", children: [" G ", _jsx(_components.span, { className: "token operator", children: "<-" }), " ", _jsx(_components.span, { className: "token punctuation", children: "(" }), "Z ", _jsx(_components.span, { className: "token percent-operator operator", children: "%*%" }), " t", _jsx(_components.span, { className: "token punctuation", children: "(" }), "Z", _jsx(_components.span, { className: "token punctuation", children: ")" }), _jsx(_components.span, { className: "token punctuation", children: ")" }), " ", _jsx(_components.span, { className: "token operator", children: "/" }), " sum2pq ", _jsx(_components.span, { className: "token operator", children: "+" }), " diag", _jsx(_components.span, { className: "token punctuation", children: "(" }), _jsx(_components.span, { className: "token number", children: "0.01" }), _jsx(_components.span, { className: "token punctuation", children: "," }), " nrow ", _jsx(_components.span, { className: "token operator", children: "=" }), " nInd", _jsx(_components.span, { className: "token punctuation", children: "," }), " ncol ", _jsx(_components.span, { className: "token operator", children: "=" }), " nInd", _jsx(_components.span, { className: "token punctuation", children: ")" }), "\n"] }), _jsxs(_components.span, { className: "code-line", children: [" IG ", _jsx(_components.span, { className: "token operator", children: "<-" }), " solve", _jsx(_components.span, { className: "token punctuation", children: "(" }), "G", _jsx(_components.span, { className: "token punctuation", children: ")" }), "\n"] }), _jsx(_components.span, { className: "code-line", children: " \n" }), _jsxs(_components.span, { className: "code-line", children: [" LHS1 ", _jsx(_components.span, { className: "token operator", children: "<-" }), " cbind", _jsx(_components.span, { className: "token punctuation", children: "(" }), "XX", _jsx(_components.span, { className: "token punctuation", children: "," }), " XW", _jsx(_components.span, { className: "token punctuation", children: ")" }), "\n"] }), _jsxs(_components.span, { className: "code-line", children: [" LHS2 ", _jsx(_components.span, { className: "token operator", children: "<-" }), " cbind", _jsx(_components.span, { className: "token punctuation", children: "(" }), "t", _jsx(_components.span, { className: "token punctuation", children: "(" }), "XW", _jsx(_components.span, { className: "token punctuation", children: ")" }), _jsx(_components.span, { className: "token punctuation", children: "," }), " ", _jsx(_components.span, { className: "token punctuation", children: "(" }), "diag", _jsx(_components.span, { className: "token punctuation", children: "(" }), _jsx(_components.span, { className: "token number", children: "1" }), _jsx(_components.span, { className: "token punctuation", children: "," }), " nrow ", _jsx(_components.span, { className: "token operator", children: "=" }), " nInd", _jsx(_components.span, { className: "token punctuation", children: "," }), " ncol ", _jsx(_components.span, { className: "token operator", children: "=" }), " nInd", _jsx(_components.span, { className: "token punctuation", children: ")" }), " ", _jsx(_components.span, { className: "token operator", children: "+" }), " lambda ", _jsx(_components.span, { className: "token operator", children: "*" }), " IG", _jsx(_components.span, { className: "token punctuation", children: ")" }), _jsx(_components.span, { className: "token punctuation", children: ")" }), "\n"] }), _jsxs(_components.span, { className: "code-line", children: [" LHS ", _jsx(_components.span, { className: "token operator", children: "<-" }), " rbind", _jsx(_components.span, { className: "token punctuation", children: "(" }), "LHS1", _jsx(_components.span, { className: "token punctuation", children: "," }), " LHS2", _jsx(_components.span, { className: "token punctuation", children: ")" }), "\n"] }), _jsx(_components.span, { className: "code-line", children: " \n" }), _jsxs(_components.span, { className: "code-line", children: [" RHS ", _jsx(_components.span, { className: "token operator", children: "<-" }), " c", _jsx(_components.span, { className: "token punctuation", children: "(" }), "Xy", _jsx(_components.span, { className: "token punctuation", children: "[" }), _jsx(_components.span, { className: "token number", children: "1" }), _jsx(_components.span, { className: "token punctuation", children: "," }), " ", _jsx(_components.span, { className: "token number", children: "1" }), _jsx(_components.span, { className: "token punctuation", children: "]" }), _jsx(_components.span, { className: "token punctuation", children: "," }), " Wy", _jsx(_components.span, { className: "token punctuation", children: ")" }), "\n"] }), _jsx(_components.span, { className: "code-line", children: " \n" }), _jsxs(_components.span, { className: "code-line", children: [" sol ", _jsx(_components.span, { className: "token operator", children: "<-" }), " solve", _jsx(_components.span, { className: "token punctuation", children: "(" }), "LHS", _jsx(_components.span, { className: "token punctuation", children: ")" }), " ", _jsx(_components.span, { className: "token percent-operator operator", children: "%*%" }), " RHS\n"] }), _jsx(_components.span, { className: "code-line", children: " \n" }), _jsxs(_components.span, { className: "code-line", children: [" return", _jsx(_components.span, { className: "token punctuation", children: "(" }), "round", _jsx(_components.span, { className: "token punctuation", children: "(" }), "sol", _jsx(_components.span, { className: "token punctuation", children: "," }), " ", _jsx(_components.span, { className: "token number", children: "3" }), _jsx(_components.span, { className: "token punctuation", children: ")" }), _jsx(_components.span, { className: "token punctuation", children: ")" }), "\n"] }), _jsxs(_components.span, { className: "code-line", children: [_jsx(_components.span, { className: "token punctuation", children: "}" }), "\n"] }), _jsx(_components.span, { className: "code-line", children: "\n" })] }) }); } function MDXContent(props = {}) { const {wrapper: MDXLayout} = { ..._provideComponents(), ...props.components }; return MDXLayout ? _jsx(MDXLayout, { ...props, children: _jsx(_createMdxContent, { ...props }) }) : _createMdxContent(props); } return { default: MDXContent };