similar to: standard error of variances and covariances of the random effects with LME

With lme4, use of mcmcsamp can be insightful. (Douglas Bates drew my attention to this function in a private exchange of emails.) The distributions of random effects are simulated on a log scale, where the distributions are much closer to symmetry than on the scale of the random effects themselves. As far as I can see, this is a straightforward use of MCMC to estimate model parameters; it is not

Using ranef() (package nlme, version 3.1-75) with an 'lme' object I can obtain random effects for intercept and slope of a certain level (say: 1) - this corresponds to (say level 1) "residuals" in MLWin. Maybe I'm mistaken here, but the results are identical. However, if I try to get the standardized random effects adding the paramter "standard=T" to the

Hello, I'm running a simulation study of a multilevel model with binary response using the binomial probit link. It is a random intercept and random slope model. GLMMPQL and lmer fail to converge on a *significant* portion of the *generated* datasets, while MlWin gives reasonable estimates on those datasets. This is unacceptable. Does anyone has similar experiences? Regards, Roel de

Dear R experts, I had a question which may not be directly relevant to R but I will be grateful if you can give me some advices. I ran a two-level multilevel model for data with repeated measurements over time, i.e. level-1 the repeated measures and level-2 subjects. I could not get convergence using lme(), so I tried MLwiN, which eventually showed the level-2 variances (random effects for

When playing around fitting unconditional growth models using R and MlwiN today, I produced two different sets of estimates that I can't reconcile and wondered if anyone here has an idea: The data is two-level repeated measures data with measures nested within child. There are two measures per child. I've fit an unconditional growth model as in Singer and Willet (2003) that allows for

Hello and thanks in advance for any help, I am using the 'multinom' function from the nnet package to calculate a multinomial logistic regression. I would like to get a matrix estimates of the estimated coefficient variances and covariances. Am I missing some easy way to extract these? Grant [[alternative HTML version deleted]]

Dear helpeRs, I would like to specify a newline command in R and pass it to latex via Sweave such that it corresponds to latex' \\ command. But that doesn't seem to be possible. If I Sweave the \n character, it just makes a new line in latex but not the \\ command. Is there a way such that the following code would result in latex in blablabla \\ blablabla on different line

Forgive me if I'm posting to the wrong place....It's my first time posting. Here's the situation: I'm using the sem package and making path diagrams using path.diagrams. Suppose I have the following code: #install.packages("ggm") require(ggm) cor = rcorr(7) nm = c("SOF", "IWF", "PWF", "FSC", "FSF", "EF",

Hi, I compare results of a simple two-level poisson estimated using lmer and those estimated using MLwiN and Stata (v.9). In R, I trype: ------------------------------------------------------------------------------------------- m2 <- lmer(.D ~ offset(log(.Y)) + (1|pcid2) + educy + agri, male, poisson) -------------------------------------------------------------------------------------------

Full_Name: Alex Galakhov Version: 2.1.0 OS: Linux (Debian) Submission from: (NULL) (195.19.131.68) After software upgrade source() does not work properly anymore. It completely ignores the chdir= parameter. This is because is.character(file) is always false after the assignment file<-file(file). The fix is (written by hand in diff-like syntax, hopefully, you'll understand it):

Hi! First of all, Im posting from Austria. Alle respondes of my server are translated to german, and if my question ist not in good english please forgive me. :-) Were runnig an asterisk und searching for a way to fix a question from one of our clients. After calling the number and no response of our client the voice-box gives response. Thats ok... but after the voice-box, which ist self-

> To sum up, I can't figure out how MLWin calculates the > standardized residuals. But I understand that this is not a > question for the R list. > Nevertheless, it would help if someone could point me to some > arguments why not to use them and stick to the results > obtainable by ranef(). Hi Dirk: Well, it is interesting that mlWin and lmer generate the same exact

Hello, has someone written by chance a function to extract the variance-covariance matrix from a lmer-object? I've noticed the VarCorr function, but it gives unhandy output. Regards, Roel de Jong

Hello I am running a net work with a NT4PDC a NT4 BDC and attent to add samba BDC. I am following the brief section 28 of the Samba-Howto collection for samba head/3.0. * Adding the samba server as a BDC from NT server manager seems OK * Getting, recording domain sid is ok (see below) > rpcclient NT4PDC -U Administrator password: > lsaquery domain IMAGERIE has sid blablabla > net

Hi there, I'm facing the decision if it would be possible to transform several more or less complex pdf files into an R Table-Format or if it has to be done manually. I think it would be a impudent to expect a complete solution, but I would be grateful if anyone could give me an advice on how the structure of such a R-program could look like, and if it's possible in general. Here

I'm depending on scp to create the .ssh directory (which is does) but it can't use the directory it just created in the current session thus the "error" message. I would like to run the command just once if possible :) Is it possible I could get a patch if this is changed? [root /root]# ls -la total 84 drwxr-x--- 9 root root 4096 Dec 17 19:18 . drwxr-xr-x 16

> > > > > > > > Hi Folks, > > > > > > > > > > > > > > > > I’m having a problem getting a the SSL cert file formatted > > > > > > > > just like icecast wants… I’m running 2.4.2 … and it > > > > > > > > doesn’t seem to want to use my combined key + cert chain > > > >

Dear list members, Where can I find code for computing the p*p variance-covariance matrix given a vector of p variances (ordered varA, varB, ..., varp) and a vector of all possible correlations (ordered corAB, corAC, ..., corp-1,p)? I know that the covariance between 2 variables is equal to the product of their correlation and their standard deviations: corAB * varA^.5 * varB^.5 and so:

Does anyone know of code to conduct hierarchical (that is, multi-level) models using R. Beyond simply requiring a nested design, I want to model explicitly the covariance between levels as is done in such multi-level modeling software as HLM or MLwin and discussed in Goldestein (1999) available online at http://www.arnoldpublishers.com/support/goldstein.htm (a nice and free resource for anyone

Dear R-list, maybe someone can help me with the following mixed-effects models problem, as I am unable to figure it out with the 'nlme-bible'. I would like to fit (in R, obviously) a so-called animal model (google e. g. "Heritability and genetic constraints of life-history" by Pettay et al.) to estimate the variance component that is due to genetic effects. The covariances of