Re: Minimization with constraint expresses as CDF[] >=
- To: mathgroup at smc.vnet.net
- Subject: [mg72490] Re: [mg72474] Minimization with constraint expresses as CDF[] >=
- From: "Chris Chiasson" <chris at chiasson.name>
- Date: Mon, 1 Jan 2007 04:15:10 -0500 (EST)
- References: <200612311007.FAA06784@smc.vnet.net>
The encoding in your message is not showing up correctly on my screen. Would you mind using only the ASCII character set? I ask because I am interested in your problem. On 12/31/06, Pratim Vakish <pratim_usc at hotmail.com> wrote: > > Hello, > > I have the following optimization problem: a minimization subject to a few linear constraints and one chance constarint formulated using CDF[] > > Needs["Statistics`MultinormalDistribution`"]; > r={{0.01,-0.016},{-0.016,0.018}}; > ndist=MultinormalDistribution[{1,3},r]; > NMinimize[{2x1+x2, > x1?0.8, > x2?2.5, > x1?0, > x2?0, > x1+x2==y, > CDF[ndist,{x1,y}]?0.9, > PrecisionGoal -> 3}, > {x1,x2,y}] > > > I keep getting the following error message. I do not understand why. Any advice is welcome. > > NMinimize::"bcons": > "The following constraints are not valid: y == x1 + x2, x1 ? 0, x2 ? 0, CDF[(MultinormalDistribution[({1, 3},({0.01,-0.016},{-0.016, 0.018}, {x1, y}] ? 0.9, > x1 ? 0.8, x2 ? 2.5}. > Constraints should be equalities, inequalities, or domain specifications involving the variables. > > > > ---------------------------------------- > > Date: Sat, 25 Nov 2006 09:53:24 -0500 > > From: hanlonr at cox.net > > To: pratim_usc at hotmail.com; mathgroup at smc.vnet.net > > Subject: Re: [mg71632] Numerical vaule of multinormal distribution > > > > Needs["Statistics`MultinormalDistribution`"]; > > r = {{0.8, 0.1}, {0.1, 0.2}}; > > ndist = MultinormalDistribution[{1, 0.4}, r]; > > CDF[ndist, {-1, 1}] > > > > 0.0124156 > > > > I would guess that you tried to use to use CDF or MultinormalDistribution before loading the package. Evaluate after starting with a fresh kernel. > > > > > > Bob Hanlon > > > > ---- Pratim Vakish <pratim_usc at hotmail.com> wrote: > > > > > > I have a basic question. > > > I want toobtain the numerical value of the cumulative normal probability distribution function at some points. > > > I have the following code to do this : > > > > > > Off[General::spell1]; > > > Needs["Statistics`MultinormalDistribution`"]; > > > r = {{0.8, 0.1}, {0.1, 0.2}}; > > > ndist = MultinormalDistribution[{1, 0.4}, r]; > > > CDF[ndist, {-1, 1}] > > > > > > > > > The output I obtain is: > > > > > > CDF[MultinormalDistribution[{1, 0.4}, {{0.8, 0.1}, {0.1, 0.2}}], {-1, 1}] > > > > > > I would like to have the numerical value of this expression. > > > How could I do? > > > > > > > > > Pratim > > > > _________________________________________________________________ > Get the Live.com Holiday Page for recipes, gift-giving ideas, and more. > www.live.com/?addtemplate=holiday > > -- http://chris.chiasson.name/