R语言逐时气温估算函数运行报argument "A"缺失错误如何解决

错误原因 #

• 你的函数定义要求传入3个独立参数Ta、A、t，但调用时仅传入了1个长度为3的向量c(28,5,7)，R会默认把这个向量赋值给第一个参数Ta，剩下的A、t均无传入值，因此触发参数缺失报错。
• 现有函数的实现逻辑和你给出的气温计算公式T(t) = Ta + Asin[w(t - to)]不匹配，计算结果会不符合预期。

修复步骤 #

1. 修正函数逻辑对齐公式 #

Toft <- function(Ta, A, t, to = 8, w = 2*pi/24){
  # Ta：日平均气温
  # A：气温正弦波动的振幅（即日最高温与平均温的差值）
  # t：待计算的时刻，单位为小时（取值范围0~23）
  # to：日最低温出现的时刻，默认设为8点可按需调整
  # w：角频率，固定为2π/24对应24小时周期
  Ta + A * sin(w * (t - to))
}

2. 正确调用函数 #

单次单时刻计算

如果要计算Ta=28、A=5、t=7对应的气温，直接传入三个独立参数即可：

Toft(Ta = 28, A = 5, t = 7)

批量逐时计算

如果需要计算一天24小时的逐时气温，直接给t传入0~23的序列即可：

# 生成逐时气温结果
hourly_temp <- Toft(Ta = 28, A = 5, t = 0:23)
# 转成数据框方便查看
temp_df <- data.frame(小时 = 0:23, 预估气温 = round(hourly_temp, 2))
head(temp_df)

向量传参适配

如果你的参数确实存储在一个向量里，可以用do.call拆分向量给对应参数：

param_vec <- c(Ta = 28, A = 5, t = 7)
do.call(Toft, as.list(param_vec))

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