similar to: Problems!

I have been trying to read in a large data set using read.table, but I've only been able to grab the first 50,871 rows of the total 122,269 rows. > f <- read.table("http://www.cs.odu.edu/~fmccown/R/Tchange_rates_crawled.dat", header=TRUE, nrows=123000, comment.char="", sep="\t") > length(f$change_rate) [1] 50871 From searching the email archives,

I was wondering if there was any way to completely remove the borders around each bar when using the barplot function. I understand how to change the color of the border (border="white"), but I need the border to be invisible so I can plot many skinny bars right next to each other with only the bar's color visible. Thanks, Frank -- Frank McCown Old Dominion University

Hi, I am using R to fit statistical models to data were the observations are means of the original data. R is used to calculate the mean before fitting the model. My problem is: When R calculates the means using tapply, the class of the means differs from the class of the original data, which gives me trouble when I want to use the original data to calculate model predictions. Here is a simple

Dear R list members, I'm encountering a problem with expand.model.frame(): Suppose that I define the following simple function (meant just to illustrate the problem): > fun <- function(model){ + expand.model.frame(model, all.vars(formula(model))) + } > and I have the following model, created with an explicit data argument: > mod Call:

I am using R 2.3.0, Is there a way to get prediction intervals using predict.gls? Any ideas would be appreciated. thanks, Spencer On 11/2/06, Frank McCown <fmccown@cs.odu.edu> wrote: > > I was wondering if anyone knows who should be contacted to add to the R > user-contributed documentation at > > http://www.r-project.org/other-docs.html > > There doesn't

I'm looking at a motherboard with the NF410 (nforce 410) chipset (AUD, VID, LAN on motherboard). Is this chipset readily know by centos 4.3? At some point nvidia changed PCI ID's and I am wonding if the NF410 PCI ID's are known by centos 4.3 before i buy? THanks, Jerry

I may be reading this wrong, but I'm trying to figure out why samba is sometimes trying to do a Unix password lookup to authenticate a samba connect. I have my samba passwords different than my Unix passwords. I'm using my own compiled version of Samba 1.9.18p10. I have an NT4 SP3 box and when I try to access a share through explorer on either of my two Linux boxes with encrypted

In preparing data sets for the lme library, I bring the data into R, convert it to the groupedData class, which is an extension of the data.frame class, and use "dump" to create an R-readable version. To cut down on the size of these files I convert the row.names attribute from expressions like c("1", "2", ..., "54") to 1:54. I have run into a curious

Can someone explain me this spec.pgram effect? Code: period.6<-c(0,0,0,0,0,10,0,0,0,0,0,10,0,0,0,0,0,10,0,0,0,0,0,10,0,0,0,0,0,10 ,0,0,0,0,0,10,0,0,0,0,0,10,0,0,0,0,0,10,0,0,0,0,0,10,0,0,0,0,0,10) period.5<-c(0,0,0,0,0,10,0,0,0,0,10,0,0,0,0,0,0,10,0,0,0,0,10,0,0,0,0,0,0,10 ,0,0,0,0,10,0,0,0,0,0,0,10,0,0,0,0,10,0,0,0,0,0,0,10,0,0,0,0,10,0) par(mfrow=c(2,1))

Send CentOS-announce mailing list submissions to centos-announce at centos.org To subscribe or unsubscribe via the World Wide Web, visit http://lists.centos.org/mailman/listinfo/centos-announce or, via email, send a message with subject or body 'help' to centos-announce-request at centos.org You can reach the person managing the list at centos-announce-owner at centos.org When

Dear R-users, I have a dataset with categories and numbers. I would like to compute and add cumulative numbers to the dataset. I do not understand the structure of by(...) or tapply(...) output enough to handle it. Here a small example -------------- d<-expand.grid(a=1:5,b=1:3,c=1:2) d$n = 10 * d$a + d$b +0.1* d$c Sn<-by(d$n,list(d$a,d$c),cumsum) str(Sn) --------- List of 10 $ : num

