Displaying 2 results from an estimated 2 matches for "auxregs".
2011 Dec 19
1
pls help to print out first row of terms(model) output in example program
...in ivnames) {
fmla <- paste( "`", i, "`", " ~ ." , sep="");
print(fmla)
maux <- lm( formula(fmla), data=dat)
results[[ i ]] <- maux
print(maux$call[2])
###fix me fix me ##
}
results
}
mcDiagnose <- function(model){
auxRegs <- lmAuxiliary(model)
auxRsq <- numeric(length=length(auxRegs))
j <- 0
for ( i in auxRegs ){
j <- j + 1
sumry <- summary(i)
auxRsq[j] <- sumry$r.squared
}
print("Drum roll please!")
print("And your R_j Squareds are (auxiliary Rsq)")
pr...
2007 Apr 10
1
When to use quasipoisson instead of poisson family
It seems that MASS suggest to judge on the basis of
sum(residuals(mode,type="pearson"))/df.residual(mode). My question: Is
there any rule of thumb of the cutpoiont value?
The paper "On the Use of Corrections for Overdispersion" suggests
overdispersion exists if the deviance is at least twice the number of
degrees of freedom.
Are there any further hints? Thanks.
--
Ronggui