version 1.19, 2004/01/07 06:53:11 |
version 1.23, 2004/01/08 06:48:32 |
Line 565 def leastsq(NormMat,ExpMat,Vars,FLAG,ID){ |
|
Line 565 def leastsq(NormMat,ExpMat,Vars,FLAG,ID){ |
|
return RET$ |
return RET$ |
} |
} |
else{ |
else{ |
RET=append(RET,[[ID,Vars,vtol(TMP1[1]*1.0)]])$ |
RET=append(RET,[[ID,Vars,vtol(TMP1[1])]])$ |
return RET$ |
return RET$ |
} |
} |
} |
} |
Line 690 def weight(PolyList,Vars,FLAG){ |
|
Line 690 def weight(PolyList,Vars,FLAG){ |
|
|
|
RET=append(RET,TMP[1])$ |
RET=append(RET,TMP[1])$ |
|
|
RET=append(RET,leastsq(TMP[0],ExpMat,Vars,FLAG,3))$ |
TMP0=leastsq(TMP[0],ExpMat,Vars,FLAG,3)$ |
|
|
|
RET=append(RET,TMP0)$ |
|
|
ExpMat=qsort(ExpMat,junban)$ |
ExpMat=qsort(ExpMat,junban)$ |
|
|
ExpMat2=[]$ |
ExpMat2=[]$ |
Line 702 def weight(PolyList,Vars,FLAG){ |
|
Line 704 def weight(PolyList,Vars,FLAG){ |
|
if(size(ExpMat)[0]!=length(ExpMat2)){ |
if(size(ExpMat)[0]!=length(ExpMat2)){ |
ExpMat=newvect(length(ExpMat2),ExpMat2)$ |
ExpMat=newvect(length(ExpMat2),ExpMat2)$ |
RET=append(RET,leastsq(0,ExpMat,Vars,FLAG,5))$ |
RET=append(RET,leastsq(0,ExpMat,Vars,FLAG,5))$ |
|
} |
|
else{ |
|
TMP0=map(ltov,TMP0)$ |
|
|
|
for(I=0;I<length(TMP0);I++) |
|
TMP0[I][0]+=2$ |
|
|
|
TMP0=map(vtol,TMP0)$ |
|
|
|
RET=append(RET,TMP0)$ |
} |
} |
|
|
return RET$ |
return RET$ |