similar to: methods on missing data

Thanks for replies, but i was not precise enough. The problem is not evaluating statistics on data with NA values. The problem is evaluation of statistics on data with length = 0. To make the problem more clear this is what i tried: This works fine: tapply(as.numeric(c(NA,2)), as.factor(c("a","b")), summary) But i need SDev, aswell, so i copied summary.default to

Dear help list, I have a timestamp in as a chron object: > (x <- chron(dates = c("12/02/11", "22/11/11"), + times = c("07:30:00", "04:00:00"), + format = c(dates = "d/m/y", times = "h:m:s"))) [1] (12/02/11 07:30:00) (22/11/11 04:00:00) Now I want to shift the timestamp by 06:30 (hh:mm) backwards, to get

Dear all, I have some problems with importing data from an Access data base via RODBC to R. The data base contains several tables, which all are imported consecutively. One table has a column with column name "NO". If I run the code attached on the bottom of the mail I get no complain, but the column name (name of the respective vector of the data.frame) is "Expr1014" instead

My first Problem is a minor one. Maybe some knows a reson. Pleas give me a hint off-line, as it seems to be something stupid. 1. I sent my first mail to r-help at lists.R-project.org and got Delivery Failure Report Your document: make error R-1.8.0 on SuSE Linux 7.3 (i386) was not delivered to: r-help at lists.R-project.org because: Delivery time expired I already had this problem

----- Original Message ----- From: "Richard A. O'Keefe" <ok at cs.otago.ac.nz> To: <paul at datavore.com> Sent: Thursday, September 04, 2003 2:56 AM Subject: Re: [R] Overlaying graphs > I do not know how to overlay the curve graphic on top of hist graphic. > > Do you know about the "add=TRUE" option for plot()? > > I am hoping to show visually

Hi I have recently been analyzing birthweight data from a clinic. The data has obvious defects in that there is digit preference on certain weights making them overrepresented. This shows as spikes in the histogram on certain well rounded weights like 2, 2.5, 3, etc. I would like to show this to government officials but can't figure out how I should present the finding in an easy to

I use knncat to make a predictive model and get misclass rate > knncat.m<-knncat(training.new,k=c(10,20),classcol=5) > knncat.m Training set misclass rate: 36.88% then I try to calculate prediction accuracy by the following: > pr.knncat.train <- predict (knncat.m,training.new,training.new,train.classcol=5,newdata.classcol=5) > tb.knncat.train <-table (pr.knncat.train,

I'm learning how to use tapply. Now I'm having a go at the following code in which dati contains almost 600 lines, Pot - numeric - are the capacities of power plants and SGruppo - text - the corresponding six technologies ("CCC", "CIC","TGC", "CSC","CPC", "TE"). .....................................................

Dear all, I have done further modifications on the histogram functions that I reported earlier this year, and I hope this can be of use and perhaps included in the distribution. I have been using this stuff a couple of months myself, and while it is nothing sophisticated, it has it's applications. :-) I did a few small modifications today to make it a bit more compact. I have modified the

Hi, I've been passing the hist() function to tapply() to quickly generate histograms based on the list of factors supplied to tapply(). However, I have not figured out how to generate titles for each of the histograms, which paste in the unique values of the list factors as part of the histogram title. I'm hoping that someone can tell me how to do this. Thanks for your time and

Dear members, I would like to pass the histogram settings to each subset of the dataframe, and generate a multiple figures graph. First, can anyone tell me how to generate a multiple figures environment? I am trying mfrow=c(2,4) and nothing appears. Secondly, I want to pass the following function in tapply() hist(x, freq=FALSE) lines(density(x), col="red") rug(x) how can I manage

This is a subject that came up in the virtio BOF session at OLS. I decided to go forward and implement something that I like, based on the latest virtio proposal at the time, which was draft III. It's not a drop-in replacement, because it's missing a host implementation. I first started my own, which is not done yet, but wanted to do one for lguest and one for emulated PCI next. It's

This is a subject that came up in the virtio BOF session at OLS. I decided to go forward and implement something that I like, based on the latest virtio proposal at the time, which was draft III. It's not a drop-in replacement, because it's missing a host implementation. I first started my own, which is not done yet, but wanted to do one for lguest and one for emulated PCI next. It's

This series adds an initial implementation for scrub. It works quite straightforward. The usermode issues an ioctl for each device in the fs. For each device, it enumerates the allocated device chunks. For each chunk, the contained extents are enumerated and the data checksums fetched. The extents are read sequentially and the checksums verified. If an error occurs (checksum or EIO), a good copy

G'Day, I want to access in a function called from tapply the current factor. In my example below, all I want to do is to write the current factor on each histogram. Needless to say my example does not work. I would be grateful for pointers in the right direction. Many thanks Bernie McConnell Sea Mammal Reserach Unit cc <- 1:10 ff <- rep(c("a","b"),5) pp<-

This patch set is a follow up to "Add a vhost RPMsg API" [1], it is marked as an RFC because firstly it depends on the RPMsg API series and secondly it is currently being reviewed on ALSA and SOF mailing lists, but any early comments from virtualisation developers would be highly appreciated too! Thanks Guennadi [1]

This patch set is a follow up to "Add a vhost RPMsg API" [1], it is marked as an RFC because firstly it depends on the RPMsg API series and secondly it is currently being reviewed on ALSA and SOF mailing lists, but any early comments from virtualisation developers would be highly appreciated too! Thanks Guennadi [1]

Dear all, I have a question about the 'sdev' value returned by the princomp function (which does principal components analysis). On the help page for princomp it says 'sdev' is 'the standard deviations of the principal components'. However, when I calculate the principal components for the USArrests data set, I don't find this to be the case: Here is how I

Dear R-users & developers, i'd like to know your opinion on the following suggestion. In short: Wouldn't it be possible to combine the properties of local regression (loess) and resistant regression (lqs) in order to cope with fitting of partially continuous data? I think the lqs-fitting could detect and eliminate the 2 outliers but perserve the discontinuity at t1 in case of the

Dear all I have a dataset that looks like: fr.wt site 1 4400 glen 2 235 glen 3 225 glen ' ' ' ' ' ' ' ' ' 82 550 glen 83 550 kom 84 550 kom ' ' ' ' ' ' ' ' ' 191 820 kom 192 2000 soet ' ' ' ' ' ' I need to do a series of histograms for each of the codes, levels or factors in