diff --git a/R/z_link-map-quap.r b/R/z_link-map-quap.r index 4eb4924..73632cc 100644 --- a/R/z_link-map-quap.r +++ b/R/z_link-map-quap.r @@ -52,7 +52,7 @@ function( fit , data , n=1000 , post , refresh=0 , replace=list() , flatten=TRUE # hacky solution -- find density function and insert whatever expression in typical link spot flik <- as.character(fit@formula[[1]][[3]][[1]]) # mu for Gaussian - if ( flik=="dnorm" ) lm <- as.character( fit@formula[[1]][[3]][[2]] ) + if ( flik=="dnorm" ) lm <- deparse( fit@formula[[1]][[3]][[2]] ) # p for binomial -- assume in third spot, after size if ( flik=="dbinom" ) lm <- as.character( fit@formula[[1]][[3]][[3]] ) # lambda for poisson