Turn on suggestions

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

Showing results for

Options

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

- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content

Posted 08-26-2010 05:07 AM
(849 views)

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

1 REPLY 1

- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content

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

**Available on demand!**

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

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.