负逆链接函数下线性预测器均值响应计算与绘图求助

背景 #

已完成Question 3中logit链接函数的计算与绘图，参考代码：

#Q3
LinearPredictor <- seq(-5,5,0.01)
head(LinearPredictor)
# compute inverse of link function
link_func <- function(x) exp(x)

# Compute mean response
MeanResponse<- as.vector(link_func(LinearPredictor))

# Plot the values
plot(LinearPredictor, MeanResponse)

当前需求：线性预测器取值范围-5到5、步长0.01，基于负逆链接函数计算并绘制均值响应与线性预测器的关系图。尝试代码link_func <- 1/-(function(x) exp(x))时，触发二元运算符错误。

错误原因 #

R中不能直接对函数对象执行算术运算（除法、取负），你写的代码把function(x) exp(x)当成数值来操作，不符合语法规则。必须把运算逻辑封装在函数内部，让函数接收数值输入后执行对应计算。

正确实现代码 #

假设负逆链接函数定义为 -1/exp(x)（若实际是1/(-exp(x))，逻辑完全一致），正确代码如下：

# 定义负逆链接函数
link_func <- function(x) {
  -1 / exp(x)
}

# 生成线性预测器序列
LinearPredictor <- seq(-5, 5, 0.01)

# 计算均值响应
MeanResponse <- link_func(LinearPredictor)

# 绘制关系图
plot(LinearPredictor, MeanResponse, 
     main = "负逆链接函数：均值响应 vs 线性预测器",
     xlab = "线性预测器", 
     ylab = "均值响应")

说明 #

函数内部的运算顺序可根据你对“负逆链接函数”的实际定义调整，核心是让函数接收LinearPredictor的每个数值元素，逐一执行算术运算，避免直接对函数对象做数值操作。

内容的提问来源于stack exchange，提问作者slycooper