calling `bnlr`or `gnlr` in a function

[guiblanchet]

I am interested in using both bnlr and gnlr within a function and I found that a chunk of code that works outside a function, when placed within a function sends an error.

To illustrate this here is a simplified example:

Code outside a function

iniParam <- list(rate = 1)

resp <- rbinom(1, n =70,  p = 0.5)
y <- cbind(resp, 1 - resp)
x <- rep(1:7, 10)

res <- gnlm::bnlr(y = y, link = "logit",
                    mu = ~ exp(-rate * x), 
                    pmu = iniParam)

rm(list = ls())

Code within a function

testFunc <- function(){
  iniParam <- list(rate = 1)

  resp <- rbinom(1, n =70,  p = 0.5)
  y <- cbind(resp, 1 - resp)
  x <- rep(1:7, 10)

  res <- gnlm::bnlr(y = y, link = "logit",
                    mu = ~ exp(-rate * x), 
                    pmu = iniParam)
  
  return(res)
}

testFunc()

I have tried to solve this problem by setting a new environment within the function or using tccov class object, but in both cases a different error arose.

I would much appreciate any insight you may have about solving this problem.

[swihart]

I'll try to look into this by Friday. Thanks!

[swihart]

How does the following suit your needs? I added three lines inside testFunc. (I spent a day looking into other solutions and concluded that solving this one might be quite involved -- I'm admittedly not a wizard with environments and scoping.)

rm(list = ls())


testFunc <- function(){
  iniParam <- list(rate = 1)
  
  resp <- rbinom(1, n =70,  p = 0.5)
  y <- cbind(resp, 1 - resp)
  x <- rep(1:7, 10)
  
  ## add these three lines:
  holdVbls <- list(y=y, x=x, iniParam=iniParam)
  attach(holdVbls)
  on.exit(detach("holdVbls"))
  
  res <- gnlm::bnlr(y = y, link = "logit",
                    mu = ~ exp(-rate * x), 
                    pmu = iniParam)
  
  return(res)
}

testFunc()

[guiblanchet]

Thanks ! This will work for me ! I will let you decide whether you want to keep this thread open or not as a reminder that you may (or not) want to play in this part of the function.