cancel
Turn on suggestions
Auto-suggest helps you quickly narrow down your search results by suggesting possible matches as you type.
Showing results for 
Search instead for 
Did you mean: 
ChristosK
‎09-19-2020
  • 141 Posts
  • 28 Likes Given
  • 0 Solutions
  • 2 Likes Received

Re: Error in generating data step

by Quartz | Level 8 in SAS Studio
‎07-22-2018 06:15 PM
‎07-22-2018 06:15 PM
I thought it had initially worked because  the error message disappeared. We now  have 24 entries. However proc freq shows that there are only 23 entries that have count data, and missing.   Hmmmm ... View more

Re: Error in generating data step

by Quartz | Level 8 in SAS Studio
‎07-22-2018 05:55 PM
‎07-22-2018 05:55 PM
How do you do that? ... View more

Re: Error in generating data step

by Quartz | Level 8 in SAS Studio
‎07-22-2018 05:38 PM
‎07-22-2018 05:38 PM
Tried onto notepad, and a word document. Found one tab which I deleted but it didnt make a difference. Whats "adding truncover"?     ... View more

Re: Error in generating data step

by Quartz | Level 8 in SAS Studio
‎07-22-2018 05:23 PM
‎07-22-2018 05:23 PM
This is interesting. Might be a bug in my SAS software? I cut and pasted your version  and still get the same error message. 1 OPTIONS NONOTES NOSTIMER NOSOURCE NOSYNTAXCHECK; 72 73 data Concordance_MAP; 74 input Measure $ DV: $15. IV: $15. CC LCL UCL Count; 75 datalines; NOTE: Invalid data for Count in line 97 1-2. RULE: ----+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8----+----9----+----0 97 LL p3min_A_MAP p3min_LL_MAP 0.4575 0.2335 0.6353 54 Measure=LA DV=p3min_A_MAPp3m IV=0.5436 CC=0.3843 LCL=0.6714 UCL=90 Count=. _ERROR_=1 _N_=21 NOTE: SAS went to a new line when INPUT statement reached past the end of a line. NOTE: The data set USER.CONCORDANCE_MAP has 23 observations and 7 variables. NOTE: DATA statement used (Total process time): ... View more

Re: Concordance correlation coefficient

