similar to: Multilevel model with lme(): Weird degrees of freedom (group level df > # of groups)

Dear R users, I have a newbie-question that I couldn't resolve after reading through several pieces of documentation and searching the archive. I have a data.frame containing experimental data from a group experiment in psychology. Each line represents a single participant, but participants were assigned to groups of three or four persons. One variable indicates each participants'

Dear R users, I still feel new to R so please apologize if I am doing something stupid here. My use of the polygon() function produces a result that I cannot comprehend: In a plot, I would like to shade the area below a normal distribution. However, I do not want the entire area to be shaded, but just the area on the right side of a vertical line that I draw through the distribution (in

Hello everyone, I am employing the waba method from the multilevel package for obtaining a within-group correlation (Description: http://bg9.imslab.co.jp/Rhelp/R-2.4.0/src/library/multilevel/man/waba.html). Does anybody know a way or a calculation for obtaining a significance value for that correlation? And another question: Does anybody know whether it is possible to save individual

Dear List, I work with multilevel data from psychological group experiments and have frequently encountered a situation for which I haven't found an elegant solution: I need to assign the value of a specific group member to all members of the group. For example, I have a group leader (identified by a binary vector) and some attribute for all group members. I want to create a new

Dear R-Users, I have a trivial problem, but extensive googling and ??'ing did not solve it: I want to obtain the sums of squares from a summary.aov() object for later use. Example: > DV <- rnorm(100) > IV1 <- as.factor(rep(c("male", "female"), times = 50)) > IV2 <- as.factor(rep(c("young", "old"), times = 50)) > >

Hello lmer() - users, A call to the lmer() function causes my installation of R (2.11.1 on Mac OS X 10.5.8) to crash and I am trying to figure out the problem. I have a data set with longitudinal data of four subsequent performance measures of 1133 individuals nested in 88 groups. The data is in long format. I hypothesize a performance increase for each individual over time and intend to

Hello, I have a problem with selecting the right type of sums of squares for an ANCOVA for my specific experimental data and hypotheses. I do have a basic understanding of the differences between Type-I, II, and III SSs, have read about the principle of marginality, and read Venable's "Exegeses on Linear Models" (http://www.stats.ox.ac.uk/pub/MASS3/Exegeses.pdf). I am pretty new to

Dear all, I was wondering if there is any plan to add authentication with OpenPGP key in addition to the DSA/RSA authentication. In the ssh.com version, you got something like that (I don't really remember the correct directive) : PgpSecreKeyfile,PgpKeyID and so on. Is it a planned feature ? Thanks. adulau -- Alexandre Dulaunoy -- http://www.foo.be/ 3B12 DCC2

================================================================ "Because things are the way they are, things will not stay the way they are." Bertolt Brecht ================================================================ Release Announcements ===================== This is the latest bugfix release of the Samba 3.0 series. Please note, that the 3.0 series will be

================================================================ "Because things are the way they are, things will not stay the way they are." Bertolt Brecht ================================================================ Release Announcements ===================== This is the latest bugfix release of the Samba 3.0 series. Please note, that the 3.0 series will be

Hi, I would find it very useful if puppet were able to replicate passwords for specified users from a master PC (the puppetmaster would do me fine, though I suspect this may not suit everyone). That would make changing passwords on my small Linux network a little easier. cheers John Dubery

I successfully switched to using mysql and all the problems relating to sqllite and stored configs, locking issues, etc. went away. However, it seems as if something new has suddenly popped up in my puppetmaster logs and the clients are failing to compile their configurations: info: Found proxy in /etc/puppet/manifests/site.pp err: Puppet::Parser::Interpreter failed with error Mysql::Error:

Hi, I''m new to puppet and trying to get started using it just in local mode. My basic question is: how can i reference a static file within the module (eg in the module''s files/ subdirectory, as described in http://reductivelabs.com/trac/puppet/wiki/ModuleOrganisation). I suppose i could try using a full.absolute path to the file, but that seems to completely defeat the purpose

I''m in a bit of a pickle and I can''t work out a safe way to do something. I currently have a define called ip that I use to define network interfaces. The built in method didn''t seem to work, so I designed this to write out the /etc/sysconfig/network-scripts/ifcfg- ethX file. However, I started thinking about what if I decide to define an eth0:0 for some purpose, and

Hi! I'm pretty raw when working with the R models (linear or not). I'm wondering has anybody worked with the NLME library and dichotomous outcomes. I have a binary outcome variable that I woul like to model in a nested (multilevel) model. I started to fit a logistic model to a NLS function, but could not suceed. I know there are better ways to do it in R with either the LRM or GLM wih

Hi R-Help. I am working on a data set with a 3-level nested structure. I have individuals nested in households and multiple observations on each individual. I assume that the individuals inside a given household are correlated and that the individuals are correlated with themselves over time. The data is not balanced. I have computed a GLMM with logit link function and two random normal

Hi, I'm trying to perform a power simulation for a simple multilevel model, using the function glmmPQL in R version 2.8.1. I want to extract the p-value for the fixed-effects portion of the regression, but I'm having trouble doing that. I can extract the coefficients (summary(fit)$coeff), and the covariance matrix (summary(fit)$varFix), but I can't grab the p-value (or the

Hello, Have someone performed a bootstrap in a multiple-mediator model? I am trying to compute a bootstrap in a multiple and multilevel mediation. Up top now, I have developed bootstraps in random coeffient models, but I am very lost concerning the mediation. Could someone to provide me some ideas about syntax in R? Thank you very much in advance, Bea

Dear all, I posted this two years ago, getting no answers or suggestions - now I am trying again, hoping something new is available in R. I am interested in an application of linear multilevel model with unequal selection probabilities at both levels. Do you know if there is an R function for multilevel pseudo-maximum likelihood estimation? Or is it possible to obtain these estimates using

Hi R folks (Dr. Bates in particular), In August 2005, Dr. Bates mentioned that the documentation for lme4 "will be in the form of a book with the working title 'Multilevel Modeling in R'" and I'm just wondering if there is an estimated date of publication or if it's still a long way off. The Rnews article does a great job of introducing the package, but I'm