MaplePrimes - answers and comments on Question, Formulate equation with different classes (categories) & subscripts

Make a procedure

Adri van der Meer — Tue, 13 Nov 2012 21:32:01 Z

restart;
para := proc(i,new)
  # new = 0 or 1
  seq(log(p[j,i]/(1-p[j,i]))=mu+tau[j]+eta[i] + new*tau[j]*eta[i],j=2..4)
end proc:
ans:=solve({para(1,1)},[seq(p[j,1],j=2..4)]);

subs( op(ans), p[2,1] );

assign(op(ans));

Taking it further

casperyc — Wed, 14 Nov 2012 02:11:03 Z

Hi Adri,

Thank you for you answer. That was very helpful.

Now, the next thing to consider is that I want to computer a 'weighted' mean.

See my worksheet. new2.mw

In the end of the file,

>par:=w1*parC1+w2*parC2+(1-w1-w2)*parC3;

is the desired output.

I have tried to write a procedure ( getpVec ) for it, but I am stuck. I think it might be a problem of declaring variables.

Since the p[i,j] are NOT declared in 'local' in getpVec.

Any help would be greatful.

Many thanks.

Casper

Try this

Adri van der Meer — Wed, 14 Nov 2012 20:49:46 Z

Casper, I'm not sure if I have all the indices right. But I think this may help you further.

restart:
para := proc(C,new)
   global p;
    # new = 0 or 1
    seq(log(p[j,C]/(1-p[j,C]))=mu+tau[j]+eta[C] + new*tau[j]*eta[C],j=2..4)
end proc:

getpVec:=proc(k,c,new)
  local i,j, parC;
  global p, wp;
  for i to c do
    solve({para(i,new)},[seq(p[j,i],j=2..(k+1))]);
    parC[i] := Vector(k, [rhs~(op(%))] );
  end do;
  add(wp[i]*parC[i],i=1..(c-1))+(1-add(wp[i],i=1..(c-1)))*parC[c];
end proc:
getpVec(3,3,1);