【发布时间】:2021-03-18 16:08:44
【问题描述】:
写了这个函数lockdown_func(beta.hat_func)。
第一件事是:我收到一个错误“参数长度为零”。 第二件事是:当我在没有日期索引的情况下计算它时,它不会改变它应该改变的值,输出向量包含每个索引的相同值。
date= c(seq(from=30, to=165))
beta.hat_func <- c(rep(x = beta.hat, times = 135))
beta.hat <- beta0[which.min(SSE)]
#implement function for modeling
lockdown_func <- function(beta.hat_func,l){
h=beta.hat_func
{
for(i in 1:length(h))
if(date[i]>60 | date[i]<110){
beta.hat_func[i]=beta.hat_func[i]*exp(-l*(date[i]-date[i-1]))
}else{
beta.hat_func[i]=beta.hat_func[i]
}
return(h)
}
}
lockdown_func(beta.hat_func,0.03)
【问题讨论】:
-
你能通过定义
date和beta.hat来让我们重现这个吗? -
当然!刚刚添加它们,beta.hat 是 SSE 估计的结果,是一个值。