datasetG <- groupedData(p ~ t | experiment, data = dataset)
mymodel <- nlme(p ~ a/(1+b * exp(-k*t)),
data = datasetG,
fixed = a + b + k ~ 1,
random = a + k ~ 1,
start = c(a = mychutes_a, b = mychutes_b , k = mychutes_k))
simulate_mymodel <- simulate_nlme(mymodel, nsim = 1000, psim = 2, value = "data.frame")
i got the sim.y values but I wanted the a, b, k values too. Does anyone know how I can get these values?
a, b, k values mixed model
