similar to: saving a data.frame to "\t" files

Greetings, I would like to plot three lines on the same figure, and I am lost. There is an answer to a similar thread… but I tried matplot and it is beyond me. An example of the data follows: Year EM IM BM 1983 9.1 16.8 -7.7 1984 12.0 18.0 -6.0 1985 13.6 19.1 -5.5 1986 12.4 17.3 -4.9 1987 14.6 20.3 -5.7 1988 20.6 23.3 -2.6 1989 25.0 27.2 -2.2 1990 28.4 30.2 -1.8 1991 33.3 31.2 2.1 1992 40.6

Greetings: I am trying to work with spdep (everything is "brand new" downloaded this morning). OS = Windows 2000 (also up to date). The code I am using follows: #example gal.county=read.geoda("lnpilnd.GAL", row.names=NULL, skip=0) summary.nb(gal.county) Error in summary.nb(gal.county) : Not a neighbours list #end The gal file works just fine in GeoDa (also up to date).

Greetings: I would like to make a table with descriptive statistics for a data.frame. I guess the question is how can I put together, in a table, the results from, say: apply(df, 2, mean, na.rm =T) apply(df, 2, median, na.rm =T) ....... Thanks, Mihai Nica Jackson State University 155 B Parkhurst Dr. Jackson, MS 39202 601 969 5423 601 914 0361 [[alternative HTML version deleted]]

Greetings: I am unsuccessful in suppressing "Given : myvariable" from a coplot. There was such a question in the past but the thread breaks down. I am sure this is a "for dummies" question :-). I tried: coplot(myvar~myvar | myvar, show.given=FALSE, xlab="....", ylab="...", main=" ") and some other variations (including without main=" ")

Greetings: I am trying to find a package/solution for panel (longitudinal) data analysis. Unfortunately it seems I don't know where to start. Could somebody offer a hint? Thanks, Mihai Nica Jackson State University 155 B Parkhurst Dr. Jackson, MS 39202 "No good deed will ever remain unpunished" [[alternative HTML version deleted]]

I just read paul murrell's new book, R graphics. now, I have always used the traditional graphics system. apparently, the new (trellis?) system is an entirely separate graphics system. after reading the book, I cannot figure out what the intrinsic capability advantage of the old graphics system is that cannot be replicated in the trellis system. if the new system's capabilities are

gReetings: Is there any way to download a (or a sample of a) crossection of stock market prices? Or is it possible to use get.hist.quote with a *wild card*? Thanks, mihai Mihai Nica 170 East Griffith St. G5 Jackson, MS 39201 601-914-0361 ____________________________________________________________________________________ 8:00? 8:25? 8:40? Find a flick in no time [[alternative HTML

Greetings: I would like to kindly ask for a little help. The rough code is: #________________________________________________________ dat=data.frame(read.delim(file="all.txt", header = TRUE, sep = "\t", quote="\"", dec=".",na.strings = "NA")) nc=read.shape("astae.shp", dbf.data=TRUE, verbose=TRUE) mappolys=Map2poly(nc) submap

Hi I am fairly new to R and I would appreciate some help to hopefully a trivial problem. I created a function: summary.aggregate <- function(y, ...) { temp.mean <- aggregate(y, FUN=mean, ...) temp.sd <- aggregate(y, FUN=sd, ...) temp.length <- aggregate(y, FUN=length, ...) temp <- data.frame(cbind(mean=temp.mean$x,stdev=temp.sd$x,n=temp.length$x)) } this outputs e.g.:

Dear r-users! I would like to cross two data frame which have the same row number but different in the number of column. Can anybody help me for this case ? Thanks a lot in advance -------------------------------------------------------------------------------- Majid Iravani PhD Student Swiss Federal Research Institute WSL Research Group of Vegetation Ecology Z?rcherstrasse 111

gReetings: I apologize for the non-programming language. I found what seems to be a strange behavior of plot(). The code follows: #_____________________________ N=3030; gn=.04; tn=1:100 n=N/(1+(N-1)*exp(-gn*tn)) N=n*(1-exp(-gn*tn))/(1-n*exp(-gn*tn)) plot(N) #strange plot N N=rep(3030, 100) plot(N) #behaves as expected ________________________________ Windows XP, R 2.4.0. or 2.4.1 Thanks,

Greetings: Using 'identify' to label points on a plot works just fine. However, when saving under 'metafile' or using the clipboard the labels dissapear. I believe it's an SDI issue. I am running last R with last Tinn-r under XP up to date. Anything I can do besides going back to MDI :-)? Thanks, Mihai Nica, ABD Jackson State University ITT Tech Instructor 170 East Griffith

Dear r-helpers! I am kind of new to R. I would like to calculate the mean of the numbers of this expression: data(USArrests) USArrests[row.names(M) == "Alabama",] class() tells me it's a ``data.frame,'' what I actually desire is to get all numbers of a row as a vector or a list to let mean() calculate the mean of the whole row. (I know this doesn't make sense

I might be sorry for asking this question :-) I have two equations and I tried to estimate them individually with "lm" and "gls", and then in a system (using systemfit) with "OLS", "WLS" and "SUR". Quite surprisingly (for myself at least) the results are identical to the last digit. Could someone (please!) give a hint as to what am I

R 1.8.1, Windows 2000. I am trying to find the "legend" for color coding and point types in a plot, which probably are standard for everybody but myself. Thanks for the tip! Mihai Nica Jackson State University "No good deed will ever remain unpunished" [[alternative HTML version deleted]]

I can't understand what is happening. This is the code and results: > agoa <- read.table(file = "C:/Users/HTPC/Documents/_Documents/Research/WithDidia/AGOAUSImports.txt", header = T, sep = "\t", dec = ".", na.strings = "NA", stringsAsFactors = T)# > str(agoa); names(agoa) 'data.frame':109 obs. of  19 variables:  $ X    : Factor w/ 39

R 1.9.1 on Win2000 or Win98SE. I am using coplot as follows: coplot(AVG~LRPI| REGION) the output seems normal but I get: "Warning message: calling par(new=) with no plot" This is the only explanation that I have for being unable to use par() with coplot for changing the way the xlab and ylab appears. From within coplot I can change the text itself but not the font, fontsize, etc,

Greetings: I would like to kindly ask help with obtaining mexval statistics (marginal explanatory value - percentage increase in SEE if the variable were left out of the regression model) for a logit (glm) model with several continuous independent variables. I believe I can do it manually for each variable, but I really hope there might be somebody who has a function already written. Writing one

Hi, I'm unable to get Kerberos4 authentication working with openssh-3.4p1. I'm getting a message that privsep is not available on my platform (Irix 6.5.15) and another message stating that compression and privsep are mutually exclusive. But, ssh decided to turn off compression, I think because of servconf.c. I think it would be more usefull to have compression enabled and disable privsep

Dears Sirs During my computational work I encountered unexpected behavior when calling "ar" function, namely # time series x<-ts(c(-0.2052083,-0.3764986,-0.3762448,0.3740089,0.2737568,2.8235722,- 1.7783313,0.2728676,-0.3273164),start=c(1978,3),frequency=4,end=c(1980,3)) # ar function res.ar<-ar(x,aic=TRUE,demean=F) # call "ar" again and ............