我想最小化建模和观察到的价差之间的均方误差(可能使用包中mse()
的)。hydroGOF
函数定义为:
KV_CDS <- function(Lambda, s, sigma_S){
KV_CDS = (Lambda * (1 + s)) / exp(-s * sigma_S) - Lambda^2)
}
目标是mse
通过在 KV_CDS 函数中为 Lambda 保留一个自由参数来最小化 KV_CDS 和 C 之间的值。
df <- data.frame(C=c(1,1,1,2,2,3,4),
Lambda=c(0.5),s=c(1:7),
sigma_S=c(0.5,0.4,0.3,0.7,0.4,0.5,0.8),
d=c(20,30,40,50,60,70,80),
sigma_B=0.3, t=5, Rec=0.5, r=0.05)