[R] ROCR package question

2009-07-24 Thread Waverley
I use ROCR to plot multiple runs' performance. Using the sample code as example: # plot ROC curves for several cross-validation runs (dotted # in grey), overlaid by the vertical average curve and boxplots # showing the vertical spread around the average. data(ROCR.xval) pred <- prediction(ROCR.xv

Re: [R] a question about replacing a column that is factor

2009-07-24 Thread milton ruser
May be you want something like: FDLFDPdata[FDLFDPdata$ID==69912008,]$tempcalved2<-'2008-03-08' FDLFDPdata$tempcalved2<-factor(FDLFDPdata$tempcalved2) ## the factor above just update the levels Case this is what you want, you can also change $ID==69912008 $ID %in% c(69912008,)

[R] how to smooth timeseries without the lagging?

2009-07-24 Thread Michael
Hi all, If I use a moving average, it will smooth the choppy time series, but it will lead to lagging... How do I smooth timeseries without the lagging effect? Thanks! __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help P

[R] a question about replacing a column that is factor

2009-07-24 Thread jlfmssm
Hi, everyone, I want to use a new value to replace a column that is a factor. However this column is replaced by , which is not what I want. I know this is because the new value is not in the levels of that column, but I don't know how to fix it. Can someone help me to figure this out? The follo

Re: [R] How to test frequency independence (in a 2 by 2 table) withmany missing values

2009-07-24 Thread Daniel Malter
Hi Tal, it is necessary to understand the method, but not for all applied researchers it is necessary, I think, to understand all the specifics or underlying processes (I know others would disagree on that). As to my suggestion: The Monte-Carlo simulations have no other purpose than showing that th

[R] regex expression to select row or column

2009-07-24 Thread Junqian Gordon Xu
I have a multidimensional data which looks like the following: "S1-a" "S2-b" "S3-c" "S4-d" "S5-a" "S6-b" "S7-c" "S8-d" "T1-A" "T1-B" "T1-C" "T1-D" "T2-A" "T2-B" "T2-C" "T2-D" I read it from csv file and would like to have 16 separate data frames like this "S1-a" "S2-a" "S1-b"

Re: [R] Making rq and bootcov play nice

