කවුරුහරි ඉන්නවද flow chart එකක් අදාල දෙන්න පුළුවන් ?
මට code එක දෙන්නම්
මේ වෙලාවේ හැටියට ලොකු උදව්වක්.
ජය වේවා.
###ETINNEE(2005) MODEL FOR HABITATS
spgen.m<-function(theta,J,m){
label<<-matrix(0,2,J) #1st row for ancestry label,2nd row for species label
if(m<1)
I=(m*(J-1))/(1-m)
a=0
s=0
for(j in 1:J){
if(m<1)
R1=I/(I+j-1)
if(m==1)
R1=1
x<-runif(1)
if(x<=R1){
a=a+1
label[1,j]<-a
R2=theta/(theta+a-1)
y<-runif(1)
if(y<=R2){
s=s+1
label[2,j]=s
}
else{
anc<-sample(1
a-1),1)
label[2,j]<-label[2,anc]
}
}
else{
ind<-sample(1
j-1),1)
label[1,j]<-label[1,ind]
label[2,j]<-label[2,ind]
}
}
return(label)
}
anc.sp<-spgen.m(50,100000,1)
මට code එක දෙන්නම්
මේ වෙලාවේ හැටියට ලොකු උදව්වක්.
ජය වේවා.
###ETINNEE(2005) MODEL FOR HABITATS
spgen.m<-function(theta,J,m){
label<<-matrix(0,2,J) #1st row for ancestry label,2nd row for species label
if(m<1)
I=(m*(J-1))/(1-m)
a=0
s=0
for(j in 1:J){
if(m<1)
R1=I/(I+j-1)
if(m==1)
R1=1
x<-runif(1)
if(x<=R1){
a=a+1
label[1,j]<-a
R2=theta/(theta+a-1)
y<-runif(1)
if(y<=R2){
s=s+1
label[2,j]=s
}
else{
anc<-sample(1
a-1),1)label[2,j]<-label[2,anc]
}
}
else{
ind<-sample(1
j-1),1)label[1,j]<-label[1,ind]
label[2,j]<-label[2,ind]
}
}
return(label)
}
anc.sp<-spgen.m(50,100000,1)