search for: cliste

Dear R help listers, I am trying to replicate results in Gelman and Hill's book (Chapter 3 in regressions and multilevel models). Below I estimated two models (chp3.1 and chp3.3 in R codes) with the same data and dependent variable but different independent variables. I have been using Stata for quite a while, and I know I can use foreach to build a loop to condense the codes (especially if I

Hi, Is there any tutorial/implementation cloud on CentOS? Or anyone with experience like to share? Regards and thanks wL

Hi All, I'm still pretty new to using R - and I was hoping I might be able to get some advice as to how to use 'apply' or a similar function instead of using nested for loops. Right now I have a script which uses nested for loops similar to this: i <- 1 for(a in Alpha) { for (b in Beta) { for (c in Gamma) { for (d in Delta) { for (e in Epsilon) { Output[i] <-

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

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

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'),

Hi, I am R beginner and am trying to figure out how to generate a complete list of species for every point, visit, and year. The code below is close but does not give me a list of species for every point, visit, and year in my data set. spplist<-unique(sumPtCt$Species) spplength<-length(spplist) Pointlist<-unique(sumPtCt$Point) Pointlength<-length(Pointlist)

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),

R users in the southeastern US (and those on the Internet-2 backbone) will be happy to hear of the new CRAN mirror at ibiblio.org aka metalab.unc.edu. ---------------------------------------------------------------------- Andrew J Perrin - http://www.unc.edu/~aperrin Assistant Professor of Sociology, U of North Carolina, Chapel Hill clists at perrin.socsci.unc.edu * andrew_perrin (at) unc.edu

Hi, I am trying to install egroupware. The installation tests have some warnings. What are the packages missing which I can use yum to install them? Checking PEAR::HTTP_WebDAV_Server is installed: False PEAR::HTTP_WebDAV_Server is needed by: WebDAV. You can install it by running: pear install HTTP_WebDAV_Server Checking PEAR::Auth_SASL is installed: False PEAR::Auth_SASL is needed by:

Hi, I have been trying to install CentOS 5.4 on a Intel SR1530SHS, Intel S3200SH mainboard.. It has a 3 x 1TB sata hotswap drives with LSI software raid onboard. I had configured the LSI to have Sata0 and Sata1 with raid 1 and the third drive as a hotspare drive. Format the harddisk and installation was a breeze. The server rebooted into a blank screen and the cursor just keep blinking. Please

Greetings- I need to extract, programatically, the standard deviations of the intercept and residuals from an lme model. These are presented by print.lme as: ... (Intercept) Residual StdDev: 1.410635 0.7800512 ... (data taken from ?lme's examples section) I can get the residuals with x$sigma where x is the fitted lme object. I can't find the intercept, though. The closest

o-c-config checks to see if local storage is configured prior to starting. If local storage is not configured then it reports an error and exits. Signed-off-by: Darryl L. Pierce <dpierce at redhat.com> --- scripts/ovirt-config-logging | 6 +++++- 1 files changed, 5 insertions(+), 1 deletions(-) diff --git a/scripts/ovirt-config-logging b/scripts/ovirt-config-logging index

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

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

I have the requirement to empty the association records and add the association records as new (as below). @task = Task.find(params[:task_id]) @task.task_users = [] params[:ids].split(", ").each do |court_user_id| task_user = TaskUser.new(:court_user_id => court_user_id) @task.task_users << task_user end @task.save This is not working. Following error is thrown.

Hi, I'm trying to match peaks between chromatographic runs. I'm able to match peaks when they are chromatographed with the same method, but not when there are different methods are used and spectra comes in to play. While searching I found the ALS package which should be usefull for my application, but I couldn't figure it out. I made some dummy chroms with R, which mimic my actual

Dear People, I have been trying to perform concatenation operations on lists without success. Consider the following example. > x <- runif(2) > y <- runif(2) > z <- runif(2) > xy <- list(x,y) > xyz <- list(x,y,z) I want to have a function foo such that foo(xy,z) is the same as xyz. This appears not to be as easy as I thought, though there must be a way. For example

Got this today. Please take a look and let me know what you think. - Dan ---------- Forwarded message ---------- From: Kevin Burge <kevin.burge at systemware.com> Date: Dec 7, 2006 6:45 AM Subject: win32-service problems with patch To: djberg96 at gmail.com Hi Daniel, Thanks for win32-service. In the process of using it for running a reliable-msg (drb) server and clients, I had some

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