Propagate critical_value through arithmetic operations

Co-Authored-By: Claude Opus 4.6 <noreply@anthropic.com>

Author: nt-williams

DESCRIPTION

@@ -1,7 +1,7 @@
 Package: ife
 Type: Package
 Title: Autodiff for Influence Function Based Estimates
-Version: 0.2.3
+Version: 0.2.4
 Authors@R: 
     c(person(given = "Nicholas",
              family = "Williams",

NEWS.md

@@ -1,3 +1,7 @@
+# ife 0.2.4
+
+* Propagate `critical_value` through arithmetic operations. Binary operations between two estimates use the maximum critical value.
+
 # ife 0.2.3
 
 * Adding `^` method for raising estimates to a power via the delta method.

R/influence_func_estimand.R

@@ -139,11 +139,11 @@ method(print, influence_func_estimate) <- function(x, ...) {
 # x + y
 method(`+`, list(influence_func_estimate, influence_func_estimate)) <- function(e1, e2) {
   check_same(e1, e2)
-  influence_func_estimate(e1@x + e2@x, e1@eif + e2@eif, e1@weights, e1@id)
+  influence_func_estimate(e1@x + e2@x, e1@eif + e2@eif, e1@weights, e1@id, max(e1@critical_value, e2@critical_value))
 }
 
 method(`+`, list(influence_func_estimate, class_numeric)) <- function(e1, e2) {
-  influence_func_estimate(e1@x + e2, e1@eif, e1@weights, e1@id)
+  influence_func_estimate(e1@x + e2, e1@eif, e1@weights, e1@id, e1@critical_value)
 }
 
 # Delegate to ife + scalar since addition is commutative
@@ -152,15 +152,15 @@ method(`+`, list(class_numeric, influence_func_estimate)) <- function(e1, e2) e2
 # x - y
 method(`-`, list(influence_func_estimate, influence_func_estimate)) <- function(e1, e2) {
   check_same(e1, e2)
-  influence_func_estimate(e1@x - e2@x, e1@eif - e2@eif, e1@weights, e1@id)
+  influence_func_estimate(e1@x - e2@x, e1@eif - e2@eif, e1@weights, e1@id, max(e1@critical_value, e2@critical_value))
 }
 
 method(`-`, list(influence_func_estimate, class_numeric)) <- function(e1, e2) {
-  influence_func_estimate(e1@x - e2, e1@eif, e1@weights, e1@id)
+  influence_func_estimate(e1@x - e2, e1@eif, e1@weights, e1@id, e1@critical_value)
 }
 
 method(`-`, list(class_numeric, influence_func_estimate)) <- function(e1, e2) {
-  influence_func_estimate(e1 - e2@x, -e2@eif, e2@weights, e2@id)
+  influence_func_estimate(e1 - e2@x, -e2@eif, e2@weights, e2@id, e2@critical_value)
 }
 
 # x / y
@@ -170,36 +170,36 @@ method(`/`, list(influence_func_estimate, influence_func_estimate)) <- function(
     stop("Division by zero: denominator estimate is zero", call. = FALSE)
   }
   eif <- (e1@eif / e2@x) - ((e2@eif / e2@x^2) * e1@x)
-  influence_func_estimate(e1@x / e2@x, eif, e1@weights, e1@id)
+  influence_func_estimate(e1@x / e2@x, eif, e1@weights, e1@id, max(e1@critical_value, e2@critical_value))
 }
 
 method(`/`, list(class_numeric, influence_func_estimate)) <- function(e1, e2) {
   if (abs(e2@x) < .Machine$double.eps) {
     stop("Division by zero: denominator estimate is zero", call. = FALSE)
   }
-  influence_func_estimate(e1 / e2@x, -e1 / e2@x^2 * e2@eif, e2@weights, e2@id)
+  influence_func_estimate(e1 / e2@x, -e1 / e2@x^2 * e2@eif, e2@weights, e2@id, e2@critical_value)
 }
 
 method(`/`, list(influence_func_estimate, class_numeric)) <- function(e1, e2) {
   if (abs(e2) < .Machine$double.eps) {
     stop("Division by zero: denominator is zero", call. = FALSE)
   }
-  influence_func_estimate(e1@x / e2, 1 / e2 * e1@eif, e1@weights, e1@id)
+  influence_func_estimate(e1@x / e2, 1 / e2 * e1@eif, e1@weights, e1@id, e1@critical_value)
 }
 
 # x * y
 method(`*`, list(influence_func_estimate, influence_func_estimate)) <- function(e1, e2) {
   check_same(e1, e2)
-  influence_func_estimate(e1@x * e2@x, e2@x * e1@eif + e1@x * e2@eif, e1@weights, e1@id)
+  influence_func_estimate(e1@x * e2@x, e2@x * e1@eif + e1@x * e2@eif, e1@weights, e1@id, max(e1@critical_value, e2@critical_value))
 }
 
 method(`*`, list(class_numeric, influence_func_estimate)) <- function(e1, e2) {
-  influence_func_estimate(e1 * e2@x, e1 * e2@eif, e2@weights, e2@id)
+  influence_func_estimate(e1 * e2@x, e1 * e2@eif, e2@weights, e2@id, e2@critical_value)
 }
 
 # x^n
 method(`^`, list(influence_func_estimate, class_numeric)) <- function(e1, e2) {
-  influence_func_estimate(e1@x^e2, e2 * e1@x^(e2 - 1) * e1@eif, e1@weights, e1@id)
+  influence_func_estimate(e1@x^e2, e2 * e1@x^(e2 - 1) * e1@eif, e1@weights, e1@id, e1@critical_value)
 }
 
 # Delegate to scalar * ife since multiplication is commutative
@@ -210,10 +210,10 @@ method(log, influence_func_estimate) <- function(x, base) {
   if (x@x <= 0) {
     stop("log() requires positive values: estimate is ", x@x, call. = FALSE)
   }
-  influence_func_estimate(log(x@x), x@eif / x@x, x@weights, x@id)
+  influence_func_estimate(log(x@x), x@eif / x@x, x@weights, x@id, x@critical_value)
 }
 
 # exp(x)
 method(exp, influence_func_estimate) <- function(x) {
-  influence_func_estimate(exp(x@x), exp(x@x) * x@eif, x@weights, x@id)
+  influence_func_estimate(exp(x@x), exp(x@x) * x@eif, x@weights, x@id, x@critical_value)
 }