cov=corrsigma[1,1]sigma[2,1]
mup=.18
I=matrix(c(1,1),nrow=2,ncol=1)
varcov=matrix(c((sigma[1,1])^2,cov,cov,(sigma[2,1])^2),nrow=2,ncol=
2,byrow=TRUE)
A=t(mu)%%solve(varcov)%%mu
B=t(mu)%%solve(varcov)%%I
C=t(I)%%solve(varcov)%%I
GG=matrix(c(A,B,B,C), nrow=2,ncol=2,byrow=TRUE)
lambda=(Cmup-B)/(AC-B^2)
delta=(A-Bmup)/(AC-B^2)
Q=mu%%lambda+I%%delta
omega=solve(varcov)%*%Q
retp=t(mu)%*%omega
sigmap=sqrt(t(omega)%%varcov%%omega)}
- Par itération du code R, on obtient :
mup 15 16 17 17 18 19 20 .....
..
25 26 26 27 28 29 29 30
sigm
ap24 25 25 26 26 27 28 .....
..35 36 38 39 40 41 43 44- Graphiquement on obtient :