<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
<html>
<head>
<meta http-equiv="content-type" content="text/html; charset=ISO-8859-1">
</head>
<body bgcolor="#ffffff" text="#000000">
I have implemented hybrid mcmc for mixed models as well as non mixed
models.<br>
It now does restarts. An example program is attached. myrun will run
the model<br>
and save the results then run -mceval to get the values for the
parameters of interest.<br>
<br>
This example comes from the bugs list archives. Just in case the link
doesn't work <br>
I have included the bugs code and the data.<br>
<br>
Johnoel says that thius new code is not available in binary form so you
will need to<br>
download the source and compile it to try it out. Johnoel is the person
to get the details <br>
from to download the source. <br>
<br>
Here is the bugs link.<br>
<br>
<a
href="https://www.jiscmail.ac.uk/cgi-bin/webadmin?A2=ind0911&L=BUGS&P=R71&1=BUGS&9=A&J=on&K=2&d=No+Match%3BMatch%3BMatches&z=4">JISCMail
- BUGS Archives</a>
<br>
<br>
<br>
<a class="moz-txt-link-freetext" href="https://www.jiscmail.ac.uk/cgi-bin/webadmin?A2=ind0911&L=BUGS&P=R71&1=BUGS&9=A&J=on&K=2&d=No+Match%3BMatch%3BMatches&z=4">https://www.jiscmail.ac.uk/cgi-bin/webadmin?A2=ind0911&L=BUGS&P=R71&1=BUGS&9=A&J=on&K=2&d=No+Match%3BMatch%3BMatches&z=4</a><br>
<br>
<br>
<br>
######## Model ###############################<br>
<br>
model<br>
{<br>
for(i in 1 : N) {<br>
<br>
Y[i] ~ dbeta(w[i],d[i])<br>
<br>
w[i] <- mu[i]*phi[i]<br>
d[i] <- phi[i] - mu[i]*phi[i]<br>
<br>
logit(mu[i]) <- b[1] + u[subject[i],1]<br>
<br>
phi[i] <- exp(-1*(b[2] + u[subject[i],2]))<br>
<br>
}<br>
<br>
# Priors for random coefficients:<br>
<br>
for(j in 1 : M) {<br>
<br>
u[j,1:2] ~ dmnorm(zero[1:2],T[1:2,1:2])<br>
<br>
}<br>
<br>
# Priors for fixed effects<br>
for (k in 1:2) { b[k] ~ dflat() }<br>
<br>
# priors for random effects:<br>
<br>
for (i in 1:2) {zero[i] <- 0}<br>
<br>
T[1:2,1:2] ~ dwish(R[,],2) <br>
<br>
V[1:2,1:2] <- inverse(T[1:2,1:2])<br>
<br>
}<br>
<br>
########## Data #####################<br>
<br>
list(N=5.00000E+02, M=1.00000E+02, R= structure(.Data= c(1.00000E+00, <br>
2.12132E-01, 2.12132E-01, 5.00000E-01), .Dim=c(2, 2)), <br>
subject=c(1.00000E+00, 1.00000E+00, 1.00000E+00, 1.00000E+00, <br>
1.00000E+00, 2.00000E+00, 2.00000E+00, 2.00000E+00, 2.00000E+00, <br>
2.00000E+00, 3.00000E+00, 3.00000E+00, 3.00000E+00, 3.00000E+00, <br>
3.00000E+00, 4.00000E+00, 4.00000E+00, 4.00000E+00, 4.00000E+00, <br>
4.00000E+00, 5.00000E+00, 5.00000E+00, 5.00000E+00, 5.00000E+00, <br>
5.00000E+00, 6.00000E+00, 6.00000E+00, 6.00000E+00, 6.00000E+00, <br>
6.00000E+00, 7.00000E+00, 7.00000E+00, 7.00000E+00, 7.00000E+00, <br>
7.00000E+00, 8.00000E+00, 8.00000E+00, 8.00000E+00, 8.00000E+00, <br>
8.00000E+00, 9.00000E+00, 9.00000E+00, 9.00000E+00, 9.00000E+00, <br>
9.00000E+00, 1.00000E+01, 1.00000E+01, 1.00000E+01, 1.00000E+01, <br>
1.00000E+01, 1.10000E+01, 1.10000E+01, 1.10000E+01, 1.10000E+01, <br>
1.10000E+01, 1.20000E+01, 1.20000E+01, 1.20000E+01, 1.20000E+01, <br>
1.20000E+01, 1.30000E+01, 1.30000E+01, 1.30000E+01, 1.30000E+01, <br>
1.30000E+01, 1.40000E+01, 1.40000E+01, 1.40000E+01, 1.40000E+01, <br>
1.40000E+01, 1.50000E+01, 1.50000E+01, 1.50000E+01, 1.50000E+01, <br>
1.50000E+01, 1.60000E+01, 1.60000E+01, 1.60000E+01, 1.60000E+01, <br>
1.60000E+01, 1.70000E+01, 1.70000E+01, 1.70000E+01, 1.70000E+01, <br>
1.70000E+01, 1.80000E+01, 1.80000E+01, 1.80000E+01, 1.80000E+01, <br>
1.80000E+01, 1.90000E+01, 1.90000E+01, 1.90000E+01, 1.90000E+01, <br>
1.90000E+01, 2.00000E+01, 2.00000E+01, 2.00000E+01, 2.00000E+01, <br>
2.00000E+01, 2.10000E+01, 2.10000E+01, 2.10000E+01, 2.10000E+01, <br>
2.10000E+01, 2.20000E+01, 2.20000E+01, 2.20000E+01, 2.20000E+01, <br>
2.20000E+01, 2.30000E+01, 2.30000E+01, 2.30000E+01, 2.30000E+01, <br>
2.30000E+01, 2.40000E+01, 2.40000E+01, 2.40000E+01, 2.40000E+01, <br>
2.40000E+01, 2.50000E+01, 2.50000E+01, 2.50000E+01, 2.50000E+01, <br>
2.50000E+01, 2.60000E+01, 2.60000E+01, 2.60000E+01, 2.60000E+01, <br>
2.60000E+01, 2.70000E+01, 2.70000E+01, 2.70000E+01, 2.70000E+01, <br>
2.70000E+01, 2.80000E+01, 2.80000E+01, 2.80000E+01, 2.80000E+01, <br>
2.80000E+01, 2.90000E+01, 2.90000E+01, 2.90000E+01, 2.90000E+01, <br>
2.90000E+01, 3.00000E+01, 3.00000E+01, 3.00000E+01, 3.00000E+01, <br>
3.00000E+01, 3.10000E+01, 3.10000E+01, 3.10000E+01, 3.10000E+01, <br>
3.10000E+01, 3.20000E+01, 3.20000E+01, 3.20000E+01, 3.20000E+01, <br>
3.20000E+01, 3.30000E+01, 3.30000E+01, 3.30000E+01, 3.30000E+01, <br>
3.30000E+01, 3.40000E+01, 3.40000E+01, 3.40000E+01, 3.40000E+01, <br>
3.40000E+01, 3.50000E+01, 3.50000E+01, 3.50000E+01, 3.50000E+01, <br>
3.50000E+01, 3.60000E+01, 3.60000E+01, 3.60000E+01, 3.60000E+01, <br>
3.60000E+01, 3.70000E+01, 3.70000E+01, 3.70000E+01, 3.70000E+01, <br>
3.70000E+01, 3.80000E+01, 3.80000E+01, 3.80000E+01, 3.80000E+01, <br>
3.80000E+01, 3.90000E+01, 3.90000E+01, 3.90000E+01, 3.90000E+01, <br>
3.90000E+01, 4.00000E+01, 4.00000E+01, 4.00000E+01, 4.00000E+01, <br>
4.00000E+01, 4.10000E+01, 4.10000E+01, 4.10000E+01, 4.10000E+01, <br>
4.10000E+01, 4.20000E+01, 4.20000E+01, 4.20000E+01, 4.20000E+01, <br>
4.20000E+01, 4.30000E+01, 4.30000E+01, 4.30000E+01, 4.30000E+01, <br>
4.30000E+01, 4.40000E+01, 4.40000E+01, 4.40000E+01, 4.40000E+01, <br>
4.40000E+01, 4.50000E+01, 4.50000E+01, 4.50000E+01, 4.50000E+01, <br>
4.50000E+01, 4.60000E+01, 4.60000E+01, 4.60000E+01, 4.60000E+01, <br>
4.60000E+01, 4.70000E+01, 4.70000E+01, 4.70000E+01, 4.70000E+01, <br>
4.70000E+01, 4.80000E+01, 4.80000E+01, 4.80000E+01, 4.80000E+01, <br>
4.80000E+01, 4.90000E+01, 4.90000E+01, 4.90000E+01, 4.90000E+01, <br>
4.90000E+01, 5.00000E+01, 5.00000E+01, 5.00000E+01, 5.00000E+01, <br>
5.00000E+01, 5.10000E+01, 5.10000E+01, 5.10000E+01, 5.10000E+01, <br>
5.10000E+01, 5.20000E+01, 5.20000E+01, 5.20000E+01, 5.20000E+01, <br>
5.20000E+01, 5.30000E+01, 5.30000E+01, 5.30000E+01, 5.30000E+01, <br>
5.30000E+01, 5.40000E+01, 5.40000E+01, 5.40000E+01, 5.40000E+01, <br>
5.40000E+01, 5.50000E+01, 5.50000E+01, 5.50000E+01, 5.50000E+01, <br>
5.50000E+01, 5.60000E+01, 5.60000E+01, 5.60000E+01, 5.60000E+01, <br>
5.60000E+01, 5.70000E+01, 5.70000E+01, 5.70000E+01, 5.70000E+01, <br>
5.70000E+01, 5.80000E+01, 5.80000E+01, 5.80000E+01, 5.80000E+01, <br>
5.80000E+01, 5.90000E+01, 5.90000E+01, 5.90000E+01, 5.90000E+01, <br>
5.90000E+01, 6.00000E+01, 6.00000E+01, 6.00000E+01, 6.00000E+01, <br>
6.00000E+01, 6.10000E+01, 6.10000E+01, 6.10000E+01, 6.10000E+01, <br>
6.10000E+01, 6.20000E+01, 6.20000E+01, 6.20000E+01, 6.20000E+01, <br>
6.20000E+01, 6.30000E+01, 6.30000E+01, 6.30000E+01, 6.30000E+01, <br>
6.30000E+01, 6.40000E+01, 6.40000E+01, 6.40000E+01, 6.40000E+01, <br>
6.40000E+01, 6.50000E+01, 6.50000E+01, 6.50000E+01, 6.50000E+01, <br>
6.50000E+01, 6.60000E+01, 6.60000E+01, 6.60000E+01, 6.60000E+01, <br>
6.60000E+01, 6.70000E+01, 6.70000E+01, 6.70000E+01, 6.70000E+01, <br>
6.70000E+01, 6.80000E+01, 6.80000E+01, 6.80000E+01, 6.80000E+01, <br>
6.80000E+01, 6.90000E+01, 6.90000E+01, 6.90000E+01, 6.90000E+01, <br>
6.90000E+01, 7.00000E+01, 7.00000E+01, 7.00000E+01, 7.00000E+01, <br>
7.00000E+01, 7.10000E+01, 7.10000E+01, 7.10000E+01, 7.10000E+01, <br>
7.10000E+01, 7.20000E+01, 7.20000E+01, 7.20000E+01, 7.20000E+01, <br>
7.20000E+01, 7.30000E+01, 7.30000E+01, 7.30000E+01, 7.30000E+01, <br>
7.30000E+01, 7.40000E+01, 7.40000E+01, 7.40000E+01, 7.40000E+01, <br>
7.40000E+01, 7.50000E+01, 7.50000E+01, 7.50000E+01, 7.50000E+01, <br>
7.50000E+01, 7.60000E+01, 7.60000E+01, 7.60000E+01, 7.60000E+01, <br>
7.60000E+01, 7.70000E+01, 7.70000E+01, 7.70000E+01, 7.70000E+01, <br>
7.70000E+01, 7.80000E+01, 7.80000E+01, 7.80000E+01, 7.80000E+01, <br>
7.80000E+01, 7.90000E+01, 7.90000E+01, 7.90000E+01, 7.90000E+01, <br>
7.90000E+01, 8.00000E+01, 8.00000E+01, 8.00000E+01, 8.00000E+01, <br>
8.00000E+01, 8.10000E+01, 8.10000E+01, 8.10000E+01, 8.10000E+01, <br>
8.10000E+01, 8.20000E+01, 8.20000E+01, 8.20000E+01, 8.20000E+01, <br>
8.20000E+01, 8.30000E+01, 8.30000E+01, 8.30000E+01, 8.30000E+01, <br>
8.30000E+01, 8.40000E+01, 8.40000E+01, 8.40000E+01, 8.40000E+01, <br>
8.40000E+01, 8.50000E+01, 8.50000E+01, 8.50000E+01, 8.50000E+01, <br>
8.50000E+01, 8.60000E+01, 8.60000E+01, 8.60000E+01, 8.60000E+01, <br>
8.60000E+01, 8.70000E+01, 8.70000E+01, 8.70000E+01, 8.70000E+01, <br>
8.70000E+01, 8.80000E+01, 8.80000E+01, 8.80000E+01, 8.80000E+01, <br>
8.80000E+01, 8.90000E+01, 8.90000E+01, 8.90000E+01, 8.90000E+01, <br>
8.90000E+01, 9.00000E+01, 9.00000E+01, 9.00000E+01, 9.00000E+01, <br>
9.00000E+01, 9.10000E+01, 9.10000E+01, 9.10000E+01, 9.10000E+01, <br>
9.10000E+01, 9.20000E+01, 9.20000E+01, 9.20000E+01, 9.20000E+01, <br>
9.20000E+01, 9.30000E+01, 9.30000E+01, 9.30000E+01, 9.30000E+01, <br>
9.30000E+01, 9.40000E+01, 9.40000E+01, 9.40000E+01, 9.40000E+01, <br>
9.40000E+01, 9.50000E+01, 9.50000E+01, 9.50000E+01, 9.50000E+01, <br>
9.50000E+01, 9.60000E+01, 9.60000E+01, 9.60000E+01, 9.60000E+01, <br>
9.60000E+01, 9.70000E+01, 9.70000E+01, 9.70000E+01, 9.70000E+01, <br>
9.70000E+01, 9.80000E+01, 9.80000E+01, 9.80000E+01, 9.80000E+01, <br>
9.80000E+01, 9.90000E+01, 9.90000E+01, 9.90000E+01, 9.90000E+01, <br>
9.90000E+01, 1.00000E+02, 1.00000E+02, 1.00000E+02, 1.00000E+02, <br>
1.00000E+02), Y=c(4.95613E-02, 7.79941E-02, 1.81326E-02, 3.09964E-02, <br>
4.36947E-02, 3.89151E-01, 2.56895E-01, 3.13902E-01, 1.95085E-01, <br>
3.75864E-01, 3.34703E-01, 2.01205E-01, 2.91234E-01, 2.04705E-01, <br>
2.15405E-01, 5.29656E-01, 4.63180E-01, 4.05392E-01, 4.98864E-01, <br>
4.18180E-01, 2.35767E-01, 3.71746E-01, 3.02582E-01, 2.78117E-01, <br>
4.09888E-01, 8.30437E-02, 1.22581E-01, 8.96370E-02, 1.24107E-01, <br>
1.30841E-01, 4.98086E-01, 3.38277E-01, 3.21507E-01, 4.27241E-01, <br>
2.89721E-01, 3.00415E-01, 2.93250E-01, 2.29410E-01, 3.25913E-01, <br>
3.53117E-01, 2.99300E-01, 3.39513E-01, 3.00704E-01, 3.27941E-01, <br>
3.43273E-01, 4.21016E-01, 2.71075E-01, 1.82676E-01, 3.58336E-01, <br>
4.11760E-01, 1.32233E-01, 1.81819E-01, 1.16145E-01, 2.39060E-01, <br>
1.07757E-01, 4.59173E-01, 2.81766E-01, 4.26113E-01, 2.88520E-01, <br>
4.34517E-01, 3.92187E-01, 3.60408E-01, 4.30432E-01, 3.91706E-01, <br>
2.83573E-01, 6.75149E-01, 6.98988E-01, 4.71903E-01, 8.07172E-01, <br>
6.88902E-01, 4.04358E-01, 5.29541E-01, 5.67765E-01, 4.61854E-01, <br>
3.70825E-01, 2.11474E-01, 5.51010E-01, 5.02223E-01, 3.29526E-01, <br>
4.16277E-01, 1.43770E-01, 2.24790E-01, 8.43102E-02, 2.33085E-01, <br>
2.51197E-01, 1.53812E-01, 1.43035E-01, 1.47347E-01, 2.07577E-01, <br>
1.42093E-01, 6.51663E-01, 5.45070E-01, 6.88870E-01, 6.14606E-01, <br>
6.17852E-01, 1.63792E-01, 7.42528E-02, 1.12385E-01, 4.23506E-02, <br>
1.56944E-01, 5.86904E-01, 2.97494E-01, 4.41823E-01, 5.33891E-01, <br>
3.59901E-01, 2.56150E-01, 3.09283E-01, 3.35121E-01, 3.11554E-01, <br>
2.51390E-01, 2.69495E-01, 2.58349E-01, 2.80504E-01, 2.77081E-01, <br>
2.41944E-01, 2.80987E-01, 1.56312E-01, 2.03176E-01, 1.94157E-01, <br>
1.93468E-01, 1.36268E-01, 1.15510E-01, 1.74921E-01, 1.37439E-01, <br>
2.34397E-01, 1.41322E-01, 1.99544E-01, 1.81958E-01, 2.01580E-01, <br>
1.87999E-01, 1.00902E-01, 1.03753E-01, 1.65586E-01, 1.57699E-01, <br>
1.11619E-01, 4.75960E-01, 5.66968E-01, 6.04824E-01, 6.14615E-01, <br>
4.40675E-01, 5.57426E-01, 5.24392E-01, 5.72538E-01, 4.17892E-01, <br>
5.75137E-01, 5.92459E-01, 5.26017E-01, 5.70187E-01, 5.13941E-01, <br>
7.23242E-01, 5.63730E-01, 4.18315E-01, 4.69322E-01, 4.71844E-01, <br>
5.31627E-01, 3.18218E-01, 2.36380E-01, 1.24041E-01, 3.27462E-01, <br>
2.89023E-01, 4.11009E-01, 4.25972E-01, 4.80736E-01, 4.48034E-01, <br>
4.37771E-01, 2.13215E-01, 2.68670E-01, 1.49843E-01, 1.39049E-01, <br>
1.88233E-01, 8.46093E-01, 7.86206E-01, 8.14992E-01, 7.79949E-01, <br>
8.11699E-01, 2.58419E-01, 1.41543E-01, 2.97552E-01, 2.50150E-01, <br>
1.89388E-01, 4.30771E-01, 5.57780E-01, 4.02770E-01, 4.16729E-01, <br>
5.38634E-01, 4.24493E-01, 5.79367E-01, 2.76611E-01, 3.99875E-01, <br>
2.95807E-01, 2.32242E-01, 1.91015E-01, 2.90013E-01, 2.31246E-01, <br>
2.94564E-01, 7.93744E-02, 7.37363E-02, 6.85535E-02, 8.20764E-02, <br>
7.57085E-02, 6.19427E-01, 5.85733E-01, 6.15442E-01, 6.54538E-01, <br>
5.37639E-01, 6.75315E-01, 5.90740E-01, 5.97983E-01, 6.57791E-01, <br>
4.64374E-01, 9.56814E-02, 1.27048E-01, 1.84311E-01, 1.70066E-01, <br>
6.21729E-02, 2.38550E-01, 2.38041E-01, 3.56752E-01, 4.35224E-01, <br>
2.87235E-01, 6.99641E-01, 6.88217E-01, 8.18052E-01, 5.80178E-01, <br>
6.69430E-01, 2.02981E-01, 1.23304E-01, 1.41922E-01, 2.26863E-01, <br>
2.40617E-01, 3.98603E-01, 2.72732E-01, 1.59025E-01, 1.90095E-01, <br>
9.75554E-02, 2.57337E-01, 3.21510E-01, 3.12309E-01, 2.47899E-01, <br>
2.88169E-01, 5.23101E-01, 3.90842E-01, 4.87574E-01, 4.63070E-01, <br>
4.29385E-01, 3.01368E-01, 4.30167E-01, 4.20292E-01, 2.30674E-01, <br>
3.13277E-01, 3.38028E-01, 2.54546E-01, 3.86777E-01, 6.32491E-01, <br>
5.86292E-01, 1.82784E-01, 1.37193E-01, 1.92350E-01, 1.56040E-01, <br>
1.54229E-01, 5.43241E-01, 4.23037E-01, 5.62903E-01, 5.57844E-01, <br>
6.36971E-01, 2.80311E-01, 4.73372E-01, 3.27982E-01, 4.58843E-01, <br>
3.80366E-01, 7.09897E-01, 6.07523E-01, 7.35242E-01, 5.61583E-01, <br>
7.96395E-01, 1.84466E-01, 1.59910E-01, 2.10035E-01, 1.94682E-01, <br>
1.54933E-01, 6.90812E-01, 3.58927E-01, 5.75018E-01, 5.31190E-01, <br>
5.11009E-01, 2.42393E-01, 1.82160E-01, 2.35204E-01, 1.93499E-01, <br>
2.07383E-01, 2.97084E-01, 1.38284E-01, 1.52242E-01, 1.67804E-01, <br>
3.84020E-01, 2.63975E-01, 3.56323E-01, 3.02182E-01, 3.11180E-01, <br>
2.98539E-01, 2.04048E-01, 4.95719E-02, 6.91643E-02, 7.46431E-02, <br>
9.82683E-02, 8.36914E-02, 9.34691E-02, 1.13242E-01, 1.27747E-01, <br>
8.94278E-02, 2.13942E-01, 2.17447E-01, 1.93899E-01, 2.25857E-01, <br>
2.20856E-01, 4.89508E-01, 7.38538E-01, 7.65385E-01, 3.71710E-01, <br>
7.12922E-01, 2.53409E-01, 2.33998E-01, 4.02738E-01, 4.81510E-01, <br>
2.58167E-01, 3.71673E-01, 4.28092E-01, 5.19156E-01, 4.12433E-01, <br>
5.58651E-01, 7.25450E-02, 1.14948E-01, 1.18412E-01, 8.38131E-02, <br>
1.30599E-01, 6.69427E-01, 4.81964E-01, 4.09181E-01, 6.40433E-01, <br>
7.56217E-01, 3.24778E-01, 3.14451E-01, 3.22245E-01, 4.59406E-01, <br>
3.99069E-01, 5.18559E-01, 4.00973E-01, 3.90888E-01, 5.28822E-01, <br>
2.64378E-01, 2.35993E-01, 1.19832E-01, 1.18974E-01, 1.30372E-01, <br>
2.66793E-01, 4.02460E-01, 4.68968E-01, 4.90212E-01, 3.73217E-01, <br>
3.79813E-01, 1.13751E-01, 1.04080E-01, 1.73163E-01, 1.97727E-01, <br>
1.66175E-01, 5.11393E-01, 4.69464E-01, 7.92517E-01, 6.07294E-01, <br>
7.95619E-01, 3.56549E-01, 4.21971E-01, 4.16195E-01, 4.98982E-01, <br>
3.26022E-01, 3.33549E-01, 2.84551E-01, 2.77705E-01, 3.72090E-01, <br>
3.57299E-01, 4.25189E-01, 4.15545E-01, 3.68058E-01, 3.33256E-01, <br>
3.92429E-01, 4.64911E-01, 5.49798E-01, 4.47110E-01, 4.72110E-01, <br>
3.74451E-01, 4.67185E-01, 4.80916E-01, 5.06014E-01, 4.49430E-01, <br>
5.16906E-01, 7.04964E-02, 1.57736E-01, 2.03266E-01, 9.98466E-02, <br>
1.40007E-01, 3.76992E-01, 4.79178E-01, 3.89255E-01, 5.03434E-01, <br>
3.71304E-01, 2.93645E-02, 5.14925E-02, 3.97149E-02, 4.76645E-02, <br>
3.61847E-02, 6.04385E-01, 6.31799E-01, 6.59257E-01, 6.12542E-01, <br>
4.94529E-01, 3.98921E-01, 5.95958E-01, 5.94514E-01, 4.37820E-01, <br>
5.78407E-01, 1.10775E-01, 1.92568E-01, 2.56543E-01, 2.16718E-01, <br>
1.02948E-01, 3.31866E-02, 3.15410E-02, 3.44858E-02, 7.84109E-02, <br>
4.72811E-02, 9.59095E-02, 1.10756E-01, 6.45108E-02, 6.18701E-02, <br>
4.99215E-02, 4.36752E-01, 3.60531E-01, 4.81942E-01, 4.53305E-01, <br>
4.36145E-01, 6.78271E-01, 5.65400E-01, 4.47615E-01, 5.75648E-01, <br>
7.97650E-01, 5.51668E-01, 5.04297E-01, 5.54140E-01, 4.48877E-01, <br>
5.64451E-01, 8.86784E-02, 1.51024E-01, 7.97259E-02, 1.36798E-01, <br>
1.24054E-01, 3.31096E-01, 1.94807E-01, 2.60465E-01, 2.66845E-01, <br>
2.41661E-01, 5.58066E-01, 6.16075E-01, 5.62075E-01, 5.47844E-01, <br>
6.02110E-01, 5.47136E-01, 4.68577E-01, 6.22336E-01, 5.53878E-01, <br>
5.08098E-01, 2.27185E-01, 1.65477E-01, 4.50822E-02, 1.01362E-01, <br>
2.41054E-01, 4.20258E-01, 4.52138E-01, 4.41439E-01, 6.33041E-01, <br>
3.75203E-01, 1.22185E-01, 9.52278E-02, 1.89553E-01, 1.12670E-01, <br>
1.41032E-01, 5.27857E-01, 5.12138E-01, 4.70416E-01, 5.54735E-01, <br>
3.86392E-01, 4.13163E-02, 2.08262E-02, 1.06521E-01, 6.14467E-03, <br>
1.20814E-02, 1.36490E-01, 3.33253E-01, 5.71203E-02, 1.53736E-01, <br>
1.19693E-01))<br>
<br>
<br>
############# Initial Values ######################<br>
<br>
list(b=c(-7.50000E-01, -4.00000E+00), T= structure(.Data=
c(1.09890E+00, <br>
-4.66224E-01, -4.66224E-01, 2.19780E+00), .Dim=c(2, 2)), u= <br>
structure(.Data= c(0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, <br>
0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, <br>
0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, <br>
0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, <br>
0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, <br>
0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, <br>
0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, <br>
0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, <br>
0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, <br>
0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, <br>
0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, <br>
0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, <br>
0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, <br>
0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, <br>
0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, <br>
0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, <br>
0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, <br>
0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, <br>
0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, <br>
0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, <br>
0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, <br>
0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, <br>
0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, <br>
0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, <br>
0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, <br>
0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, <br>
0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, <br>
0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, <br>
0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, <br>
0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, <br>
0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, <br>
0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, <br>
0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, <br>
0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, <br>
0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, <br>
0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, <br>
0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, <br>
0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, <br>
0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, <br>
0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, <br>
0.00000E+00), .Dim=c(100, 2)))<br>
<br>
-------------------------------------------------------------------<br>
<br>
</body>
</html>