Re: [R] how to save plots automatically

Great! Thanks, Yihui and David. On Wed, Sep 23, 2009 at 1:04 PM, Yihui Xie wrote: > Something like > > for(i in 1:20){ > pdf(sprintf("myplot%d.pdf",i)) > # your calculation and plots here > dev.off() > } > > If you have 100 plots in total, then you should have 5 graphs in each > single pdf. >

[R] synchronisation of time series data using interpolation

Readers, I have data with different time stamps that I wish to plot (for example): data set 1 time(hh:mm:ss),datum 01:00:00,500 01:00:15,600 01:00:30,750 01:00:45,720 01:01:00,700 01:01:15,725 01:01:30,640 01:01:45,710 data set 2 time,datum 01:00:12,20 01:01:01,55 01:01:55,22 The time interval i

[R] Odp: Re tain current graphs in figure

Hi r-help-boun...@r-project.org napsal dne 25.09.2009 05:31:58: > > I want to know, how do I retain the current plot and axes properties such > that subsequent graphing commands add to the existing graph. Are you looking for ?lines, ?points Regards Petr > > Thank you very much!! > > --

[R] Re tain current graphs in figure

I want to know, how do I retain the current plot and axes properties such that subsequent graphing commands add to the existing graph. Thank you very much!! -- View this message in context: http://www.nabble.com/Retain-current-graphs-in-figure-tp25606069p25606069.html Sent from the R help mail

Re: [R] [R-SIG-Finance] Does anybody know how to connect to KDB from within R?

Michael wrote: Does anybody know how to connect to KDB from within R? Please give me some pointers... Thanks a lot! Michael, KX distributes R glue code with kdb. http://kx.com/q/interfaces/r/ Talk to them if you need support on using it, your firm presumably paid enough for the license.

Re: [R] [R-SIG-Finance] Does anybody know how to connect to KDB from within R?

Dear all, I need to compute yield corporate bond yields. I have coupon, market price and maturity. Could anybody give an advice how to compute this. I found the package termstr in R. Is this one which I need? Best, Sven -Original Message- From: r-sig-finance-boun...@stat.math.ethz.ch [ma

Re: [R] How to download from github

Hmm, clicking on the 'Download' button and then on either the 'TAR' or 'ZIP' icons is working fine for me. It might take a while for the actual download to start-- GitHub has to compress the files which can take a half a minute or more. Also, GitHub appears to be preparing for a move to a new set o

Re: [R] Penalized Logistic Regression - Query

On Thu, 24 Sep 2009, Axel Urbiz wrote: Dear R users, Is there any package that I could use to perform Penalized Logistic Regression (i.e. Ridge/Lasso regularization) including also an offset term in the model (i.e. a variable with a known coefficient of 1 rather than an estimated coefficient)?

[R] How to download from github

Hi: Is my first attempt to try to download from github. Nothing happens by clicking on the 'download' button. Could anyone give me a hint on how to get all the files from the link below? Thanks http://github.com/hadley/ggplot2-bayarea/tree/0a8bf71dea38cfbf2d928eb713d24dfd928359fc Felipe D.

Re: [R] superimposing xyplots on same scale

Sorry, doubleYScale is not appropriate, since you specifically want a common y scale. I think Baptiste was suggesting to use layer(), rather than as.layer(): something like xyplot(pct_compl ~ date, col="red", ...) + layer(panel.xyplot(date, time_pct, col = "blue")) But again, I would def

[R] simulating a model

Dear useRs, I have written an ecological model, based on the epidemiology SIR model. I've been trying to simulate it in R. However, I can't simulate it properly. Two guesses: my script isn't right; I'm not setting the parameters properly I have uploaded an image to the model here: http://img24.im

[R] Penalized Logistic Regression - Query

Dear R users, Is there any package that I could use to perform Penalized Logistic Regression (i.e. Ridge/Lasso regularization) including also an offset term in the model (i.e. a variable with a known coefficient of 1 rather than an estimated coefficient)? I couldn't find any package that would all

Re: [R] packGrob and dynamic resizing

Hi baptiste auguie wrote: Dear all, I'm trying to follow an old document to use Grid frames, Creating Tables of Text Using grid Paul Murrell July 9, 2003 As a minimal example, I wrote this, gf <- grid.frame(layout = grid.layout(1, 1), draw = TRUE) label1 <- textGrob("test", x = 0, just =

[R] Does anybody know how to connect to SAS from within R?

And what might be the benefit doing that? Thanks a lot! __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-c

Re: [R] superimposing xyplots on same scale

Um, no, please use doubleYScale() in latticeExtra, rather than as.layer(). However, in your case, why don't you just do a standard lattice plot like xyplot(time_pct + pct_compl ~ date, par.settings = simpleTheme(col = c("blue", "red")), type="b", pch=15, ...) If you gave a reproducib

Re: [R] Date formats in as.Date

Try this: > library(zoo) > as.yearmon("Sep-1981", "%b-%Y") [1] "Sep 1981" > as.Date(as.yearmon("Sep-1981", "%b-%Y")) [1] "1981-09-01" > as.Date(paste(1, "Sep-1981"), "%d %b-%Y") [1] "1981-09-01" On Thu, Sep 24, 2009 at 7:15 PM, Worik R wrote: > I have trouble with this: > > as.Date("Sep-1981"

[R] Adding to tick marks to Scatter3d

Is it possible to add tick marks to a scatter3d chart?? The initial chart is created in R Commander from the Excel add-in and I have been able to make numerous changes directly to the chart but I can't figure out how to produce tick marks.? Thanks for your help. [[alternative HTML versi

Re: [R] Creating grid graphics with grid.layout

Paul Murrell wrote: Hi If using integer indices is awkward, you could push a viewport tree with the relevant viewports given useful names. Then you can downViewport() using the names (if those plots are lattice plots, then you might even be able to use the 'draw.in' argument). Is that the s

Re: [R] keeping all rows with the same values, and not only unique ones

test[which(test[,"total"] %in% needed),] --- On Fri, 25/9/09, Dimitri Liakhovitski wrote: > From: Dimitri Liakhovitski > Subject: [R] keeping all rows with the same values, and not only unique ones > To: "R-Help List" > Received: Friday, 25 September, 2009, 8:52 AM > Dear R-ers, > > I have a

[R] subsetting from a vector or matrix

I realize this should be simple but I'm having trouble subsetting vectors and matrices, for example extracting all values meeting a certain criterion, from a vector. Cannot seem to figure out the correct syntax and help page not very helpful. Or should I be using some other function than subset.

Re: [R] Maximum likelihood estimation of parameters make no biological sense

Hi, Your results are do to using an unstable parameterization of the Von Bertalanffy growth curve, combined with the unreliable optimization methods supplied with R. I coded up your model in AD Model Builder which supplies exact derivatives through AD. I used your starting values and ran the mo

Re: [R] Date formats in as.Date

Worik, You need a day! as in: as.Date("1-Sep-1981", format="%d-%b-%Y") ## first of the month HTH, Jim Porzak Ancestry.com San Francisco, CA www.linkedin.com/in/jimporzak use R! Group SF: www.meetup.com/R-Users/ On Thu, Sep 24, 2009 at 4:15 PM, Worik R wrote: > I have trouble with this: > >

Re: [R] graphics mailing list?

Hi baptiste.auguie wrote: (Sorry about the double post earlier, googlemail is having hiccups today) 2009/9/24 Romain Francois : Why just grid ? why not a list for all kind of graphics ? I figured that a good share of the traffic on r-help might be considered "graphics"-related, while I was

Re: [R] Creating grid graphics with grid.layout

Hi If using integer indices is awkward, you could push a viewport tree with the relevant viewports given useful names. Then you can downViewport() using the names (if those plots are lattice plots, then you might even be able to use the 'draw.in' argument). Is that the sort of thing you are

[R] Date formats in as.Date

I have trouble with this: as.Date("Sep-1981", format="%b-%Y") Returns "NA" >From documentation for strftime '%b' Abbreviated month name in the current locale. (Also matches full name on input.) '%Y' Year with century. What am I doing wrong? cheers Worik [[alternative H

Re: [R] keeping all rows with the same values, and not only unique ones

Dimitri - Use %in% instead of ==: test[test$total %in% needed,] x y total 1 1 2 7 2 2 3 7 5 5 6 9 6 6 7 9 - Phil Spector Statistical Computing Facility

[R] keeping all rows with the same values, and not only unique ones

Dear R-ers, I have a data frame "test": test<-data.frame(x=c(1,2,3,4,5,6,7,8),y=c(2,3,4,5,6,7,8,9),total=c(7,7,8,8,9,9,10,10)) test I have a vector "needed": needed<-c(7,9) needed I need the result to look like this: 1 2 7 2 3 7 5 6 9 6 7 9 When I do the following: result<-test[test["total"]==n

[R] Creating grid graphics with grid.layout

Hi, I recently created a function which uses grid with a viewport constructed with grid.layout() to position four plots on a plot device. My question is what's the best way to programatically traverse each element of the viewport? The grid is 5x5 and the four plots occupy positions (2,2), (4

Re: [R] [R-SIG-Finance] Does anybody know how to connect to KDB from within R?

Cool! I've seen the following files in that folder: Index of /q/interfaces/r/w32 NameLast modified SizeDescription Parent Directory - R.dll 12-Feb-2009 18:58 2.7M c.o 12-Feb-2009 18:5823K k.h 12-Feb-2009 18

Re: [R] Access to conditioning variables (lattice)

> The problem I am facing is that I do not know how to generically access > the conditioning variables within the panel function. In this simple > case, I can achieve what I want to do with the following call : > >    do_something_with(Type[which.packet()]) > I just continued testing and I realize

Re: [R] Color of the plot which correspond to the group of the observations

On 09/25/2009 04:28 AM, FMH wrote: Dear All, Let: dp: depth of the river tp: temperature with respect to depth These pair of observations are in 3 different groups i.e: Obs. 1,3,5,7 from the first group Obs. 2,4 and 10 from second group Obs 6,8 and 9 from third group. We can have a simple

Re: [R] Non-parametric test for location with two unpaired setsof data measured on ordinal scale.

Greg, I used the term location because I did not want to use the terms mean or median for the exact reason that you gave; these to values can be different in a given distribution. I want to test the null hypothesis that the data come from a single distribution. This is often done by comparing a

Re: [R] save txt file

a<-rbinom(4,10,0.8) b<-rbinom(2,6,0.7) # See ?write.table for many options write.table(t(a), file="D:/myfile.txt") # t() to give row of data. write.table(t(b), file="D:/myfile.txt", append=TRUE) --- On Wed, 9/23/09, Eiger wrote: > From: Eiger > Subject: [R] save txt file > To: r-help@r-pr

Re: [R] Non-parametric test for location with two unpaired sets of data measured on ordinal scale.

What do you mean by location? I can think of examples where 2 distributions have the same median but different means, or the same means but different medians. Are you willing to assume that the distributions are exactly the same under the null hypothesis? (not just the same 'center/location')

Re: [R] Downloading currency data from from Yahoo

you simply had the wrong symbol. a search on yahoo.com suggests the symbol might be "GBP=x", getSymbols('GBP=x',src='yahoo') str(get("GBP=X")) An ‘xts’ object from 2007-01-03 to 2009-09-23 containing: Data: num [1:707, 1:6] 0.51 0.51 0.52 0.52 0.52 0.52 0.51 0.51 0.51 0.51 ... - attr(*, "dim

Re: [R] col headers in read.table()

Try this; mat <- read.table('test.txt',header=T,row.names=1,sep='\t', check.names = FALSE) On Thu, Sep 24, 2009 at 2:25 PM, Tim Smith wrote: > Hi, > > I was trying to read in a file test.txt, which has the following data: >  norm   norm    norm    class   class   class > a       1       2      

[R] col headers in read.table()

Hi, I was trying to read in a file test.txt, which has the following data: norm normnormclass class class a 1 2 3 4 5 6 b 3 4 5 6 7 8 c 5 6 7 8 9 10 in my R code, I do the f

[R] Non-parametric test for location with two unpaired sets of data measured on ordinal scale.

Please forgive a stats question. I have to sets of data (unpaired) measured on an ordinal scale. I want to test to see if the two sets are different (i.e. do they have the same location): set1: 1,3,2,2,4,3,3,2,2 set: 4,4,4,3,3,5,4,4 What is the most appropriate non-parametric test to test lo

Re: [R] Color of the plot which correspond to the group of the observations

Try these three options, dp <- c(1,4,3,2,5,7,9,8,9,2) tp <- 1:10 group <- factor(c(1, 2, 1, 2, 1, 3, 1, 3, 3, 2), label=letters[1:3]) plot(tp,dp, type= 'p', col = group) d <- data.frame(dp=dp, tp=tp, group=group) library(lattice) xyplot(dp~tp, data=d, groups=group, auto.key=TRUE) library(

Re: [R] Color of the plot which correspond to the group of the observations

It is easier in lattice dp <- c(1,4,3,2,5,7,9,8,9,2) tp <- 1:10 gg <- rep(1:3, c(3,3,4)) ddff <- data.frame(dp=dp, tp=tp, gg=gg) xyplot(dp ~ tp, groups=gg, data=ddff) __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLE

Re: [R] more strange behavior of Revolution R 1.3.0

The best place for questions specific to REvolution R is the REvolution forums: http://forums.revolution-computing.com/ (Note: due to anti-spam moderation, it may take a little while for your post to appear.) In this particular case, it looks to me like you're getting a significant speedup after

[R] Color of the plot which correspond to the group of the observations

Dear All,   Let: dp: depth of the river tp: temperature with respect to depth These pair of observations are in 3 different groups i.e: Obs. 1,3,5,7 from the first group Obs. 2,4 and 10 from second group Obs 6,8 and 9 from third group. We can have a simple scatter plot, between depth as y-axis a

Re: [R] panel.text question

Try this: xyplot(y ~ x | a, panel=function(x, y, subscripts, ...){ panel.loess(x, y) panel.text(0, 2, label=c('best','better','bad','worst')[tail(subscripts, 1)/100]) }) On Thu, Sep 24, 2009 at 2:45 PM, Osman Al-Radi wrote: > Dear R-help, > > I would like to add text to each

[R] multicomp plotting

I am trying to plot my multiple comparison data. Can anyone give me some input of the error I am getting. The data and code is appended below. Thanks ../Murli library(multcomp) sig.data<-structure(list(X = 1:63, Cell.lines = structure(c(1L, 6L, 13L, 25L, 33L, 42L, 2L, 7L, 14L, 26L, 34L, 43L

[R] error in amer function

Dear all, I am trying to reproduce the exemple in the vignette "Using lme4 to fit Generalized Additive Mixed Models" with my dataset. But... > mod <- amer(pasvig ~ -1 + harvf + tp(dias,by=harvf) + (1 | pac), data=exemplo) Erro em if (from == to) rep.int(from, length.out) else as.vector(c(fro

[R] panel.text question

Dear R-help, I would like to add text to each of four panels in a plot generated by xyplot in lattice library. A sample code is given below, the plot generated has the first label repeated in all panels! How can I get the labels to be different in each panel? library(lattice) x <- rnorm(400) y <

[R] Time Series Format Error

I have rec'd the following error: P34annual <- read.table("A:\\Data\\Output\\Sparrow\\Hydro_Data\\P34_Annual.txt", header=TRUE, sep=",", stringsAsFactors= FALSE, skip=1) P34annual$GS <- rep(1.86, dim(P34annual)[1]) P34annual$Depth <- as.numeric(P34annual$P34_stage) - as.numeric(P34annual$G

Re: [R] multinormial runs tests?

You can do this by simulation: Generate data from a multinomial of the same length as your data (the sample function can help) using either theoretical or observed probabilities. Measure the length of the longest run, or the number of runs (the rle function can help). Repeat this a bunch of ti

[R] Does anybody know how to connect to KDB from within R?

Please give me some pointers... Thanks a lot! __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, r

Re: [R] RODBC problem

Schalk, This worked. Thanks for the hint. Walt Schalk Heunis wrote: Walt I get the same message using R2.9.2 on Vista. Using sqlFetch(con,'Sheet1') seems to however. HTH Schalk Heunis On Thu, Sep 24, 2009 at 6:23 PM, Data Analytics Corp. wrote: Hi, I'm attempting to use the RODBC

[R] October-November R /S Courses: Nationwide (1) R/S+ Fundamentals and (2) R/S-Plus Advanced Programming. in San Francisco, New York City, Boston, Washington DC, Seattle and Salt Lake City

XLSolutions Corporation is proud to announce our October-November R /S course schedule in New York City, San Francisco, Boston, Washington DC, Seattle and Salt Lake City. Taught by top R/S+ gurus! http://www.xlsolutions-corp.com/rplus.asp (1-a) R/S-PLUS: An Introduction to R and S October-Nov

Re: [R] RODBC problem

Walt I get the same message using R2.9.2 on Vista. Using sqlFetch(con,'Sheet1') seems to however. HTH Schalk Heunis On Thu, Sep 24, 2009 at 6:23 PM, Data Analytics Corp. wrote: > Hi, > > I'm attempting to use the RODBC package on Windows Vista to import an excel > spreadsheet.  The spreadsh

Re: [R] RODBC problem

Try it without the '$' in the table name, that has worked for me in the past. -- Gregory (Greg) L. Snow Ph.D. Statistical Data Center Intermountain Healthcare greg.s...@imail.org 801.408.8111 > -Original Message- > From: r-help-boun...@r-project.org [mailto:r-help-boun...@r- > project.o

Re: [R] Access to conditioning variables (lattice)

> -Original Message- > From: r-help-boun...@r-project.org > [mailto:r-help-boun...@r-project.org] On Behalf Of Martin D. Lepage > Sent: Thursday, September 24, 2009 7:38 AM > To: r-help@r-project.org > Subject: [R] Access to conditioning variables (lattice) > > [using R version 2.8.1 (200

Re: [R] Post-Hoc tests for Friedman Test?

On Sep 20, 2009, at 9:35 PM, David Winsemius wrote: On Sep 20, 2009, at 9:05 PM, j...@terraspark.com wrote: Hi there all, This is my first post to the list and I'll first say a few things: - R is great! - The archives of this list have helped me solve all of my problems/ questions so far

Re: [R] Problem with xtabs(), exclude=NULL, and counting NA's

>> Also see addNA. Works great. Sometimes R drives me crazy, but Hadley, you make it much easier for me > That is nice. The addNA function does not exactly jump off the page for the > (too) casual reader. In the context of the OP's original problem, these > lines of code are illustrative: He

Re: [R] more strange behavior of Revolution R 1.3.0

Please discuss with REvolution support. Most of us do not have a version of REvolution R. Uwe Ligges Jason Liao wrote: It runs more than twice as slowly using 8 core than using a single core in inverting large matrix. Tested on 8 core Windows XP 64 machine. n = 1000 n.simu = 100

Re: [R] how to make a function recognize the name of an object/vector given as argument

Try this, testFun <- function(x,y) plot(x,y, main=paste("plot of",deparse(substitute(x)),"and", deparse(substitute(y))) ) a1 <- 5:8 b1 <- 9:6 testFun(a1,b1) ?deparse HTH, baptiste 2009/9/24 Wolfgang Raffelsberger : > Dear guRus, > I'd like to learn how to make a function recognize the name

[R] how to make a function recognize the name of an object/vector given as argument

Dear guRus, I'd like to learn how to make a function recognize the name of an object/vector given as argument If I have : testFun <- function(x,y) plot(x,y, main=paste("plot of",names(x),"and",names(y)) ) # this just a simple example ... a1 <- 5:8 b1 <- 9:6 testFun(a1,b1) # Returns the pl

[R] multinormial runs tests?

Dear R users, I would like to test the randomness in a series of N values (N>=2). I know that runs.test works for dichotomous factor only: x <- rep(c(1,2), 50) runs.test(factor(x)) However it doesn't work for series that can take any N values (N>2): x <- rep(c(1,2,5,4),50) runs.test(factor(x))

[R] Access to conditioning variables (lattice)

[using R version 2.8.1 (2008-12-22)] Hello, I'm trying to access the conditioning variables of an xyplot within a 'panel' function but I have not been able to figure out how to do so. Here is a simple example that describes what I wish to do (the problem lies with the commented line): dataset <-

Re: [R] scaled Schoenfeld residuals

hi thanks, I see that cox.zph is plotting and smoothing the "scaled Schoenfeld" residuals as generated by R, but since the term is already in the literature with a formula, maybe the help should clarify the offset. I found it confusing anyway. thanks for help greg Thomas Lumley tlumley at u.w

[R] more strange behavior of Revolution R 1.3.0

It runs more than twice as slowly using 8 core than using a single core in inverting large matrix. Tested on 8 core Windows XP 64 machine. > n = 1000 > n.simu = 100 > func1 = function() + { + x = rnorm(n*n) + dim(x)=c(n,n) + y = solve(x) + } > > setMKLthreads(1) > system.t

[R] lmer() vs. "fixed effects" regression

Hi, First, some quick terminology I am using: Fixed effects = model with unit dummy variables Random effects = model without unit dummy variables, integrating unit-level variance out of likelihood I am confused about the difference between the multilevel modeling framework of lmer() and a "fixed

[R] basic cubic spline smoothing

Hello, I come from a non statistics background, but R is available to me, and I needed to test an implementation of smoothing spline that I have written in c++, so I would like to match the results with R (for my unit tests) I am following http://www.nabble.com/file/p25569553/SPLINES.PDF SPLINE

[R] RODBC problem

Hi, I'm attempting to use the RODBC package on Windows Vista to import an excel spreadsheet. The spreadsheet has three worksheets the last of which is blank. Following an example in Phil Spector's book (p. 34), after creating a connection named con I did the following: > con RODBC Connecti

Re: [R] graphics mailing list?

(Sorry about the double post earlier, googlemail is having hiccups today) 2009/9/24 Romain Francois : > Why just grid ? why not a list for all kind of graphics ? I figured that a good share of the traffic on r-help might be considered "graphics"-related, while I was aiming at discussing less docu

[R] basic cubic spline smoothing (resending because not sure about pending)

Hello, I come from a non statistics background, but R is available to me, and I needed to test an implementation of smoothing spline that I have written in c++, so I would like to match the results with R (for my unit tests). I am following Smoothing Splines, D.G. Pollock (available online)

Re: [R] dotchart to barplots

I decided to use your tip and plot the bars using different shades of grey as follows barplot(t(as.matrix(intersect.data[,2:5])),col=c('black','grey40','darkgrey','white'), beside = T, horiz = T, legend.text = names(intersect.data)[-1], axes=TRUE, border=TRUE,plot.grid=F,cex=2,

Re: [R] Downloading data from from internet

Thanks for explaining this, Charlie. Just for completeness and to make things a little easier, the XML package has a function named readHTMLTable() and you can call it with a URL and it will attempt to read all the tables in the page. tbls = readHTMLTable('http://www.rateinflation.com/consumer

Re: [R] making R print on screen

Thanks a lot, everyone! On Thu, Sep 24, 2009 at 11:41 AM, Mario Valle wrote: > ?flush.console > > Ciao! >        mario > ld7631 wrote: >> Hello! >> >> I am running a "for" loop. In the loop I am producing some >> intermediary results and asking R to print it (of the type below). >> However, I not

Re: [R] making R print on screen

?flush.console Ciao! mario ld7631 wrote: > Hello! > > I am running a "for" loop. In the loop I am producing some > intermediary results and asking R to print it (of the type below). > However, I noticed - when the task is complicated and takes a lot of > time, R does not print those inter

Re: [R] making R print on screen

On Sep 24, 2009, at 10:26 AM, ld7631 wrote: Hello! I am running a "for" loop. In the loop I am producing some intermediary results and asking R to print it (of the type below). However, I noticed - when the task is complicated and takes a lot of time, R does not print those intermediary results

Re: [R] making R print on screen

On 09/24/2009 05:26 PM, ld7631 wrote: Hello! I am running a "for" loop. In the loop I am producing some intermediary results and asking R to print it (of the type below). However, I noticed - when the task is complicated and takes a lot of time, R does not print those intermediary results immed

Re: [R] xtable - print - suppress output

On Sep 24, 2009, at 10:20 AM, Thomas Lumley wrote: On Mon, 21 Sep 2009, David Winsemius wrote: On Sep 21, 2009, at 5:52 PM, Martin Batholdy wrote: I use xtable to convert data.frames to html tables. But when I use the print-command I always get the whole output printed even if I just wan

[R] making R print on screen

Hello! I am running a "for" loop. In the loop I am producing some intermediary results and asking R to print it (of the type below). However, I noticed - when the task is complicated and takes a lot of time, R does not print those intermediary results immediately, but prints them in batches - or d

Re: [R] Statistical analysis

Since todays ground water may be influenced by yesterdays rainfall, you may want to look at the dynlm package and possibly lag.plot and the zoo package. -- Gregory (Greg) L. Snow Ph.D. Statistical Data Center Intermountain Healthcare greg.s...@imail.org 801.408.8111 > -Original Message

Re: [R] How to show number in the %f format?

There is also the formatC function, whose description is Formatting numbers individually and flexibly, using 'C' style format specifications. -Don At 2:28 AM -0400 9/24/09, David Winsemius wrote: On Sep 23, 2009, at 6:42 PM, Peng Yu wrote: On Wed, Sep 23, 2009 at 5:16 PM, David

Re: [R] graphics mailing list?

Why just grid ? why not a list for all kind of graphics ? On 09/24/2009 04:34 PM, baptiste.auguie wrote: Dear all, Would it make sense to have a separate mailing list (special interest group*) for Grid graphics? (or is there one already?) I don't feel comfortable asking questions about the de

Re: [R] graphics mailing list?

On 9/24/2009 10:34 AM, baptiste.auguie wrote: Dear all, Would it make sense to have a separate mailing list (special interest group*) for Grid graphics? (or is there one already?) I don't feel comfortable asking questions about the design of new a new grid class in R-help where I'm guessing mos

[R] graphics mailing list?

Dear all, Would it make sense to have a separate mailing list (special interest group*) for Grid graphics? (or is there one already?) I don't feel comfortable asking questions about the design of new a new grid class in R-help where I'm guessing most people won't be interested. Of course having y

Re: [R] Problem with xtabs(), exclude=NULL, and counting NA's

On Sep 24, 2009, at 10:06 AM, hadley wickham wrote: wtf <- factor(x, levels(c(levels(wtf), NA), exclude=NULL) xtabs (~ wtf, exclude=NULL, na.action=na.pass) Also see addNA. That is nice. The addNA function does not exactly jump off the page for the (too) casual reader. In the context of

[R] graphics mailing list?

Dear all, Would it make sense to have a separate mailing list (special interest group*) for Grid graphics? (or is there one already?) I don't feel comfortable asking questions about the design of new a new grid class in R-help where I'm guessing most people won't be interested. Of course having y

Re: [R] Downloading data from from internet

Bogaso wrote: > > Hi all, > > I want to download data from those two different sources, directly into R > : > > http://www.rateinflation.com/consumer-price-index/usa-cpi.php > http://eaindustry.nic.in/asp2/list_d.asp > > First one is CPI of US and 2nd one is WPI of India. Can anyone please g

Re: [R] pipe data from plot(). was: ROCR.plot methods, cross validation averaging

Yes, that's exactly what I am after. Thank you for clarifying my problem for me! I'll try to dive into the plot.performance function. Best, Tim >>> Tobias Sing 9/24/2009 9:57 AM >>> Tim, if I understand correctly, you are trying to get the numerical values of averaged cross-validation curves.

Re: [R] xtable - print - suppress output

On Mon, 21 Sep 2009, David Winsemius wrote: On Sep 21, 2009, at 5:52 PM, Martin Batholdy wrote: I use xtable to convert data.frames to html tables. But when I use the print-command I always get the whole output printed even if I just want to save the html table into a variable; table <- prin

Re: [R] scaled Schoenfeld residuals

On Wed, 23 Sep 2009, Greg Dropkin wrote: hi sorry if this has been discussed before, but I'm wondering why the scaled Schoenfeld residuals do not follow the defining formula for obtaining them from the ordinary Schoenfeld residuals, but are instead offset by the estimated parameter values. B

Re: [R] pipe data from plot(). was: ROCR.plot methods, cross validation averaging

David, Thank you for your reply. Yes, I can access the y-values slot with p...@y-values but, note that in the cross-validation example (ROCR.xval), the plot function averages across the list of ten vectors in the y-values slot. I might be able to create a function to average across these ten ve

[R] Modelling

Dear R-users, Suppose I have the followin g sample of data, 0 1 2 4 3 1 2 1 3 1 1 3 3 4 1 0 1 2 1 2 1 4 1 4 2 1 2 2 1 1 The first variable is the response variable where 0 is defective and 1 normal. The other four factors( x1,x2,x3,x4) that influence th

Re: [R] Problem with xtabs(), exclude=NULL, and counting NA's

> wtf <- factor(x, levels(c(levels(wtf), NA), exclude=NULL) > xtabs (~ wtf, exclude=NULL, na.action=na.pass) Also see addNA. Hadley -- http://had.co.nz/ __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read

Re: [R] Bug

You have found a bug. It would be best to use dput(test1) to show unambiguously display what is in test1 but in the absence of that I will assume that its as in test1 shown below. > library(sqldf) > test1 <- data.frame(sale_date = as.Date(c("2008-08-01", "2031-01-09", + "1990-01-03", "2007-02-03"

Re: [R] Statistical analysis

Rainfall data is widely accepted as Random walk process and hence it is non-stationary. Therefore if correlation or regression coef. is measured on raw data then you may land in the world of spurious measures. I would suggest you to check whether unit root is there in your data or not first. If it

Re: [R] P-value and R-squared variable selection criteria

Don't throw out the baby with the bath water just yet. Note that even though your first model is insignificant, the R-squared is very high. This is because you fit the whole model with intercept and three coefficients on 1 degree of freedom. You need to first import the data, then run the model, an

Re: [R] pipe data from plot(). was: ROCR.plot methods, cross validation averaging

Tim, if I understand correctly, you are trying to get the numerical values of averaged cross-validation curves. Unfortunately the plot function of ROCR does not return anything in the current version (it's a good suggestion to change this). If you want a quick fix, you could change the plot.perfo

Re: [R] superimposing xyplots on same scale

Hi, try ?as.layer in the latticeExtra package. HTH, baptiste 2009/9/24 Larry White : > I have two xyplots that i want to superimpose (code below).  By default they > are displayed on slightly different y scales (one runs from 10 to 25, the > other from 10 to 30). I would like to force them both

Re: [R] any advice on web interfaces to R?

Thanks! On Wed, Sep 23, 2009 at 9:42 PM, jverzani wrote: > Mitchell Maltenfort gmail.com> writes: > >> >> I saw http://cran.r-project.org/doc/FAQ/R-FAQ.html#R-Web-Interfaces >> and I'm still not sure yet which platform (Linux, Windows, etc.) I'll >> be working on -- and no, it's not under my

Re: [R] pipe data from plot(). was: ROCR.plot methods, cross validation averaging

On Sep 24, 2009, at 9:09 AM, Tim Howard wrote: All, I'm trying again with a slightly more generic version of my first question. I can extract the plotted values from hist(), boxplot(), and even plot.randomForest(). Observe: # get some data dat <- rnorm(100) # grab histogram data hdat <- h

Re: [R] P-value and R-squared variable selection criteria

Lucas This problem is very old --- older than keypunches. There are several methods for selecting variables (forward, backwards, both, all subsets) using a variety of criteria (p-values, R^2, adjusted R^2, Cp, AIC, BIC, and more). Be sure you understand the methods, especially the tendency to

[R] Downloading data from from internet

Hi all, I want to download data from those two different sources, directly into R : http://www.rateinflation.com/consumer-price-index/usa-cpi.php http://eaindustry.nic.in/asp2/list_d.asp First one is CPI of US and 2nd one is WPI of India. Can anyone please give any clue how to download them dir

Re: [R] pipe data from plot(). was: ROCR.plot methods, cross validation averaging

Whoops, sorry. Here is the full set with the missing lines: library(ROCR) data(ROCR.xval) pred <- prediction(ROCR.xval$predictions, ROCR.xval$labels) perf <- performance(pred,"tpr","fpr") RCdat <- plot(perf, avg="threshold") RCdat Thanks. Tim >>> David Winsemius 9/24/2009 9:25 AM >>> On Sep 24,

  1   2   >