search for: mcpar

This fairly harmless looking piece of code, which worked in 0.61 (and in S-PLUS) fails in 0.62.1. R> "mcpar<-" <- function(x,mcpar) {attr(x,"mcpar") <- mcpar; x} R> mcpar(x) <- c(1,100,1) Error in mcpar<-(*tmp*, value = c(1, 100, 1)) : unused argument to function The error message gives a hint about how to work around the problem - just replace the argument "mcpar&q...

Hi all, I have managed to get the r-scale usage algorithm to work, but I need to obtain the final results from this. As I understand it, this code is designed to generate a matrix after processing and store it somewhere????? Here is the code..... I get this part of the code, it all makes sense. ## if(nchar(Sys.getenv("LONG_TEST")) != 0) {R=1000} else {R=5} { data(customerSat)

...s. For example: > set.seed (1) > group <- rep (1:5,10) > a <- rnorm (5,-3,3) > y <- rnorm (50, a[group], 2) > fit <- lmer (y ~ 1 + (1 | group)) > mcmcsamp(fit) (Intercept) log(sigma^2) log(grop.(In)) [1,] -2.771979 0.6909418 1.750876 attr(,"mcpar") [1] 1 1 1 attr(,"class") [1] "mcmc" I want "grop" in the above output to be "group", and I want "(In)" to be "(Intercept)". For my purposes (fitting models in general settings) it is important to have descriptive variable names...

...102 2.755102 2 0.616483 3.640738 I can't figure out why everything is rounded for the first model but not for the second. Also, mcmcsamp() works for M1 but not for M2: > mcmcsamp(M1) (Intercept) log(sigma^2) log(grop.(In)) [1,] 9.099073 0.5711817 3.246981 attr(,"mcpar") [1] 1 1 1 attr(,"class") [1] "mcmc" > mcmcsamp(M2) Error: inconsistent degrees of freedom and dimension Error in t(.Call("mer_MCMCsamp", object, saveb, n, trans, PACKAGE = "Matrix")) : unable to find the argument 'x' in selecting...

Dear All, With kind help from several friends on the list, I am getting close. Now here are something interesting I just realized: for random effects, lmer reports standard deviation instead of standard error! Is there a hidden option that tells lmer to report standard error of random effects, like most other multilevel or mixed modeling software, so that we can say something like "randome