thanks alot ... View more

If I want to save sas data set in a specific sheet in the excel file suppose in sheet3, what can I do? ... View more

the rest of code2: con c1: sum {i in ISET}(col1**2 *lambda/(lambda+k)**2)<=10; con c2: sum {i in ISET}(col2**2 *lambda/(lambda+k)**2)<=10; con c3: sum {i in ISET}(col3**2 *lambda/(lambda+k)**2)<=10; con c4: (max_lamb+k)/(min_lamb+k)<= max_lamb/min_lamb; solve with nlpc / tech =qne; print k ; ... View more

the rest of code 2: con co1: sum {i in ISET}(col1**2 *lambda/(lambda+k)**2)<=10; con co2: sum {i in ISET}(col2**2 *lambda/(lambda+k)**2)<=10; con co3: sum {i in ISET}(col3**2 *lambda/(lambda+k)**2)<=10; con co4: (max_lamb+k)/(min_lamb+k)<= max_lamb/min_lamb; solve with nlpc / tech =qne; print k ; ... View more

I want my code using proc iml (code1) to be the same as the one using proc optmodel (code2), Is they are the same ? as they give me different results for K. If they are not the same, What is missing in code 1 to be the same as code 2? moreover, in code 2 we can tell SAS that K is positive (var k >=0;) so, I want to put K as positive variable in code 1 as well, How can I do so? code1: Proc iml; start fun(K); sumf=0; do i = 1 to 3; sumf = sumf + (estvar*lambda)+ (g* k)**2 /(lambda+ k)**2; end; return (sumf); finish fun; start con(k); c=j(4,1,0); sumc1=0; do i = 1 to 3; sumc1 = sumc1 + (col1**2 *lambda/(lambda+k)**2); end; c[1]= 10-sumc1; sumc2=0; do i= 1 to 3; sumc2 = sumc2 + (col2**2 *lambda/(lambda+k)**2); end; c[2]= 10-sumc2; sumc3=0; do i= 1 to 3; sumc3 = sumc3 + (col3**2 *lambda/(lambda+k)**2); end; c[3]= 10-sumc3; c[4]=(max(lambda)/min(lambda))-(max(lambda)+k)/(min(lambda)+k); return (c); finish con; optn=j(1,11,.); optn[4]=1; optn[6]=2; optn[11]=0; optn[10]=4; call nlpqn (rc, kres, "fun", k, optn) nlc="con"; code2: proc optmodel; var k >=0; set ISET; num lambda {ISET}; num g {ISET}; num col1 {ISET}; num col2 {ISET}; num col3 {ISET}; num min_lamb = min {i in ISET} lambda; num max_lamb = max {i in ISET} lambda; read data sasdata.sxlamb2 into ISET = [_n_] lambda g; read data work.datavif into ISET = [_n_] col1 col2 col3; num estvar ; read data sasdata.m3data into estvar; min f = sum {i in ISET}((estvar*lambda)+(k**2 *g **2)/(lambda+k)**2); con co1: sum {i in ISET}(col1**2 *lambda/(lambda+k)**2)<=10; con co2: sum {i in ISET}(col2**2 *lambda/(lambda+k)**2)<=10; con co3: sum {i in ISET}(col3**2 *lambda/(lambda+k)**2)<=10; con co4: (max_lamb+k)/(min_lamb+k)<= max_lamb/min_lamb; solve with nlpc / tech =qne; print k ; ... View more

I want my code using proc iml (code1) to be the same as the one using proc optmodel (code2), Is they are the same ? as they give me different results for K. If they are not the same, What is missing in code 1 to be the same as code 2? moreover, in code 2 we can tell SAS that K is positive (var k >=0;) so, I want to put K as positive variable in code 1 as well, How can I do so? code1: Proc iml; start fun(K); sumf=0; do i = 1 to 3; sumf = sumf + (estvar*lambda)+ (g* k)**2 /(lambda+ k)**2; end; return (sumf); finish fun; start con(k); c=j(4,1,0); sumc1=0; do i = 1 to 3; sumc1 = sumc1 + (col1**2 *lambda/(lambda+k)**2); end; c[1]= 10-sumc1; sumc2=0; do i= 1 to 3; sumc2 = sumc2 + (col2**2 *lambda/(lambda+k)**2); end; c[2]= 10-sumc2; sumc3=0; do i= 1 to 3; sumc3 = sumc3 + (col3**2 *lambda/(lambda+k)**2); end; c[3]= 10-sumc3; c[4]=(max(lambda)/min(lambda))-(max(lambda)+k)/(min(lambda)+k); return (c); finish con; optn=j(1,11,.); optn[4]=1; optn[6]=2; optn[11]=0; optn[10]=4; call nlpqn (rc, kres, "fun", k, optn) nlc="con"; code2: proc optmodel; var k >=0; set ISET; num lambda {ISET}; num g {ISET}; num col1 {ISET}; num col2 {ISET}; num col3 {ISET}; num min_lamb = min {i in ISET} lambda; num max_lamb = max {i in ISET} lambda; read data sasdata.sxlamb2 into ISET = [_n_] lambda g; read data work.datavif into ISET = [_n_] col1 col2 col3; num estvar ; read data sasdata.m3data into estvar; min f = sum {i in ISET}((estvar*lambda)+(k**2 *g **2)/(lambda+k)**2); con co1: sum {i in ISET}(col1**2 *lambda/(lambda+k)**2)<=10; con co2: sum {i in ISET}(col2**2 *lambda/(lambda+k)**2)<=10; con co3: sum {i in ISET}(col3**2 *lambda/(lambda+k)**2)<=10; con co4: (max_lamb+k)/(min_lamb+k)<= max_lamb/min_lamb; solve with nlpc / tech =qne; print k ; ... View more

Another question, I want c[1], c[2], c[3], and c[4] which are my 4 constraints to be greater than or equal zero, is this correct in my code? ... View more

Is there any option to increase the decimal places of k optimal inside call statement? ... View more

what is the optimal k value? is it kres in my code? ... View more