consistency | R Documentation |
Given a hyper2
object, calculate the maximum likelihood point in
two ways and plot one against the other to check for consistency.
consistency(H, plot=TRUE, ...)
H |
A |
plot |
If |
... |
Further arguments, passed to |
Given a hyper2
object, calculate the maximum likelihood
estimate of the players' strengths using maxp()
; then reverse
the pnames
attribute and calculate the players' strengths
again. These two estimates should be identical but small differences
highlight numerical problems. Typically, the differences are small if
there are fewer than about 25 players.
Reversing the pnames()
is cosmetic in theory but is a
non-trivial operation: for example, it changes the identity of the
fillup from the last player to the first.
Returns a named three-row matrix with first row being the direct evaluate, second row being the reverse of the reversed evaluate, and the third being the difference
Robin K. S. Hankin
ordertrans
# consistency(icons)
x <- icons
y <- icons
pnames(y) <- rev(pnames(y))
gradient(x,indep(equalp(x)))
gradient(y,indep(equalp(y)))
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.