> expect4.ec<-array(rep(0,4^w), rep(4,w))
> for(i in 1:4){
for(j in 1:4){
for(k in 1:4){
for(m in 1:4){
expect4[i,j,k,m]<-
+ 48*49*prob.ec[i]*prob.ec[j]*prob.ec[k]*prob.ec[m]
} } } }