Displaying 1 result from an estimated 1 matches for "squareproduct".
2006 Apr 08
1
add lines to a plot with a loop without erase the last one
...{
data[j,i]<-((alpha*s[j]) / (s[j] + beta)) + (epsilon[j]*(s[j]/(s[j] +
beta)))
b<-data[,i]
}
#print(b)
X<-b[-1]/s[-1]
#print(X)
#print(X)
#print(length(b[-1]))
#print(length(X))
meanX<-mean(X)
#print(meanX)
meanV<-mean(b)
#print(meanV)
betachap[i]<-((meanX*squareproduct(b[-1],b[-1]))-
(meanV*squareproduct(X,b[-1])))/((meanV*squareproduct(X,X))-
(meanX*squareproduct(X,b[-1])))
alphachap[i]<-meanV + betachap[i]*meanX
varchap[i]<-((squareproduct(b[-1],b[-1])) +
(2*betachap[i]*squareproduct(X,b[-1])) +
(((betachap[i]^2)*squareproduct(X,X))))/(length(s)-2)
#u&l...