How to calculate the (odd ,even) condition in the following Poisson distribution for p_z parameter? i need some help in the code bellow

when q=[0.1,0.2,0.3,0.4,0.5] and p =1-q, z=6

[Torsten]

By "Solve the equation" you mean: "Calculate P_z" ?

[Walter Roberson]

Solve it for what parameter ?

[Sofy]

clc;
 format long
 close all;
 clear all;
z=7;
x=1:z-1;
q=0.1:0.1:0.5;
%q=0.1 ;
lamda=2.456349;
for i=1:z-1
    y(i)=poisspdf(x(i),lamda);
    y1(i)=1-sum(y(i));
end
for k=1:5
    v(k)=abs((q(k)/(1-q(k))))
end
   % v=abs(q/(1-q));
syms n
    y2=symprod(v,n,1,6);
Pz=1-y1(i)*(1-y2)

i did this code but i think the result is not correct ,coz i don't know how to find the value of lamda by using if condition

[Walter Roberson]

Sofy comments to Torsten

yeah calculate p_z when q=[0.1,0.2,0.3,0.4,0.5] p=1-q

[Torsten]

And what are q_1, p_1, q_2, p_2 ? In your case q_1 = 0.1, p_1 = 0.9, q_2 = 0.2, p_2 = 0.8 ?

[Walter Roberson]

Sofy comments

yeah

Answer 1

[Torsten]

function main
  p1 = 0.9;
  p2 = 0.8;
  q1 = 1-p1;
  q2 = 1-p2;
  z = 6;
  vec = zeros(1,z);
  vec(1:2:end) = q1/p1;
  vec(2:2:end) = q2/p2;
  lambda = sum(vec);
  terms_right = [1-fliplr(cumprod(vec)), 0];
  terms_left = lambda.^(0:z)./factorial(0:z);
  P_z = 1 - exp(-lambda)*sum(terms_left.*terms_right)
end 

[Rik]

Comment posted as flag by Sofy:

thx bro appreciate your help