我在RStudio中使用Shiny工作,我一直在尝试运行这段代码。
# Vamos a simular modelos poisson compuestos con diferentes
# distribuciones de severidad.
library(actuar)
library(shiny)
# Define UI for application that draws a histogram
ui <- fluidPage(
# Número de simulaciones
headerPanel('Número de simulaciones'),
numericInput(inputId = "n",label = NULL,
value = 10,min = 1,max = 20000),
# Gráfica S exponencial.
headerPanel('Exponencial'),
sidebarPanel(
sliderInput(inputId = 'lambda1',label = 'Lambda 1',value = 7,
min = 0, max = 15),
sliderInput(inputId = 'rate',label = 'rate',value = 5,
min = 0, max = 15),
),
mainPanel(
plotOutput('plot1')
)
)
# Define server logic required to draw a histogram
server <- function(input, output) {
output$plot1 <- renderPlot({
#Preparativos para los gráficos:
set.seed(20)
n <- input$n
lambda1 <- input$lambda1
rate <- input$rate
S1 <- rcompound(n = n, #Genera n
model.freq = rpois(lambda1), #N~Poi(lambda1)
model.sev = rexp(rate = 2)) #Y~Exp(rate)
MASS::truehist(S1,
col=rainbow(125, start = 0.5, 1),
main = "exp",nbins = 125)
abline(h=0,v=0,col="black",lwd=2)
})
}
# Run the application
shinyApp(ui = ui, server = server)
我的问题是在服务器函数中,不管我怎么做,它一直说找不到lambda1,但我在sliderInput中定义了lambda1!另一方面,这段代码可以完美地工作。
n = 1234
#Exponencial
lambda1 <- 7 ; rate <- 5
S1 <- actuar::rcompound(n = n, #Genera n
model.freq = rpois(lambda = lambda1), #N~Poi(lambda1)
model.sev = rexp(rate = rate)) #Y~Exp(rate)
这到底是怎么回事?我已经试了好几天了!
我看了一下下面的代码 rcompound
但我没有找到如何解决这个问题。
一种可能性是将 lambda1
到全球环境。
lambda1 <<- input$lambda1
但这是在全球环境中写的... ...
我只想用我自己的 rcompound
函数,而不是使用一个包。这并不复杂。
rcompound <- function(n, lambda, rate){
N <- rpois(n, lambda)
vapply(N, function(k) sum(rexp(k, rate = rate)), numeric(1))
}