blupf90으로 다형질 모형(same model and no missing record)의 육종가 구하기
R. A. Mrode, Linear Models for the prediction of Animal Breeding Values, 2nd Edition. p85 Example 5.1
data
4 1 4.5 6.8
5 2 2.9 5.0
6 2 3.9 6.8
7 1 3.5 6.0
8 1 5.0 7.5
animal, sex, trait1, trait2
data.txt로 저장
pedigree
1 0 0
2 0 0
3 0 0
4 1 0
5 3 2
6 1 2
7 4 5
8 3 6
pedi.txt로 저장
renumf90을 위한 parameter 파일 작성
# Parameter file for program renf90; it is translated to parameter
# file for BLUPF90 family programs.
DATAFILE
data.txt
TRAITS
3 4
FIELDS_PASSED TO OUTPUT
WEIGHT(S)
RESIDUAL_VARIANCE
40.0 11.0
11.0 30.0
EFFECT
2 2 cross numer
EFFECT
1 1 cross numer
RANDOM
animal
FILE
pedi.txt
FILE_POS
1 2 3
PED_DEPTH
0
(CO)VARIANCES
20.0 18.0
18.0 40.0
설명
DATAFILE
data.txt
자료 파일 이름
TRAITS
3 4
자료 파일에서 관측치의 위치(컬럼)
FIELDS_PASSED TO OUTPUT
WEIGHT(S)
RESIDUAL_VARIANCE
40.0 11.0
11.0 30.0
잔차 분산-공분산 행렬
EFFECT
2 2 cross numer
둘째 컬럼이 첫째 형질과 둘째의 분류 효과. 숫자로 되어 있음
EFFECT
1 1 cross numer
첫째 컬럼이 첫째 형질과 둘째의 분류 효과. 숫자로 되어 있음
RANDOM
animal
임의 개체 효과
FILE
pedi.txt
혈통 파일 이름
FILE_POS
1 2 3
혈통 파일은 animal, sire, dam
PED_DEPTH
0
끝까지 혈통 추적
(CO)VARIANCES
20.0 18.0
18.0 40.0
개체 효과의 분산-공분산 행렬
실행 화면
생성된 파일
renf90.tables
Effect group 1 of column 1 with 2 levels, effect # 1
Value # consecutive number
1 3 1
2 2 2
고정효과 : 원래 번호, 개수, 새로운 번호
renf90.dat
4.5 6.8 1 2
2.9 5.0 2 5
3.9 6.8 2 3
3.5 6.0 1 1
5.0 7.5 1 4
trait1, trait2, 고정효과, 개체효과
renadd03.ped
1 2 5 1 0 2 1 0 0 7
7 0 0 3 0 0 0 0 2 2
2 6 0 2 0 1 1 1 0 4
3 6 7 1 0 2 1 0 1 6
6 0 0 3 0 0 0 2 0 1
4 8 3 1 0 2 1 0 0 8
8 0 0 3 0 0 0 2 0 3
5 8 7 1 0 2 1 0 1 5
renumbered 된 혈통. 자세한 설명은 single trait animal model 참조
renf90.par
# BLUPF90 parameter file created by RENF90
DATAFILE
renf90.dat
NUMBER_OF_TRAITS
2
NUMBER_OF_EFFECTS
2
OBSERVATION(S)
1 2
WEIGHT(S)
EFFECTS: POSITIONS_IN_DATAFILE NUMBER_OF_LEVELS TYPE_OF_EFFECT[EFFECT NESTED]
3 3 2 cross
4 4 8 cross
RANDOM_RESIDUAL VALUES
40.00000 11.00000
11.00000 30.00000
RANDOM_GROUP
2
RANDOM_TYPE
add_animal
FILE
renadd02.ped
(CO)VARIANCES
20.00000 18.00000
18.00000 40.00000
설명
DATAFILE
renf90.dat
자료 파일의 이름
NUMBER_OF_TRAITS
2
형질의 수
NUMBER_OF_EFFECTS
2
효과의 수(sex, animal)
OBSERVATION(S)
1 2
관측치의 위치
WEIGHT(S)
EFFECTS: POSITIONS_IN_DATAFILE NUMBER_OF_LEVELS TYPE_OF_EFFECT[EFFECT NESTED]
3 3 2 cross
4 4 8 cross
셋째 컬럼이 trait1과 trait2의 효과, 레벨 개수는 2, 분류 효과
넷째 컬럼이 trait1과 trait2의 효과, 레벨 개수는 8, 분류 효과
RANDOM_RESIDUAL VALUES
40.00000 11.00000
11.00000 30.00000
잔차 효과의 분산-공분산 행렬
RANDOM_GROUP
2
효과 중 둘째 효과가 임의 효과
RANDOM_TYPE
add_animal
additive genetic animal effect
FILE
renadd02.ped
혈통 파일의 이름
(CO)VARIANCES
20.00000 18.00000
18.00000 40.00000
개체 효과의 분산-공분산 효과
blupf90 실행 화면
solutions 결과 파일
trait/effect level solution
1 1 1 4.36086700
2 1 1 6.79989768
1 1 2 3.39726169
2 1 2 5.88029605
1 2 1 -0.31611757
2 2 1 -0.47898372
1 2 2 -0.01023894
2 2 2 -0.01267064
1 2 3 0.27580827
2 2 3 0.51723842
1 2 4 0.24375551
2 2 4 0.39196151
1 2 5 -0.27033146
2 2 5 -0.47783034
1 2 6 0.15091558
2 2 6 0.27959802
1 2 7 -0.01539252
2 2 7 -0.00761009
1 2 8 -0.07839191
2 2 8 -0.17034150
형질 2개, 효과 2개
1번 animal은 원래 7번 개체, trait1의 육종가는 –0.31611757,
trait2의 육종가는 –0.47898372
'Animal Breeding > BLUPF90' 카테고리의 다른 글
blupf90으로 다형질 모형(different model)의 육종가 구하기 (0) | 2014.03.12 |
---|---|
blupf90으로 다형질 모형(same model with missing record)의 육종가 구하기 (0) | 2014.03.12 |
blupf90으로 common environmental effect 모형의 육종가 구하기 (0) | 2014.02.25 |
blupf90으로 repeatability 모델의 육종가 구하기 (0) | 2014.02.25 |
blupf90으로 단형질 개체 모형(Single Trait Animal Model)의 육종가(Breeding Value) 구하기 (0) | 2014.01.28 |