similar to: [ibiblio.org #1674] Mirroring request (fwd)

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

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

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

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

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

Greetings- Using R 1.5.0 under linux and the latest MASS and nlme, I am trying to develop a three-level (two levels of nesting) model with a dichotomous oucome variable. The unconditional model is thus: > doubt1.pql<-glmmPQL(fixed = r.info.doubt ~ 1, random = ~1 | groupid/participantid, + family = binomial, data = fgdata.10statements.df) iteration 1 iteration 2 iteration 3 iteration 4

Can anyone shed any light on this? > doubt.demographic.pql<-glmmPQL(random = ~ 1 | groupid/participantid, + fixed = r.info.doubt ~ + realage + minority + female + education + income + scenario, + data = fgdata.df[coded.resource,], + na.action=na.omit, +

Recently I've been having this happen frequently: aperrin at perrin:~/microcultures$ R R : Copyright 2002, The R Development Core Team Version 1.5.0 (2002-04-29) 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 project with many

I'm grappling with a problem and would appreciate any thoughts on it. I'm revising a paper for resubmission to a journal. For the paper, I've coded each "turn" in a series of conversations with several binary codes. (A turn is one package of statements made by one speaker, starting with the beginning of the speech and ending when the speaker stops or is interrupted.) The

Greetings- a reviewer for a paper of mine noted an anomaly in some models I ran using glmmPQL (from the MASS package). Specifically, the models are three-level hierarchical probit models estimated using PQL under R. The anomaly is that the log-likelihoods decrease (or, alternatively -2logLik increases) as variables are added to the null model. This is unusual, and I'm trying to figure out

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. I'm running some models under R using glmmPQL from MASS. These are three-level models (two grouped levels and the individual level) with dichotomous outcomes. There are several statistics of interest; for the moment, I have two specific questions: 1.) This question refers to the following model (I present first the call, then the output of summary():

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

In the R-2.3.0 directory, what do you get if you type bin/R? Does R starts up OK, or do you get an error? What were the last few lines of output from running configure? Andy _____ From: Pramod Anugu [mailto:pramod.r.anugu@jsums.edu] Sent: Wednesday, May 31, 2006 11:23 AM To: Liaw, Andy; 'Andrew Perrin' Cc: r-help@stat.math.ethz.ch Subject: RE: [R] Running R [Broadcast]

Hi! Meanwhile also the location is fixed: it is happening at metalab (http://metalab.at/) - a place for geeks. Thus, we meet there at Thursday, 20.9.2007, 19:00 CEST (=local Vienna time). Metalab is located next to the city hall: http://metalab.at/wiki/Lage Metalab is no pub/restaurant. Thus, don't come hungry! Nevertheless liquid food (drinks) is available. We meet in the library (in

Is this really caused by a buggy server, or is this an interoperability problem? It seems to work ok when I specify -o "protocol 1" on the command line. Thanks, Greg [gleblanc at grego1 gleblanc]$ ssh -v login.metalab.unc.edu SSH Version OpenSSH_2.2.0p1, protocol versions 1.5/2.0. Compiled with SSL (0x0090581f). debug: Reading configuration data /etc/ssh/ssh_config debug: Applying

Greetings- Running R 2.4.0 under Debian Linux, I am getting a memory error trying to read a very large file: > library(foreign) > oldgrades.df <- read.spss('Individual grades with AI (Nov 7 2006).sav',to.data.frame=TRUE) Error: cannot allocate vector of size 10826 Kb This file is, granted, quite large: aperrin at perrin:/data0/grading$ ls -l total 630304 -r-xr-xr-x 1 aperrin

I've installed the openss* rpm's from metalab.unc.edu/pub/Linux/distributions/redhat/contrib/libc6/i386 on a redhat 6.2 system. sshd is running but refuses all connections from all hosts including localhost. The client reports debug: Connection established. ssh_exchange_identification: Connection closed by remote host debug: Calling cleanup 0x8056160(0x0)

Hi, I am trying to install qvalue, however its giving installation error -> Error : package 'tcltk' does not have a name space ERROR: lazy loading failed for package ?qvalue? * removing ?/home/sbw/R/x86_64-unknown-linux-gnu-library/2.12/qvalue? The downloaded packages are in ?/tmp/RtmpKnS1X4/downloaded_packages? Warning message: In install.packages(pkgs = pkgs, repos = repos, ...)