Hello, I am a graduate student (in ecology, so not a great background in mathematics or coding languages) attempting an optimization for a part of my thesis.  Currently I have the coding below in an attempt to get the "program" to work with pseudo-data.  This has been constructed based on publications, which use other programming languages, but the functions are "...local minimizer for smooth non-linear functions subject to bound-constrained parameters, and uses a quasi-newton method". Thus I believe the correct call is NLPQN, however when running this I do not get the expected optimization. Proc IML; start qfasa(y);     prey={17 41 5, 40 30 25}; /*pseudo-profile of three prey items with 2 variables each*/     predator={13.4, 29.5}; /*pseudo-profile of one predator with 2 variables that are the proportional addition of the prey items*/     predic=(prey#y);     diet=predic[,+];     KL1=((predator-diet)#log10(predator/diet));  /*Kulback-Liebler Distance*/     KL=KL1[+,];      Return(KL);  /*currently must be a scalar, possible to get this for each predator column when more are added? (do loops?)*/ finish qfasa; y={.5 .25 .25};  /* Initial guess at proportions*/ c={0 0 0, 1 1 1}; /*Constraints on parameters*/ optn={0 2}; call nlpqn(rc,xr,"qfasa",y,optn,c) ; run; */Answer should be { .25 .15 .60}*/ Thanks for any help! ... View more

Hello- I am attempting to divide the values in rows of matrix A by the values in rows of matrix B. the columns of both are the same variables, the rows are different though and the number of rows are different. I have create a scaled down version below of what I am trying to do: A{ 1 3, 4 5} B{ 5 9, 6 2, 7 1} C{ .2 .33, .166 1.5, .14 3, .8 .55, .66 2.5, .57 5} Simply using C=A/B doesn't work, But I am wondering what does? Or if I need to use something other than IML? ... View more

Hey guys, I can't seem to find this anywhere online, but R directions are out there so i figured there must exist a way for SAS to do it too. Basically I have 2 quantitative variables and one classification variable. I am looking to be able to make a plot of the means of both Q-variables and their respective confidence intervals for each class. Is there something else these are called?  Its basically a scatter plot but I cant seem to get it to work correctly (Below is what I am trying to make) Any direction is appreciated. Thanks, Austin ... View more

Hello, I am trying to do a simple merge of two data sets. The first one has all of the data needed, the second one has a few missing areas that I would like to have filled with zeros Say this is data set 1: input id spp count; cards; 1w 250 16 1w 322 30 1w 99 2 2r 10 1 2r 99 1 ; And this is the second data set: input id prey amount; cards; 1w 250 10 1w 322 15 2r 10 1 ; What I need it to look like is: input id spp prey count amount; cards; 1w 250 250 16  10 1w 322 322 30  15 1w  99   99   2    0 2r  10   10   1    1 2r  99   99   1    0 ; The key being where the item is not present in the second set is is given a zero as later I will be obataining the difference between count and amount and continuing on with my manipulation.  I just cant seem to get my head around how to get this to work the way I need to. Thanks for your time! ... View more