索赔金额模拟每个合约的不同lambda

Question

我正在模拟R中的索赔数据，我之前使用了下面的代码，我假设索赔的数量是Poisson分布参数lambda，我也假设这个lambda从一个契约变化到另一个，因此它跟随了Gamma分布;

[<-

现在不使用变化的lambda（rgamma（1,8）），我希望使用下面的向量为每个合约使用不同的lambda;

install.packages("actuar")
library("actuar")
portfolio.1<-simul(nodes=list(contract=10,year=5),model.freq=expression(contract=rgamma(1,8),year=rpois(contract)),model.sev=expression(contract = NULL,year = rlnorm(6,1)))
sev.1 <- severity(portfolio.1)
freq.1 <- frequency(portfolio.1)
agg.1 <- round(aggregate(portfolio.1),2)

所以每个lambda对应于每个契约，但我没有在R中编写代码。 愿你协助。