Displaying 1 result from an estimated 1 matches for "poolednum".
2010 Nov 08
1
try (nls stops unexpectedly because of chol2inv error
...n 1:length(D))
{
print(paste(p1, " [D] of ",length(D))); flush.console();
print(D[p1]);
myReject.pooled = myReject.pooled.1 = MAX.pooled = rep(-1,nsim);
gsim = 0; ## good simulations
for(i in 1:nsim)
{
doubleBreak = F;
print(paste(i, " of ",nsim)); flush.console();
tData = NULL;
pooledNum = matrix(0,nrow=k,ncol=k); ##numerator as weighted sum AS
(n_k-1)cov.scaled
pooledDen = 0; ##denominator as correction AS N-k
#Sigma_pooled = ((omit.1-1)*summary.nls.1$cov.scaled +
(omit.2-1)*summary.nls.2$cov.scaled +
(omit.L-1)*summary.nls.L$cov.scaled)/(sum(omit.1,omit.2,omit.L)-L);
for(j in...