beefpack - khon kaen university · interface which allows you to create blup report within 5 steps....
Post on 09-Feb-2021
7 Views
Preview:
TRANSCRIPT
-
BeefPack Version 2.5, 2006
Monchai Duangjinda Department of Animal Science
Khon Kaen University, Khon Kaen, Thailand 40002
Ignacy Misztal Department of Animal and Dairy Science
The University of Georgia, Athens, GA30602
Shogo Tsuruta Department of Animal and Dairy Science
The University of Georgia, Athens, GA30602
-
2
© 2006
Department of Animal Science Khon Kaen University
& Department of Animal and Dairy Science
The University of Georgia
-
3
Contents
1. Introduction 4
A. Background B. Specification C. Conditions of use D. Program download E. Online registration
2. Installation 6 A. Installation B. How to start
3. Genetic Evaluation At A Glance 8 A. General view B. Report view C. Advanced options
4. Model Descriptions 11 A. Basic animal model B. Animal with PE model C. Maternal Model D. Multi-trait model
5. How Program Functions? 14 A. Program and file components B. The way program work C. Solutions at convergence
6. Genetic Evaluation by Examples 16 A. Setup directory and preferences B. Data and pedigree files C. Single trait analysis example D. Maternal with PE model example E. Multi-trait analysis example
7. Rapid Analysis from Previous Parameters 20 A. Re-analysis by old parameters B. Creating BLUP report from previous solutions
8. Parameter file example 21 A. Repeatability model B. Multi-trait model C. Random regression test-day model
9. Trouble shooting 24 10. Who’s who? 25
-
4
I. Introduction
A. Background
“BLUPF90” and related programs were developed in the lab of Ignacy Misztal with
the purpose of providing comprehensive computing capabilities to problems related to mixed models in animal breeding. See http://nce.ads.uga.edu/~ignacy for details and documentation. These programs are mostly written in Fortran 90 and have a line-mode interface. “BeefPack” is a set of programs branched from “PCPack” with graphical user interface for simply used in Windows. BeefPack was developed with specific purpose for beef genetic evaluation, therefore parts of programs in the BLUPF90 family are available, which are: BLUPF90 : BLUP estimation using PCG. REMLF90 : Variance estimation using REML by EM algorithm. AIREMLF90 : Variance estimation using REML by AI Algorithm.
RENUMMAT : Renumber program for creating data file and additive pedigree with animal ID in order number.
ACCF90 : Approximate accuracy for BLUP solutions for direct and maternal models
All programs were compiled separately under Microsoft Windows using Visual Fortran version 5.1. Users can run all programs separately using MS-DOS prompt or from the menu in BeefPack. It allows users to create all parameter files and creating BLUP reports with a point-and-click interface written in Visual BASIC application for MS-EXCEL.
B. Overall Features
BeefPack performs under EXCEL. It manages genetic evaluation using wizard interface which allows you to create BLUP report within 5 steps. Only original data and pedigree are requires. It will automatically renumber the animals and count the number of effects in the model. After variance estimation or BLUP analysis were performed, it will create BV report with the original animal ID. Users can keep data files and program files separately, however specified directories are required before the analysis. BV with accuracy and genetic trend report for direct genetics up to 4 traits and for maternal genetics up to 2 traits are available.
B. Specification
BLUPF90-BeefPack 2 requires windows 95/98/ME/XP/2000 environment to install and Excel 98/2000/XP/2002/2003 for running applications. It also requires at least 32 MB for memory and 5 MB of disk space for storing programs.
For BLUP and variance estimation, BLUPF90 and REMLF90 support general single and multiple trait models, i.e. sire-maternal grandsire, animal model, maternal model, repeatability model, and dominance model, with missing values and different models for each trait. Random regression models are also supported. AIREMLF90 may not support some models and some particular structure of data. ACCF90 is an approximation that works with repeatability and maternal models.
With BeefPack wizard interface, however, the single trait will support only simple animal, animal with PE model and direct-maternal with PE model. Multiple traits will support up to 4 traits for direct genetics and up to 2 traits with maternal genetics. Variance estimation and BV report with accuracy are available for both single and multiple traits.
http://nce.ads.uga.edu/~ignacy
-
5
C. Conditions of Use BLUPF90-BeefPack is distributed free of charge for academic and scientific use
under the conditions that it remains copyright. The use of any applications and compiled programs from BeefPack needs to be credited in any publications. For commercial or grant project, personal communication for further agreements is required with any of the authors. There is no guarantee for its correctness and there is no-service for user purpose. However, specific questions, criticism and bug reports are invited. Please email to monchai@kku.ac.th
D. Program Download BLUPF90-BeefPack has been made available on the CD. However, the updated
version is available at BLUPF90 homepage, http://agserver.kku.ac.th/monchai/BlupF90. The complete package provides program files, manuals and examples.
E. Online Registration
Online registration is requested for further breeding and genetic group connection and
update version information. The registration page for BLUPF90-BeefPack is also available at BLUPF90 homepage, http://agserver.kku.ac.th/monchai/BlupF90.
http://agserver.kku.ac.th/monchai/BlupF90http://agserver.kku.ac.th/monchai/BlupF90
-
6
II. Installation A. Installation
BLUPF90-BeefPack is stored in one installation file named “BeefPack.Exe.” To
install the programs, do the following steps:
1. Run BeefPack.EXE.
2. Setup Wizards will tell you through the
installation steps.
3. BeefPack will install all F90 programs to the
default directory “C:\BLUPF90-BP”. To modify the directory name, click BROWSE button.
4. Click FINISH button to complete the installation.
B. How to Start DOS environment
Each program can be run directly from directory \BLUPF90. The following are executable programs that can be called at DOS prompt. For example, to run BLUPF90 program, type the following at the command prompt: C:\BLUPF90-BP\BLUPF90
The programs that can be called from the dos prompt are: BLUPF90 : BLUP estimation using PCG. REMLF90 : Variance estimation using REML by EM algorithm. AIREMLF90 : Variance estimation using REML by AI Algorithm. RENUMMAT : Renumber program for creating data file and additive pedigree with
animal ID in order number. ACCF90 : Approximate accuracy for BLUP solutions for direct and maternal models
-
7
Note: All programs require specific parameter files, which need to be created before calling the
programs. For parameter file examples, see section 6.
Windows environment
1. At the windows task bar, click on Start
Menu > Programs > BLUPF90 > BeefPack
2. Click BlupF90 BeefPack icon.
3. Excel program will and BeefPack will
ask for macro enables. Click Enable Macros button. (Note: BLUPF90-BeefPack cannot execute without macros for Visual Basic Applications. If there is a problem, enable macros by clicking menu Tools > Macro > Security in Excel. Then, set security level to Medium)
4. Click OK button at splash dialog to start
the program.
5. The main menu should look as below:
-
8
III. Genetic Evaluation At A Glance
A. General View
1. BeefPack is a set of programs in BLUPF90
family with a specific purpose for dairy cattle evaluation.
2. BeefPack performs variance estimation using
REML and BLUP analysis with wizard interface.
3. BLUP analysis can be done directly with true
variance or variance ratio. In addition, re-estimate variance components with REML or AIREML can also apply.
4. If multiple trait model is interested, Direct and
maternal genetic effects can be evaluated.
5. Original data and pedigree files are used in the
analysis. Therefore, all effects in the model are simply specified by column number in data file. For multi-trait model, different contemporary groups for each trait can be specified. For maternal direct or PE model, dam column is not required in data file (program will automatically creat from pedigreed).
-
9
B. Report View
6. BLUP report with accuracy is created with
original ID in Excel format. Therefore, sorting, filtering can simply perform using Excel functions.
7. If REML or AIREML is performed. New
variance estimates are kept in separate file.
8. Genetic trend report is also created if require. All
graphic properties can be modified using general Excel features.
-
10
C. Advanced Options
9. Previous renumbered data and pedigree file with
parameter file can be modified and reanalysis.
10. BV reports for single trait and multiple trait from
previous solutions from BeefPack or another BLUP family can be recreated with specified options.
11. Multiple trait BV reports for
direct and maternal can be created up to 4 traits.
12. Report with accuracy can also
be created.
-
11
IV. Model Descriptions
The main objective of BLUPF90-BeefPack is to accomplish various model generally
used in beef cattle evaluation with friendly graphic interface for PC and Windows users. Using powerful features from programs of BLUPF90, BeefPack can accomplish a wide range of genetic evaluation. BeefPack can estimate variance components with REML and perform BLUP analysis from linear mixed model which includes random effect of animal as additive genetic or permanent environment effect. The following will describe more details of model that can be used in the analysis.
A. Basic Animal Model
BeefPack provide basic animal model which allows animals in the data and animals in the pedigree to be included in the analysis so that all known relationships can be taken into account. Other effects, fixed and random, can be including for comprehensive use of mixed model. Fixed effects used in the model can be fitted as cross-classified variable and covariate. Combination of fixed effects such as herd-year-season can be performed during the analysis, therefore, no adding step of data preparation is required. Normally, all traits used in BeefPack should be continuous rather than ordinal scale for proper use in the analysis of linear mixed model. Model with single record per animal such as daily gain yield can be analyzed by basic animal model as follow:
εZaXβy ++= , and , ⎥⎥⎦
⎤
⎢⎢⎣
⎡=⎥
⎦
⎤⎢⎣
⎡2e
2aV
σσ
I00A
εa
where y is vector of response variable, β is vector of fixed effects, a is vector of
random additive genetic effects, ε is vector of random residual, X and Z are incident matrices related to fixed and random effects, A is numerator relationship matrix, is
additive genetic variance, and is residual variance.
2aσ
2eσ
To perform BLUP, Henderson’s MME can be written as:
⎥⎦
⎤⎢⎣
⎡′′
=⎥⎦
⎤⎢⎣
⎡⎥⎦
⎤⎢⎣
⎡+′′
′′− yZ
yXaβ
AZZXZZXXX
1α, where
2a
2e
σσ
α =
B. Animal with PE Model
If multiple parity records are available, the maternal permanent environment effect due to dam need to be taken into account. Fitting maternal permanent environment effect as uncorrelated random effects is generally used in genetic evaluation. The model for analysis is:
εWcZaXβy +++= , and , ⎥⎥⎥
⎦
⎤
⎢⎢⎢
⎣
⎡
=⎥⎥⎥
⎦
⎤
⎢⎢⎢
⎣
⎡
2e
2c
2a
Vσ
σσ
I000I000A
εca
where y is vector of response variable, β is vector of fixed effects, a is vector of
random additive genetic effects, c is vector of random permanent environment effects, ε is vector of random residual, X , and W Z are incident matrices related to fixed and random effects, A is numerator relationship matrix, is additive genetic variance, is additive
genetic variance, and is residual variance.
2aσ
2cσ
2eσ
-
12
To perform BLUP, Henderson’s MME can be written as:
⎥⎥⎥
⎦
⎤
⎢⎢⎢
⎣
⎡
′′′
=⎥⎥⎥
⎦
⎤
⎢⎢⎢
⎣
⎡
⎥⎥⎥
⎦
⎤
⎢⎢⎢
⎣
⎡
+′′′′+′′′′′
−
yWyZyX
caβ
IWWZWXWWZAZZXZWXZXXX
1
γα , where
2a
2e
σσ
α = , 2c
2e
σσ
γ =
C. Maternal Model
If multiple parity records are available, also such trait could be effected by maternal genetic, therefore the maternal permanent environment effect and maternal genetic due to dam need to be taken into account. Fitting maternal permanent environment effect as uncorrelated random effects and maternal genetic effect as correlated random effects are generally used in genetic evaluation. The model for analysis is:
εWcmZaZXβy 21 ++++= , and ,
⎥⎥⎥⎥⎥
⎦
⎤
⎢⎢⎢⎢⎢
⎣
⎡
=
⎥⎥⎥⎥
⎦
⎤
⎢⎢⎢⎢
⎣
⎡
2
2
2
2
e
c
mam
ama
V
σσ
σσσσ
I0000I0000AA00AA
εcma
where y is vector of response variable, β is vector of fixed effects, a is vector of
random additive direct genetic effects, is vector of random additive maternal genetic effects, is vector of random permanent environment effects,
mc ε is vector of random residual,
X , , and are incident matrices related to fixed and random effects, W 1Z 2Z A is
numerator relationship matrix, is additive direct genetic variance, is additive maternal
genetic variance, , is direct-maternal genetic covariance, is maternal permanent
environmental variance, and is residual variance.
2aσ
2mσ
amσ2cσ
2eσ
To perform BLUP, Henderson’s MME can be written as:
⎥⎥⎥⎥
⎦
⎤
⎢⎢⎢⎢
⎣
⎡
′′′′
=
⎥⎥⎥⎥
⎦
⎤
⎢⎢⎢⎢
⎣
⎡
⎥⎥⎥⎥
⎦
⎤
⎢⎢⎢⎢
⎣
⎡
+′′′′′+′+′′′+′+′′′′′′
−−
−−
yWyZyZyX
cmaβ
IWWZWZWXWWZAZZAZZXZWZAZZAZZXZWXZXZXXX
2
1
1
21
221
122
11
211
111
21
γαααα
2
2221
1211
,
where , 1
2
22
2221
1211 −
⎥⎥⎦
⎤
⎢⎢⎣
⎡=
⎥⎥⎦
⎤
⎢⎢⎣
⎡
mam
amae σσ
σσσαααα
2c
2e
σσ
γ =
D. Multi-trait Model
BeefPack can also perform multi-trait analysis. Estimation of genetic correlations among traits and multivariate BLUP analysis can be accomplished. However, for graphic user interface, not greater than 4 traits is available. To perform beyond this, parameter editing is required and do the analysis from menu BLUP>Use old parameters. Multi-trait analysis can perform similar model for all traits, missing records for some traits, different model for each trait, also different contemporary group for each trait. For example, the following traits for birth weight, weaning weight and yearling weight can be analyzed together.
-
13
⎥⎥⎥
⎦
⎤
⎢⎢⎢
⎣
⎡+
⎥⎥⎥
⎦
⎤
⎢⎢⎢
⎣
⎡
⎥⎥⎥
⎦
⎤
⎢⎢⎢
⎣
⎡+
⎥⎥⎥
⎦
⎤
⎢⎢⎢
⎣
⎡
⎥⎥⎥
⎦
⎤
⎢⎢⎢
⎣
⎡+
⎥⎥⎥
⎦
⎤
⎢⎢⎢
⎣
⎡
⎥⎥⎥
⎦
⎤
⎢⎢⎢
⎣
⎡+
⎥⎥⎥
⎦
⎤
⎢⎢⎢
⎣
⎡
⎥⎥⎥
⎦
⎤
⎢⎢⎢
⎣
⎡=
⎥⎥⎥
⎦
⎤
⎢⎢⎢
⎣
⎡
3
2
1
2
1
2
1
2
1
2
1
3
2
1
3
2
1
3
2
1
3
2
1
3
2
1
εεε
0cc
0000W000W
0mm
0000M000M
aaa
Z000Z000Z
βββ
X000X000X
yyy
Where , AGma
⊗=⎥⎦
⎤⎢⎣
⎡V [ ] IWc ⊗=V , [ ] IRε ⊗=V
where is vector of response variable for trait 1, 2 and 3; is
vector of fixed effects; is vector of random additive direct genetic effects; is vector of random additive maternal genetic effects; is vector of random maternal permanent environmental effects; is vector of random residual,
321 yy,y , 321 β,β,β
321 a,a,a 21 m,m
21 c,c
321 ε,ε,ε X , , W and M Z are incident matrices related to fixed and random effects, A is numerator relationship matrix, G is matrix of direct-maternal genetic variance-covariance for trait 1, 2 and 3; is matrix of maternal permanent environmental variance-covariance for trait 1 and 2;
WR is matrix of
residual variance-covariance for trait 1,2 and 2.
-
14
V. How Program Functions?
A. Program and file components After installation, programs are stored in main or user-specified directory, ie.
C:\BLUPF90-BP. Two sub-directories of examples and helps are also created. Each directory will find the following programs and files.
Main directory:
Name Type DescriptionBeefPack.xls XLS - Main graphic user interface.
- Creating BLUP and REML parameter files. - Creating BLUP and VCE report in excel format. - Creating genetic trend. - Computing BV and h2 for test day model.
BLUPF90.EXE PROG - Computing BLUP solutions. REMLF90.EXE PROG - Estimating variance components using REML with EM
algorithm. AIREMLF90.EXE PROG - Estimating variance components using REML with AI
algorithm. ACCF90.EXE PROG - Computing approximate accuracy for BLUP solutions. RENUMMAT.EXE PROG - Renumbering animal in data and pedigree file in
consecutive order number. Examples directory:
Name Type DescriptionBFDAT.PRN TXT - Data file for analysis with single trait and multi-trait
with maternal effects and repeated records. BFPED.PRN TXT - Pedigree file for analysis with BFDAT. BFDAT.FMT TXT - Describe column number format for BFDAT.PRN
Helps directory:
Name Type DescriptionWHOSWHO.TXT TXT - Accredit for key persons involved in BLUPF90 family. Manual-BP.PDF PDF - Manual for BLUPF90 BeefPack
B. The way program works
When performed the analysis with wizard interface in BeefPack. All parameters entered in the form will be kept in particular Excel sheets. They will be written with corrected format as parameter file for RENUM and BLUPF90 using visual basic. Batch file to call the program with the parameter is also need to be created. VB in Excel has specific function to operate EXE file in this batch without closing the Excel program. Pedigree and solutions from the analysis are read to Excel sheets to join back the original animal ID, and also genetic trend will be created using chart function if user required.
Data and Pedigree files used in BeefPack must be ASCII or TEXT file. All data must
be in number, except for animal ID, Sire, and Dam that can be alpha-numeric format. If create from Excel, Save as PRN file (Text file delimited with space) is preferable than Tab delimited or comma delimited.
-
15
If analyzing data is kept in different directory of programs. All execute programs will be copied to the data directory. After analysis, a few file will be created, which can be copied to new name if need.
Suppose the original files for the analysis are BFDAT.PRN and BFPED.PRN, some additional files after analysis, which might be useful for later analysis, are:
Name Type DescriptionREBFDAT.PRN TXT - Renumbered data file. REBFPED.PRN TXT - Renumbered Pedigree file. RENUM.PAR TXT - Parameter file for program RENUMMAT. RENUM.MSG TXT - Log file from renumbering. Describe the levels of fixed
and random effects after renum. RENUM.PRN TXT - The details from renumbering. Describe how effects
are combined and replications of each effects. BLUP.PAR TXT - Parameter files for REML and BLUP analysis SOLUTIONS TXT - BLUP solutions file. SOLUTIONS_VCE TXT - This file keeps variance estimation if performed.
B. Solutions at Convergence
The default for convergence is 1d-08, however, users can choose their own by selecting menu Options > Set Options.
-
16
VI. Genetic Evaluation by Examples
A. Setup Directory and Preferences
1. The correct directory for program files and data
files need to be checked or set up before starting the analysis. Select menu Directory >Setup Directory.
2. Other options like convergence levels and
accuracy in the reports can be chosen. Select menu Options > Set Options.
B. Data and Pedigree files
Single trait analysis The example will show how to predict breeding value for milk production adjusted for 305 day. BFDAT.PRN is data file and BFPED.PRN is pedigree file. Multi-trait analysis The example will show how to estimate genetic correlation and predict multivariate breeding value for milk305, fat and protein yield. BFDAT.PRN is data file and BFPED.PRN is pedigree file. File format Format for BFDAT.PRN is stored in BFMDAT.FMT. There are 11 columns of:
#1 ID #2 Sex #3 Herd #4 Parity #5 Age of dam (month) #6 Month at birth #7 Year at birth #8 Birth weight #9 Month at weaning #10 Year at weaning #11 Weaning weight
Format for pedigree is the same. There are 4 columns of: #1 Animal ID #2 Sire ID #3 Dam ID #4 Year of birth.
Generally, year of birth can be omit from pedigree file, however genetic trend cannot
to be created. For maternal direct or PE model, dam column is not required in data file (program will automatically creat from pedigreed).
Note: all files are stored in C:\BLUPF90-BP\EXAMPLES.
-
17
Start
Finish
VCE report
C. Single Trait Analysis Example
Start Select menu BLUP > Single trait. BeefPack wizard will show up, click Next button to start analysis. 6. Input data file (BFdat.prn) and pedigree file
(BFped.prn). Note: User can click View button to check if there are correct files in the directory. Click at Path button to change path name if need.
7. Input either true variance or variance ratio for effects in analysis model. In this example, “0.4”, “0.1” and “0.5” are entered for heritability, maternal permanent environment and error ratio for birth weight. Note: To perform BLUP only, do not check on Re-estimate of variance check box. If this box is checked, all variances entered will be used as prior values for REML, and BLUP solutions will created with final variance estimates.
8. Input column number for fixed effects, animal and traits. In this examples, “3 6 7” are input to combine herd-month-year at birth as contemporary group, “2” is column for sex, “4” is column for parity, animal id and birth weight are in column “1” and “8”, respectively.
9. Enter filename to keep BLUP solutions and create BV report.
10. Wizard will ask for creating genetic trend. Click No button if only BLUP report is needed. Click Yes with another filename to keep genetic trend.
Finish Click View variance estimation button to view variance components estimated by REML (if Re-estimates in is checked). Click View EBV / Genetic Trend button to see reports as bellows.
EBV report
Genetic trend
-
18
Start
Finish
VCE report
D. Maternal with PE Model Example
Start Select menu BLUP > Single trait. BeefPack wizard will show up, click Next button to start analysis. 1. Input data file (BFdat.prn) and pedigree file
(BFped.prn). Note: User can click View button to check if there are correct files in the directory. Click at Path button to change path name if need.
2. Input either true variance or variance ratio for effects in analysis model. In this example, “0.30”, “0.15” and “-0.05” are entered for direct, maternal and covariance ratio; “0.10” and “0.50” for maternal permanent environment and error ratio for weaning weight. Note: To perform BLUP only, do not check on Re-estimate of variance check box. If this box is checked, all variances entered will be used as prior values for REML, and BLUP solutions will created with final variance estimates.
3. Input column number for fixed effects, animal and traits. In this examples, “3 9 10” are input to combine herd-month-year at birth as contemporary group, “2” is column for sex, “4” is column for parity, animal id and birth weight are in column “1” and “11”, respectively.
4. Enter filename to keep BLUP solutions and create BV report.
5. Wizard will ask for creating genetic trend. Click No button if only BLUP report is needed. Click Yes with another filename to keep genetic trend.
Finish Click View variance estimation button to view variance components estimated by REML (if Re-estimates in is checked). Click View EBV / Genetic Trend button to see reports as bellows.
EBV report
Genetic trend
-
19
E. Multi-trait Analysis Example
VCE report
Finish
EBV report
Start Select menu BLUP > Multiple trait. BeefPack wizard will show up, click Next button to start analysis. 1. Input data file (LMdat.prn) and pedigree file
(LMped.prn). 2. Input either true variance or variance ratio for
effects in analysis model. In this example, variance-covariances are entered for additive direct, permanent environment and error for trait birth, fat and protein yield. Note: To perform BLUP only, do not check on Re-estimate of variance check box. If this box is checked, all variances entered will be used as prior values for REML, and BLUP solutions will created with final variance estimates.
3. Input column number for fixed effects, animal and traits. In this examples, “3 6 7” and “3 9 10”are input to combine herd-month-year at birth and herd-month-year at wean as contemporary group for each trait, “2” is column for sex, “4” is column for parity, animal id and birth weight are in column “1” and “11”, respectively . Birth and weaning weight are in column “8” and “11”, respectively.
4. Enter filename to keep BLUP solutions and create BV report.
5. Wizard will ask for creating genetic trend. Click No button if only BLUP report is needed. Click Yes with another filename to keep genetic trend.
Finish Click View variance estimation button to view variance components estimated by REML (if Re-estimates in is checked). BeefPack also reports genetic, environment and phenotypic correlations.
Click View EBV / Genetic Trend button to see multivariate BLUP reports with accuracy as bellows.
Start
-
20
VII. Rapid Analysis from Previous Parameters A. Re-analysis By Old Parameters
Without using wizards, which always start from the first step, BeefPack has a useful feature to help user perform BLUP from previous parameter files. The parameters can be modified for complicated model or other models not available in wizards.
B. Creating BLUP Report From Previous Solutions All solutions from the analysis of BLUPF90 family programs. Such as BLUPF90,
REMLF90/AIREMLF90, BLUPF90-PCPACK, can be brought into BeefPack to create BV reports. This step requires solutions file and renum-pedigree file from previous analysis.
1. Select menu BLUP > Use old parameters. Enter parameter filename. The default from the latest analysis is blup.par.
2. User can click at View/Edit button to check or modify the parameter before the analysis.
3. If Re-estimates variance is “YES”, option for REML or AIREML needs to be selected.
4. Click “OK” button at Start Computing section. 5. Click View button to see BLUP solutions and
Estimated variance (if required).
1. Select menu BLUP > Create BV reports. Enter solution and renum-
pedigree filenames. 2. User can click at View button to check
the correct file name. Click Path to change path.
3. If solutions are from multi-trait BLUP
analysis, specify number of trait. 4. Enter number of animal effect or
maternal effect in solutions file.
5. Enter BV filename to save. Click “OK”
button at Start Computing section.
6. Click View button to see BLUP report.
-
21
VIII. Parameter File Examples A. Animal with PE Model
Data File: BFdat.prn Id sex hd par age mb yb bwt mw yw wwt 39001CP 2 89 2 59 2 1996 27 9 1996 200 39001LO 1 94 3 64 2 1996 28 9 1996 214 39001TP 1 141 4 95 2 1996 29 10 1996 209 39001UB 2 151 2 58 2 1996 29 9 1996 0 39001UD 2 178 2 65 2 1996 0 9 1996 228 39001YL 1 344 1 43 1 1996 26 7 1996 211 39002KA 2 90 2 59 2 1996 20 8 1996 0 39002LO 1 94 2 60 2 1996 27 9 1996 199 39002MK 2 96 2 78 3 1996 22 10 1996 209 39002TP 1 141 8 137 2 1996 27 10 1996 208 Pedigree File: BFped.prn
Id sire dam yob 39001CP 31131MK 34075MK 1996 39001KA 35003SS 32287AWD 1996 39001LO 35010MK 33059LO 1996 39001MK 28831AI 27134MK 1996 39001TP 31397UDT 31542UTV 1996 39001UB 35014TP 34162UMA 1996 39001UD 29044MK 33053TP 1996 39001UT 35007UB 34710ULO 1996 39001YL 35006LO 35021USB 1996
Renum Data File: ReBFdat.prn Cg sex par id bwt dam
5 2 2 1 27 2283 5 1 3 2 28 2284 5 1 4 3 29 2285 5 2 2 4 29 2286 5 2 2 5 0 2287 1 1 1 6 26 2288 5 2 2 7 20 2289 5 1 2 8 27 2290 9 2 2 9 22 2291 5 1 8 10 27 2292 Renum Pedigree File: ReBFped.prn
Id sire dam x YOB x x x orig ID 1 3498 2283 1 1996 2 1 0 39001CP 2 3499 2284 1 1996 2 1 0 39001LO 3 3500 2285 1 1996 2 1 0 39001TP 4 3501 2286 1 1996 2 1 0 39001UB 5 3502 2287 1 1996 2 1 0 39001UD 6 3503 2288 1 1996 2 1 0 39001YL 7 3504 2289 1 1996 2 1 0 39002KA 8 3499 2290 1 1996 2 1 0 39002LO 9 3505 2291 1 1996 2 1 0 39002MK 10 3506 2292 1 1996 2 1 0 39002TP
Note: *Cg is a combination of herd-month-year at birth
Parameter File: BLUP.PAR
DATAFILE ReBFDat.prn NUMBER_OF_TRAITS 1 NUMBER_OF_EFFECTS 5 OBSERVATION(S) 5 WEIGHT(S) EFFECTS: POSITIONS_IN_DATAFILE NUMBER_OF_LEVELS TYPE_OF_EFFECTS 1 48 cross 2 2 cross 3 11 cross 4 3768 cross 6 3768 cross RANDOM_RESIDUAL VALUES 0.5 RANDOM_GROUP 4 RANDOM_TYPE add_animal FILE ReBFped.prn (CO)VARIANCES 0.4 RANDOM_GROUP 5 RANDOM_TYPE diagonal FILE (CO)VARIANCES 0.1
Note: Effect of contemporary groups in column #1 has 48 levels.
Effect of sex in column #2 has 2 levels. Effect of parity in column #3 has 11 levels. Effect of direct genetic from animal in column #4 has 3768 levels.
Effect of maternal PE from dam in column #6 has 3768 levels.
Effect number 4 is random effect with add_animal type. This is genetic effects which require pedigree file.
Effect number 5 is random effect with diagonal type for PE, which is uncorrelated.
Analysis of single trait. Trait is in column #5 which is birth weight.
-
22
B. Matenal with PE Model
Data File: BFdat.prn Id sex hd par age mb yb bwt mw yw wwt 39001CP 2 89 2 59 2 1996 27 9 1996 200 39001LO 1 94 3 64 2 1996 28 9 1996 214 39001TP 1 141 4 95 2 1996 29 10 1996 209 39001UB 2 151 2 58 2 1996 29 9 1996 0 39001UD 2 178 2 65 2 1996 0 9 1996 228 39001YL 1 344 1 43 1 1996 26 7 1996 211 39002KA 2 90 2 59 2 1996 20 8 1996 0 39002LO 1 94 2 60 2 1996 27 9 1996 199 39002MK 2 96 2 78 3 1996 22 10 1996 209 39002TP 1 141 8 137 2 1996 27 10 1996 208 Pedigree File: BFped.prn
Id sire dam yob 39001CP 31131MK 34075MK 1996 39001KA 35003SS 32287AWD 1996 39001LO 35010MK 33059LO 1996 39001MK 28831AI 27134MK 1996 39001TP 31397UDT 31542UTV 1996 39001UB 35014TP 34162UMA 1996 39001UD 29044MK 33053TP 1996 39001UT 35007UB 34710ULO 1996 39001YL 35006LO 35021USB 1996
Renum Data File: ReBFdat.prn Cg sex par id wwt dam
7 2 2 1 200 2283 47 1 3 2 214 2284 129 1 4 3 209 2285 151 2 2 4 0 2286 177 2 2 5 228 2287 220 1 1 6 211 2288 27 2 2 7 0 2289 47 1 2 8 199 2290 97 2 2 9 209 2291 129 1 8 10 208 2292 Renum Pedigree File: ReBFped.prn
Id sire dam x YOB x x x orig ID 1 3498 2283 1 1996 2 1 0 39001CP 2 3499 2284 1 1996 2 1 0 39001LO 3 3500 2285 1 1996 2 1 0 39001TP 4 3501 2286 1 1996 2 1 0 39001UB 5 3502 2287 1 1996 2 1 0 39001UD 6 3503 2288 1 1996 2 1 0 39001YL 7 3504 2289 1 1996 2 1 0 39002KA 8 3499 2290 1 1996 2 1 0 39002LO 9 3505 2291 1 1996 2 1 0 39002MK 10 3506 2292 1 1996 2 1 0 39002TP
Note: *Cg is a combination of herd-month-year at wean
Parameter File: BLUP.PAR
DATAFILE ReBFDat.prn NUMBER_OF_TRAITS 1 NUMBER_OF_EFFECTS 6 OBSERVATION(S) 5 WEIGHT(S) EFFECTS: POSITIONS_IN_DATAFILE NUMBER_OF_LEVELS TYPE_OF_EFFECTS 1 233 cross 2 2 cross 3 11 cross 4 3768 cross 6 3768 cross 6 3768 cross RANDOM_RESIDUAL VALUES 0.50 RANDOM_GROUP 4 5 RANDOM_TYPE add_animal FILE ReBFped.prn (CO)VARIANCES
Note: Effect of contemporary groups in column #1 has 48 levels.
Effect of sex in column #2 has 2 levels. Effect of parity in column #3 has 11 levels. Effect of direct genetic from animal in column #4 has 3768 levels.
Effect of maternal genetic from dam in column #6 has 3768 levels.
Effect of maternal PE from dam in column #6 has 3768 levels.
Effect number 4 and 5 is random effect with add_animal type. This is genetic effects which require pedigree file. Note: covariance structure required direct, maternal and covariances.
Effect number 5 is random effect with diagonal type for PE, which is uncorrelated.
Analysis of single trait. Trait is in column #5 which is weaning weight.
0.30 -0.05 -0.05 0.15 RANDOM_GROUP 6 RANDOM_TYPE diagonal FILE (CO)VARIANCES 0.10
-
23
C. Multi-trait Model
Data File: BFdat.prn Id sex hd par age mb yb bwt mw yw wwt 39001CP 2 89 2 59 2 1996 27 9 1996 200 39001LO 1 94 3 64 2 1996 28 9 1996 214 39001TP 1 141 4 95 2 1996 29 10 1996 209 39001UB 2 151 2 58 2 1996 29 9 1996 0 39001UD 2 178 2 65 2 1996 0 9 1996 228 39001YL 1 344 1 43 1 1996 26 7 1996 211 39002KA 2 90 2 59 2 1996 20 8 1996 0 39002LO 1 94 2 60 2 1996 27 9 1996 199 39002MK 2 96 2 78 3 1996 22 10 1996 209 39002TP 1 141 8 137 2 1996 27 10 1996 208
Pedigree File: BFped.prn
Id sire dam yob 39001CP 31131MK 34075MK 1996 39001KA 35003SS 32287AWD 1996 39001LO 35010MK 33059LO 1996 39001MK 28831AI 27134MK 1996 39001TP 31397UDT 31542UTV 1996 39001UB 35014TP 34162UMA 1996 39001UD 29044MK 33053TP 1996 39001UT 35007UB 34710ULO 1996 39001YL 35006LO 35021USB 1996
Renum Data File: ReLMdat.prn Cg1 Cg2 sex par id bwt wwt dam 4 7 2 2 1 27 200 2283 44 47 1 3 2 28 214 2284 160 129 1 4 3 29 209 2285 199 151 2 2 4 29 0 2286 237 177 2 2 5 0 228 2287 303 220 1 1 6 26 211 2288 14 27 2 2 7 20 0 2289 44 47 1 2 8 27 199 2290 112 97 2 2 9 22 209 2291 160 129 1 8 10 27 208 2292
Renum Pedigree File: ReBFped.prn
Id sire dam x YOB x x x orig ID 1 3498 2283 1 1996 2 1 0 39001CP 2 3499 2284 1 1996 2 1 0 39001LO 3 3500 2285 1 1996 2 1 0 39001TP 4 3501 2286 1 1996 2 1 0 39001UB 5 3502 2287 1 1996 2 1 0 39001UD 6 3503 2288 1 1996 2 1 0 39001YL 7 3504 2289 1 1996 2 1 0 39002KA 8 3499 2290 1 1996 2 1 0 39002LO 9 3505 2291 1 1996 2 1 0 39002MK 10 3506 2292 1 1996 2 1 0 39002TP
Note: *Cg1, Cg2 is a combination of herd-month-year at birth and wean, respectively.
Parameter File: BLUP.PAR
DATAFILE ReBFDat.prn NUMBER_OF_TRAITS 2 NUMBER_OF_EFFECTS 6 OBSERVATION(S) 6 7 WEIGHT(S) EFFECTS: POSITIONS_IN_DATAFILE NUMBER_OF_LEVELS TYPE_OF_EFFECTS 1 2 233 cross 3 3 2 cross 4 4 11 cross 5 5 3768 cross 8 8 3768 cross 8 8 3768 cross RANDOM_RESIDUAL VALUES 10 5 5 100 RANDOM_GROUP 4 5 RANDOM_TYPE add_animal FILE ReBFped.prn (CO)VARIANCES 44.0 85.5 -9.8 0 85.5 668 0 0 -9.8 0 13.6 0 0 0 0 20 RANDOM_GROUP 6 RANDOM_TYPE diagonal FILE (CO)VARIANCES 10 5 5 10
Note: Effect of contemporary groups for trait 1 and 2 in column #1 and #2 has maximum at 233 levels.
Effect of sex in column #3 has 2 levels. Effect of parity in column #4 has 11 levels. Effect of direct genetic from animal in column #5 has 3768 levels.
Effect of maternal genetic from dam in column #8 has 3768 levels.
Effect of maternal PE from dam in column #8 has 3768 levels.
Effect number 4 and 5 is random effect with add_animal type. This is genetic effects which require pedigree file. Note on covariance structure required direct1, direct2, covariance between direct; maternal1, maternal2, covariance between maternal; and covariance between each direct and each maternal.
Effect number 5 is random effect with diagonal type for PE. Note on covariance structure.
Analysis of bivariate trait. Trait is in column #6 and #7 which are birth weight and weaning weight.
-
24
IX. Trouble Shooting Q: How to get accuracy in BV reports?
This option has to modify from main menu. Click menu Option > Set option, then change the Add accuracy in reports option to “YES”.
Q: All menus disappear and all buttons do not work properly.
BeefPack might be opened in disable macro. To use the interfaces in BeefPack, macro always need to be set to “ENABLE”. If no dialog appears for setting the macro, read details on page 7.
Q: How to decrease or increase convergence of program?.
This option has to modify from main menu, like adding accuracy. Click menu Option > Set option, then change the Set convergence to option to “1d-06”, “1d-08”, “1d-10”, etc.
Q: External programs of BLUPF90 seem not run?.
Generally, BeefPack calls the external programs of BLUPF90 to do the analysis in windows environment. However, if nothing seems to happen, check the correct directory, data and pedigree file. Check variance structure in model (multi-trait or random regression analysis), i.e. non-symmetric, not-positive definite, missing element.
Q: How to analyze multi-trait more than 4 triats, or random regression with other test day function?.
BeefPack does not support all models. For very complicated model, BLUPF90-PCPACK is more flexible. However, PCPACK is not easy-clicked handle, background in animal model and BLUP is important.
Q: Cannot find the BV report.
BV report and genetic trend reports (in Excel format) are kept in DATA directory, not in Program directory.
-
25
X. Who’s who
========================== Who's who in BLUPF90-PCPAK .... ========================== < Program >-------------------- BLUPF90 : Ignacy Misztal REMLF90 : Ignacy Misztal, Shogo Tsurata AIREML90 : Shogo Tsurata, Ignacy Misztal, Tom Druet RENUMMAT : Ignacy Misztal RENDOMN : Nicholus Gengler, Ignacy Misztal Windows Interface : Monchai Duangjinda PC-Compilation : Monchai Duangjinda < Libraries >-------------------- DENSEOP : Tomasz Strabel, Ignacy Misztal FSPAK : Miguel Perez-Enciso, Ignacy Misztal, Mauricio Elzo FSPAK90 : Ignacy Misztal Ginv : Rohan Fernando IOUNF : Ignacy Misztal LAPACK90 : UNI-C, Denmark; Univ. of Tennessee, USA; NAG Ltd., UK Univ. of California Berkeley, Courant Institute, Argonne National Lab, and Rice University, Alan Miller, Jack Dongarra, Sven Hammarling SPARSEM : Ignacy Misztal PROB : Luis Varona, Ignacy Misztal
Monchai DuangjindaDepartment of Animal ScienceKhon Kaen University, Khon Kaen, Thailand 40002Ignacy MisztalDepartment of Animal and Dairy ScienceThe University of Georgia, Athens, GA30602Shogo TsurutaDepartment of Animal and Dairy ScienceThe University of Georgia, Athens, GA30602( 2006Department of Animal ScienceKhon Kaen UniversityDepartment of Animal and Dairy ScienceThe University of GeorgiaContentsI. IntroductionA. BackgroundB. Overall FeaturesB. SpecificationC. Conditions of UseD. Program DownloadE. Online Registration
II. InstallationA. InstallationB. How to StartDOS environmentWindows environment
III. Genetic Evaluation At A GlanceA. General ViewB. Report ViewC. Advanced Options
IV. Model DescriptionsA. Basic Animal ModelB. Animal with PE ModelC. Maternal ModelD. Multi-trait Model
V. How Program Functions?A. Program and file componentsB. The way program worksB. Solutions at Convergence
VI. Genetic Evaluation by ExamplesA. Setup Directory and PreferencesB. Data and Pedigree filesC. Single Trait Analysis ExampleD. Maternal with PE Model ExampleE. Multi-trait Analysis Example
VII. Rapid Analysis from Previous ParametersA. Re-analysis By Old ParametersB. Creating BLUP Report From Previous Solutions
VIII. Parameter File ExamplesIX. Trouble ShootingX. Who’s who
top related