Description Usage Arguments Author(s) References See Also Examples
summary
methods for the mdscore
objects
1 2 |
object |
object resulting from a run of the |
... |
not currently used |
Damiao N. da Silva damiao@ccet.ufrn.br
da Silva-Junior AHM, da Silva DN, Ferrari SLP (2014). mdscore: An R Package to Compute Improved Score Tests in Generalized Linear Models. Journal of Statistical Software, 61(2), 1-16., http://www.jstatsoft.org/v61/c02/
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 | library(Sleuth3)
d <- transform(case1102, TLrat = Brain/Liver, Ltime = log(Time),
Lwrat = log((Weight + Loss)/Weight),
Treat = factor(Treatment == "BD",
labels=c("NS", "BD"))
)
fitf <- glm(TLrat ~ Ltime * Treat + Days + Sex + Lwrat
+ Tumor + Treat*Lwrat, data = d,
family = Gamma("log")
)
X <- model.matrix(fitf)
fit0 <- glm(TLrat ~ Ltime * Treat + Sex + Lwrat + Tumor + Days,
data=d, family=Gamma("log"))
test <- mdscore(fit0, X1=X[,9], phi=NULL)
summary(test)
|
Loading required package: MASS
Df Value P-value
Score 1 4.38 0.0363
Modified score 1 3.72 0.0537
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.