Hi everybody! I''m new to this wondeful world that Ruby On Rails is. I''m trying to pass multiple parameters using link_to, but i don''t even know if it''s possible. Here''s the code in where the link_to is used : -- Posted via http://www.ruby-forum.com/. -- You received this message because you are subscribed to the Google Groups "Ruby on Rails:

Suppose I have a small dataframe > aa Target Eaten ID 50 TPP 0 1 51 TPP 1 2 52 TPP 3 3 53 TPP 1 4 54 TPP 2 5 50.1 GPA 9 1 51.1 GPA 11 2 52.1 GPA 8 3 53.1 GPA 8 4 54.1 GPA 10 5 And I want to reshape it into ID TPP GPA 1 1 0 9 2 2 1 11 3 3 3 8 4 4 1 8 5 5 2 10 I realise that

Hi all, I have this figure: http://docs.google.com/Doc?id=df5zfsj4_103rjt2v4d5 created with the following steps: > x [1] 90.4 57.8 77.0 103.7 55.4 217.5 68.1 85.3 152.0 113.0 97.1 89.9 [13] 68.1 83.7 77.4 34.5 104.9 170.3 88.6 88.1 108.8 77.4 85.6 82.7 [25] 81.3 108.0 49.5 71.0 85.7 99.3 203.5 275.9 51.1 84.8 16.5 72.6 [37] 160.5 158.3 136.7 140.0 98.4 116.1

Dear all, can somebody please help me how to calculate "McFadden R-square" for a LOGIT model? Corresponding definition can be found here: http://publib.boulder.ibm.com/infocenter/spssstat/v20r0m0/index.jsp?topic=%2Fcom.ibm.spss.statistics.help%2Falg_plum_statistics_rsq_mcfadden.htm Here is my data: Data <- structure(c(1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 1, 1, 1, 0, 0, 1, 1,

Hi, When I run this code a part of my x-axis and y-axis changes color. Can somebody tell me what is wrong? Also, is there a way to control the color of the average lines? monthplot(AirPassengers+500, ylim=c(min(AirPassengers), max(AirPassengers+500)), ylab="") par(new=T) monthplot(AirPassengers, col="blue", ylim=c(min(AirPassengers), max(AirPassengers+500)),

Dear colleagues, I am relatively new to R and time series and so I am experiencing difficulties in interpreting the output of "arima" in MA models (but not in AR models). I cannot make sense of the 1st innovations returned by "arima". In an AR(1) model I expect data[t]=phi1*data[t-1]+a[t] and in a MA(1) model data[t]=a[t]+theta1*a[t-1]. My interpretation from R-help is

Dear colleagues, ?arima? returns directly the 1-step ahead errors but I am interested in obtaining other h-step ahead errors for several ARIMA models I have fitted. Is there any way I can obtain this with R? Any help would be appreciated. Sincerely, Nuno Prista _________________________ CO - FCUL, Lisboa, Portugal CQFE - ODU, Norfolk, USA

Dear colleagues, Can anyone of you tell me how to write a "nabla" symbol in an R graph? Thanks in advance, Nuno ______________________________________________ Centro de Oceanografia - IO-FCUL, Portugal Center for Quantitative Fisheries Ecology - ODU, USA [[alternative HTML version deleted]]

  Hi r-users,   I would like to use gev and my data (annual rainfall ) is as follows:   > head(dat,20) A B C D E F G H I J 1 45.1 41.5 58.5 50.1 46.0 49.1 37.7 49.1 59.8 54.0 2 50.3 39.8 49.4 56.4 49.4 48.8 42.1 49.8 49.4 58.3 3 41.7 39.3 44.6 39.1 35.7 41.5 40.8 40.8 38.5 45.6 4 50.7 33.9 48.4 28.2 35.5 39.1 61.4 17.0 30.7 38.3 5 39.3 30.6 46.9 23.8 25.8