| ²é¿´: 323 | »Ø¸´: 1 | ||
xuhaibocn¾èÖú¹ó±ö (СÓÐÃûÆø)
|
[ÇóÖú]
Çó°ïæ½â ÊÍËľä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))); |
» ²ÂÄãϲ»¶
085801 ×Ü·Ö275 ±¾¿ÆÐÂÄÜÔ´ Çóµ÷¼Á
ÒѾÓÐ22È˻ظ´
298·Ö 070300Çóµ÷¼Á
ÒѾÓÐ3È˻ظ´
²ÄÁÏÓ뻯¹¤306·ÖÕÒµ÷¼Á
ÒѾÓÐ13È˻ظ´
070305¸ß·Ö×Ó»¯Ñ§ÓëÎïÀí 304·ÖÇóµ÷¼Á
ÒѾÓÐ15È˻ظ´
260Çóµ÷¼Á
ÒѾÓÐ4È˻ظ´
299Çóµ÷¼Á
ÒѾÓÐ5È˻ظ´
Çóµ÷¼ÁÇóµ÷¼Á
ÒѾÓÐ7È˻ظ´
265Çóµ÷¼Á
ÒѾÓÐ15È˻ظ´
070300»¯Ñ§Çóµ÷¼Á
ÒѾÓÐ5È˻ظ´
296²ÄÁÏר˶Çóµ÷¼Á
ÒѾÓÐ13È˻ظ´
xuhaibocn
¾èÖú¹ó±ö (СÓÐÃûÆø)
- Ó¦Öú: 0 (Ó×¶ùÔ°)
- ½ð±Ò: 32.6
- Ìû×Ó: 60
- ÔÚÏß: 24.4Сʱ
- ³æºÅ: 690185
- ×¢²á: 2009-01-09
- רҵ: Êг¡ÓªÏú
|
Õû¸ö³ÌÐò %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














»Ø¸´´ËÂ¥