millironx/beefblup
1
0
Fork 0
mirror of https://github.com/MillironX/beefblup.git synced 2024-12-22 09:08:16 +00:00
Code Issues Releases Packages Activity

Removed null observations from incidence matrix

Browse source
This commit is contained in:
Thomas A. Christensen II 2018-09-28 09:12:24 -06:00
parent bd44916781
commit dfddabde76
1 changed files with 4 additions and 0 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

4
MATLAB/beefblup.m
View file

@ -187,6 +187,10 @@ Y = cell2mat(data(:, 5));
% The identity matrix for random effects % The identity matrix for random effects
Z = eye(numanimals, numanimals); Z = eye(numanimals, numanimals);
% Remove items where there is no data
nullobs = find(isnan(Y));
Z(nullobs, nullobs) = 0;
% Prompt for heritability % Prompt for heritability
h2 = input('What is the heritability for this trait? >> '); h2 = input('What is the heritability for this trait? >> ');
lambda = (1-h2)/h2; lambda = (1-h2)/h2;