similar to: lme How to validate a model with a validation set and a t est set

Dear all, I am looking at habituation of dogs trotting on a treadmill. I record the ground reaction force and I analyze it with several discrete variables (maximum, minimum,...) For each variable, I get between 40 and 50 data per sample. I record data at time 1 min, 2 min, and 4 min a day, and I have 4 days of measurement (one day a week). That means I have 12 samples : Day1_Min1, Day1_Min2,

Dear all, I am looking at habituation of dogs trotting on a treadmill. Each record is made of 40 to 50 data for the same variable (for example Peak). I get each record at several minutes (1, 2 and 4) for each session. And I have 4 sessions of training (one session a week). The aim is to study the effect of the factor (Minute) on the Peak variable, and to study the changes of this effect

I'm trying to find a clever way to re-map data from a database query into a data.frame. Querying a database often returns a table (data.frame) like this: GeneID MethodID Value 6 1 123 6 2 456 6 3 987 7 1 234 7 3 432 8 2 190 8 3 34 8 1 864 Note that GeneID=7 doesn't have a value for MethodID=2. Note that GeneID=8 doesn't have the

Hi, Has anyone seen this problem before, and does anyone have a solution? When I query a database using DBI and ROracle to retrieve a date in a format other than the system default, I get a core dump. System: Redhat Linux 7.3 uname: Linux 2.4.20-20.7smp #1 SMP Mon Aug 18 14:46:14 EDT 2003 i686 Oracle version 8i R --version: R 1.7.1 (2003-06-16). DBI_0.1-6, ROracle_0.5-0 Start up R, load

I might be missing something but I thought that AUC was a measure for comparing ROC curves, so there is nothing else needed to "compare" them. The larger AUC is the higher correlation of 2 variables compared. No other measures or calculations are needed. Jarek Tuszynski -----Original Message----- From: r-help-bounces at stat.math.ethz.ch [mailto:r-help-bounces at stat.math.ethz.ch] On

Dear all, I try to compare the performances of several parameters to diagnose lameness in dogs. I have several ROC curves from the same dataset. I plotted the ROC curves and calculated AUC with the ROCR package. I would like to compare the AUC. I used the following program I found on R-help archives : From: Bernardo Rangel Tura Date: Thu 16 Dec 2004 - 07:30:37 EST

When plotting using the cex.axis argument to boxplot(), the size of the plotting symbols beyond the whiskers of the boxplot are being changes. Example: par(mfrow=c(2,1)) boxplot(c(rnorm(10), 10), horizontal=TRUE, main="Test", las=2, cex.axis=2) boxplot(c(rnorm(10), 10), horizontal=TRUE, main="Test", las=2, cex.axis=1) This is from the following line in bxp()

hi, I'd like to extract data from very large files (ca 1-2 Mio lines), and I already know which lines of these files I need. Is there a way to do it without filling the memory with the whole file content? The lines I need are randomly distributed within the files. Thanks, Marc -----Original Message----- From: r-help-bounces at stat.math.ethz.ch [mailto:r-help-bounces at

Hi, I have some very noisy, relatively sparse data; a biological response of roughly ~8 subjects at ~8 times points). I've been following the data trend using a lowess line, over-plotted with several values of bandwidth, 'f <- seq(0.3, 0.9, by=0.1)'. At this point, we have no models for these data. I wonder if there is any way under R to assign some sort of confidence interval to

I have been using gnls with the weights argument (and varConstPower) to specify a variance structure for curve fits. In attempting to extract the parameters for the variance model I am seeing results I don't understand. When I simply display the model (or use "summary" on the model), I get what seem like reasonable values for both "power" and "const". When I

I have installed the binary packages for Debian GNU/Linux release 2.3 (woody) in the U.S. mirror of the CRAN archive. They should propagate to the main CRAN archive within a day and to the other mirrors within two days. These packages have been compiled with gcc-3.0 and g77-3.0. I believe the testing distribution currently provides only a snapshot of gcc-3.0, not the latest released version, so

I have installed the binary packages for Debian GNU/Linux release 2.3 (woody) in the U.S. mirror of the CRAN archive. They should propagate to the main CRAN archive within a day and to the other mirrors within two days. These packages have been compiled with gcc-3.0 and g77-3.0. I believe the testing distribution currently provides only a snapshot of gcc-3.0, not the latest released version, so

All, I was having trouble trying to create a new class of data and pass it on to splom (in the lattice library). I mentioned this to Martin Morgan after a talk he gave. Following is not so much a question, but rather an answer from Morgan that might be useful to others. Here is the edited part of an email conversation with him: On Thursday, August 26, 2010 1:36 PM, Martin Morgan wrote:

Works perfectly. Thank you. -Hugh Rand -----Original Message----- From: Spencer Graves [mailto:spencer.graves at pdf.com] Sent: Sunday, January 15, 2006 6:41 PM To: Rand, Hugh Cc: 'r-help at lists.R-project.org' Subject: Re: [R] trouble with extraction/interpretation of variance structure para meters from a model built using gnls and varConstPower How about this: >

Hi, Is it possible to add labels to variables in R (so as to have a better description of what the variables represent)? Thanks, Neil

Dear R-gurus, I have an interpretation problem regarding lme models. I am currently working on dog locomotion, particularly on some variation factors. I try to figure out which limb out of 2 generated more dispersed data. I record a value called Peak, around 20 times for each limb with a record. I repeat the records during a single day, and on several days. I tried to build two models, one

We have an aging Linux Cluster here, running Red Hat 7.3. We have some business reasons for not upgrading the OS version. I don't recall the exact hardware (dual Pentium III, 4 Gbyte RAM, 1 GHz clock?), but it was pretty good in 2001 or so. We recently tried to build and install R, ver 2.6.0, for this cluster. It built and apparently ran correctly, but it failed "make check".

Hello there, I am doing analysis on survival data. How do I pick out a probability of survival at a chosen landmark time(for example, 3 years, 4 years) from the result of "survfit"? or any other functions? As I know, the result of 'survfit' only have the probabilities for all event or cencor time. Thank you very much Lisa Wang Princess Margaret Hospital Toronto, Ca tel:

Hello How do I use function 'MERGE" to combine the FILE A and FILE B below to make FILE C? Thank you FILE A 140 151 167 30.1 11.4 40 FILE B 140 167 5.7 30.3 FILE C 140 151 167 30.1 11.4 40 5.7 NA 30.3

Event: 2005 Insightful User Conference Dates: Oct 26-27, 2005 Location: Princeton, NJ URL: http://www.insightful.com/news_events/2005uc/ for details on pricing, hotel accommodations and to register for this event. The Insightful 2005 User Conference is being held October 26th-27th in Princeton, NJ. This year's conference focuses on the techniques and methodologies pivotal to the