similar to: Verifying lm results

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- 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'd like to get summary statistics (really just a mean would be fine) for a vector in a data frame, but split based on the value of another vector. That is, I have a data frame (hcd.df) with variables datecat (which is always 1 or 2) and auth.sum (-8..+8). I've used xtabs to get chi-square comparisons, but what I need now is a simple mean of auth.sum where datecat is 1 and another where

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?

Greetings- I am in the process of constructing a nonlinear model using nlme. The model is attempting to fit a nested data structure from some public-opinion data (the data are from individuals nested within organizations). The question I have is fairly simple (I hope). The data were collected in two stages: a 15,000-subject randomly-sampled telephone interview (the SCREEN), with 2,517 subjects

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.

Greetings. I've been experimenting with some algorithms for document classification (specifically, a Naive Bayes classifier and a kNN classifier) and I would now like to calculate some inter-rater reliability scores. I have the data in a PostgreSQL database, such that for each document, each measure (there are 9) has three variables: ap_(measure), nb_(measure), and knn_(measure). ap is me

My apologies if this is obvious - I assume it is but still can't seem to make it work. I've got a data frame acps.df; I want to run an analysis on a subset of that data frame, defined as the subset for which the value of OGMCAT (a column in the data frame) is anything OTHER than 9, 11, or 12. So I'd like to define a new data frame that consists of the same variables as in acps.df,

Greetings- I have a strange problem displaying a graph from R (1.3.1, linux) in a LaTeX document of documentclass seminar. I'm using graphicx to include the file: \usepackage{graphicx} ... \resizebox{\textwidth}{\textheight}{\includegraphics{crime.eps}} When I do this, the entire slide (including the page number) is rotated 180 degrees. Any ideas why this happens? The graph was created

Greetings. I'd like to generate a graph that displays two distributions as side-by-side bar graphs, then plots a smoothed line of the distribution on top. The idea is to be able to visually compare the distributions. I've done the following, but the axes don't line up. Any suggestions? barplot(table(hcd.df$datecat, hcd.df$auth.sum), beside=TRUE,

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

Earlier today I posted a problem importing an R graph into a LaTeX file of seminar class: specifically, the graphic was showing up rotated 180 degrees, along with the rest of the page it was on. In a real victory for open-source software, I got lots of responses with three distinct approaches, each of which appears to solve the problem. Try getting fast, correct help from Microsoft on a Sunday

I'd like to produce a series of simple line graphs for my methods class that show the three questions used on a repeated survey to make up a particular index. The data frame is: > efficacy.df year complicated havesay dontcare 1 1952 71 68 63 2 1954 NA NA NA 3 1956 64 71 71 4 1958 NA NA NA 5 1960

Greetings. I'm trying to replicate an analysis I did a few years ago, then in SPSS, using the SPSS GLM command: GLM n_diffpt WITH age_i inc_i join_i work_i educ_i give_i cs_i eff_i age_a inc_a join_a work_a educ_a give_a cs_a eff_a /METHOD = SSTYPE(3) /INTERCEPT = INCLUDE /PRINT = PARAMETER ETASQ RSSCP GEF /CRITERIA = ALPHA(.05) /DESIGN = age_i*age_a inc_i*inc_a

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

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

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