Re: [R] Count between interval

2010-02-22 Thread Bosken
Hi all, Thank you for the reactions! It works. Gr. Bosken -- View this message in context: http://n4.nabble.com/Count-between-interval-tp1564537p1564995.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing

Re: [R] how to make R plot under Linux

2010-02-22 Thread Vojtěch Zeisek
Dne Po 22. února 2010 19:33:45 Kevin E. Thorpe napsal(a): K. Elo wrote: Hi! 22.02.2010 19:53, xin wei wrote: hi, Kevin and K.Elo: thank you for the suggestion. Can you be more specific on these? (like how exactly get into x-switch or man ssh). I am totally ignorant about linux and

Re: [R] how to make R plot under Linux

2010-02-22 Thread Cedrick W. Johnson
I've managed to successfully use R (based on a remote Linux server) and have th e graphics piped back to me via SSH on a windows machine.. Take a look at XMing on the windows side, along with PuTTY. On 2/22/2010 2:41 PM, Vojtěch Zeisek wrote: Dne Po 22. února 2010 19:33:45 Kevin E. Thorpe

Re: [R] Extract information from S4 object

2010-02-22 Thread Martin Morgan
On 02/22/2010 08:46 AM, Jay wrote: Tahnk you. But, when I try the command you both suggested I get a NULL as the results. names(object1 @ x.values) NULL Where did I go wrong? On Feb 22, 4:34 pm, David Winsemius dwinsem...@comcast.net wrote: On Feb 22, 2010, at 8:05 AM, Jay wrote:

[R] Re-assigning variables stored as character strings in another variable

2010-02-22 Thread Anthony Damico
Is there any way to get the last line of this code to double the contents of a and b without naming them directly? #create variables a and b a-5 b-10 #store variable names a and b in variables c and d c-a d-b e-c(c,d) #loop through both variables for (i in e){ #print the numbers five and ten

Re: [R] gsub patterns from vector elements w/out loop?

2010-02-22 Thread Bert Gunter
I don't now whether or not this is less kludgy than the explicit loop (loops are implicitly used anyway), but Reduce(function(x1,x2)gsub(x2,something else,x1),x, init=y) does it using gsub. Note, however, that this will **not** work if the something else regular expression itself contained

[R] env() for lme4

2010-02-22 Thread Jerosch-Herold, Michael
I have been looking at a draft of Doug Bate's book on mixed effects models (http://lme4.r-forge.r-project.org/book), and in chapter one (section 1.4.2 Matrices and Vectors in the Fitted Model Object, to be exact) he refers to the function env() for checking the environment of a fitted lmer

Re: [R] Re-assigning variables stored as character strings in another variable

2010-02-22 Thread Henrique Dallazuanna
See ?get and ?assign On Mon, Feb 22, 2010 at 4:49 PM, Anthony Damico ajdam...@gmail.com wrote: Is there any way to get the last line of this code to double the contents of a and b without naming them directly? #create variables a and b a-5 b-10 #store variable names a and b in variables c

Re: [R] Re-assigning variables stored as character strings in anothervariable

2010-02-22 Thread William Dunlap
Bill Dunlap Spotfire, TIBCO Software wdunlap tibco.com -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Anthony Damico Sent: Monday, February 22, 2010 11:50 AM To: r-help@r-project.org Subject: [R] Re-assigning variables

[R] Siegel-Tukey test for equal variability (code)

2010-02-22 Thread Daniel Malter
Hi, I recently ran into the problem that I needed a Siegel-Tukey test for equal variability based on ranks. Maybe there is a package that has it implemented, but I could not find it. So I programmed an R function to do it. The Siegel-Tukey test requires to recode the ranks so that they express

Re: [R] env() for lme4

2010-02-22 Thread Marianne Promberger
Hi, I have been looking at a draft of Doug Bate's book on mixed effects models (http://lme4.r-forge.r-project.org/book), and in chapter one (section 1.4.2 Matrices and Vectors in the Fitted Model Object, to be exact) he refers to the function env() for checking the environment of a fitted

[R] relative file path

2010-02-22 Thread Rob Forler
Hello, Is there a way to find where a script is located within a script? getwd() doesn't do what I want because it depends on where R was called from. I want something like source(randomFile) and within randomFile there is a function called whereAmI() which returns c:\blah\blah2\randomFile.R In

Re: [R] how to make R plot under Linux

