Skip to main content

Half-whiskers plot of difference in treatment effect

Below an alternative graphical representation of difference in treatment effects.




data ests;
missing z;
infile datalines delimiter=',' MISSOVER;
input behandling time VASe VASs lowere uppere lowers uppers;
datalines;
1,0,100,.z,100,100,.z,.z
1,1,78.88888889,.z,78.88888889,120,.z,.z
1,2,53.23552736,.z,53.23552736,100,.z,.z
1,3,27.7394636,.z,27.7394636,72.34042553,.z,.z
1,4,13.67521368,.z,13.67521368,78.72340426,.z,.z
1,5,7.948717949,.z,7.948717949,76.59574468,.z,.z
1,6,7.948717949,.z,7.948717949,77.65957447,.z,.z
1,7,2.948717949,.z,2.948717949,82.9787234,.z,.z
1,8,2.948717949,.z,2.948717949,81.91489362,.z,.z
2,0,.z,100,.z,.z,100,100
2,1,.z,64.40092166,.z,.z,11.11111111,64.40092166
2,2,.z,39.42105263,.z,.z,0,39.42105263
2,3,.z,22,.z,.z,0,22
2,4,.z,17.44444444,.z,.z,0,17.44444444
2,5,.z,9.545454545,.z,.z,0,9.545454545
2,6,.z,6.666666667,.z,.z,0,6.666666667
2,7,.z,5.409356725,.z,.z,0,5.409356725
2,8,.z,2.631578947,.z,.z,0,2.631578947
;;
run;

ods graphics on /noborder;
title 'Treatment effects';
proc sgplot data=ests noautolegend;
highlow x=time low=lowere high=uppere/HIGHCAP=Serif LINEATTRS=(color=cx003399); 
highlow x=time low=lowers high=uppers/LOWCAP=Serif LINEATTRS=(color=cxA23A2E); 
vbox VASe/category=time NOOUTLIERS EXTREME NOMEAN NOCAPS BOXWIDTH=0 CONNECT=MEAN meanattrs=(size=9) CONNECTATTRS=(color=cx003399);
vbox VASs/category=time NOOUTLIERS EXTREME NOMEAN NOCAPS BOXWIDTH=0 CONNECT=MEAN meanattrs=(size=9) CONNECTATTRS=(color=cxA23A2E);
scatter x=time y=VASe/MARKERATTRS=(symbol=SquareFilled color=cx003399);
scatter x=time y=VASs/MARKERATTRS=(symbol=CircleFilled color=cxA23A2E);
yaxis label="Treatment";
xaxis label="Time";
run;

Comments

Popular posts from this blog

Alder/korrekt århundrede udfra cpr nummer

De fleste, der arbejder med registre eller databaser, står ofte med problemstillingen, at alder er uoplyst, medens cpr-nummer er kendt. Hvordan regner man den ud? Følgende regel er gældende: Hvis syvende ciffer er 0, 1, 2 eller 3 er man født i det 20. århunderede (1900-tallet) Ligeledes, hvis syvende ciffer er 4 eller 9, og årstallet (femte og sjette ciffer) er større end eller lig 37. Endelig er man født i det 19. århundrede (1800-tallet) hvis syvende ciffer er 5, 6, 7 eller 8 og årstallet er større end eller lig 58. Nedenfor finder du eksempel i SAS kode: En lille makro, der udover fødselsdato også udregner køn samt den præcise alder givet datovariabel. Kilde: Opbygning af CPR nummeret, cpr.dk proc format library=work; value gender 0="Female" 1="Male" ; run; %macro agefromCPR(cpr,datevar=inddto,birthvar=birth,agevar=age); dy_temp=input(substrn(&cpr,1,2),2.); mt_temp=input(substrn(&cpr,3,2),2.); yr_temp=input(substrn(&cpr,5,2),

Comorbidity indexes in SQL

Generating Elixhauser comorbidity index from Danish National Health Register as relational database. ( ICD 10 Coding  in SAS) A lookup-table based version of Charlson comorbidity index I made in SQL. A similar approach can be applied to Elixhauser. SELECT V_CPR, MAX(EI1)+MAX(EI2)+MAX(EI3)+MAX(EI4)+MAX(EI5)+ MAX(EI6)+MAX(EI7)+MAX(EI8)+MAX(EI9)+MAX(EI10)+ MAX(EI11)+MAX(EI12)+MAX(EI13)+MAX(EI14)+MAX(EI15)+ MAX(EI16)+MAX(EI17)+MAX(EI18)+MAX(EI19)+MAX(EI20)+ MAX(EI21)+MAX(EI22)+MAX(EI23)+MAX(EI24)+MAX(EI25)+ MAX(EI26)+MAX(EI27)+MAX(EI28)+MAX(EI29)+MAX(EI30)+MAX(EI31) AS Elixhauser FROM (SELECT V_CPR, -- Congestive Heart Failure CASE WHEN DIAG LIKE 'DI099%' OR DIAG LIKE 'DI110%' OR DIAG LIKE 'DI130%' OR DIAG LIKE 'DI132%' OR DIAG LIKE 'DI255%' OR DIAG LIKE 'DI420%' OR DIAG LIKE 'DI425%' OR DIAG LIKE 'DI426%' OR DIAG LIKE 'DI427%' OR DIAG LIKE 'DI428%' OR DIAG LIKE 'DI429%' OR D

HackRF on Windows 8

This technical note is based on  an extract from thread . I have made several changes and added recommendations. I have experienced lot of latency using GnuRadio and HackRF on Pentoo Linux, so I wanted to try out GnuRadio on Windows. HackRF One is a transceiver, so besides SDR capabilities, it can also transmit signals, inkluding sweeping a given range, uniform and Gaussian signals. Pentoo Linux provides the most direct access to HackRF and toolboxes. Install Pentoo Linux on a separate drive, then you can use osmocom_siggen from a terminal to transmit signals such as near-field GSM bursts, which will only be detectable within a meter. Installation of MGWin and cmake: Download and install the following packages: - MinGW Setup (Go to the Installer directory and download setup file) - CMake (I am using CMake 3.2.2 and I installed it in C:\CMake, this path is important in the commands we must send in the MinGW shell) Download and extract the packages