Перевести код SAS на Python графика расхождения KL —

#python-3.x #data-visualization #translate #information-theory #sas-iml

#python-3.x #визуализация данных #перевести #информация-теория #sas-iml

Вопрос:

Как выглядит приведенный ниже код SAS при переводе на python?

 /* Plot the K-L div versus lambda for a sequence of Poisson(lambda) models */
lambda = do(4, 16, 0.1);
KL = j(1, ncol(lambda), .);
do i = 1 to ncol(lambda);
   poisPDF = pdf("Poisson", Levels, lambda[i]);
   KL[i] = KLDiv(Proportion, poisPDF); 
end;
 
title "K-L Divergence from Poisson(lambda)";
call series(lambda, KL) grid={x y} xvalues=4:16 label={'x' 'K-L Divergence'};
  

Источник

Комментарии:

1. Это не базовый код SAS. Вместо этого это матричный код, написанный для PROC IML.

2. можно ли его перевести на python