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: 
Miracle
‎02-26-2024
  • 315 Posts
  • 136 Likes Given
  • 13 Solutions
  • 50 Likes Received

Errors from Proc iml

by Barite | Level 11 in SAS Programming
‎03-02-2015 12:53 AM
‎03-02-2015 12:53 AM
Dear all, How are you? I'm trying to run the code from this paper http://support.sas.com/resources/papers/proceedings14/2042-2014.pdf with my own data(8 5-level Likert items) but there is errors. Variable 'sumcor' seems like causing the problem. According to the paper, "Note that ∑ ∑ is just the sum of all values in the correlation matrix, since in a correlation matrix the diagonal elements are all 1, they sum to be the number of elements and the double sum adds up all off diagnonal elements (the correlations), note that each one would be repeated twice given the symmetry in the matrix." I'm a novice user of proc iml and your help is greatly appreciated.  I have printed out a partial log output FYI. LOGIC(CALCALPHA.VARLIST):  %DO loop index variable I is now 8; loop will iterate again. SYMBOLGEN:  Macro variable I resolves to 8 MLOGIC(CALCALPHA.VARLIST):  %DO loop index variable I is now 9; loop will not iterate again. MPRINT(CALCALPHA.VARLIST):  rAttitudes01 rAttitudes02 rAttitudes03 rAttitudes04 rAttitudes05 rAttitudes06 rAttitudes07 rAttitudes08 MLOGIC(CALCALPHA.VARLIST):  Ending execution. MPRINT(CALCALPHA):   } into itemcor; SYMBOLGEN:  Macro variable QS resolves to rAttitudes MPRINT(CALCALPHA):   section="rAttitudes"; SYMBOLGEN:  Macro variable CORDAT resolves to corrmatrix MPRINT(CALCALPHA):   datused="corrmatrix"; MPRINT(CALCALPHA):   numbitem=trace(itemcor); MPRINT(CALCALPHA):   one=j(numbitem,1); MPRINT(CALCALPHA):   print 'The Input Correlation'; MPRINT(CALCALPHA):   print itemcor; MPRINT(CALCALPHA):   print 'Number of items' numbitem; MPRINT(CALCALPHA):   sumcor=(one)*itemcor*one; ERROR: (execution) Matrices do not conform to the operation. operation : * at line 4010 column 1 operands  : one, itemcor one      7 rows      1 col     (numeric)          1          1          1          1          1          1          1 itemcor      8 rows      8 cols    (numeric)   0.999999 0.7418163 0.6005298 0.6328579 0.2149593 -0.085045 0.1591052 0.1087784 0.7418163  0.999999 0.6696118  0.639722 0.1674687 -0.127974 0.0821654 0.0697031 0.6005298 0.6696118  0.999999 0.8357702 0.2290657 -0.118532  0.136779    0.1264 0.6328579  0.639722 0.8357702  0.999999 0.2232193 -0.166345 0.1338478  0.137883 0.2149593 0.1674687 0.2290657 0.2232193  0.999999 -0.136473 -0.047612 0.1070489 -0.085045 -0.127974 -0.118532 -0.166345 -0.136473  0.999999 0.0824758 -0.076128 0.1591052 0.0821654  0.136779 0.1338478 -0.047612 0.0824758  0.999999  0.080517 0.1087784 0.0697031    0.1264  0.137883 0.1070489 -0.076128  0.080517  0.999999 statement : ASSIGN at line 4010 column 1 MPRINT(CALCALPHA):   alpha=(numbitem/(numbitem-1))*(1-(numbitem / sumcor)); ERROR: (execution) Matrix has not been set to a value. operation : / at line 4010 column 1 operands  : numbitem, sumcor numbitem      1 row       1 col     (numeric)   7.999992 sumcor      0 row       0 col     (type ?, size 0) statement : ASSIGN at line 4010 column 1 MPRINT(CALCALPHA):   print alpha; ERROR: Matrix alpha has not been set to a value. statement : PRINT at line 4010 column 1 MPRINT(CALCALPHA):   L = eigval(itemcor); MPRINT(CALCALPHA):   theta=(numbitem/(numbitem-1))*(1-1/L[1]); MPRINT(CALCALPHA):   print 'Theta'; MPRINT(CALCALPHA):   print theta; SYMBOLGEN:  Macro variable OUT resolves to alpha MPRINT(CALCALPHA):   create alpha var {section datused alpha theta}; MPRINT(CALCALPHA):   append; SYMBOLGEN:  Macro variable OUT resolves to alpha MPRINT(CALCALPHA):   close alpha; NOTE: The data set WORK.ALPHA has 1 observations and 4 variables. MPRINT(CALCALPHA):   quit; NOTE: Exiting IML. NOTE: The SAS System stopped processing this step because of errors. NOTE: PROCEDURE IML used (Total process time):       real time           0.06 seconds       cpu time            0.06 seconds MPRINT(CALCALPHA):   run; MLOGIC(CALCALPHA):  Ending execution. ... View more

