Re: problems with Nmaximization
- To: mathgroup at smc.vnet.net
- Subject: [mg67828] Re: [mg67815] problems with Nmaximization
- From: Bob Hanlon <hanlonr at cox.net>
- Date: Sun, 9 Jul 2006 04:51:03 -0400 (EDT)
- Reply-to: hanlonr at cox.net
- Sender: owner-wri-mathgroup at wolfram.com
m={75,36,12,4,2,0,0,0,0}; n={571-75,338-36,173-12,93-4,41-2,21,10,6,4}; t=Range[0,Length[m]-1]; s=Simplify[Total[ m*Log[(h*p*(1-p)^t)/(1-h+h*(1-p)^t)]- n*Log[1-(h*p*(1-p)^t)/(1-h+h*(1-p)^t)]]] eps=10^-6; Plot3D[s, {h, eps, 1-eps}, {p, eps, 1-eps}]; Limit[Limit[s, h -> 1], p -> 1] Infinity Bob Hanlon ---- narsi <narasimhan.sowmyanarayanan at gmail.com> wrote: > I am trying to use nmax to evaluate a function using Nmaximize. The > function is as follows. > > \!\(m[1] = 75\[IndentingNewLine] > m[2] = 36\[IndentingNewLine] > m[3] = 12\[IndentingNewLine] > m[4] = 4\[IndentingNewLine] > m[5] = 2\[IndentingNewLine] > m[6] = 0\[IndentingNewLine] > m[7] = 0\[IndentingNewLine] > m[8] = 0\[IndentingNewLine] > m[9] = 0\[IndentingNewLine] > n[1] = 571 - 75\[IndentingNewLine] > n[2] = 338 - 36\[IndentingNewLine] > n[3] = 173 - 12\[IndentingNewLine] > n[4] = 93 - 4\[IndentingNewLine] > n[5] = 41 - 2\[IndentingNewLine] > n[6] = 21\[IndentingNewLine] > n[7] = 10\[IndentingNewLine] > n[8] = 6\[IndentingNewLine] > n[9] = 4\[IndentingNewLine]\[IndentingNewLine] > NMaximize[{âË?â??\+\(t = 1\)\%9\((m[ > t]*Log[\((\(h\ *p*\((1 - > p)\)\^\((t - 1)\)\)\/\(1 - h + h*\((1 - p)\)\^\((t - 1)\)\))\)] - > n[t]*Log[1 - \((\(h\ *p*\((1 - p)\)\^\((t - 1)\)\)\/\(1 - h + > h*\((1 - p)\)\^\((t - 1)\)\))\)])\), 1 ââ?°Â¥ h ââ?°Â¥ 0, 1 ââ?°Â¥ p ââ?°Â¥ 0}, > {h, p}]\) > > I cannot get this to converge. It gives me the following message > > \!\(ower::"infy" \(\(:\)\(\ \)\) " > Infinite expression \!\(1\/0.`\) encountered"\) > > > Can someone please help with this error in mathematica. I dont know > what the problem is. > > Thanks > > narsi >