search for: bvar

...the ones plotted on page 14 of the document available at http://www.eric.ed.gov/ERICDocs/data/ericdocs2/content_storage_01/0000000b/80/2b/b3/94.pdf). Harold Doran mentioned (http://tolstoy.newcastle.edu.au/~rking/R/help/05/08/10638.html) that the posterior modes' variances can be found in the bVar slot of lmer objects. However, when I fit e.g. this model: lmertest1<-lmer(mathtot~1+(m_escs_c|schoolid),hlmframe) then lmertest1 at bVar$schoolid is a three-dimensional array with dimensions (2,2,28). The factor schoolid has 28 levels, and there are random effects for the intercept and m_escs...

...the following would produce confidence intervals for the estimated BLUPs from a linear mixed effect model: OrthoFem<-Orthodont[Orthodont$Sex=="Female",] fm1OrthF. <- lmer(distance~age+(age|Subject), data=OrthoFem) fm1.s <- coef(fm1OrthF.)$Subject fm1.s.var <- fm1OrthF. at bVar$Subject fm1.s0.s <- sqrt(fm1.s.var[1,1,]) fm1.s0.a <- sqrt(fm1.s.var[2,2,]) fm1.s[,1]+outer(fm1.s0.s, c(-2,0,2)) fm1.s[,2]+outer(fm1.s0.a, c(-2,0,2)) [,1] [,2] [,3] [1,] 11.08578 14.48469 17.88359 [2,] 13.86631 17.26521 20.66412 [3,] 13.37444 16.77334 20.17224 [4,] 13.557...

...er of variables increases. Is there a way I can do this kind of processing using a loop? I tried defining a vector to hold the names for the "c variables" (e.g. c1,c2, ... cn) and creating new variables in a loop using code like: avars<-c("a1","a2","a3") bvars<-c("b1","b2","b3") cvars<-c("c1","c2","c3") for(i in 1:3){ df$cvars[i]<-df$avars[i]/df$bvars[i] } But the variable references don't resolve properly with this particular syntax. Any help would be much appreciated. Che...

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

Using lmer in the lme4 package, you can compute the conditional variance-covariance matrix of the random effects using the bVar slot: bVar: A list of the diagonal inner blocks (upper triangles only) of the positive-definite matrices on the diagonal of the inverse of ZtZ+Omega. With the appropriate scale factor (and conversion to a symmetric matrix) these are the conditional variance-covariance matrices of the random effect...

I want to compute confidence intervals for the random effect estimates for each subject. From checking on postings, this is what I cobbled together using Orthodont data.frame as an example. There was some discussion of how to properly access lmer slots and bVar, but I'm not sure I understood. Is the approach shown below correct? Rick B. # Orthodont is from nlme (can't have both nlme and lme4 loaded at same time!) # OrthoFem<-Orthodont[Orthodont$Sex=="Female",] # http://tolstoy.newcastle.edu.au/R/help/06/03/23758.html library(lme4)...

..., Harold: Thanks for the clarification. I thought I had read the original post. Obviously, I had misread it. Thanks again. spencer graves Doran, Harold wrote: > Yes, it is a different issue. ranef() extracts the empirical Bayes > estimates, which are the empirical posterior modes. The bVar slot > holds the corresponding posterior variances of these modes. > > Technically, (according to D. Bates) the values in the bVar slot are > the the diagonal elements of (Z'Z+\Omega)^{-1}. > > The original post was asking how to test and compare a specific random > ef...

I have not been able to find any programs for running vector autoregressions with R. I am interested in running Bayesian VARs and also running VARs that run all combinations of variables in the vector. Is anyone currently developing this? -Nirav Mehta

...&& BufPtr[2] == 0xc0 && BufPtr[3] == 0xde; } Second, I change ParseBitcodeInto as this: original: http://llvm.org/docs/doxygen/html/BitcodeReader_8cpp_source.html#l01971 ... errs() << "parsebitcodeinto sniff the signature\n"; uint32_t bvar = Stream.Read(8); errs() << "B :" << bvar << "\n"; if (bvar != 'B') { errs() << "B :" << bvar << "\n"; return Error(InvalidBitcodeSignature); } if (Stream.Read(8) != 'C') { errs() << &quot...

...; BufPtr[3] == 0xde; > } > > > Second, I change ParseBitcodeInto as this: > original: http://llvm.org/docs/doxygen/html/BitcodeReader_8cpp_source.html#l01971 > ... > errs() << "parsebitcodeinto sniff the signature\n"; > uint32_t bvar = Stream.Read(8); > errs() << "B :" << bvar << "\n"; > if (bvar != 'B') { > errs() << "B :" << bvar << "\n"; > return Error(InvalidBitcodeS...

...>> } >> >> >> Second, I change ParseBitcodeInto as this: >> original: http://llvm.org/docs/doxygen/html/BitcodeReader_8cpp_source.html#l01971 >> ... >> errs() << "parsebitcodeinto sniff the signature\n"; >> uint32_t bvar = Stream.Read(8); >> errs() << "B :" << bvar << "\n"; >> if (bvar != 'B') { >> errs() << "B :" << bvar << "\n"; >> return Error...

Hello, I having a weird problem while writing a bytecode module to a string, and after read/parse it for unsing on a jit. I write a pass to export function to module, and put this module inside a global variable. I use WriteBitcodeToFile for this. For debuging, after this write, I try to load the exported module with parseBitcodeFile. This two step works. After, while the compiled program is

...> >>>> <http://llvm.org/docs/doxygen/html/BitcodeReader_8cpp_ >>>> source.html#l01971> >>>> ... >>>> errs() << "parsebitcodeinto sniff the signature\n"; >>>> uint32_t bvar = Stream.Read(8); >>>> errs() << "B :" << bvar << "\n"; >>>> if (bvar != 'B') { >>>> errs() << "B :" << bvar << &quot...

attempting to do multivariate modelling in R with known future conditions (in this case variable 'b') using MSBVAR and hc.forecast. The sample code (a paired down representation) does not give anywhere near the expected results - I am assuming that a forecast 8 steps out would approximate 'a' as the sequence 1.1,2.1,3.1,100.1 corresponding to the input set. I have varied the input length to be longer...

Hi all, I have the following code, When I run the code, it never terminate this is because of the while loop i am using. In general, if you need a loop for which you don't know in advance how many iterations there will be, you can use the `while' statement so here too i don't know the number how many iterations are there. So Can some one suggest me whats going on? I am using the

I can't find a way to retrieve z values and p-values from the output from lmer in the lme4 package. How is this done? Rick B.

Hi all, I got another VAR question here and really appreciate if somebody would help me out :) I have five time series, say A,B,C,D,E. My objective is to predict the series A using the rest, that is, B, C, D and E. A Vector Autoregression Model should work here. But first of all, I should select which series of B, C, D and E to be include in the VAR model, as well as the number of lags. I wonder

Hi, I have been having a tedious issue with lmer models with lots of factors and lots of levels. In order to get the basic information at the beginning of the print out I also have to generate these enormous tables as well. Is there a method command to leave off all of the effects and correlations? Or, do I have to go to string commands?

Hello All, I would like to retrieve some of the results from the lmer(...) function in library lme4. If I run a model, say fm.1 <- lmer(y ~ 1 + (1 | x), data = dog) and try names(fm.1), I get NULL. Is there anyway to retrieve the information? Thanks

Hi, I would like to know whether it is possible to obtain a value of significance for random effects when aplying the lme or related functions. The default output in R is just a variance and standard deviation measurement. I feel it would be possible to obtain the significance of these random effects by comparing models with and without these effects. However, I'm not used to perform