How to join the points?

by Barite | Level 11 in Graphics Programming
‎02-22-2015 07:50 PM
‎02-22-2015 07:50 PM
Dear everyone, How are you? Can I please ask how do I join the 2 points (amount in 2000 & 2001) of the same individual in each level of the x-axis? Thank you for your help. FYI here is the code. ods graphics / antialias=off; proc sgplot data=tmp;   scatter x=item01 y=amt / group=year groupdisplay=cluster clusterwidth=0.5;   xaxis type=discrete; run; ods graphics off; ... View more

Re: Removing select cases from a dataset

by Barite | Level 11 in SAS Procedures
‎02-04-2015 10:51 PM
3 Likes
‎02-04-2015 10:51 PM
3 Likes
Try this, proc sort data=a; by ID Incident_number; run; proc sort data=b; by ID Incident_number; run; data want; merge a(in=a) b(in=b); by ID Incident_number; if a=b; drop VarA VarB; run; proc print noobs; run; ... View more

Proportion test on repeated measures data.

by Barite | Level 11 in Statistical Procedures
‎11-04-2014 12:20 AM
‎11-04-2014 12:20 AM
Dear all. Wish you all well. I wish to find out if the proportion of yes to problem is the same between time 0 and 1 given that they had used(used=1). The data structure is as id   time    used    problem 1    0        0          1 1    1        1          0 2    0        .           0 2    1        1          0 ... etc My initial thought was proc genmod data=tmp descending; class id time used; model problem=time used time*used / dist=bin link=identity; repeated sub=id / type=ind; lsmeans time*used / ilink cl diff; run; Is this the right approach? Does lsmean give me the right proportion and the proportion difference? :smileyconfused: Because I was reading this 46997 - Estimating the risk (proportion) difference for matched pairs data with binary response. and it says "Unlike the default logit link function, the identity link does not ensure that the model produces valid probability estimates. Errors may be result when fitting such models depending on the model and the data." Thank you for your help ... View more

Re: loop or array?

by Barite | Level 11 in SAS Programming
‎11-03-2014 11:37 PM
‎11-03-2014 11:37 PM
Thanks everyone for your valuable input ... View more

loop or array?

