Displaying 1 result from an estimated 1 matches for "sapirstein".
2018 May 05
1
error in chol.default((value + t(value))/2) : , the leading minor of order 1 is not positive definite
...of call
to nlme - the start values are taken from fitted values in brms
library(ggplot2)
windows(record=TRUE)
#generate 3*10? rats - add fixed effects to the four parameters
according to the three groups - add random effects pr each rat - add
residual random effect
#Parameter values taken from Sapirstein AJP 181:330-6, 1955
set.seed(1234)
Time <- seq(1,60,by=1)
A <- 275; B <-? 140;? g1 <- 0.1105; g2 <- .0161
N <- 30
AA <- rep(A,30)+rnorm(30,0,30);BB <- rep(B,30)+rnorm(30,0,15) ;
gg1 <- rep(g1,30)+rnorm(30,0,0.01); gg2 <- rep(g2,30)+rnorm(30,0,0.001)
TRT <- gl(3...