Description Usage Arguments Details Value Note See Also Examples
Add (co)variance estimates to a PIM
1 2 | add.varianceestimate(object, varianceestimator = varianceestimator.sandwich(),
estimator = estimator.nleqslv(), verbosity = 0)
|
object |
|
varianceestimator |
Function like the result of |
estimator |
Function like the result of |
verbosity |
The higher this value, the more levels of progress and debug information is displayed (note: in R for Windows, turn off buffered output) |
The result should be the same as you would have gotten if you passed along the
varianceestimator
to the original pim
call.
Important: this function will only work for PIMs where you specified keep.data=TRUE
.
An object of class pim
. The covariance related items are replaced
in object
, or added. Typically, these are vcov
and morevarfitinfo
Most variance estimators do not actually use the estimator
. The only one know
at the time of this writing is varianceestimator.bootstrap
. So in practice
this parameter can probably be mostly ignored.
estimator.nleqslv
varianceestimator.sandwich
pim
1 2 3 4 5 | set.seed(1)
iris$out<-factor(sample(2, nrow(iris), replace=TRUE))
iris$xord<-as.ordered(iris$Species)
pima<-pim(out~Sepal.Length, data=iris, keep.data=TRUE, link="identity")
pimb<-add.varianceestimate(pima, varianceestimator=varianceestimator.H0())
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.