by Barite | Level 11 in SAS Programming
‎10-28-2014 02:16 AM
‎10-28-2014 02:16 AM
Dear all, Wishing you all well. Can I please ask how do I rewrite my sas code below so I don't have to repeat the if-else statements 32 times? :smileyconfused: Thank you very much   q71l_type1=0;   if q71l_a_q1_1=1 or q71l_a_q1_1=1 or q71l_b_q1_1=1 or q71l_c_q1_1=1 or q71l_aa_q1_1=1 or q71l_d_q1_1=1 or  q71l_e_q1_1=1 or           q71l_f_q1_1=1 or q71l_g_q1_1=1 or q71l_h_q1_1=1 or q71l_i_q1_1=1 or q71l_j_q1_1=1 or q71l_k_q1_1=1 or q71l_cc_q1_1=1 or     q71l_n_q1_1=1 or q71l_o_q1_1=1 or q71l_dd_q1_1=1 or  q71l_p_q1_1=1 or q71l_q_q1_1=1 or q71l_r_q1_1=1 or q71l_s_q1_1=1 or     q71l_bb_q1_1=1 then q71l_type1=1;   if q71l_a_q1_1=. and q71l_a_q1_1=. and q71l_b_q1_1=. and  q71l_c_q1_1=. and q71l_aa_q1_1=. and q71l_d_q1_1=. and q71l_e_q1_1=.     and q71l_f_q1_1=. and q71l_g_q1_1=. and q71l_h_q1_1=. and q71l_i_q1_1=. and q71l_j_q1_1=. and q71l_k_q1_1=. and q71l_cc_q1_1=.     and q71l_n_q1_1=. and q71l_o_q1_1=. and q71l_dd_q1_1=. and  q71l_p_q1_1=. and q71l_q_q1_1=. and q71l_r_q1_1=. and                       q71l_s_q1_1=. and q71l_bb_q1_1=. then q71l_type1=.;       up to   q71l_type32=0;   if q71l_a_q1_32=1 or q71l_a_q1_32=1 or q71l_b_q1_32=1 or q71l_c_q1_32=1 or q71l_aa_q1_32=1 or q71l_d_q1_32=1 or                              q71l_e_q1_32=1 or q71l_f_q1_32=1 or q71l_g_q1_32=1 or q71l_h_q1_32=1 or q71l_i_q1_32=1 or q71l_j_q1_32=1 or q71l_k_q1_32=1 or      q71l_cc_q1_32=1 or  q71l_n_q1_32=1 or q71l_o_q1_32=1 or q71l_dd_q1_32=1 or  q71l_p_q1_32=1 or q71l_q_q1_32=1 or                          q71l_r_q1_32=1 or q71l_s_q1_32=1 or q71l_bb_q1_32=1 then q71l_type32=1;   if q71l_a_q1_32=. and q71l_a_q1_32=. and q71l_b_q1_32=. and  q71l_c_q1_32=. and q71l_aa_q1_32=. and q71l_d_q1_32=. and                  q71l_e_q1_32=.    and q71l_f_q1_32=. and q71l_g_q1_32=. and q71l_h_q1_32=. and q71l_i_q1_32=. and q71l_j_q1_32=. and                    q71l_k_q1_32=. and q71l_cc_q1_32=.    and q71l_n_q1_32=. and q71l_o_q1_32=. and q71l_dd_q1_32=. and  q71l_p_q1_32=. and              q71l_q_q1_32=. and q71l_r_q1_32=. and  q71l_s_q1_32=. and q71l_bb_q1_32=. then q71l_type32=.; ... View more

Re: difference in proportions

by Barite | Level 11 in Statistical Procedures
‎10-12-2014 08:48 AM
‎10-12-2014 08:48 AM
Hi Keshan. Thanks for your reply. By the way, I thought I've already moved this post to SAS Statistical Procedures. Have I not? :smileyconfused: ... View more

Re: difference in proportions

by Barite | Level 11 in Statistical Procedures
‎10-11-2014 10:33 PM
‎10-11-2014 10:33 PM
Hi Keshan. Thanks for your reply. Actually, maybe I need to do a McNemar's due to the correlated proportions of before and after from individuals. I thought I need to create dummy variables for the two 10-level categorical variables in order to test for the before and after proportion for each level. If I just did a Fisher's exact test between the two 10-level freq_before and 10-level freq_after variable, it will not tell me which level of before and after proportions are different. Am I making any sense here? :smileyconfused: Any comment and insight is greatly appreciated ... View more

Re: difference in proportions