2009-07-24 Thread John Gardner
Amazing! Thanks Frank. - John On Fri, Jul 24, 2009 at 7:23 PM, Frank E Harrell Jr wrote: > Quantile regression is now supported in the Design package and bootcov has > been updated accordingly.  To get the new code run > > require(Design) > source('http://biostat.mc.vanderbilt.edu/cgi-bin/viewvc.

[R] R programmer

2009-07-24 Thread Mantel, Michael
Good afternoon: I am a Senior Staffing Specialist and the Technical Recruiter for mount Sinai Medical Center on the Upper East side of Manhattan NYC. I am NOT an agent. I am recruiting for an R programmer who also has some JAVA experience. If you are interested or know anyone who is please ge

[R] Lists

2009-07-24 Thread Alberto Lora M
Hi Everybody I have the following problem suppose that we a<-c("uno","dos","tres") I am working with a while cycle and the idea is in each iteration adding an item to a list In the first iteration the resultshould be: [[1]] [1] "uno" In the second [[1]] [1] "uno" [[2]] [1] "dos" And the fina

[R] How to convert a matrix to a dataframe and makes its first row be the column name

2009-07-24 Thread yabado
Hi I want to convert a matrix to a dataframe. The matrix is like this Site A B C D X 1 2 3 4 Y 5 6 7 8 Z 9 1011 12 The converted dataframe need to have the first row as column names and the values of each column (1, 2, 3,

Re: [R] Aggregate, max and time of max

2009-07-24 Thread Carlos J. Gil Bellosta
Hello, I believe that by( data.ex, data.ex[,c(3,4)], function(x) x[which.max(x[,1]),] ) does what you want. Then, do.call( rbind, by( data.ex, data.ex[,c(3,4)], function(x) x[which.max(x[,1]),] ) ) looks somewhat nicer. Best regards, Carlos J. Gil Bellosta http://www.datanalytics.com On Fr

Re: [R] [Rd] How to create a permanent dataset in R.

2009-07-24 Thread S Ellison
>>> Liviu Andronic 07/24/09 7:06 PM >>> >.. use save.image() to save R's workspace and re-load it when > re-opening R via load(). Or you could just use save.image in your default directory for R and R will then open it automatically. I'd use that sparingly though; it is not always a blessing to h

Re: [R] Making rq and bootcov play nice

2009-07-24 Thread Frank E Harrell Jr
Quantile regression is now supported in the Design package and bootcov has been updated accordingly. To get the new code run require(Design) source('http://biostat.mc.vanderbilt.edu/cgi-bin/viewvc.cgi/*checkout*/Design/trunk/R/bootcov.s') source('http://biostat.mc.vanderbilt.edu/cgi-bin/viewvc.

Re: [R] GGPlot 2 and odfweave

2009-07-24 Thread Felipe Carrillo
Hi John: Use print before qplot or ggplot. Something like this: <>= print(qplot(subscriber,data=TV09,geom="bar",fill=subscriber, xlab=NULL,ylab=NULL)+coord_flip()+ opts(legend.position="none")) @ Felipe D. Carrillo  Supervisory Fishery Biologist  Department of the Interior  US Fi

Re: [R] Computer Modern

2009-07-24 Thread cameron.bracken
Mr Derik wrote: > > I am trying to use computer modern fonts in postscript files for a latex > document. Ultimately I want to automate this through sweave. > Hi There- I admit this is shameless self promotion but might I suggest the tikzDevice package. I am one of the developers and we jus

Re: [R] How to test frequency independence (in a 2 by 2 table) withmany missing values

2009-07-24 Thread Tal Galili
Hello dear Robert and Daniel, I read your replies with much interest, thank you very much for them! Dear Daniel, Your approach is, as Robert said, "whopped me alongside the head" as well. So much so that I fear it will be much more then my current level of understanding can take for the project I

Re: [R] Navigate to Index page of a package from R command prompt

2009-07-24 Thread Steven McKinney
Okay, I've worked out much of the details on my PC and Mac. What I'd like to be able to do is enter the command > help(package = "survival") or > ?survival:: and get the usual hyperlinked help page displayed instead of the static "text only" display or an error message. The help() function re

Re: [R] Most elegant way to use formals() in building functions

2009-07-24 Thread Duncan Murdoch
On 24/07/2009 6:35 PM, Frank E Harrell Jr wrote: Dear Group: I want to create a function having a ... argument and to have the default arguments evaluated, as thus: g <- function(a, b, ...) a+b formals(g) <- alist(a=,b=2+3,...=) g function (a, b = 2 + 3, ...) a + b But I want the default arg

Re: [R] single row/column-indexing on matrices

2009-07-24 Thread Murat Tasan
awesome! thanks much! On Jul 24, 6:41 pm, Duncan Murdoch wrote: > On 24/07/2009 6:34 PM, Murat Tasan wrote: > > > > > > > hi all - quick question: > > > i have a matrix m, say nrow=5, ncol=4. > > in a function i'd like to retrieve certain rows or columns from m, but > > which rows/cols are not k

Re: [R] single row/column-indexing on matrices

2009-07-24 Thread Duncan Murdoch
On 24/07/2009 6:34 PM, Murat Tasan wrote: hi all - quick question: i have a matrix m, say nrow=5, ncol=4. in a function i'd like to retrieve certain rows or columns from m, but which rows/cols are not known ahead of time. the function should return a sub-matrix (i.e. still of class 'matrix'). w

[R] Most elegant way to use formals() in building functions

2009-07-24 Thread Frank E Harrell Jr
Dear Group: I want to create a function having a ... argument and to have the default arguments evaluated, as thus: g <- function(a, b, ...) a+b formals(g) <- alist(a=,b=2+3,...=) g function (a, b = 2 + 3, ...) a + b But I want the default argument for b to be evaluated as 5. How can this b

[R] single row/column-indexing on matrices

2009-07-24 Thread Murat Tasan
hi all - quick question: i have a matrix m, say nrow=5, ncol=4. in a function i'd like to retrieve certain rows or columns from m, but which rows/cols are not known ahead of time. the function should return a sub-matrix (i.e. still of class 'matrix'). when selecting a single column (or row), the

[R] Specifying package dependencies

2009-07-24 Thread Christopher Bare
Hi, I'm trying to put together an R package. My library has dependencies on three other libraries: RSQLite, gaggle (a component of Bioconductor), and another package that is available as a download from the author's website. Is it possible to specify each of those dependencies in such a way that

[R] Replace rgb with different colors for image display

2009-07-24 Thread Bob Meglen
I have created pseudo-color images from principal component scores derived from multispectral data. I am using BiOps package stmt such as, imgGetRGBFromBands(PC_1,PC_2,PC_3) where PC_x are grey scale images. I would like to display these 3 pc scores as rgb images where all three r g and b la

[R] suggestion for paired t-tests

2009-07-24 Thread Jack Tanner
There's a funny inconsistency in how t.test handles paired=T or paired=F. If x and y parameters are lists, paired=F works, but paired=T doesn't. > lg=read.csv("my.csv") > a = subset(lg, condition=="a")["score"] > b = subset(lg, condition=="b")["score"] > t.test(a,b) > t.test(a,b, paired=TRUE) Erro

Re: [R] How to pass a character argument which contains expressions to arg.names in barplot?

2009-07-24 Thread roger koenker
Surely a fortune, because fortuna can be cruel. url:www.econ.uiuc.edu/~rogerRoger Koenker emailrkoen...@uiuc.eduDepartment of Economics vox: 217-333-4558University of Illinois fax: 217-244-6678Urbana, IL 61801 On Jul 24

Re: [R] How to test frequency independence (in a 2 by 2 table) withmany missing values

2009-07-24 Thread Robert A LaBudde
Daniel's response about using a mixed model "whopped me alongside the head", so now I'm thinking more clearly. Here are a few more comments: 1. The LRT based on ML or the resampling alternative is probably the most powerful among reasonable algorithms to do a test on your problem. You might

Re: [R] [Rd] How to create a permanent dataset in R.

2009-07-24 Thread Ben Bolker
Liviu Andronic wrote: > > (redirecting to r-help; it seems more appropriate for such a question) > Hello, > > On Fri, Jul 24, 2009 at 8:11 AM, Albert EINstEIN > wrote: >> Actually, we know that If we create a dataset in R ,after closing the >> session the dataset automatically is closed. I tri

Re: [R] How to pass a character argument which contains expressions to arg.names in barplot?

2009-07-24 Thread jcano
Dear Uwe Sorry for not included the data, but even as that, it worked! Thank you very much for your unvaluable time and effort Cheers Javier Cano -- View this message in context: http://www.nabble.com/How-to-pass-a-character-argument-which-contains-expressions-to-names.arg-in-barplot--t

Re: [R] How to test frequency independence (in a 2 by 2 table) withmany missing values

2009-07-24 Thread Daniel Malter
Hi Tal, you can use the lme4 library and use a random effects model. I will walk you through the example below (though, generally you should ask a statistician at your school about this). At the very end, I will include a loop for Monte-Carlo simulations that shows that the estimation of the fixe

Re: [R] Include files?

2009-07-24 Thread Mark Knecht
On Fri, Jul 24, 2009 at 10:22 AM, Duncan Murdoch wrote: > On 7/24/2009 12:25 PM, Mark Knecht wrote: >> >> Hi, >>   I have 15 or 20 functions I've written to convert the sort of data >> I'm working with. They are currently in their own R file which I load >> by hand in Rgui before loading and runnin

Re: [R] Duplicated date values aren't duplicates

2009-07-24 Thread Tim Clark
Don and Jim, Thanks! I got it! Duplicated is only returning one of the two duplicated dates (the second date). It all makes sense now! Tim Tim Clark Department of Zoology University of Hawaii --- On Fri, 7/24/09, Don MacQueen wrote: > From: Don MacQueen > Subject: Re: [R] Duplicated

Re: [R] Colors column don't give correct color with plot()

2009-07-24 Thread Vivek Ayer
Problem solved! Thanks! I'm still a noob at learning R, so eventually I'll figure out all these different datatypes. Vivek On Fri, Jul 24, 2009 at 11:35 AM, Henrique Dallazuanna wrote: > Isn't a bug, > > try this > > plot(d[,1],d[,2],col=as.character(d[,3])) > > its because d[,3] is a factor. > >

Re: [R] Colors column don't give correct color with plot()

2009-07-24 Thread Henrique Dallazuanna
Isn't a bug, try this plot(d[,1],d[,2],col=as.character(d[,3])) its because d[,3] is a factor. Or in read.csv: d <- read.csv("blah.csv",header=FALSE,sep="", stringsAsFactor = FALSE) On Fri, Jul 24, 2009 at 3:23 PM, Vivek Ayer wrote: > Hey guys, > > I'm having trouble getting the correct col

[R] Colors column don't give correct color with plot()

2009-07-24 Thread Vivek Ayer
Hey guys, I'm having trouble getting the correct colors, when they are read from a csv file. Here's an example: 1 1 black 1 2 green 1 3 green 1 4 black 1 5 peachpuff Call it blah.csv. Then I ran: data <- read.csv("blah.csv",header=FALSE,sep="") and then to plot the data, I ran: plot(data[,1],

Re: [R] Colors column don't give correct color with plot()

2009-07-24 Thread Sarah Goslee
Not a bug, just another failure to understand what R is doing when it imports data or to actually look at your data. By default, R converts strings to factors during the import process, and as you found that completely changes the results (converts them to numbers, and uses those to specify colors

Re: [R] Colors column don't give correct color with plot()

2009-07-24 Thread Erik Iverson
My guess is that > class(data[,3]) gives "factor" as the result, try casting to character if so. -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Vivek Ayer Sent: Friday, July 24, 2009 1:24 PM To: r-help@r-project.org Subject: [

[R] Colors column don't give correct color with plot()

2009-07-24 Thread Vivek Ayer
Hey guys, I'm having trouble getting the correct colors, when they are read from a csv file. Here's an example: 1 1 black 1 2 green 1 3 green 1 4 black 1 5 peachpuff Call it blah.csv. Then I ran: data <- read.csv("blah.csv",header=FALSE,sep="") and then to plot the data, I ran: plot(data[,1],

Re: [R] How to test frequency independence (in a 2 by 2 table) with many missing values

2009-07-24 Thread Robert A LaBudde
1. Drop all subjects where both test results are not present. These are uninformative on the difference. 2. Perform McNemar test on the remaining table. (Only the different pairs are informative.) This will give you a p-value on the data that actually contrasts the two tests. (For your examp

Re: [R] metafor

2009-07-24 Thread Chuck Cleland
On 7/24/2009 11:40 AM, Frank Pearson wrote: > I had found the author's (Wolfgang Viechtbauer) earlier meta-analytic code in > R, MiMa, useful. so I have been exploring metafor using an example dataset > from MiMa. metafor provides a lot more. However, MiMa provided parameter > estimates, stand

Re: [R] CI wiskers

2009-07-24 Thread Daniel Malter
Here it is... http://www.nabble.com/plotting-confidence-intervals-td24482119.html Marc's answer is probably the way to go Daniel - cuncta stricte discussurus - -Ursprüngliche Nachricht- Von: Daniel Malter [mailto:dan...@umd.edu] Gesend

Re: [R] CI wiskers

2009-07-24 Thread Daniel Malter
I recently answered this (or, rather, basically the same) question in a thread with an example in it. I am trying to find it back. Just a sec - cuncta stricte discussurus - -Ursprüngliche Nachricht- Von: r-help-boun...@r-project.org [mailto

Re: [R] [Rd] How to create a permanent dataset in R.

2009-07-24 Thread Liviu Andronic
(redirecting to r-help; it seems more appropriate for such a question) Hello, On Fri, Jul 24, 2009 at 8:11 AM, Albert EINstEIN wrote: > Actually, we know that If we create a dataset in R ,after closing the > session the dataset automatically is closed. I tried for creating dataset > permanently in

[R] CI wiskers

2009-07-24 Thread mfreidin
I have a matrix containing means and CIs (lower and upper in two columns, so three columns for every data point) for several points. I have to build a graph of these means accompained by the CIs (as wiskers). No problems with making the graph of means, but I don't know how to introduce CIs. Can

[R] How to test frequency independence (in a 2 by 2 table) with many missing values

2009-07-24 Thread Tal Galili
Hello dear R help group. My question is statistical and not R specific, yet I hope some of you might be willing to help. *Experiment settings*: We have a list of subjects. each of them went through two tests with the answer to each can be either 0 or 1. *Goal*: We want to know if the two experim

Re: [R] Trouble with If statements

2009-07-24 Thread PDXRugger
Alright guys i figured it out, no need to reply. Sorry for any bother... PDXRugger wrote: > > Hey everyone, > Consider the below. I would like to look up each of the items in > "data" and store the result in "BinStore". In this example it isnt > storing any value but the last. In my

[R] Trouble with If statements

2009-07-24 Thread PDXRugger
Hey everyone, Consider the below. I would like to look up each of the items in "data" and store the result in "BinStore". In this example it isnt storing any value but the last. In my actual code i am getting another error altogether, i get a 1:In if (VacAcresVals.CandTaz <= 4) (BinNumbe

[R] metafor

2009-07-24 Thread Frank Pearson
I had found the author's (Wolfgang Viechtbauer) earlier meta-analytic code in R, MiMa, useful. so I have been exploring metafor using an example dataset from MiMa. metafor provides a lot more. However, MiMa provided parameter estimates, standard errors, z values, etc. for individual moderators

Re: [R] Include files?

2009-07-24 Thread Duncan Murdoch
On 7/24/2009 12:25 PM, Mark Knecht wrote: Hi, I have 15 or 20 functions I've written to convert the sort of data I'm working with. They are currently in their own R file which I load by hand in Rgui before loading and running my main programs. Is there any way to have this file included in

Re: [R] ggplot2 : commands on one line vs two lines.

2009-07-24 Thread John Kane
--- On Fri, 7/24/09, Kingsford Jones wrote: > From: Kingsford Jones > Subject: Re: [R] ggplot2 : commands on one line vs two lines. > To: "John Kane" > Cc: r-h...@stat.math.ethz.ch > Received: Friday, July 24, 2009, 2:26 AM > Hi John, > > Thank you for the easily cut-and-pastable example. >

Re: [R] Stting non-overlaping text positions

2009-07-24 Thread Greg Snow
A general solution is not simple, some options include: The thigmophobe.labels function in the plotrix package will place labels on the side of a point furthest from the closest point (works well in some cases, but does not guarantee non-overlap). The dynIdentify and TkIdentify functions in the

Re: [R] Include files?

2009-07-24 Thread cls59
Mark Knecht wrote: > > Hi, >I have 15 or 20 functions I've written to convert the sort of data > I'm working with. They are currently in their own R file which I load > by hand in Rgui before loading and running my main programs. > >Is there any way to have this file included in my R pr

Re: [R] Stting non-overlaping text positions

2009-07-24 Thread John Kane
You might want to have a look at thigmophobe.labels in the plotrix package. --- On Fri, 7/24/09, Ana Conesa wrote: > From: Ana Conesa > Subject: [R] Stting non-overlaping text positions > To: r-help@r-project.org > Received: Friday, July 24, 2009, 8:29 AM > Dear List > > I am having problems

[R] Aggregate, max and time of max

2009-07-24 Thread Afshartous, David
All, For data consisting of serial measurements on subjects, one may use the aggregate function to say compute the peak response for each subject for each design condition. Is there a way to alter this or another one-liner to also retain the time at which the peak occurred and thus avoid writing

Re: [R] How to pass a character argument which contains expressions to arg.names in barplot?

2009-07-24 Thread Uwe Ligges
jcano wrote: Hi all Can anybody help me with this? I am trying to include in an automatic way the argument in arg.names in a barplot. I generate the labels I want to appear below the bars with a for loop, and they contain subscripts, so I need to use expression anch<-0.05 esp<-4 for (i

Re: [R] CI wiskers

2009-07-24 Thread John Kane
Or a simple approach is to use arrows to draw the CI whiskers Example == # Uses arrows to produce confidence intervals for a set of values. low <- c(312.9460, 312.9419, 312.9422, 312.9380 ) mass <- c(312.9476, 312.9435, 312.9438

Re: [R] computing the radius of an arc

2009-07-24 Thread Greg Snow
There is a function "rsaga.local.morphometry" in the RSAGA package that says it computes curvature (among other things). It looks like that function was designed for a different type of data than yours is, but it may work, or if not, then you may be able to adapt some of the code to work with y

Re: [R] Include files?

2009-07-24 Thread Mark Knecht
Thanks. I think source is what I was looking for. Cheers, Mark On Fri, Jul 24, 2009 at 9:30 AM, Erik Iverson wrote: > Several solutions, easiest might be to use the source function in your > program.  Also, you could create a .Rprofile if you use them every time you > use R. > > See ?source and

Re: [R] Fwd: Making rq and bootcov play nice

2009-07-24 Thread Frank E Harrell Jr
roger koenker wrote: John, You can make a local version of bootcov which either: deletes these arguments from the call to fitter, or modify the switch statement to include rq.fit, the latter would need to also modify rq() to return a fitFunction component, so the first option is simpl

Re: [R] Include files?

2009-07-24 Thread jim holtman
?Startup ?.Rprofile On Fri, Jul 24, 2009 at 12:25 PM, Mark Knecht wrote: > Hi, >   I have 15 or 20 functions I've written to convert the sort of data > I'm working with. They are currently in their own R file which I load > by hand in Rgui before loading and running my main programs. > >   Is ther

Re: [R] Include files?

2009-07-24 Thread Erik Iverson
Several solutions, easiest might be to use the source function in your program. Also, you could create a .Rprofile if you use them every time you use R. See ?source and ?Startup Erik -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Beha

[R] Include files?

2009-07-24 Thread Mark Knecht
Hi, I have 15 or 20 functions I've written to convert the sort of data I'm working with. They are currently in their own R file which I load by hand in Rgui before loading and running my main programs. Is there any way to have this file included in my R program like #include might in C?

Re: [R] how to calculate growth rate of a variable

2009-07-24 Thread Albyn Jones
Your data will have all sorts of patterns (diurnal, seasonal) in addition to long term trend. I'd start by smoothing out the cyclic patterns with loess or gam, then use a secant approximation to the slope on the smoothed series. albyn On Fri, Jul 24, 2009 at 06:13:00PM +0530, Yogesh Tiwari wro

[R] Fwd: Making rq and bootcov play nice

2009-07-24 Thread roger koenker
John, You can make a local version of bootcov which either: deletes these arguments from the call to fitter, or modify the switch statement to include rq.fit, the latter would need to also modify rq() to return a fitFunction component, so the first option is simpler. One of t

Re: [R] CI wiskers

2009-07-24 Thread Steve Lianoglou
Hi, On Jul 24, 2009, at 11:43 AM, Dieter Menne wrote: mfreidin wrote: I have a matrix containing means and CIs (lower and upper in two columns, so three columns for every data point) for several points. I have to build a graph of these means accompained by the CIs (as wiskers). No probl

Re: [R] String manipulations

2009-07-24 Thread Marc Schwartz
On Jul 24, 2009, at 10:35 AM, Peng Yu wrote: Hi, I am looking for a reference page that gives me a complete list of R string manipulation operations. At this moment, I am looking for something like basename and dirname. The following webpage doesn't give me any information. Can somebody point m

Re: [R] String manipulations

2009-07-24 Thread Gabor Grothendieck
Try this from within R: help.search(keyword = "character", package = "base") On Fri, Jul 24, 2009 at 11:35 AM, Peng Yu wrote: > Hi, > > I am looking for a reference page that gives me a complete list of R > string manipulation operations. At this moment, I am looking for > something like basenam

Re: [R] CI wiskers

2009-07-24 Thread Dieter Menne
mfreidin wrote: > > I have a matrix containing means and CIs (lower and upper in two columns, > so three columns for every data point) for several points. I have to build > a graph of these means accompained by the CIs (as wiskers). No problems > with making the graph of means, but I don't know

Re: [R] EOF revisited

2009-07-24 Thread Steve Lianoglou
Hi, On Jul 24, 2009, at 11:12 AM, wrote: I learnt from this forum to test for EOF reached with fiunction readLines as follows: con <- file("MyFle.txt","r") repeat { line <- readLines(con,n=1) if (length(line) == 0) break } It works fine if I read one line at a time. Since now I hav

[R] String manipulations

2009-07-24 Thread Peng Yu
Hi, I am looking for a reference page that gives me a complete list of R string manipulation operations. At this moment, I am looking for something like basename and dirname. The following webpage doesn't give me any information. Can somebody point me a better resource? http://wiki.r-project.org/

[R] EOF revisited

2009-07-24 Thread mauede
I learnt from this forum to test for EOF reached with fiunction readLines as follows: con <- file("MyFle.txt","r") repeat { line <- readLines(con,n=1) if (length(line) == 0) break } It works fine if I read one line at a time. Since now I have a huge file so that it would take forever

Re: [R] Lag representation in ccf() while zoo object is used?

2009-07-24 Thread Gabor Grothendieck
On Fri, Jul 24, 2009 at 10:31 AM, Keith wrote: > Dear All, > > I have 2 time-series data sets and would like to check the cross > correlation. These data sets were set as a zoo object, called data, > and in general look like: > >                                        V1             V2 > 2007-01-01

Re: [R] computing the radius of an arc

2009-07-24 Thread Nair, Murlidharan T
Thanks, for all the suggestions. Indeed I am interested in computing curvature (http://en.wikipedia.org/wiki/Curvature) of the curve that is given as a discrete set of points. I have the curve in 3D but I can certainly write it in 2D as well. Code is given at the bottom. I will try with the sugg

[R] Stting non-overlaping text positions

2009-07-24 Thread Ana Conesa
Dear List I am having problems to set text labels in a scatter plot. At some areas dots are close and labels overlap and texts cannot be read. I have too many dots to adjust this manually. Is there any function that will calculate the label offsets to avoid this overlapping problem? Thanks

Re: [R] str(data.frame) after subsetting reflects original structure, not subsetted structure?

2009-07-24 Thread Marc Schwartz
On Jul 24, 2009, at 8:17 AM, Bryan Hanson wrote: I find that after subsetting (you may prefer "conditional selection") a data frame and assigning it to a new object, the str(new object) reflects the original data frame, not the new one: A <- rnorm(20) B <- factor(rep(c("t", "g"), 10)) C <-

[R] Lag representation in ccf() while zoo object is used?

2009-07-24 Thread Keith
Dear All, I have 2 time-series data sets and would like to check the cross correlation. These data sets were set as a zoo object, called data, and in general look like: V1 V2 2007-01-01 00:00:00 0.0 0.176083 2007-01-01 01:00:00 0.0

[R] how to calculate growth rate of a variable

2009-07-24 Thread Yogesh Tiwari
Dear R Users, If a variable, say CO2(ppm), is varying with time. Then how to calculate CO2 (ppm) growth rate /a-1 I have CO2 time series (1991-2000), as: time, year, month, day, hour, min, sec, lat, long, height, CO2 1991.476722 1991 6 24 0 5 0 -38.93 145.15 4270 353.680 1991.476741 1991 6 24 0 15

Re: [R] roclet

2009-07-24 Thread hadley wickham
On Fri, Jul 24, 2009 at 9:08 AM, Steve Lianoglou wrote: > Hi, > > On Jul 24, 2009, at 8:49 AM, > wrote: > >> Hi, >>  can anyone tell me what a "roclet" is exactly? >>  Is this terminology specific to R, or is it used in other languages >> aswell? > > Never heard of it before and google isn't retu

Re: [R] str(data.frame) after subsetting reflects original structure, not subsetted structure?

2009-07-24 Thread Bryan Hanson
Thanks Marc and Ben... Your answers were most helpful. I suspected something had been written about it, but was having trouble formulating a reasonable search query. I was looking in the help page for str(), which was sort of a dead end. Bryan * Bryan Hanson Professor of Chemistry

Re: [R] roclet

2009-07-24 Thread Steve Lianoglou
Hi, On Jul 24, 2009, at 8:49 AM, > wrote: Hi, can anyone tell me what a "roclet" is exactly? Is this terminology specific to R, or is it used in other languages aswell? Never heard of it before and google isn't returning anything interesting or related to R. Can you put it into co

Re: [R] Duplicated date values aren't duplicates

2009-07-24 Thread Don MacQueen
Look at results of table( mydata$DateTime ) and I think you will see that some are duplicated. Specifically, the two in your dupes object. -Don At 5:50 PM -0700 7/23/09, Tim Clark wrote: Dear list, I just had a function (as.ltraj in Adehabitat) give me the following error: "Error in as.

[R] Does roxygen work with S4 classes

2009-07-24 Thread Renaud Gaujoux
Hi, I'm trying to run roxygen to generate Rdoc for code that define S4 classes. Apparently the files generated/updated by function 'roxygenize' don't follow the naming convention of package.skeleton myClass-class.Rd and myGeneric-methods.Rd. I've got duplicate files for each class and method,

Re: [R] str(data.frame) after subsetting reflects original structure, not subsetted structure?

2009-07-24 Thread Ben Bolker
Bryan Hanson wrote: > > I find that after subsetting (you may prefer "conditional selection") a > data > frame and assigning it to a new object, the str(new object) reflects the > original data frame, not the new one: > > A <- rnorm(20) > B <- factor(rep(c("t", "g"), 10)) > C <- factor(rep(c("

Re: [R] multiple error bars when NA present

2009-07-24 Thread William Revelle
At 1:38 AM -0500 7/24/09, Junqian Gordon Xu wrote: package: psych function: error.bars() R version: 2.9.1 OS: both linux and windows a<-c(1,2,3,4,5) b<-c(1,2,3,4,4) c<-c(1,2,3,NA,5) data<-data.frame(a,b,c) error.bars(data[,1:2],ylim=c(0,6)) looks fine however error.bars(data,ylim=c(0.6)) s

[R] str(data.frame) after subsetting reflects original structure, not subsetted structure?

2009-07-24 Thread Bryan Hanson
I find that after subsetting (you may prefer "conditional selection") a data frame and assigning it to a new object, the str(new object) reflects the original data frame, not the new one: A <- rnorm(20) B <- factor(rep(c("t", "g"), 10)) C <- factor(rep(c("h", "l"), 10)) D <- data.frame(A, B, C) s

[R] roclet

2009-07-24 Thread news
Hi, can anyone tell me what a "roclet" is exactly? Is this terminology specific to R, or is it used in other languages aswell? -- aleblanc __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting

Re: [R] legend title in qplot

2009-07-24 Thread RON70
PLEASE IGNORE MY PREVIOUS MAIL RON70 wrote: > > Thank you for your reply. Yes I am reading that document and ofcourse > simultaniously trying to apply in my current problem as well. Yes still I > am learning. > > Here I have my code : > > library(ggplot2) > library(plyr) > dat = rnorm(1000);

Re: [R] Updating DataFrame with Function

2009-07-24 Thread Ben Bolker
BostonR wrote: > > This should be very simple but I am wrestling with updating a data frame > from within a function. Here is a simple example: > ### SET UP DATA FRAME >> tFrame <- data.frame(T=c(1:5)) >> >> tFrame > T > 1 1 > 2 2 > 3 3 > 4 4 > 5 5 >> > # Simple function that updates w

Re: [R] rows missing after dataset loaded to R

2009-07-24 Thread Rnewbie
Many thanks to all of you :handshake: Gabor Grothendieck wrote: > > Another situation would be if you have comment characters in strings that > are intended to be content. > > On Thu, Jul 23, 2009 at 5:35 PM, Greg Snow wrote: >> Some programs quote everything to be "safe", others only quote whe

[R] Updating DataFrame with Function

2009-07-24 Thread BostonR
This should be very simple but I am wrestling with updating a data frame from within a function. Here is a simple example: ### SET UP DATA FRAME > tFrame <- data.frame(T=c(1:5)) > > tFrame T 1 1 2 2 3 3 4 4 5 5 > # Simple function that updates with 1000 > tUpdate<- function (x){ + tFrame$

Re: [R] legend title in qplot

2009-07-24 Thread RON70
Thank you for your reply. Yes I am reading that document and ofcourse simultaniously trying to apply in my current problem as well. Yes still I am learning. Here I have my code : library(ggplot2) library(plyr) dat = rnorm(1000); variable = rep(c("Variable:1", "Variable:2"), each=500); coll = rep

[R] binary digit

2009-07-24 Thread Угодай n/a
Does R has package for providing work for binary digit: arithmetic operation, convert to/from decimal digit, etc? I not found it, but think that CRAN contain it. [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://

Re: [R] Conditional sorting

2009-07-24 Thread Dimitris Rizopoulos
try this: d <- data.frame(RUN = rep(1:3, each = 3), ID = 1:9, AUC = runif(9,1,100)) d[order(d$RUN, d$AUC), ] I hope it helps. Best, Dimitris Riley, Steve wrote: Greetings! I am trying to figure out how to order a data frame by one variable conditioned on another. Here is an example of

[R] Conditional sorting

2009-07-24 Thread Riley, Steve
Greetings! I am trying to figure out how to order a data frame by one variable conditioned on another. Here is an example of what I have: d <- data.frame(RUN = rep(1:3, each = 3), ID = 1:9, AUC = runif(9,1,100)) > d RUN ID AUC 1170.2 1286.5 1320.1 24

Re: [R] downsampling

2009-07-24 Thread Philipp Pagel
On Fri, Jul 24, 2009 at 03:16:58AM -0600, Warren Young wrote: > Michael Knudsen wrote: > >On Fri, Jul 24, 2009 at 9:32 AM, Jan Wiener > >wrote: > > > >>x=sample(1:5, 115, replace=TRUE) > >> > >>How do I downsample this vector to 100 entries? Are there any R > >>functions or packages that provide s

Re: [R] Random # generator accuracy

2009-07-24 Thread Duncan Murdoch
On 23/07/2009 2:48 PM, Jim Bouldin wrote: Thanks Greg, that most definitely was it. So apparently the default is sampling without replacement. Fine, but this brings up a question I've had for a bit now, which is, how do you know what the default settings are for the arguments of any given funct

Re: [R] adjusting grid on Xaxis ticks

2009-07-24 Thread jim holtman
is this what you want: plot(0, xlim=as.POSIXct(c('2009-01-01', '2009-12-31')), xaxt='n') x <- seq(as.POSIXct('2009-01-01'), as.POSIXct('2009-12-13'), by='1 month') axis(1, at=x, labels=format(x, "%b")) abline(v=x, lty=2) On Thu, Jul 23, 2009 at 11:37 PM, wapita wapita wrote: > > Hello, > I have

Re: [R] SweaveOpts(eval=false) not working - CORRECTION

2009-07-24 Thread Friedrich Leisch
> On Thu, 23 Jul 2009 15:09:47 +0200, > Rainer M Krug (RMK) wrote: > On Thu, Jul 23, 2009 at 3:04 PM, Duncan Murdoch wrote: >> On 23/07/2009 8:58 AM, Rainer M Krug wrote: >>> >>> Hi >>> >>> I hope this is the right mailing list - if not, could you please refer >>> me to a m

Re: [R] ggplot question

2009-07-24 Thread ONKELINX, Thierry
Have a look at the colour argument of geom_line() http://had.co.nz/ggplot2/geom_line.html You need colour = A_Factor_With_Names_Of_Timeseries I recommend to you read the first chapters of Hadley's book on ggplot2 (it's on the ggplot2 website). That will answer much of your basic ggplot2 questions

Re: [R] computing the radius of an arc

2009-07-24 Thread Hans W Borchers
Nair, Murlidharan T iusb.edu> writes: > > Hi!! > > I am interesting in computing the radius of an arc that best approximates > a curve. Is there an R function that I can use to draw an arc? > Nothing useful came up when I searched help.search. Does anyone have any > suggestion to do this? >

  1   2   >