作者blizzard258 (鯛魚燒)
看板Statistics
標題[程式] 有關使用SAS建立ACE模型的問題
時間Wed May 24 01:15:14 2017
TITLE請註名軟體類別 如[程式] stata
發文前請把灰色的註解刪除(可使用ctrl+y刪除一整行)
[軟體程式類別]:
SAS
[程式問題]:
ACE Model
(家族聚集性模型 變異數P = 變異數A + 變異數C + 變異數E
變異數P=總表現型變異 total phenotypic variance
變異數A=附加遺傳因素變異 additive genetic variance
變異數C=共享環境因素變異 common environmental variance
變異數E=特異環境因素變異 non-shared enviromental variance)
[軟體熟悉度]:
新手
[問題敘述]:
目前我取得了一筆資料是有關家族性疾病的資料
約略格式如下(該筆資料約有1百多萬筆以上)
罹病狀況相關的欄位中 有罹病為1 無罹病為0 沒有資料則為.
ID age 觀察時間(年) 自身罹病 配偶罹病 手足數目 手足罹病 父親罹病 母親罹病
A 30 7 0 . 2 0 1 .
B 60 6 1 1 3 1 0 1
C 70 10 0 0 . . . 1
D 45 5 1 0 2 1 1 1
E 37 8 1 1 . . 0 0
F 51 2 0 . . . . .
目前我已用SAS的cox regression跑出雙親、手足與配偶罹病對自身罹病的相對風險,
也計算出該疾病的分年齡層分性別盛行率,
想請教大家如果想要知道該疾病的附加遺傳因素、共享環境影響與特異環境影響程度,
也就是建立ACE model的話,
程式碼的部分應該如何編寫才比較好,
我已上網和翻閱許多書籍,
發現可以用配偶(遺傳度ar=0)以及手足(遺傳度ar=0.5)來計算ACE,
以及好像是用 proc mixed來建立,
但是有關如何實際使用SAS來建立ACE medel這部分的文章較少,
希望各位大大們能夠提供一些編寫上的建議,
如果有不用使用SAS也能估算ACE的其他方法,
也希望大大們能不吝指教,
真的非常謝謝大家!!!
[程式範例]:
--
※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 163.25.119.94
※ 文章網址: https://webptt.com/m.aspx?n=bbs/Statistics/M.1495559716.A.9AC.html