| 查看: 593 | 回复: 4 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
帮忙解释加下划线的四句code，谢谢
|
||
|
Sample Text %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% Sample Text %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% Sample Text |
回复此楼» 猜你喜欢
各位大神，目前国内有哪些比较好用的逆合成软件？
已经有14人回复
-
国社科系统bug了，是不是要放榜了？
已经有10人回复
-
现代”学阀”该如何界定
已经有9人回复
-
课题组招2027级博士 上海工程技术大学 激光智能制造方向
已经有4人回复
-
上海工程技术大学 激光智能制造课题组 2027级博士研究生招生
已经有3人回复
-
上海工程技术大学激光智能制造课题组｜2027级博士研究生招生公告
已经有11人回复
-
上海工程技术大学激光智能制造课题组-博士研究生
已经有10人回复
-
申博发邮件
已经有11人回复
-
我的奶奶
已经有3人回复
5楼2014-02-09 13:06:26
|
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))); %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 03:55:51
3楼2014-02-09 03:57:07