| 查看: 334 | 回复: 1 | ||
[求助]
求帮忙解 释四句code
|
|
就这四句 ndr(:,i)=tqp*de+(sqrt(dsv)*sqrt(de)).*dr(:,i); prp2=reshape(prp2',n*nd,6); prp2=mean(reshape((prod(prp2',1))',nd,n),1); f=-sum(sum(log(prp2))); |
回复此楼» 猜你喜欢
提交了我也来说说感想
已经有11人回复
-
评审有感
已经有30人回复
-
青B发送上会通知了吗
已经有9人回复
-
西安交大新媒学院副院长用撤稿论文结题
已经有6人回复
-
论文撤稿了
已经有8人回复
-
博士申请
已经有7人回复
-
化学专业申博
已经有4人回复
-
河北省自然科学基金
已经有9人回复
-
某211大学教师把个人教师官方主页改成:我跑了我跑了我跑了!官宣跑路！
已经有5人回复
-
26/27申博自荐
已经有9人回复
|
整个程序 %Procedure for the likelihood% function f=lpr(x) global de dr pres uptop nd n prp prp2 ndr qp; % Initial Quality Perception; % True Quality Perception; % Detailing Signal Variance; % Initial Quality Variance iqp=x(1); tqp=x(2); dsv=exp(x(3)); iqv=1; ndr=zeros(n*6,nd); qp=ndr; for i=1:nd %Distribution of Signals for each time period% ndr(:,i)=tqp*de+(sqrt(dsv)*sqrt(de)).*dr(:,i); %Mean quality perceptions for each time period% qp(:,i)=((iqp/iqv)+(ndr(:,i)/dsv))./((1/iqv)+(de/dsv)); end %Prescription probability% prp=exp(qp)./(1+exp(qp)); %No prescription probability% nprp=1./(1+exp(qp)); %Choose appropriate probability or 1 depending upon behavior and number of periods physician active in data% prp2=zeros(n*6,nd); for i=1:nd prp2(:,i)=((prp(:,i).*pres)+(nprp(:,i).*(1-pres))).^uptop; end %Reshaping to compute likelihood% prp2=reshape(prp2',n*nd,6); prp2=mean(reshape((prod(prp2',1))',nd,n),1); f=-sum(sum(log(prp2))); |
2楼2014-02-09 14:16:38