2010-02-22 Thread xin wei
hi, Guys: thank you so much for all the suggestion. Now I seem to be able to set up x11 forwarding in PUTTY. however, I still could not get plot and I get the following error msg: Error in function (display = , width, height, pointsize, gamma, bg, : X11 I/O error while opening X11 connection

[R] Creating regularly spaced time series from irregular one

2010-02-22 Thread Ivan Kalafatic
Hello, I have a series of intraday (high-frequency) price data in the form of POSIX timestamp followed by the value. I sucesfuly loaded that into its package object. I would like to create from it a regularly spaced time series of prices (for example 1min, 5min, etc apart) so i could calcualte

Re: [R] relative file path

2010-02-22 Thread Henrik Bengtsson
Not directly what you're asking for, but sourceTo() of R.utils has an argument 'chdir=FALSE' that when TRUE will change the working directory to the directory of the script being evaluated. /Henrik On Mon, Feb 22, 2010 at 9:44 PM, Rob Forler rfor...@uchicago.edu wrote: Hello, Is there a way

Re: [R] [R-SIG-Finance] Creating regularly spaced time series from irregular one

2010-02-22 Thread Jeff Ryan
Firstly, don't cross-post. Second, take a look at the archives on both these lists for answers to your questions. 'its' is rather old, and not where you want to be looking. Take a look at xts for fast time-series manipulation like you need, specifically to.period, endpoints, and align.time.

Re: [R] relative file path

2010-02-22 Thread Hrishi Mittal
There was a brief discussion about this on stackoverflow - http://stackoverflow.com/questions/1815606/rscript-determine-path-of-the-executing-script. - Try http://prettygraph.com Pretty Graph , the easiest way to make R-powered graphs on the web. -- View this message in context:

Re: [R] relative file path

2010-02-22 Thread Duncan Murdoch
On 22/02/2010 3:44 PM, Rob Forler wrote: Hello, Is there a way to find where a script is located within a script? getwd() doesn't do what I want because it depends on where R was called from. I want something like source(randomFile) and within randomFile there is a function called whereAmI()

Re: [R] (Somewhat) broken EPS files produced

2010-02-22 Thread Wartan Hachaturow
On Mon, Feb 22, 2010 at 4:10 PM, Marc Schwartz marc_schwa...@me.com wrote: It would be helpful to attach the .Rnw SWeave file (if small) or post it someplace that we can look at it, along with any data to replicate your process. Sure, you can grab both at

Re: [R] relative file path

2010-02-22 Thread William Dunlap
If you restructure your code as a bunch of functions calling functions instead of scripts sourcing scripts then the problem would largely go away. Put the functions into a package with a namespace and it becomes more reliable. Bill Dunlap Spotfire, TIBCO Software wdunlap tibco.com

Re: [R] Problem with installing genetics package

2010-02-22 Thread Juliet Hannah
I just installed it, and it worked fine. sessionInfo() R version 2.10.1 (2009-12-14) i386-pc-mingw32 locale: [1] LC_COLLATE=English_United States.1252 LC_CTYPE=English_United States.1252 [3] LC_MONETARY=English_United States.1252 LC_NUMERIC=C [5] LC_TIME=English_United States.1252 attached

[R] big panel: filehash, bigmemory or other

2010-02-22 Thread Eric Fail
Dear R-list I'm on my way to start a new project on a rather big panel, consisting of approximately 8 million observations in 30 waves of data and about 15 variables. I have a similar data set that is approximately 7 gigabytes in size. Until now I have done my data management in SAS, and

Re: [R] how to make R plot under Linux

2010-02-22 Thread Duncan Murdoch
On 22/02/2010 3:46 PM, xin wei wrote: hi, Guys: thank you so much for all the suggestion. Now I seem to be able to set up x11 forwarding in PUTTY. however, I still could not get plot and I get the following error msg: Error in function (display = , width, height, pointsize, gamma, bg, : X11

Re: [R] how to make R plot under Linux

2010-02-22 Thread Don MacQueen
The classic test for a properly set up X11 is to issue the command xclock at the prompt (on the linux box). (this assumes the linux box has xclock installed, but that is highly likely) If it works, the X11 forwarding is set up. If not, then, evidently not, but in the latter case at least

[R] tickmarks with time in dotplots

2010-02-22 Thread Mikkel Grum
Hi, I'm drawing lattice dotplots with time along the x-axis as in: C -data.frame(c(A, B),Sys.time()+ rnorm(50)*3600) names(C) - c(Name, Time) dotplot(Name ~ Time, data = C, horizontal = TRUE) On my display, the x-axis shows tick marks every two hours. I would like to show something more

Re: [R] relative file path

2010-02-22 Thread Rob Forler
I appreciate everyone's response. I will look at the options found in stack overflow. My issue is that i'm trying to load a file (.csv) relative to the R file. I pass the directory structure around but that directory structure may be put into something I'm unaware of. If packages can solve this

[R] Boxplots with similar number of classes....

2010-02-22 Thread kupz
I have four boxplots, stacked on top of one another (vertically). There are five classes that the continuous variable is set against. The problem I have is the data for the second boxplot only contains values for four of the five classes. I would like to include the space in this graph so there

Re: [R] how to make R plot under Linux

2010-02-22 Thread Saeed Abu Nimeh
Try to install xming in your windows box http://www.straightrunning.com/XmingNotes/. Make sure to run xming before plotting. Saeed On Mon, Feb 22, 2010 at 12:46 PM, xin wei xin...@stat.psu.edu wrote: hi, Guys: thank you so much for all the suggestion. Now I seem to be able to set up x11

[R] lmom: plotting log Pearson Type III

2010-02-22 Thread Hutchinson,David [PYR]
Can anyone show me how to add a log Pearson type III plot using the evdistq() command to an extreme value plot using the lmom package? Attached sample code below... Thanks in advance, Dave library(lmom) # annual maximum daily streamflows Mackenzie River mackenzieRiver = c(26600, 30300,

Re: [R] legend mathematical annotation problem

2010-02-22 Thread Gavin Simpson
On Sat, 2010-02-13 at 13:57 -0500, David Winsemius wrote: On Feb 13, 2010, at 1:45 PM, David Winsemius wrote: On Feb 13, 2010, at 1:13 PM, Mark Heckmann wrote: 1) I need to plot a legend containing the mathematical symbol greater-than-or-equal sign. And I want the text to start

Re: [R] NMDS ordination

2010-02-22 Thread Gavin Simpson
On Sat, 2010-02-13 at 11:33 -0800, Aisyah wrote: Hi Im currently trying to plot my NMDS data together with fitted variables (envfit funct) on an ordination plot. The plot function shows two displays=sites and sp. I was wondering how to plot it so that the sites come up as different points

Re: [R] tickmarks with time in dotplots

2010-02-22 Thread Phil Spector
Mikkel - Suppose you want tick marks at every half hour: ticks = seq(as.POSIXct('2010-02-22 12:00:00 PST'), as.POSIXct('2010-02-22 17:00:00 PST'),by='30 min') dotplot(Name ~ Time, data = C, horizontal = TRUE,scales=list(x=list(at=ticks))) Hope this helps.

Re: [R] tickmarks with time in dotplots

2010-02-22 Thread Hrishi Mittal
Thanks Phil, I was also trying to get at something like that but couldn't quite figure it out as quickly. Mikkel, I guess for your specific example you would need to define ticks like this: ticks = seq(min(C$Time),max(C$Time),by='30 min') - Try http://prettygraph.com Pretty Graph , the

Re: [R] how to make R plot under Linux

2010-02-22 Thread David Scott
Saeed Abu Nimeh wrote: Try to install xming in your windows box http://www.straightrunning.com/XmingNotes/. Make sure to run xming before plotting. Saeed On Mon, Feb 22, 2010 at 12:46 PM, xin wei xin...@stat.psu.edu wrote: hi, Guys: thank you so much for all the suggestion. Now I seem to be

Re: [R] Boxplots with similar number of classes....

2010-02-22 Thread David Winsemius
On Feb 22, 2010, at 5:11 PM, kupz wrote: I have four boxplots, stacked on top of one another (vertically). There are five classes that the continuous variable is set against. The problem I have is the data for the second boxplot only contains values for four of the five classes. I

Re: [R] [R-SIG-Finance] Creating regularly spaced time series from irregular one

2010-02-22 Thread Gabor Grothendieck
You probably want to ask the author of the its package directly about that. Note that as.zoo has an its method and as.its has a zoo method and they will allow you to convert back and forth between its and zoo so you can effectively use the functionality of both packages: # converting back and

[R] How to change a venn command into a named object that can be plotted like a lattice object

2010-02-22 Thread George Chen
Hello, I am plotting data as a venn diagram but would like to be able to control how it is plotted like a lattice object. Right now, it plots right away. I would like to name it and then plot at will. I thought to convert the whole thing to a PostScript file then get it back into R via

Re: [R] How to change a venn command into a named object that can be plotted like a lattice object

2010-02-22 Thread Hrishi Mittal
George, Unless, Venn Diagrams are produced as lattice objects, I don't think you can save them to modify or update later on. However, if you are just looking for a shortcut to avoid calling the plotting function again and again you could use the recordPlot() and replayPlot() functions. Add this

Re: [R] Lists into matrices within lists...again

2010-02-22 Thread David Winsemius
On Feb 22, 2010, at 8:27 PM, ewaters wrote: Related questions to this have been asked before, but I have tried all options they gave me unsuccessfully (do.call and unlist). I start with three lists of summary statistics, 100 elements each, which I bind together: None of this

[R] matching on two criteria

2010-02-22 Thread LCOG1
Howdy y'all, I am looking to use the match function to match a data attribute to another data frame but i need it to do so on two criteria to ensure an appropriate match. The following matches incorrectly. I know the example looks pointless but its just an example. Perhaps there is another

[R] (quite possibly OT) Re: how to make R running on a Linux server display a plot on a Windows machine

2010-02-22 Thread Emmanuel Charpentier
Dear Xin, Le lundi 22 février 2010 à 09:53 -0800, xin wei a écrit : hi, Kevin and K.Elo: thank you for the suggestion. Can you be more specific on these? (like how exactly get into x-switch or man ssh). I am totally ignorant about linux and SSH:( Memory limitation forces me to switch

Re: [R] how to make R plot under Linux

2010-02-22 Thread Zhuanshi He
HI, Please try to setup environment variable DISPLAY on your remote machine. ex, csh: setenv DISPLAY localhostIPAddress:20.0 or bash export DISPLAY=localhostIPAddress:20.0 where localhostIPAddress is your IP address of your current working machine (local machine) Usually, I just use

[R] Normal distribution (Lillie.test())

2010-02-22 Thread Bosken
Hi all, I have a dataset of 2000 numbers ( it's noise measured with a scoop ) Now i want to know of my data is normal distributed (Gaussian distribution). I did already: - 68-95-99.7 test - Q-Q-plot and now i used nortest library and the Lilli.test() However i don't understad the output?

[R] Use 2 ifelse to sort data

2010-02-22 Thread Chunhao
Dear R users, I have a question how to use 2 ifelse to sort my data. Such as from 11 to 20 assign to A; 6 to 10 assign to B, and the rest of them assign to C a-1:20 tt-ifelse(a10, A,no=ifelse( 5 a =10, B, C)) Many Thanks Chunhao -- View this message in context:

Re: [R] Use 2 ifelse to sort data

2010-02-22 Thread Bill.Venables
a - 1:20 tt - ifelse(a 10, A, ifelse(a 5, B, C)) The problem you have is with '5 a = 10'. Such double sided inequalities are used in mathematics but not in R. Here you only need the first part, but if you did need both you would need to write 5 a a = 10 Look carefully and spot the

Re: [R] Lists into matrices within lists...again

2010-02-22 Thread Peter Alspach
Tena koe Edward It is difficult to know the best approach from the information supplied, but using unlist() at the appropriate place in your code, maybe unlist(stats[[i]]$means), will probably get you where you want to be. HTH . Peter Alspach -Original Message- From:

Re: [R] Use 2 ifelse to sort data

2010-02-22 Thread David Winsemius
On Feb 22, 2010, at 8:13 PM, Chunhao wrote: Dear R users, I have a question how to use 2 ifelse to sort my data. Such as from 11 to 20 assign to A; 6 to 10 assign to B, and the rest of them assign to C a-1:20 tt-ifelse(a10, A,no=ifelse( 5 a =10, B, C)) Two way comparisons are not

[R] Short cut for documenting arguments in Rd file

2010-02-22 Thread Frank E Harrell Jr
The following construct in Rd files has always been accepted by R: \arguments{ \item{formula}{} \item{family}{} \item{data}{} \item{weights}{} \item{subset}{} \item{na.action}{} \item{start}{} \item{offset}{} \item{control}{} \item{model}{} \item{method}{} \item{x}{}

Re: [R] new to R Project

2010-02-22 Thread spencerg
Hi, Ivan: Are you trying to read standard text or csv files? If yes, then try ?read.table at a command prompt. Or are you trying to connect directly from R to a database system? If yes, which database system? This can be done from R, but as far as I know, finding the functions and

Re: [R] Goodness of fit test for count data

2010-02-22 Thread Moshe Olshansky
You can compute the conditional probability that your variable equals k given that it is non-zero. For example, if X has poisson distribution with parameter lambda then P(X=k/X!=0) = P(X=k)/(1-P(X=0)) = (exp(-lambda)/(1-exp(-lambda))*lambda^k/k! Now you can find lambda for which the sum of

Re: [R] new to R Project

2010-02-22 Thread spencerg
Hi, Ivan: p.s. If you'd like to generate reports using LaTeX, I suggest you also try (LaTeX - ???LaTeX). ### Are you trying to read standard text or csv files? If yes, then try ?read.table at a command prompt. Or are you trying to connect directly from R to

Re: [R] Short cut for documenting arguments in Rd file

2010-02-22 Thread RICHARD M. HEIBERGER
Frank, you can collect multiple arguments into a single \item{} You example simplified becomes \item{x, y, contrasts}{see \code{\link{glm}}} Rich [[alternative HTML version deleted]] __ R-help@r-project.org mailing list

Re: [R] Error with KNN

2010-02-22 Thread Nancy Adam
Hi Uwe, Thank you for your email. I need to compute one value indicating to the performance of knn but I don’t know anything other than MSE. Can you please tell me how I do that? Many thanks, Nancy Date: Mon, 22 Feb 2010 17:19:48 +0100 From: lig...@statistik.tu-dortmund.de To:

[R] Reports generation

2010-02-22 Thread chinna
Hi friends, Thanks for the answers.. now i am connecting to database but i dont know how to generate reports with that data can anyone please help me. if possible give the documents for report generation. thanks in advance. -- View this message in context:

Re: [R] Quadprog help

2010-02-22 Thread Moshe Olshansky
Hi Sergio, Having singular Dmat is certainly a problem. I can see two possibilities: 1) try to eliminate X1,...,X9, so that you are left with P1,...,P6 only. 2) if you can not do this, add eps*X1^+...+eps*X9^2 to your matrix Dmat so that it is positive definite (eps is a small positive number).

Re: [R] Fortune candidate! Re: new to R Project

2010-02-22 Thread spencerg
I have a concern with the sentiment expressed: When one is new to a language, it is often not easy to know where to start in the available documentation. It's hard to RTFM when one does not know which FM2R. I currently subscribe to email help for three other open source

Re: [R] matching on two criteria

2010-02-22 Thread jim holtman
Exactly what results are you expecting? Your 'match' has the same variables in both locations and will therefore return the first match: match(TazProperties..$Props,TazProperties..$Props) [1] 1 1 1 4 4 4 4 8 8 8 8 8 so you are getting the results that you ask for. On Mon, Feb 22, 2010 at 5:42

[R] Thanks Friends

2010-02-22 Thread chinna
Hi Bart, Thanks for the answer. connecting to a db, extracting the data you want, analyse this and put in a report, with graphs isn't that hard, but if you are starting out, you have to accept that you will have to learn the basics. I am learning R .after some reserch now i am

[R] generating reports from database through R

2010-02-22 Thread durga chennu
Hi, after some research now i am connecting to the database .but i am not getting any reports can u please tell me any suggestions or ideas. Regards chinna. [[alternative HTML version deleted]] __ R-help@r-project.org mailing list

Re: [R] Thanks Friends

2010-02-22 Thread Bart Joosen
Chinna, glad to hear that you are connected to your db. About your reports: we don't know what you are looking for, so how on earth can we help you. As already pointed out: List Requirements: PLEASE do read the posting guide [1]http://www.R-project.org/posting-guide.html and provide

Re: [R] forecasted Results

2010-02-22 Thread Dieter Menne
chinna wrote: i am connecting to the database and i am getting data then how can i get forecasted results of this using R project. A tourist in Berlin asks: How do I get to the Philharmony? Berliner: Pratice! Dieter -- View this message in context:

Re: [R] Alternatives to linear regression with multiple variables

2010-02-22 Thread Dieter Menne
Guy Green wrote: I wonder if someone can give some pointers on alternatives to linear regression (e.g. Loess) when dealing with multiple variables. For two variables, there is also interp.loess in package tcp. It can be rather slow depending on the parameters, so I fear a generalization

<    1   2