Enforce the zero power sum condition

Sometimes a hyper2 object is unbalanced in the sense that its powers do not sum to zero. This is rectified by balance(), which modifies the power of the bracket corresponding to the sum of all pnames accordingly.

balance(H)

Arguments

H

object of class hyper2 or hyper3

Details

This is just a convenience function, all it does is


    H[pnames(H)] <- 0
    H[pnames(H)] <- -sum(pnames(H))
    H

(the first line ensures that H[pnames(H)] is over-written correctly by the second). Package vignette zeropower discusses the zero power sum condition.

Value

Returns a balanced hyper2 object

Author

Robin K. S. Hankin

See also

print.hyper2

Examples



H <- hyper2()
H["a"] <- 6
H["b"] <- 3
H[c("a","c")] <- 7
H <- balance(H)
maxp(H)
#>            a            b            c 
#> 8.124981e-01 1.875009e-01 1.002659e-06