作者bcs (= ="frailty..gggg XD)
看板Statistics
标题[问题] matlab转至R的问题
时间Mon May 11 13:10:21 2009
刚试了从matlab的程式转到R，但是出来的结果很不相同
检查许多遍觉得应该会是相同(rowsum应该都是1)，
但可能有细微的地方没弄清，请问哪个转码的过程弄错了？
谢谢
%Matlab
function pr = markov(n,rho);
pr=zeros(n,n);
p1=pr;p2=pr;p3=pr;p4=pr;
p=0.5*(1+rho);
pr(1:2,1:2)=[p 1-p;1-p p];
for i = 3:n;
p1(1:i-1,1:i-1)=pr(1:i-1,1:i-1);
p2(1:i-1,2:i )=pr(1:i-1,1:i-1);
p3(2:i ,1:i-1)=pr(1:i-1,1:i-1);
p4(2:i ,2:i )=pr(1:i-1,1:i-1);
pr=p*(p1+p4)+(1-p)*(p2+p3);
pr(2:i-1,1:i)=0.5*pr(2:i-1,1:i);
end;
a=markov(4,.2)
sum(a,2)
#R
markov<-function(n,rho){
pr=matrix(0,n,n);
p1=p2=p3=p4=pr;
p=(1+rho)*0.5;
pr[1:2,1:2]=matrix(c(p,1-p),2,2);
for (i in 3:n){
p1[1:i-1,1:i-1]=pr[1:i-1,1:i-1];
p2[1:i-1,2:i ]=pr[1:i-1,1:i-1];
p3[2:i ,1:i-1]=pr[1:i-1,1:i-1];
p4[2:i ,2:i ]=pr[1:i-1,1:i-1];
pr=p*(p1+p4)+(1-p)*(p2+p3);
pr[2:i-1,1:i]=0.5*pr[2:i-1,1:i];
}
return(pr);
}
a=markov(4,.2);
rowSums(a);
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 140.112.86.139
※ 编辑: bcs 来自: 140.112.86.139 (05/11 13:12)
※ 编辑: bcs 来自: 140.112.86.139 (05/11 13:12)
1F：→ clickhere:1:(i-1) 05/11 15:28
2F：推 lin15:b大可以问一下要怎麽转吗... 05/11 16:31
3F：→ bcs:用手工调的 05/13 10:50