keep_age <- which(rowSums(portfolio_LTC$ec) > 5e2)
keep_duration <- which(colSums(portfolio_LTC$ec) > 1e3)
d <- portfolio_LTC$d[keep_age, keep_duration]
ec <- portfolio_LTC$ec[keep_age, keep_duration]
fit <- WH_2d(d, ec)
plot(fit)
plot(fit, "std_y_hat")
Run the code above in your browser using DataLab