如何用S编程语言进行紧急医疗数据转换？
|
您将如何在R中执行此操作？
（这是一个数据准备任务）
从不良事件数据集中，导出治疗紧急不良事件
事件数据集：
对于每个身体系统和首选术语，每个患者一行
要么：
该不良事件发生在基准后时期，而不是基准时期，或者
即使事件确实发生在基准期间，
基线后发生的严重性高于
基线
变量：
severity = 1 , 2 , 3 (integer code for mild moderate severe)
patid visit bodysys prefterm
Baseline rows are rows such that visit<=2
Post baseline rows are rows such that visit>2
data base1_dset ;
set ae_dset ;
if visit<=2 ;
proc sort data=base1_dset ;
by patid bodysys prefterm severity ;
data base2_dset ;
set base1_dset ;
by patid bodysys prefterm severity ;
if last.prefterm ;
data post1_dset ;
set ae_dset ;
if visit> 2 ;
proc sort data=post1_dset ;
by patid bodysys prefterm severity ;
data post2_dset ;
set post1_dset ;
by patid bodysys prefterm severity ;
if last.prefterm ;
rename severity = severity2 ;
data new_ae_dset ;
merge base2_dset post2_dset ;
by patid bodysys prefterm ;
if severity2>severity or severity==. ;
inlist ae_dset ;
if not ( visit<=2 ) deleterow ;
select severity=max(severity) by patid bodysys prefterm ;
sendoff(base2_dset) patid bodysys prefterm severity ;
inlist ae_dset ;
if not ( visit>2 ) deleterow ;
select severity2=max(severity) by patid bodysys prefterm ;
sendoff(post2_dset) patid bodysys prefterm severity2 ;
inlist base2_dset post2_dset ;
mergeby patid bodysys prefterm ;
if not ( severity2>severity or severity is null ) deleterow ;
sendoff(new_ae_dset) patid bodysys prefterm severity2 ;
没有找到相关结果
已邀请:
1 个回复
赣借
smallvisit <- ae_dset[ ae_dset$visit <= 2, ] bigvisit <- ae_dset[ ae_dset$visit > 2, ] nams <- c(\"patid\", \"bodysys\", \"prefterm\") smallvisitsorted <- smallvisit[ do.call( order, smallvisit[nams] ), ] smallvisitsplit <- split( smallvisitsorted, smallvisitsorted[nams], drop=TRUE ) last <- function(a){ tail( a, 1 ) } smallvisitlast <- as.data.frame( t( sapply( smallvisitsplit, last ) ) ) mergedvisit <- merge( bigvisit, smallvisitlast, by=nams, all.x=TRUE ) new_ae_dset <- mergedvisit[ mergedvisit$severity.x > mergedvisit$severity.y | is.na( mergedvisit$severity.y ) , ]看起来像然后，使用此代码，将看起来