SAS Optimization, and SAS Simulation Studio

turn on suggestions

Auto-suggest helps you quickly narrow down your search results by suggesting possible matches as you type.

Showing results for

Find a Community

Topic Options

- Subscribe to RSS Feed
- Mark Topic as New
- Mark Topic as Read
- Float this Topic for Current User
- Bookmark
- Subscribe
- Printer Friendly Page

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Report Inappropriate Content

08-26-2010 05:07 AM

Hi everyone,

I have SAS 9.2.

This is more of a math side question hope someone has an answer.

I'm working with acceptance sampling plans specifically double sampling by variables.

So I have nonlinear function binomial pdf and which is iterated and calculated according to a array of probabilities [0..1].

I would like to equalize (I mean not exact match, equivalent) the result of my binomial pdf to bivariate normal distribution to find my 3 decision variables which are used to calculate bivariate function. As an example:

do p=0.0001 to 0.9999 by 0.0001;

Paatt = probbnml(p,n,c);

Zp = probit(1-p);

Pa1interup=sqrt(n1)*(Zp-k1);

Pa2interinup=sqrt(n1)*(Zp-k2);

Pa2interindwn=sqrt(n1)*(Zp-k1);

Pa2interoutup=sqrt(n1+n2)*(Zp-k3);

PA1 = probnorm(sqrt(n1)*(Zp-k1));

PAa = probbnrm(sqrt(n1+n2)*(Zp-k3), sqrt(n1)*(Zp-k2), n1/(n1+n2));

PAb = probbnrm(sqrt(n1+n2)*(Zp-k3), sqrt(n1)*(Zp-k1), n1/(n1+n2));

PA2 = PAa - PAb;

PA = PA1 + PA2;

end;

I would like to make the*PA* equivalent to *Pattr* by changing *k1,k2,k3* for the whole nonlinear function. So its more of a approximation of k values from the result of *Pattr* vector. Any information may be helpful.

Thanks

I have SAS 9.2.

This is more of a math side question hope someone has an answer.

I'm working with acceptance sampling plans specifically double sampling by variables.

So I have nonlinear function binomial pdf and which is iterated and calculated according to a array of probabilities [0..1].

I would like to equalize (I mean not exact match, equivalent) the result of my binomial pdf to bivariate normal distribution to find my 3 decision variables which are used to calculate bivariate function. As an example:

do p=0.0001 to 0.9999 by 0.0001;

Paatt = probbnml(p,n,c);

Zp = probit(1-p);

Pa1interup=sqrt(n1)*(Zp-k1);

Pa2interinup=sqrt(n1)*(Zp-k2);

Pa2interindwn=sqrt(n1)*(Zp-k1);

Pa2interoutup=sqrt(n1+n2)*(Zp-k3);

PA1 = probnorm(sqrt(n1)*(Zp-k1));

PAa = probbnrm(sqrt(n1+n2)*(Zp-k3), sqrt(n1)*(Zp-k2), n1/(n1+n2));

PAb = probbnrm(sqrt(n1+n2)*(Zp-k3), sqrt(n1)*(Zp-k1), n1/(n1+n2));

PA2 = PAa - PAb;

PA = PA1 + PA2;

end;

I would like to make the

Thanks

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Report Inappropriate Content

Posted in reply to bermur

09-03-2010 06:17 PM

This forum is for Operations Research and Mathematical Optimization.

Maybe more relevant information can be found in the Statistics forum at the following link:

http://support.sas.com/forums/forum.jspa?forumID=46

Maybe more relevant information can be found in the Statistics forum at the following link:

http://support.sas.com/forums/forum.jspa?forumID=46