by Barite | Level 11 in Statistical Procedures
‎10-10-2014 09:30 PM
‎10-10-2014 09:30 PM
Hi Keshan. Thanks for your reply. My data looks like this. The sample size is small with only about 100 individuals. id freq_before freq_after 1 10 1 2 5 2 3 6 8 . . . . . . . . . 100 2 3 where 10=everyday, 9=once every few days, ...., 1=never for both freq_before and freq_after variable. I was wondering if I should create dummy variables for the two 10-level categorical variables and do a Chi-square test/Fisher's exact test like this? :smileyconfused: proc freq data=tmp; tables before_level1*after_level1 / chisq; run; . . . proc freq data=tmp; tables before_level10*after_level10 / chisq; run; Any comments and insight is greatly appreciated ... View more

difference in proportions

by Barite | Level 11 in Statistical Procedures
‎10-09-2014 09:32 PM
‎10-09-2014 09:32 PM
Dear all, I wish you all well. I have two(freq_b4 and freq_after) 10-level category variables which I would like to test if the proportions of each level from the 2 variables differ. Am I doing it right? Your insight is greatly appreciated. My SAS code is as below:  proc freq data=tmp compress; exact binomial; tables freq_b4 / binomial(level=1 p=0.05 ); . . . tables freq_b4 / binomial(level=10 p=0.1 ); run; ... View more

Re: Longitudinal weighting

by Barite | Level 11 in Statistical Procedures
‎10-04-2014 08:16 AM
‎10-04-2014 08:16 AM
Dear Keshan and Steve, Thanks for your reply. I wish to analyse only the 2013 data and to create weighting for the 2013 only participants so it will be a nationally representative of the population. However, I'm not sure how do I go about doing the weighting given that there are only drinkers in the 2013 data. Thank you. ... View more

rake-trim sas macro

by Barite | Level 11 in SAS Programming
‎09-27-2014 10:27 AM
‎09-27-2014 10:27 AM
Hi, wish you all well. Can I please ask where can I download the rake-trim sas macro by Louis Rizzo? Or the %RAKE_AND_TRIMM SAS MACRO as in Extreme Survey Weight Adjustment as a Component of Sample Balancing (a.k.a. Raking) by David Izrael and et al? Thank you very much. ... View more

Longitudinal weighting

by Barite | Level 11 in Statistical Procedures
‎09-23-2014 01:39 AM
‎09-23-2014 01:39 AM
Dear all, Wish you all well.  My apology that this is not SAS-related. I have to do weighting for an alcohol longitudinal study(2011-2013). How do I go about getting the final weights? Some people stopped responding as the survey progressed leaving only the drinkers responded in the end(2013). Does it sound statistically correct by multiplying the sampling weights by postratification weights(weighted on census by gender, age and ethnicity) then by attrition weights to get the final weights? Any input will be greatly appreciated ... View more

Re: Proc Logistic - Class statement

by Barite | Level 11 in SAS Procedures
‎06-30-2014 01:21 AM
‎06-30-2014 01:21 AM
Just a suggestion, maybe you could change your proc format like this so you don't need to do the if-else statement. proc format; Value pkyrcat3f 0='never cig smk'                 1-19 ='<20 pkyrs'                 20-39='20-<40'                 40-high='>=40'; proc logistic data=pan2.pkyr; /*1=case, 2= control*/ format packy pkyrcat3f.; class age10yg(ref=first)sex race educat (ref=last) bmiad3cat (ref=first) etoh(ref=first) diabetes(ref=last) pancreat(ref=last) packy (ref='never cig smk') / param=ref; model status= age10yg sex race educat bmiad3cat etoh diabetes pancreat hisp gb_stones packy / selection = forward include=11 ctable pprob= (0 to 1 by .1) outroc=roc; run; quit; ... View more

Re: Compare SAS programs

by Barite | Level 11 in SAS Programming
‎06-05-2014 11:44 PM
‎06-05-2014 11:44 PM
Thanks everyone for your valuable contribution i shall try them all next time. In the mean time, I'm using MsWord to do the comparison. ... View more
Contact Me
Online Status
Offline
Date Last Visited
‎02-26-2024 06:33 PM
Likes given to
Likes from
Quick links