similar to: multilevel (hierarchical) statistical analysis

Hi, does anyone out there have a recommendation for multilevel / random effects and longitudinal analysis? My dream book would be something that's both accessible to a non-statistician but rigorous (because I seem to be slowly turning into a statistician) and ideally would use R. Peter

After upgrading to 1.7.0 under debian linux, I can't get e1071 working properly. The first problem I had was that g++-3.0 was the standard compiler but wasn't installed, so I installed it. e1071 then installed correctly, but I get the following: aperrin at perrin:~/afshome/papers/authoritarian/R$ R R : Copyright 2003, The R Development Core Team Version 1.7.0 (2003-04-16) R is free

Sorry if this is obvious, but my S skills aren't great and I haven't been able to find it documented anywhere. I want to write a new function for use with lme objects; the function will simply calculate an ICC (aka "rho") for each level of a mixed-effects model. What I need for this is pretty simple: (c(var1..varn, residual)) / sum(c(var1..varn, residual)) where var1..varn

A while ago (April 2002) there was a short thread on software for generalized linear mixed models. Since that time, has anyone written or found R code to use a glmm to analyze binomial data? I study CWD in white-tailed deer, and I'd like to do a similar analysis as Kleinschmidt et al. (2001, Am. J. Epidemiology 153: 1213-1221) used to assess control for spatial structure in malaria

Greetings- I'm working with data that are multilevel in nature and have a dichotomous outcome variable (presence or absence of an attribute). As far as I can tell from reading archives of the R and S lists, as well as Pinheiro and Bates and Venables and Ripley, - nlme does not have the facility to do what amounts to a mixed-effects logistic regression. - The canonical alternative is

Greetings- An odd situation has developed. I use the following code to create .eps files of two very similar graphs: postscript(file='resources.bygt.eps', onefile=FALSE, horizontal=TRUE) barplot(resources.bygt.matrix, beside = TRUE, legend.text=c('narrative','doubt'),

Returning after the long weekend, I get the following: aperrin at perrin:~/afshome/papers/microcultures/R$ R R : Copyright 2003, The R Development Core Team Version 1.7.0 (2003-04-16) R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type `license()' or `licence()' for distribution details. R is a collaborative

Please forgive my ignorance on grapics. I'm trying to make a relatively simple plot with two line plots, same axes, by mean over a series of dates. I can make the plot well like this: plot(sort(tapply(first.anti.auth.sum,date,mean), partial=1), type="l", col="yellow",ylim=c(0,2.0)) par(new=TRUE) plot(sort(tapply(first.pro.auth.sum,date,mean), partial=1),

Hi, Can R calculate an intraclass correlation coefficient for clustered data, when the outcome variable is dichotomous? By now I calculate it by hand, estimating between- and intracluster variance by one-way ANOVA - however I don't feel very comfortable about this, since the distributional assumptions are not really met.... Maybe anyone can help me? Best regards and many many thanks,

Greetings- I have some data on which I need to do something like a Hierarchical Linear Model (please bear with me, I'm only learning the technique so I don't know yet if my language is correct). Essentially I'm analyzing data at two levels simultaneously; data are about individuals an organizations of which they are members. Can someone point me toward an appropriate package in R?

In the non linear mixed effects package a groupedData object can be created to facilitate modeling. The gD object includes a formula of the form 'response variable' ~ 'primary covariate' | 'grouping factor'. In experiments creating response surfaces there are 2 or more primary covariates. Is there any way to use the groupedData() function to include 2 primary

Sorry for the off-topic question, but I know there are some talented LaTeX users out there. Which bibliography style gives only the year in text citations (e.g "for further details, see Anderson (1992)" )? Thanks Jason -- Indigo Industrial Controls Ltd. http://www.indigoindustrial.co.nz 64-21-343-545 jasont at indigoindustrial.co.nz

Can R conduct content analysis? I search "content analysis" in mail archives and only find the concord package can compute Krippendorff's alpha. If R cannot cope with content analysis, does anyone can tell me some other softwares? Thank you very much!

Is anyone using R with postgres? I'd like to do so, but cannot seem to find any reasonable explanation of how to do so. Joshua Gramlich Chicago, IL

Any ideas on this? nujoma:~> R R : Copyright 2002, The R Development Core Team Version 1.5.0 (2002-04-29) ... Error: couldn't find function ".Alias" Error: couldn't find function "attach" [Previously saved workspace restored] > q() Error: couldn't find function "q" > Suspended nujoma:~> kill %1 This is on a brand-new installation of R

Hi Folks, Could anyone point me to a good reference on linear regression models? Specifically, I am trying to gain an intuitive feel for how the standard error values are calculated for the parameter estimates. My understanding is that these are computed using the variance-covariance matrix computed from the input data matrix. Although I think I understand the math, I still don't have a

I'm asking on this list because the question will be more easily understood than on a Perl list. Lots of talented people using R also know how to use Perl, so it's easier to ask them. Running a Perl script with the -d switch can do some pretty neat things, but as far as I can tell, every line has to be done individually. Is there a way to specify where to stop in the way browser() does

Hello R community, I would like to test for learning effects by subjects in my experiment. Each subject participates in six consecutive auction rounds of the same treatment. The response variable is the efficiency of an auction outcome measured by a real number. Since the efficiency increases over the six rounds, I suppose that subjects learn about the rules of the auction institution, but I

Dear list In HLM, one can specify a covariate at one of the "levels". For example, if the data structure are repeated observations nested within students nested within schools, school size might be a covariate that is used at level 3, but not at the other levels. In HLM this is rather easy to do. However, how can one specify a covariate in R for only one of the levels? I have a

---------- Forwarded message ---------- Hello, I am analyzing a weighted, stratified, clustered survey data set with approximately 1 million observations and 50 variables. I am new to R (I'm a Stata user), and so far couldn't find any documentation on how to handle survey data. In other words, is there a specific package to handle a combination of weigths, clusters and strata. I am also