powered by
Converts a kinship matrix to a Genetic Relation Matrix (GRM)
kinship2grm(x)
a square matrix containing kinship coefficients
a GRM matrix
for off-diagonal elements, \(GRM=2 \times x_{ij}\); for diagonal elements, \(GRM=1+ x_{ii}\)
# NOT RUN { dos<-matrix(sample(0:2,replace=TRUE,size=1000),nrow=10) #dosage matrix for 10 inds at 100 loci ks<-beta.dosage(dos) # kinship matrix kinship2grm(ks) # }
Run the code above in your browser using DataLab