作者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/cn.aspx?n=bbs/Statistics/M.1495559716.A.9AC.html