by Quartz | Level 8 in Statistical Procedures
‎07-22-2018 05:08 PM
‎07-22-2018 05:08 PM
Rick:   I have been using your macro to generate concordance coefficients but in one example the lower confidence interval is negative. Is there a problem somewhere here?  Whilst I understand that the direction of effect of any correlation coefficient just describes the direction of the relationship, I am puzzled in this example as this is the only model where  I have found a negative value .   Christos    title "Concordance Correlation Coefficient"; proc sgplot data=Schumann.NumBP noautolegend; scatter y=i3min_LL_DIA x= i3min_A_DIA; lineparm x=0 y=0 slope=1 / clip; yaxis label='i3min_LL_DIA' min=28 max=108; xaxis label='i3min_A_DIA' min=38 max=111; run; proc iml; ******************************************************************************* * Enter the appropriate SAS data set name in the use statement and enter the * * appropriate variable names in the read statements. * *******************************************************************************; use Schumann.NumBP; read all var {i3min_LL_DIA} into var1; read all var {i3min_A_DIA} into var2; start concor(x1, x2);nonmiss=loc(x1#x2^=.);var1=x1[nonmiss];var2=x2[nonmiss];free nonmiss;n=nrow(var1); mu1=mean(var1);mu2=mean(var2);Sigma = cov( var1||var2 );sigma11=Sigma[1,1];sigma22=Sigma[2,2];sigma12=Sigma[1,2]; lshift=(mu1-mu2)/((sigma11#sigma22)##0.25);rho=sigma12/sqrt(sigma11#sigma22);z=log((1+rho)/(1-rho))/2;se_z=1/sqrt(n-3); t=tinv(0.975,n-3);z_low=z-(se_z#t);z_upp=z+(se_z#t);rho_low=(exp(2#z_low)-1)/(exp(2#z_low)+1);rho_upp=(exp(2#z_upp)-1)/(exp(2#z_upp)+1);crho=(2#sigma12)/((sigma11+sigma22)+((mu1-mu2)##2)); z=log((1+crho)/(1-crho))/2;if sigma12^=0 then do;t1=((1-(rho##2))#(crho##2))/((1-(crho##2))#(rho##2));t2=(2#(crho##3)#(1-crho)#(lshift##2))/(rho#((1-(crho##2))##2));t3=((crho##4)#(lshift##4))/(2#(rho##2)#((1-(crho##2))##2));se_z=sqrt((t1+t2-t3)/(n-2));end; else se_z=sqrt(2#sigma11#sigma22)/((sigma11+sigma22+((mu1-mu2)##2))#(n-2));t=tinv(0.975,n-2); z_low=z-(se_z#t);z_upp=z+(se_z#t);crho_low=(exp(2#z_low)-1)/(exp(2#z_low)+1);crho_upp=(exp(2#z_upp)-1)/(exp(2#z_upp)+1); Results=n//mu1//mu2//sigma11//sigma22//sigma12//rho_low//rho//rho_upp//crho_low//crho//crho_upp;r_name={'SampleSize' 'Mean_1' 'Mean_2' 'Variance_1' 'Variance_2' 'Covariance' 'Corr LowerCL' 'Corr' 'Corr UpperCL' 'ConcCorr LowerCL' 'ConcCorr' 'ConcCorr UpperCL'};print Results [rowname=r_name format=best6. L='Estimated Correlation and Concordance Correlation'];finish concor; run concor(var1, var2); ... View more

Error in generating data step

by Quartz | Level 8 in SAS Studio
‎07-22-2018 04:52 PM
‎07-22-2018 04:52 PM
Experiencing a strange problem:   Its with the 3rd line from the bottom:   LL p3min_A_MAP p3min_LL_MAP 0.4575 0.2335 0.6353 54  Prior to entering the count values , I have 24 entries in my dataset. I have changed the spacing, variable name, treated it as a categorical value, and tried entering this one line as the only value in the dataset (which actually works). I asked a colleague to try it and he had no problems with this code. Have also restarted SAS Studio, and rebooted the virtual machine .   Any solutions?         data Concordance_MAP; input Measure $ DV: $15. IV: $15. CC LCL UCL Count; datalines; LA BI_A_MAP BI_LA_MAP 0.5813 0.4376 0.6960 89 LL BI_A_MAP BI_LL_MAP 0.3992 0.1926 0.5719 58 Nex BI_A_MAP BI_Nex_MAP 0.6511 0.5149 0.7553 89 UA BI_A_MAP BI_UA_MAP 0.5274 0.356 0.6644 77 LA i3min_A_MAP i3min_LA_MAP 0.5942 0.4548 0.7051 89 LL i3min_A_MAP i3min_LL_MAP 0.2979 0.0579 0.5054 54 Nex i3min_A_MAP i3min_Nex_MAP 0.8376 0.764 0.8897 90 UA i3min_A_MAP i3min_UA_MAP 0.2369 0.049 0.4086 75 LA i15min_A_MAP i15min_LA_MAP 0.6234 0.4911 0.7275 90 LL i15min_A_MAP i15min_LL_MAP 0.5584 0.3475 0.7155 54 Nex i15min_A_MAP i15min_Nex_MAP 0.6498 0.5169 0.7521 90 UA i15min_A_MAP i15min_UA_MAP 0.5681 0.4265 0.6825 74 LA i30min_A_MAP i30min_LA_MAP 0.6298 0.4908 0.7376 90 LL i30min_A_MAP i30min_LL_MAP 0.5181 0.3028 0.6831 54 Nex i30min_A_MAP i30min_Nex_MAP 0.6895 0.5624 0.7848 90 UA i30min_A_MAP i30min_UA_MAP 0.6076 0.4677 0.7178 74 LA i45min_A_MAP i45min_LA_MAP 0.6877 0.5604 0.7832 89 LL i45min_A_MAP i45min_LL_MAP 0.4402 0.2084 0.6252 52 Nex i45min_A_MAP i45min_Nex_MAP 0.7509 0.6422 0.8300 89 UA i45min_A_MAP i45min_UA_MAP 0.5898 0.4348 0.7108 73 LA p3min_A_MAP p3min_LA_MAP 0.5436 0.3843 0.6714 90 LL p3min_A_MAP p3min_LL_MAP 0.4575 0.2335 0.6353 54 Nex p3min_A_MAP p3min_Nex_MAP 0.6559 0.5297 0.7538 90 UA p3min_A_MAP p3min_UA_MAP 0.424 0.2609 0.5636 72 ; run; 1 OPTIONS NONOTES NOSTIMER NOSOURCE NOSYNTAXCHECK; 72 73 data Concordance_MAP; 74 input Measure $ DV: $15. IV: $15. CC LCL UCL Count; 75 datalines; NOTE: Invalid data for Count in line 102 1-2. RULE: ----+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8----+----9----+----0 102 LL p3min_A_MAP p3min_LL_MAP 0.4575 0.2335 0.6353 54 Measure=LA DV=p3min_A_MAPp3m IV=0.5436 CC=0.3843 LCL=0.6714 UCL=90 Count=. _ERROR_=1 _N_=21 NOTE: SAS went to a new line when INPUT statement reached past the end of a line. NOTE: The data set USER.CONCORDANCE_MAP has 23 observations and 7 variables. NOTE: DATA statement used (Total process time): real time 0.02 seconds cpu time 0.02 seconds         ... View more

Re: specifying a keylegend...

by Quartz | Level 8 in Graphics Programming
‎07-21-2018 08:59 PM
‎07-21-2018 08:59 PM
Thank you Jeff!   Worked perfectly and I was able to position it where I needed.   Christos  ... View more

Re: specifying a keylegend...

by Quartz | Level 8 in Graphics Programming
‎07-21-2018 07:23 PM
‎07-21-2018 07:23 PM
Ok so I was able to get a keylegend by specifying "Measures" instead of "IV":   But I can't control the appearance or location of this keylegend with a keylegend statement.     proc sgplot data=Concordance_MAP_lab nowall noborder ; title 'Concordance coefficient by time'; styleattrs datacontrastcolors=(CX016381 CX1BC9FF CXFF7F00 CX814101 ) ; scatter x=DV y=CC / group=Measure yerrorlower=LCL yerrorupper=UCL noerrorcaps markerattrs=(symbol=circlefilled size =5 ) groupdisplay=cluster clusterwidth =0.2; yaxis display =(noticks nolabel noline) VALUEATTRS=(Color=Black Family=Arial Size=8 Style=Normal Weight=Normal) ; xaxis display =(noticks nolabel noline) VALUEATTRS=(Color=Black Family=Arial Size=8 Style=Normal Weight=Normal); refline 0.75 / axis=y discretethickness=1 lineattrs=(color=cxd0e0f0) LABEL = ('0.75'); /* keylegend "Measure"/location=outside; */ run;   ... View more

specifying a keylegend...

by Quartz | Level 8 in Graphics Programming
‎07-21-2018 06:29 PM
‎07-21-2018 06:29 PM
data Concordance_MAP; input Measure $ DV: $15. IV: $15. CC LCL UCL ; datalines; LA BI_A_MAP BI_LA_MAP 0.5813 0.4376 0.696 LL BI_A_MAP BI_LL_MAP 0.3992 0.1926 0.5719 Nex BI_A_MAP BI_Nex_MAP 0.6511 0.5149 0.7553 UA BI_A_MAP BI_UA_MAP 0.5274 0.356 0.6644 LA i3min_A_MAP i3min_LA_MAP 0.5942 0.4548 0.7051 LL i3min_A_MAP i3min_LL_MAP 0.2979 0.0579 0.5054 Nex i3min_A_MAP i3min_Nex_MAP 0.8376 0.764 0.8897 UA i3min_A_MAP i3min_UA_MAP 0.2369 0.049 0.4086 LA i15min_A_MAP i15min_LA_MAP 0.6234 0.4911 0.7275 LL i15min_A_MAP i15min_LL_MAP 0.5584 0.3475 0.7155 Nex i15min_A_MAP i15min_Nex_MAP 0.6498 0.5169 0.7521 UA i15min_A_MAP i15min_UA_MAP 0.5681 0.4265 0.6825 LA i30min_A_MAP i30min_LA_MAP 0.6298 0.4908 0.7376 LL i30min_A_MAP i30min_LL_MAP 0.5181 0.3028 0.6831 Nex i30min_A_MAP i30min_Nex_MAP 0.6895 0.5624 0.7848 UA i30min_A_MAP i30min_UA_MAP 0.6076 0.4677 0.7178 LA i45min_A_MAP i45min_LA_MAP 0.6877 0.5604 0.7832 LL i45min_A_MAP i45min_LL_MAP 0.4402 0.2084 0.6252 Nex i45min_A_MAP i45min_Nex_MAP 0.7509 0.6422 0.83 UA i45min_A_MAP i45min_UA_MAP 0.5898 0.4348 0.7108 LA p3min_A_MAP p3min_LA_MAP 0.5436 0.3843 0.6714 LL p3min_A_MAP p3min_LL_MAP 0.4575 0.2335 0.6353 Nex p3min_A_MAP p3min_Nex_MAP 0.6559 0.5297 0.7538 UA p3min_A_MAP p3min_UA_MAP 0.424 0.2609 0.5636 ; run; /* ods listing style=listing; */ proc template; define style mystyle; parent=styles.sasweb; class graphwalls / frameborder=off; class graphbackground / color=white; end; run; ods listing close; ods html style=mystyle path='.' file='newgraph.html'; ods graphics / reset=all border=off width=600 height=480; proc sgplot data=Concordance_MAP nowall noborder ; title 'Concordance coefficient by time'; styleattrs datacontrastcolors=(CX016381 CX1BC9FF CXFF7F00 CX814101 ) ; scatter x=DV y=CC / group=IV yerrorlower=LCL yerrorupper=UCL noerrorcaps markerattrs=(symbol=circlefilled size =5 ) groupdisplay=cluster clusterwidth =0.2; yaxis display =(noticks nolabel noline) VALUEATTRS=(Color=Black Family=Arial Size=8 Style=Normal Weight=Normal) ; xaxis display =(noticks nolabel noline) VALUEATTRS=(Color=Black Family=Arial Size=8 Style=Normal Weight=Normal); refline 0.75 / axis=y discretethickness=1 lineattrs=(color=cxd0e0f0) LABEL = ('0.75'); run; I tried entering a keylegend statement before the "run" line, specifying "Measure" as the categorical variable I wanted to be described in the kelegend using the colours : CX016381 CX1BC9FF CXFF7F00 CX814101.   It is not making a difference so don't know whet the error is. Without the keylegend statement  I just get a decryption for "IV"   Also I'm getting this error:   NOTE: Writing HTML Body file: newgraph.html ERROR: Insufficient authorization to access /opt/sasinside/SASConfig/Lev1/SASApp/newgraph.html. ERROR: No body file. HTML output will not be created. ... View more

Re: Changing the length of a character variable in a data input statem...

by Quartz | Level 8 in SAS Studio
‎07-20-2018 08:51 PM
‎07-20-2018 08:51 PM
Thank you so much!   I never would have thought of this approach on my own. Was trying to do this with a length statement and format / informat statements. ... View more

Changing the length of a character variable in a data input statement....

by Quartz | Level 8 in SAS Studio
‎07-20-2018 07:51 PM
‎07-20-2018 07:51 PM
I have values within some character variables that are longer than 8, which is the default.   How do I lengthen these values so that the data input step  can read these in?   data Concordance_MAP; input DV $ IV $ CC LCL UCL ; datalines; BI_A_MAP BI_LA_MAP 0.5813 0.4376 0.696 etc..... ; run; ... View more

Re: SAS Studio won't connect to servers

by Quartz | Level 8 in SAS Studio
‎07-08-2018 08:57 PM
‎07-08-2018 08:57 PM
Didn't have to update the browser. All I had to do was manually force every programme off.  Previous attempts are rebooting were not enough. Had to force shut down everything following the latest SAS update.   Thanks again.   Christos ... View more

Re: Concordance correlation coefficient

by Quartz | Level 8 in Statistical Procedures
‎07-08-2018 08:46 PM
‎07-08-2018 08:46 PM
Most of this seems to work,  however I'm getting an error message with this part:     proc gplot uniform data=dice_baselineall; plot BI_Nex_MAP*BI_A_MAP var1*var2/overlay vaxis=axis1 haxis=axis2 nolegend frame; axis1 label=(a=90 'BI Nex MAP') minor=none; axis2 label=('BI A line MAP') minor=none; symbol1 value=star color=black interpol=none; symbol2 value=none color=black interpol=join; title "Concordance Correlation Coefficient"; run;     1 OPTIONS NONOTES NOSTIMER NOSOURCE NOSYNTAXCHECK; 72 73 74 proc gplot uniform data=dice_baselineall; ERROR: Procedure GPLOT not found. 75 plot BI_Nex_MAP*BI_A_MAP var1*var2/overlay vaxis=axis1 haxis=axis2 nolegend frame; 76 axis1 label=(a=90 'BI Nex MAP') minor=none; 77 axis2 label=('BI A line MAP') minor=none; 78 symbol1 value=star color=black interpol=none; 79 symbol2 value=none color=black interpol=join; 80 title "Concordance Correlation Coefficient"; 81 run;   NOTE: The SAS System stopped processing this step because of errors. NOTE: PROCEDURE GPLOT used (Total process time): real time 0.00 seconds cpu time 0.01 seconds      The proc IML seems to work just fine though. ... View more

SAS Studio won't connect to servers

by Quartz | Level 8 in SAS Studio
‎07-08-2018 08:10 AM
‎07-08-2018 08:10 AM
Following the update, l cannot connect to SAS Studio: Here is the error message:       "The application could not log on to the server "localhost:8591". No server is available at that port on that machine. Connection refused"   Any ideas? ... View more

Re: Non breaking spaces and formatting for gaps on an axis

by Quartz | Level 8 in ODS and Base Reporting
‎03-07-2018 04:41 PM
‎03-07-2018 04:41 PM
Thanks. I'm working with SAS Studio University Edition, and it is up to date: First problem is with Formatting: proc format; value $Namefmt "Age" ="Age" "BMI" ="BMI" "Sex_num" ="Sex (Male)" "CABG_i" ="CABG isolated" "Valve_i" ="Valve isolated" "CABG_V" ="CABG / Valve" "OHT" ="Heart Transplant" "Redo1" ="Redo" "IABP_ever1" ="Hx of IABP" "Hyperlip" ="Hyperlipidemia" "NIDDM1" ="NIDDM" "COPD" ="COPD" "CHF" ="CHF" "Afib" ="Afib" "CM" ="Cardiomyopathy" "CAD" ="CAD" "MI" ="Hx of MI" "AICD" ="AICD" "CRF" ="Renal Failure" "PVD" ="PVD" "ASA" ="Aspirin" "BB" ="Beta Blockers" "Steroids" ="Steroids" "Amio" ="Amiodarone" "LMWH" ="LMWH" "Digoxin" ="Digoxin" "AA1" ="Alpha Antagonists" "Statin1" ="Statins" "Nitrates1" ="Nitrates " "Diurectics1" ="Diuretics" "pre_Na" ="preop Na" "pre_BUN" ="preop BUN" "pre_Cr" ="preop Cr" "pre_Gluc" ="preop Gluc" "pre_Neutro" ="preop Neutrophils" "pre_Lymph" ="preop Lymphocytes" "pre_Hgb" ="preop Hb" "pRBC" ="PRBC" "FFP" ="FFP" "PLT" ="Platelets" "Cryo" ="Cryoprecipitate" "CPBmin" ="Bypass time" "AoXCmin" ="Cross Clamp time" "Ind_Temp1" ="induction Temp(C)" "Ind_HR1" ="ind HR" "Ind_SBP1" ="ind SBP" "Ind_MAP1" ="ind MAP" "Ind_DBP1" ="ind DBP" "Ind_PAS1" ="ind PAS" "Ind_mPAP1" ="ind mPAP" "Ind_PAD1" ="ind PAD" "pulse_p" ="ind pulse pressure" "A" ="." "B" ="." "C" ="."; run; data new; set HEERDT.VP_Table2b2; format Name $Namefmt.; run; I crated the missing values for A,B,C  so they would come out as empty lines in the subsequent tables. I also tried this with a data step to see whether it made any difference. The "Induction BP's" did not get formatted consistently. Some were continuous words, and some had a space after "Ind".   This was the code before attempting the formatted approach: proc template; define style mystyle; parent=styles.sasweb; class graphwalls / frameborder=off; class graphbackground / color=white; end; run; ods listing close; ods html style=mystyle path='.' file='newgraph.html'; ods graphics / reset=all border=off width=800 height=900; /* sd OR's */ /* alternate colorbands */ proc sgplot data =HEERDT.VP_Table2b2 nowall noborder ; Title "OR (95% CI) for covariates across Vasoplegia Severity" ; styleattrs datacontrastcolors=(CX18577A CX30CC6E ) ; scatter x=sdOR y=Name / xerrorupper=sdUCL xerrorlower=sdLCL noerrorcaps markerattrs=(symbol=circlefilled size =5 ) group=Severity groupdisplay=cluster clusterwidth=0.6; yaxis DISCRETEORDER= DATA reverse ; xaxis display =(noticks nolabel noline) VALUEATTRS=(Color=Black Family=Arial Size=8 Style=Normal Weight=Normal)type =log logbase =10; yaxis display =(noticks nolabel noline) VALUEATTRS=(Color=Black Family=Arial Size=8 Style=Normal Weight=Normal) valueshalign=left colorbands=odd colorbandsattrs=(transparency=0.7); refline 1.0 / axis=x discretethickness=1 lineattrs=(color=cxd0e0f0); yaxistable pdp2_char sdORCI /class =Severity classdisplay =stack position=right location=outside VALUEATTRS=(Color=Black Family=Arial Size=8 Weight=Normal ) VALUEHALIGN=CENTER VALUEJUSTIFY= CENTER CLASSORDER=DATA COLORGROUP=Severity; xaxis display=(nolabel) values=(0.4 1.0 1.5 2.0 2.5 5 ); run;  I added this in my first attempt:   proc sgplot data =HEERDT.VP_Table2b2 nowall noborder ; Title "OR (95% CI) across Vasoplegia Severity, units=1 s.d." ; styleattrs datacontrastcolors=(CX18577A CX30CC6E ) ; format Name $Namefmt.; .......(rest is the same as above).   Then I tried this valuesformat approach: proc sgplot data =Heerdt.Vasoplegia_space nowall noborder ; Title "OR (95% CI) for covariates across Vasoplegia Severity" ; styleattrs datacontrastcolors=(CX18577A CX30CC6E ) ; scatter x=sdOR y=Name / xerrorupper=sdUCL xerrorlower=sdLCL noerrorcaps markerattrs=(symbol=circlefilled size =5 ) group=Severity groupdisplay=cluster clusterwidth=0.6; yaxis DISCRETEORDER= DATA reverse valuesformat= $Namefmt.; (the rest being the same)   I do get the spaces, but the formatting steps do not result in an empty value on the y axis for "A,B, and C". Those letters do appear on the y axis, and there is an empty space across that level. So its frustrating that these letters are still there, and that the formatted values on that Y axis do not appear as requested.   The only difference  from the examples that Sanjay gave us were that these non breaking spaces had only one measure at that categorical axis. My values have two because im comparing OR's actoss two levels of severity of an outcome. Might just be that we cant do this in this locked down version of SAS. Sanjay Kindly agreed to look at my data and the code, so it will be interesting to see his solution.   Much appreciated, Christos      ... View more
Contact Me
Online Status
Offline
Date Last Visited
‎09-19-2020 03:15 PM
Likes given to
Likes from
View All
Quick links