楼主: 雪荆
7062 12

[原创博文] 发病人年数sas程序 [推广有奖]

  • 2关注
  • 0粉丝

本科生

54%

还不是VIP/贵宾

-

威望
0
论坛币
5 个
通用积分
0.0004
学术水平
0 点
热心指数
0 点
信用等级
0 点
经验
831 点
帖子
60
精华
0
在线时间
103 小时
注册时间
2011-7-27
最后登录
2020-2-20

+2 论坛币
k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群
赵安豆老师微信:zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

感谢您参与论坛问题回答

经管之家送您两个论坛币!

+2 论坛币
请求各位高手帮帮忙,我现在想要分析HIV的发病密度,但是想了好久没有想到正确的程序,急求各位高手解答!
要分析的是一个随访数据,每一个观察对象因脱失,其随访年数不同。R1-R5表示第一次-第五次检测结果,T1-T5表示每一次检测结果对应的时间。数据有几种类型:1、(见obs1)第一次检测结果是阳性,那么他得随访人年数是0;2、(见OBS2、OBS3)检测都为0,那么他的随访人年数是T3-T1(存在缺失)和T5-T1(无缺失)。4、最后一种是可能有时候因为检测的假阴性率的原因导致之前检测为阳性的,之后检测结果却为阴性(该疾病不能转阴),但是计算的时候就要选择第一次阳转的时间距最开始的时间作为发病人年数。数据处理比较困难,麻烦各位高手多多帮忙!O(∩_∩)O谢谢!

OBS

R1

R2

R3

R4

R5

T1

T2

T3

T4

T5

1

1

.

.

.

.

2006/11/1

 

 

 

 

2

0

0

0

.

.

2006/10/25

2007/12/12

2010/9/15

 

 

3

0

0

0

0

0

2006/11/6

2007/11/2

2008/6/2

2008/6/6

2009/6/16

4

0

0

1

1

1

2006/11/7

2009/11/26

2010/6/21

2010/8/24

2011/6/28

5

0

1

0

.

.

2006/11/13

2007/11/1

2009/6/17

 

 

二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

关键词:sas程序 检测结果 数据处理 obs HIV 程序 检测

回帖推荐

hopewell 发表于2楼  查看完整内容

本帖被以下文库推荐

沙发
hopewell 发表于 2011-12-18 18:22:04 |只看作者 |坛友微信交流群
  1. data raw;
  2.     infile datalines truncover;
  3.     input (R1-R5 T1-T5) (5*:8. 5*:yymmdd10.);
  4.     format t1-t5 yymmdd10.;
  5. datalines;
  6. 1 . . . . 2006/11/1
  7. 0 0 0 . . 2006/10/25 2007/12/12 2010/9/15
  8. 0 0 0 0 0 2006/11/6 2007/11/2 2008/6/2 2008/6/6 2009/6/16
  9. 0 0 1 1 1 2006/11/7 2009/11/26 2010/6/21 2010/8/24 2011/6/28
  10. 0 1 0 . . 2006/11/13 2007/11/1 2009/6/17
  11. ;
  12. data out(drop=i index);
  13.     set raw;
  14.     array rarr{*} r:;
  15.     do i=1 to dim(rarr) until(rarr(i) ne 0);
  16.         index=i-(missing(rarr(i)));
  17.     end;
  18.     year=intck('year',t1,choosen(index,of t:));
  19. run;
复制代码
已有 2 人评分经验 论坛币 学术水平 热心指数 信用等级 收起 理由
bakoll + 3 + 3 精彩帖子
雪荆 + 1 + 1 + 1 热心帮助其他会员

总评分: 经验 + 3  论坛币 + 3  学术水平 + 1  热心指数 + 1  信用等级 + 1   查看全部评分

观钓颇逾垂钓趣 种花何问看花谁

使用道具

藤椅
雪荆 发表于 2011-12-18 20:53:34 |只看作者 |坛友微信交流群
嗯,谢谢谢谢!真是感激不尽,纠结了好久,多谢您的指点啦!

使用道具

板凳
webgu 发表于 2011-12-18 21:10:35 |只看作者 |坛友微信交流群
发病密度一般是的定义形式是:出现的病例数/观察人年数,楼立直接定义发病人年数是否合适?
SAS资源
1. SAS 微信:StatsThinking
2. SAS QQ群:348941365

使用道具

报纸
webgu 发表于 2011-12-18 21:18:23 |只看作者 |坛友微信交流群
hopewell 发表于 2011-12-18 18:22
对于前几次检测结果缺失和中间有假阳性的还是不好解决。
SAS资源
1. SAS 微信:StatsThinking
2. SAS QQ群:348941365

使用道具

地板
雪荆 发表于 2011-12-18 21:53:40 |只看作者 |坛友微信交流群
其实是我不是太会用SAS程序,不知道怎么用sas直接计算出来。好像我用到数据里面有观测序号ID的时候就出错了。提示Invalid argument to function INTCK at line 1594 column 14.
现在好像还是不对!望您再多多指点一下。我尝试过用宏来编,但是不太会。

使用道具

7
雪荆 发表于 2011-12-18 21:53:51 |只看作者 |坛友微信交流群
其实是我不是太会用SAS程序,不知道怎么用sas直接计算出来。好像我用到数据里面有观测序号ID的时候就出错了。提示Invalid argument to function INTCK at line 1594 column 14.
现在好像还是不对!望您再多多指点一下。我尝试过用宏来编,但是不太会。

使用道具

8
hopewell 发表于 2011-12-18 22:04:57 |只看作者 |坛友微信交流群
水平太差,看不懂,啥叫阳转,假阴性,1是阳性发病的意思吗,缺失值有什么特征?
观钓颇逾垂钓趣 种花何问看花谁

使用道具

9
雪荆 发表于 2011-12-18 22:12:15 |只看作者 |坛友微信交流群
嗯,1表示检测结果为阳性,如果第一次检测结果是阴性的,第二次检测为阳性,就是发病(阳转)。对于不能转阴的疾病,如果已经检测出阳性结果,之后再检测出现结果为阴性,说明这个结果可能是由于检测技术或其他原因导致的错误的结果,阴性结果是不可靠的,所以是假阴性。您能重新帮帮忙看看用什么语句能够解决这个问题呢?其实是很希望能够直接计算出发病密度的。谢谢啦!O(∩_∩)O~

使用道具

10
雪荆 发表于 2011-12-18 22:12:21 |只看作者 |坛友微信交流群
嗯,1表示检测结果为阳性,如果第一次检测结果是阴性的,第二次检测为阳性,就是发病(阳转)。对于不能转阴的疾病,如果已经检测出阳性结果,之后再检测出现结果为阴性,说明这个结果可能是由于检测技术或其他原因导致的错误的结果,阴性结果是不可靠的,所以是假阴性。您能重新帮帮忙看看用什么语句能够解决这个问题呢?其实是很希望能够直接计算出发病密度的。谢谢啦!O(∩_∩)O~

使用道具

您需要登录后才可以回帖 登录 | 我要注册

本版微信群
加好友,备注cda
拉您进交流群

京ICP备16021002-2号 京B2-20170662号 京公网安备 11010802022788号 论坛法律顾问:王进律师 知识产权保护声明   免责及隐私声明

GMT+8, 2024-4-26 22:06