similar to: lme4 / HLM question

Dear R listers, I know I'm breaking the rules by asking a "homework" related question-- I hope you'll forgive me. I am a social psychology graduate student, and the only one in my department who uses R. I successfully completed my multiple regression and structural equation modeling courses using R (John Fox's car and sem packages were a big help, as was his book).

Dear R-users, I'm running a lmer model using the lme4 package. My dependent variable is dichotomous and I'm using the "binomial" family. The results are slightly different from the HLM results based on a Bernoulli distribution. I read that a Bernoulli distribution is an extension of a binomial distribution. Is that right? If so, how can I adapt my R model to a Bernoulli

Hi R users! I am trying to learn some multilevel analysis, but unfortunately i am now very confused. The reason: http://www.ats.ucla.edu/stat/hlm/seminars/hlm_mlm/mlm_hlm_seminar.htm http://www.ats.ucla.edu/stat/sas/seminars/sas_mlm/mlm_sas_seminar.htm and MlmSoftRev. pdf from mlmRev package. >From what i see, the first two links seem to declare the level one variable as a random part (i

Hi I am trying to convert SAS codes to R, but some of the result are quite different from SAS. When I ran proc mixed, I have an option ddfm=bw followed by the model. How can I show this method in R?(I am thinking that this maybe the reason that I can't get the similar results) below is my SAS codes: proc mixed data=test covtest empirical; class pair grade team school; model score = trt

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

Hello dear uesRs, I am working on modeling both level one and level two heteroscedasticity in HLM. In my model, both error variance and variance of random intercept / random slope are affected by some level two variables. I found that nlme is able to model heteroscedasticity. I learned how to use it for level one heteroscedasticity but don't know how to use it to model the level

If you don't get a response it is because you did not read the Posting Guide which indicates that the R-sig-ME mailing list is where this question would have been on-topic. -- Sent from my phone. Please excuse my brevity. On August 16, 2017 6:17:03 AM PDT, b88207001 at ntu.edu.tw wrote: >Hello dear uesRs, > >I am working on modeling both level one and level two

A better place for this post would be on R's mixed models list: r-sig-mixed-models . Cheers, Bert Bert Gunter "The trouble with having an open mind is that people keep coming along and sticking things into it." -- Opus (aka Berkeley Breathed in his "Bloom County" comic strip ) On Wed, Aug 16, 2017 at 6:17 AM, <b88207001 at ntu.edu.tw> wrote: > Hello dear

Dear all, I have a question regarding setIs and method dispatch in S4 classes: Consider the following example: ##################################################### ## A02 "contains" A01 by setIs: setClass("A01", representation(a="numeric",b="numeric")) setClass("A02", representation(a="numeric",d="numeric"))

I seem to come to this problem alot, and I can find my way out of it with a loop, but I wish, and wonder if there is a better way. Here's an example (lmer1-5 are a series of lmer objects): bs=data.frame(bic=BIC(lmer1,lmer2,lmer3,lmer4,lmer5)$BIC) rownames(bs)=c('lmer1','lmer2','lmer3','lmer4','lmer5') best=rownames(bs)[bs==min(bs)] > best [1]

Hi, I was wondering if I might be able to ask some advice about doing residual plots for the lmer function in the lme4 package. Our group's aim is to find if the expression staining of a particular gene in a sample (or "core") is related to the pathology of the core. To do this, we used the lmer function to perform a logistic mixed model below. I apologise in advance

Hi everyone, I have posted a similar question to this list, but I don't get a reply. I really want to solve this problem, so I post it again... I am trying to use R to fit some mixed-effects models for a nested data. The data is a simulated data with 111 subjects. Each subject has 6 waves' data. Below are the first two subjects' data : > simu1[1:12,] Grouped Data: gf ~ age |

Dear list members, the following hlm was constructed: hlm <- groupedData(laut ~ design | grpzugeh, data = imp.not.I) the grouped data object is located at and can be downloaded: www.anicca-vijja.de/lg/hlm_example.Rdata The following works: library(nlme) summary( fitlme <- lme(hlm) ) with output: ... AIC BIC logLik 425.3768 465.6087 -197.6884 Random effects:

Suppose I have a mixed-effects model where yij is the jth sample for the ith subject: yij= beta0 + beta1(age) + beta2(age^2) + beta3(age^3) + beta4(IQ) + beta5(IQ^2) + beta6(age*IQ) + beta7(age^2*IQ) + beta8(age^3 *IQ) +random intercepti + eij In R how can I get an F test against the null hypothesis of beta6=beta7=beta8=0? In SAS I can run something like contrast age*IQ 1,

Hello, This is a follow up question to my previous one http://tolstoy.newcastle.edu.au/R/e4/help/08/02/3600.html I am attempting to model relationship satisfaction (MAT) scores (measurements at 5 time points), using participant (spouseID) and couple id (ID) as grouping variables, and time (years) and conflict (MCI.c) as predictors. I have been instructed to include random effects for the

Full_Name: Zivan Karaman Version: 2.7.2 OS: Windows XP Submission from: (NULL) (195.6.68.214) I'm puzzled by the readRegistry function. Shouldn't the "hive" argument be something like c("HLM", "HCR", "HCU", "HU", "HCC", "HPD") rather than c("HLM", "HCR", "HCU", "HU",

data from multilevel units,first sample the class ,and then the student in calss.following is the 2-level model. and the level-1 model deals with the student,and the level-2 model deals with the class level the students belong to. Level-1 Model Y = B0 + B1*(ZLEAD) + B2*(ZBUL) + B3*(ZSHY) + R Level-2 Model B0 = G00 + U0 B1 = G10 + G11*(ZWARMT) + U1 B2 = G20 + G21*(ZWARMT) + G22*(ZABLET) +

Hi all, Thanks to all who offered advice on using F95 in R. Now I'm trying to compile a test package using gfortran, Linux 2.4.21 and R 2.1.0. I was able to successfully compile and use a test F95 routine by setting my environment variables as follows in bash: export PATH=~/bin/:$PATH export F77=gfortran export LD_LIBRARY_PATH=~/bin/irun/lib export GFORTRAN_STDIN_UNIT=-1 Now I'm

Dear list, I am currently working with a rather large data set on body temperature regulation in wintering birds. My original model contains quite a few dependent variables, but I do not (of course) wish to keep them all in my final model. I've fitted the following model to the data: >

> On Tue, 12 Sep 2000, Magill, Brett wrote: > > > 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) > > The lme()