similar to: Lots of data problems

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

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

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

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

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

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

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

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

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

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():

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

Greetings - is there any update on a PostgreSQL driver for the DBI package? If not, what's the currently-preferred method of creating a link from a PostgreSQL database and R? Thanks. ---------------------------------------------------------------------- Andrew J Perrin - andrew_perrin (at) unc.edu - http://www.unc.edu/~aperrin Assistant Professor of Sociology; Book Review Editor, _Social

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]

Greetings- I'm venturing into a new (for me) area of using R: some presentation graphics based on R data. I'd like to compare two distributions of ordinal variables in relatively small (N=203 and 207) subsets of an R dataset. I can get a reasonably good picture of the distributions' differences with: x11() plot(table(hcd.df$auth.sum[hcd.df$datecat==1]),

Greetings- A relatively small request, but here goes. I have some crosstabs, generated with table(), in R: > table(pg.df$grouptype,pg.df$age) 0 1 2 3 4 B 1 3 6 8 3 C 1 3 4 10 7 P 0 2 9 9 8 S 2 9 4 5 16 U 0 15 15 3 0 X 0 0 0 0 0 I would like to export these in latex format using xtable(). But it appears there's no appropriate method for

Greetings- I'm working through Pinheiro and Bates' _Mixed Effects Models in S and S-Plus_ using R (1.3.1 for linux). On page 13 (okay, so I haven't got that far :)) is: plot.design( ergoStool) which returns on my system: > plot.design(ergoStool) Error: couldn't find function "plot.design" any ideas? Thanks.