search for: fratio

...model<-lm(yi~xi) stands[st]<-Stand[i-1] intercepts[st]<-model$coefficients[1] slopes[st]<-model$coefficients[2] mses[st]<-sum(resid(model)^2)/(length(xi)-2) ssr<-var(yi)*(length(xi)-1)-sum(resid(model)^2) rsquares[st]<-ssr/(var(yi)*(length(xi)-1)) fRatio<-ssr/mses[st] pValues[st]<-1-pf(fRatio,1,length(xi)-2) #Increment the stand number, zero the within stand collections, #and start again st<-st+1 xi<-numeric(0) yi<-numeric(0) xi[1]=X[i] yi[1]=Y[i] } } #Make your data set standEstimates<-data.frame...

...(cov)))-1 probp = pchisq(Gp, dfp, lower.tail=F) #total G Gt = sum(G) dft = sum(dfs) probt = pchisq(Gt, dft, lower.tail=F) #heterogeneity G Gh = Gt-Gp dfh = dft - dfp probh = pchisq(Gh, dfh, lower.tail=F) #F ratio Fratio = (Gp/dfp)/(Gh/dfh) probf = pf(Fratio, dfp, dfh, lower.tail=F) res = list(TotG=c(Gt, dft, probt), PooledG=c(Gp, dfp, probp), HetG=c(Gh, dfh, probh), Fratio=c(Fratio, dfp, dfh, probf)) res } #end hetg ______________________________________________________________...