Re: Here is the problem with DiscreteMarkovProcess...
- To: mathgroup at smc.vnet.net
- Subject: [mg129338] Re: Here is the problem with DiscreteMarkovProcess...
- From: bhuvaneshbhatt at gmail.com
- Date: Sat, 5 Jan 2013 02:20:25 -0500 (EST)
- Delivered-to: l-mathgroup@mail-archive0.wolfram.com
- Delivered-to: l-mathgroup@wolfram.com
- Delivered-to: mathgroup-newout@smc.vnet.net
- Delivered-to: mathgroup-newsend@smc.vnet.net
- References: <kc2phn$md5$1@smc.vnet.net>
`
On Wednesday, January 2, 2013 8:14:15 PM UTC-6, Richard Palmer wrote:
> The code
>
>
>
> len=3;
>
> v=RandomInteger[1000,len];v=v/Total[v];v
>
> Total[v]
>
> tm = RandomReal[1, {len, len}];
>
> tm = #/Total[#] & /@ tm;
>
> dmp = DiscreteMarkovProcess[v, tm]
>
> Graph[dmp]
>
> MarkovProcessProperties[dmp]
>
>
>
> works as advertised.
>
>
>
> Change len=7 and the code does not work.
>
> On my system, len=4 works, but len=5 does not.
>
>
>
> What is going on?
>
This was a bug, and it has been fixed and should work correctly in an upcoming release. In the meantime, here is a workaround:
Unprotect[DiscreteMarkovProcess];
ProcessParameterQ[dmp] ^= True;
Protect[DiscreteMarkovProcess];
I apologize for the inconvenience.
Bhuvanesh
Wolfram Research