Description Usage Arguments Details Value Note Author(s) Examples
Performs a global test based on the cgOneFactorFit object, to assess whether there are any significant differences amongst levels of the factor, i.e. amongst the groups. A cgOneFactorGlobalTest class object is created.
1 2 | ## S4 method for signature 'cgOneFactorFit'
globalTest(fit, display="print", ...)
|
fit |
A fit object of class |
display |
One of three valid values:
|
... |
Additional arguments. Only one is currently valid:
For other possible |
The notion of a global F test, or equivalently, of R^2,
for resistant & robust linear models is
murky, as no clear theoretical analogue to the ordinary classical
least squares approach exists. The approach taken here is ad-hoc,
which is essentially to re-fit a linear model with lm()
and weights
from the resistant & robust fit. This ad-hoc approach is taken when
there are 3 or more groups.
If there are only 2 groups, then the comparisonsTable.cgOneFactorFit
method is used with the rlm()
model component.
Creates an object of class cgOneFactorGlobalTest
, with the
following slots:
ols.gpval
The p-value of a global F test applied
to the olsfit
component of the cgOneFactorFit
object, unless model="rronly"
is specified. Will not be appropriate in
the case where a valid aftfit
component is present in the
cgOneFactorFit
object.
rr.gpval
The p-value of an ad-hoc global test applied
to the rrfit
component of the cgOneFactorFit
object, if a valid resistant & robust fit object is present.
See the Details section
above. If rrfit
is a simple character value of
"No fit was selected."
, or model="olsonly"
was
specified, then the value is NULL
.
aft.gpval
The p-value of a global chi-square test applied
to the aftfit
component of the cgOneFactorFit
object if a valid accelerated failure time fit object is present.
If aftfit
is a simple character value of
"No fit was selected."
, then the value is NULL
.
uv.gpval
The p-value of a global F test applied
to the uvfit
component of the cgOneFactorFit
object if a valid unequal variances fit object is present.
Contact cg@billpikounis.net for bug reports, questions, concerns, and comments.
Bill Pikounis [aut, cre, cph], John Oleynick [aut], Eva Ye [ctb]
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 | data(canine)
canine.data <- prepareCGOneFactorData(canine, format="groupcolumns",
analysisname="Canine",
endptname="Prostate Volume",
endptunits=expression(plain(cm)^3),
digits=1, logscale=TRUE, refgrp="CC")
canine.fit <- fit(canine.data)
canine.globalTest <- globalTest(canine.fit)
globalTest(canine.fit, model="both")
globalTest(canine.fit, model="olsonly")
globalTest(canine.fit, model="rronly")
data(gmcsfcens)
gmcsfcens.data <- prepareCGOneFactorData(gmcsfcens, format="groupcolumns",
analysisname="cytokine",
endptname="GM-CSF (pg/ml)",
logscale=TRUE)
gmcsfcens.fit <- fit(gmcsfcens.data, type="aft")
globalTest(gmcsfcens.fit)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.