Displaying 1 result from an estimated 1 matches for "temppzu".
Did you mean:
temppzi
2006 Jan 23
1
will vectorization help in this case?
...)
{
flst=list(Pzu, Pzi, array(Puzu[,as.character(r[[l,"u"]])],
dim=dim(Puzu)[1], dimnames=dimnames(Puzu)[1]), array(Pizi[,as.character(r[[l,
"i"]])], dim=dim(Pizi)[1], dimnames=dimnames(Pizi)[1]) );
tempa = 1;
for(j in 1:length(flst)) tempa = tempa %X% flst[[j]];
tempPzu = apply(tempa, 1, sum);
tPzu = tPzu + tempPzu;
tPuzu[,as.character(r[[l,"u"]])] = tPuzu[,as.character(r[[l,"u"]])] +
tempPzu;
tempPzi = apply(tempa, 2, sum);
tPzi = tPzi + tempPzi;
tPizi[,as.character(r[[l,"i"]])] = tPizi[,as.character(r[[l,"i...