similar to: Trends for many units

************* Content Filter Notification ************** The following mail was blocked since it contains sensitive content. Source mailbox: <s-luppescu at uchicago.edu> Destination mailbox(es): <r-help at stat.math.ethz.ch> Policy: Dirty Words(shit) Action: Quarantine Recipient, Content filter has detected a sensitive e-mail. ******************* End of message *******************

Is the width= parameter in barplot() supposed to work? I couldn't get it to work in my plot, and even in the example, data(VADeaths, package = "base") barplot(VADeaths, width=rep(0.1, 4)) the plot looks identical regardless of what I put in for the width. I looked at the source for barplot() and it looks like it SHOULD work (but what do I know). Am I doing something wrong?

Is there an efficient way to do linear model selection by choosing the model with the highest BIC from all possible models? ______________________________________________________________________ Stuart Luppescu -=-=- University of Chicago $(B:MJ8$HCRF`H~$NIc(B -=-=- s-luppescu at uchicago.edu http://www.consortium-chicago.org/people/sl.html http://musuko.uchicago.edu/pubkey.asc

I would like to make my horizontal barplot vertical axis labels perpendicular to the axis. I tried las=1, srt=90; I even tried yaxt='n' thinking I'd put the labels in using mtext, but the axis continued to be drawn. Can anyone help me with this? My barplot() statement looks like this: barplot(height=foo$rebint[o1], names=foo$Unit[o1], horiz=TRUE, col=mycolors[foo$type[o1]],

I want to estimate the center of a distribution with lots of outliers in one tail, and thought I would use a function such as S-plus's location.m() with psi.fun=bisquare (as per MASS 3 p. 131). However, R seems not have such a function, so my questions are: 1) Is there an R equivalent to location.m()? 2) Would huber() give me results that are similar (i.e., close enough)? Thanks.

I would like to make the print format of the axis value labels appear as ``nn%'' instead of just the numbers. Is there a way to do this? Also, the x axis value labels on my horizontal barplot are rotated 90 degrees (perpendicular to the axis) no matter what I do. I tried las=1 but that didn't work. Is this doable? Thanks. $platform [1] "i686-pc-linux-gnu" $arch [1]

I seem to remember this coming up before, but I can't find it any messages I've saved or in the archives (searching by subject). I want to rotate mtext so that it's perpendicular to the right side. I tried srt=90 and lots of other values, but it seems to be ignored. Is there a way to do this? ______________________________________________________________________ Stuart Luppescu

I have plain-text data in lower triangular form that I want to read in. Does anyone know of an easy way to do this? ______________________________________________________________________ Stuart Luppescu -=-=- University of Chicago $(B:MJ8$HCRF`H~$NIc(B -=-=- s-luppescu at uchicago.edu http://www.consortium-chicago.org/people/sl.html http://musuko.uchicago.edu/pubkey.asc for PGP

Hello, I'm trying to do a little rlm of some data that looks like this: UNIT COHORT perdo adjodds 1010 96 0.39890 1.06894 1010 97 0.48113 1.57500 1010 98 0.36328 1.21498 1010 99 0.44391 1.38608 It works fine like this: rlm(perdo ~ COHORT, psi=psisquare) But the problem is that I have about 100 UNITs, and I want to do a

Hello, I'm trying to read data in from a file using scan(). The last field is a character string that contains blanks. I had read it in in S-Plus using this code: ifile <- list(entry=0,measure=0,st=0,count=0,score=0, error=0,inmsq=0,instd=0,outms=0,outstd=0,displc=0,ptbis=0,a=0, r="",name="") if.widths <- c(1, 5, 8, 3, 6, 7, 7, 7, 7, 7, 7, 7, 7, 2,

I'm stumped with this. When I execute the lines in the function singly, they run fine, but when I run the function, I get this error on the read.table() line: Error in count.fields(file, sep, quote, skip) : can't open file fspci1.dat Can anyone tell my why this should be so? Here is the program: library(rpart) wait <- function(str="Press a key when ready...")

I want to leave a little extra room at the bottom for a legend in a box plot, but I find that the boxplot function ignores ylim. Is there any way around this? Thanks. ______________________________________________________________________ Stuart Luppescu -=-=- University of Chicago $(B:MJ8$HCRF`H~$NIc(B -=-=- s-luppescu at uchicago.edu

Has the loess.smooth() function changed? It used to work, but now it causes R to abort with a segmentation fault. I stole the function points.lines() from V&R 1st ed. pp. 67--68, but now it only works if I remove the line with loess.smooth. Here's the function I'm using: points.lines <- function(x, y, ...) { cor1 <-round(cor(x, y, use="pairwise"), digits=2)

I need to calculate an index of predictive association, similar to Goodman and Kuskal's (1954) for a 2x2x2 array. I have two questions: 1) Is there an easy way to do this in R? (Of course, the answer to this kind of question is nearly always ``yes.'' Perhaps I should have asked, ``How can I do this in R?'') 2) Is it reasonable to try this with a 3-dimensional array? The way I

Hello, As you can see from my signature in this message, I use the R fortune function to generate a fortune, which is then fed to the signature program, which constructs a named pipe containing the fortune-bearing sig, which is then included in mail messages. The problem is that it's got extraneous junk in it and I can't figure out how to get rid of it. This is the command that generates

This is really a question of how to program this *BETTER*. It works as I have done it, but is quite ugly. I want to do a 3d scatterplot of the upper triangle of a matrix, where the z-values are the values in the matrix, and the row and column indices are the y- and x-values. The complete (11 by 11) matrix is mmtop94.2. Here is my awkward code: mmtop94.2[lower.tri(mmtop94.2)] <- NA # Here i

I'm trying to install packages in Gentoo Linux using install.packages(). The packages download fine, but when it comes to installation, I get this message: WARNING: ignoring environment value of R_HOME /var/tmp/portage/R-1.5.1/image//usr/lib/R/bin/Rcmd: /var/tmp/portage/R-1.5.1/image//usr/lib/R/bin/Rcmd: No such file or directory I have the R_HOME set to the correct value (/usr/lib/R), but I

I can't remember the programming term -- ``implicit declaration''? Has this changed between 1.4.1 and 1.5? This code used to run as is: for (i in 1:length(diff.time)) { if (diff.time[i] < 16) { total.time[i] <- ifelse(i==1, 7.5, total.time[i-1]) + diff.time[i] last.time[i] <- 0 } else { total.time[i] <- 7.5 last.time[i-1] <- 1 } } Now in 1.5 I

Do MD5 digests exist somewhere for the source packages? If not, is there some other way to verify the integrity of the files we've downloaded? Thanks. -- Stuart Luppescu -=- s-luppescu at uchicago.edu University of Chicago -=- CCSR $B:MJ8$HCRF`H~$NIc(B -=- Kernel 2.4.19-xfs-r1 There *is* no such thing as a civil engineer. -------------- next part

As of a moment ago (9 a.m., Monday, Jan. 8, 2001), the S Programming book by Venables and Ripley was once again selling for 17.98 American dollars at Amazon.com. Alas, I notice that "Mixed Effects Models in S and S-Plus" by Pinheiro and Bates now lists at Amazon for 97.25 American dollars, as opposed to about $70 a little while ago.