BookmarkSubscribeRSS Feed
jpunch
Calcite | Level 5
I am trying to maximize a likelihood for independent founders using OPTMODEL (I have already written the code to perform a simple grid search but want something quicker). I have a dataset QTp with only one variable (also called QTp). Is the following code valid?

proc optmodel;

set S;

number t{S};

var q GE 0 init 0.5;

var muaa init 0;

var muab init 0;

var mubb init 0;

var sigma>=0.000001 init 0.05;

max logl= sum{i in S} LOG(1/(SQRT(2*3.14)*sigma))+
sum{i in S} LOG(((1-q)**2)*(EXP(-((t-muaa)**2)/(2*((sigma)**2)))+ (2*q*(1-q))*EXP(-((t-muab)**2)/(2*((sigma)**2)))+
((q)**2)*EXP(-((t-mubb)**2)/(2*((sigma)**2))));

read data QTp into S=[_n_] t=QTp;

solve with nlpc / tech=cgr printfreq=1 maxiter=5000;

print q sigma muaa muab mubb;

quit;

Message was edited by: jpunch

I also want q LE 1 but can't put it in the text window. Thanks.

Message was edited by: jpunch Message was edited by: jpunch
1 REPLY 1
There were several parentheses that were missing.
If you try the following code you should not get any syntax errors:

proc optmodel;

set S;

number t{S};

var q >= 0 <= 1 init 0.5;

var muaa init 0;

var muab init 0;

var mubb init 0;

var sigma>=0.000001 init 0.05;

max log= sum{i in S}
LOG(1/(SQRT(2*3.14)*sigma)
)
+
sum{i in S}
LOG(
((1-q)**2)*
(EXP(
-((t-muaa)**2)/(2*((sigma)**2)))
+
(2*q*(1-q))*EXP(
-((t-muab)**2)/(2*((sigma)**2))
)
+
((q)**2)*EXP(
-((t-mubb)**2)/(2*((sigma)**2))
)
)
);


read data QTp into S=[_n_] t=QTp;

solve with nlpc / tech=cgr printfreq=1 maxiter=5000;

print q sigma muaa muab mubb;

quit;

sas-innovate-2024.png

Available on demand!

Missed SAS Innovate Las Vegas? Watch all the action for free! View the keynotes, general sessions and 22 breakouts on demand.

 

Register now!

Multiple Linear Regression in SAS

Learn how to run multiple linear regression models with and without interactions, presented by SAS user Alex Chaplin.

Find more tutorials on the SAS Users YouTube channel.

Discussion stats
  • 1 reply
  • 818 views
  • 0 likes
  • 2 in conversation