-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathDataJagsPrM.R
33 lines (27 loc) · 1.68 KB
/
DataJagsPrM.R
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
###############################################################################################################
################################# Data needed for the jags model ####################################
###############################################################################################################
#dummy variables for nr.relapses.2y.prior.study
SMSCdataC$nr.relapses.2y.prior.study1<-SMSCdataC$nr.relapses.2y.prior.study
SMSCdataC$nr.relapses.2y.prior.study2<-SMSCdataC$nr.relapses.2y.prior.study
SMSCdataC$nr.relapses.2y.prior.study1[which(SMSCdataC$nr.relapses.2y.prior.study1==2)]<-0
SMSCdataC$nr.relapses.2y.prior.study2[which(SMSCdataC$nr.relapses.2y.prior.study2==1)]<-0
SMSCdataC$nr.relapses.2y.prior.study2[which(SMSCdataC$nr.relapses.2y.prior.study2==2)]<-1
#give the data
jagsdataSMSC <- list(
Nobservations=nrow(SMSCdataC),
outcome=as.numeric(factor(SMSCdataC$outcome))-1,
npf=10,
npid=length(unique(SMSCdataC$patient.id)),
subj=as.integer(factor(SMSCdataC$patient.id)),
age=SMSCdataC$age-mean(SMSCdataC$age),
disease.duration=SMSCdataC$disease.duration-mean(SMSCdataC$disease.duration),
edss=SMSCdataC$edss-mean(SMSCdataC$edss),
nr.Gd.enhanced.lesions=as.factor(SMSCdataC$nr.Gd.enhanced.lesions),
nr.relapses.2y.prior.study1=as.factor(SMSCdataC$nr.relapses.2y.prior.study1),
nr.relapses.2y.prior.study2=as.factor(SMSCdataC$nr.relapses.2y.prior.study2),
months.since.last.relapse=SMSCdataC$months.since.last.relapse-mean(SMSCdataC$months.since.last.relapse),
treatment.naive.prior.visit=as.factor(SMSCdataC$treatment.naive.prior.visit),
gender=as.factor(SMSCdataC$gender),
treatment.during.cycle=as.factor(SMSCdataC$treatment.during.cycle)
)