# Compute the power-transformed identification function.
df <- data.frame(
y = rep(x = 2, times = 3),
x = 1:3,
a = c(1, 2, 3)
)
df$meanpower_if <- meanpower_if(x = df$x, y = df$y, a = df$a)
print(df)
Run the code above in your browser using DataLab