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 generaliza

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: http://n4.nabble.com/Gene

[R] forecasted Results

2010-02-22 Thread chinna
I have a table having the data as follows country Revenue us $68967 aus $60087 newz $78965 china $67846 i am connecting to the database and i am getting data then how can i get forecasted results of this using R project. -- View this message in contex

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 commente

[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 https:

[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 connectin

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

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 collabo

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). Y

[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: http://n4.nabble.com/new-to-R-P

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: nan

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 https://stat.ethz.ch/

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

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 square

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 an

[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}{} \it

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(a>10, "A",no=ifelse( 5< a <=10, "B", "C")) Two way comparisons a

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

2010-02-22 Thread Moshe Olshansky
Hi, As far as I understand, D is the value of (Kolmogorov-Smirnov) statistic and p-value is the probability to get that (or greater) value for normally distributed variables (so in your case you would most probably reject the hypothesis that your data is normal). --- On Tue, 23/2/10, Bosken w

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: r-help-bou

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

[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(a>10, "A",no=ifelse( 5< a <=10, "B", "C")) Many Thanks Chunhao -- View this message in context: http://n4.nabbl

[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?

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 localhost

[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 wa

[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 swi

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 represents

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 l

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

2010-02-22 Thread ewaters
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: statslist <- as.data.frame(cbind (means, vars, mcrs)) I then take 100 sampl

[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 grImport

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 for

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

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

[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, 34000,

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 wrote: hi, Guys: thank you so much for all the suggestion. Now I seem to be able to set up x11 for

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 ea

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 woul

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] 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 wrote: > > hi, Guys: > thank you so much for all the suggestion. Now I seem to be able to set up > x11 forwarding in PUTTY.

[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 i

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] 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 mo

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 yo

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, : X1

[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] 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 ba

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 > -Origi

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

2010-02-22 Thread Wartan Hachaturow
On Mon, Feb 22, 2010 at 4:10 PM, Marc Schwartz 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 http://wart.highsecure.ru/report.tar.gz (Be warned, th

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() w

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: http://n4.

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. The

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 wrote: > Hello, > > Is there a way to find where a scr

[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

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 connectio

[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] 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 fi

[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 vari

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 variab

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 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 and d >

[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 objec

[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

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 wrote: >> On Feb 22, 2010, at 8:05 AM, Jay wrote: >> >> >> >>>

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 n

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 abou

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

[R] useR! 2010, Submission/early registration deadline: March 1

2010-02-22 Thread Katharine Mullen
Dear R-help members, This message is to remind you that the submission deadline for abstracts for the R User Conference, useR! 2010, is one week away: Submission deadline: Monday, March 1, 2010 The deadline for early registration is also Monday, March 1, 2010. I encourage you all to su

[R] Goodness of fit test for count data

2010-02-22 Thread pinusan
Dear all, I am trying to test goodness of fit. I assume that a data follow Poisson or Negative binomial distribution. I can test the goodness of fit in case of no truncated data. However, I could not find any good function or packages when a data is truncated. For example, a frequency table

[R] First Appearance of value in splitted matrix

2010-02-22 Thread jorgusch
Hello, I have a large matrix (foo), I have to split by weeks and in which measures start at different points of time. Now I am looking for the number of the week, in which a specific value appears. I do not need the date or any other information from a row, but (due to structure of result table)

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

2010-02-22 Thread Zhuanshi He
Hi, You just need to install X11 windows software on windows, such as xming, and start xming before ssh -X use...@remotehost Hope this helps. On 2/22/10, xin wei wrote: > > > thank you for reply. I just type: hist(x) from SSH terminal, expecting a > histogram to pop up like what i got under

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

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

2010-02-22 Thread Gabor Grothendieck
Here are a few possibilities using gsubfn in the gsubfn package x <- c("one","two") y <- paste(rep(x,2),"blah") library(gsubfn) # 1 gsubfn("\\w+", w ~ if (w %in% x) "something else" else w, y) # 2 gsubfn("\\w+", list(one = "something else", two = "something else"), y) # 3 L <- sapply(x, functio

Re: [R] adding infrequent date labels to x-axis

2010-02-22 Thread emorway
Hello Hrishi, That worked great, and in the process I learned some new ways of going about writing R code. Thank you very much for helping me out! Eric -- View this message in context: http://n4.nabble.com/adding-infrequent-date-labels-to-x-axis-tp1564804p1564943.html Sent from the R help mai

Re: [R] adding infrequent date labels to x-axis

2010-02-22 Thread Hrishi Mittal
Hi emorway, You are right. I tried the command only with one year's data, in which case it plots month names as labels. But for multiple years, the labels become years. I found this old post useful - http://n4.nabble.com/Month-tick-marks-on-a-plot-td879121.html#a879121. Picking the code from the

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

2010-02-22 Thread Kevin E. Thorpe
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 SSH:( Memory limitation forces me to switch from windows to Lin

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

2010-02-22 Thread Liaw, Andy
You can try the locfit package, which I believe can handle up to 5 variables. E.g., R> library(locfit) Loading required package: akima Loading required package: lattice locfit 1.5-6 2010-01-20 R> x <- matrix(runif(1000 * 3), 1000, 3) R> y <- rnorm(1000) R> mydata <- data.frame(x, y) R> str(m

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

2010-02-22 Thread Vojtěch Zeisek
Dne Po 22. února 2010 18:53:55 xin wei napsal(a): > 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 Hello, in Linux You can plot as in Windows, for example using hist(x). Bu

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

2010-02-22 Thread K. Elo
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 > SSH:( Memory limitation forces me to switch from windows to Linux > c

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

2010-02-22 Thread Christian Raschke
Marianne, The function substring2 from the Hmisc packages works nicely for me to do this (without explicit gsub, though): x<- c("one","two") y<- paste(rep(x,2),"blah") y [1] "one blah" "two blah" "one blah" "two blah" substring2(y,x)<- "something else" y [1] "something else blah" "something

Re: [R] SEM error

2010-02-22 Thread John Fox
Dear Jarret, Uwe, and Dan, Sorry -- I missed the initial question. What's a bit odd here is that the singularity occurs only in the computation of the modification indices. It might help to look at the conditioning of the covariance matrix of the parameter estimates (i.e., the eigenvalues or si

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

2010-02-22 Thread xin wei
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 from windows to Linux cluster. Xin -- View this message in context:

Re: [R] adding infrequent date labels to x-axis

2010-02-22 Thread Hrishi Mittal
Hi emorway, It seems to me that all you need is one command: plot(as.Date(cropped.cast1$date,"%Y-%m-%d"),cropped.cast1$Frac_ET_Satsfd_mean) - Try http://prettygraph.com Pretty Graph , the easiest way to make R-powered graphs on the web. -- View this message in context: http://n4.nabble.c

Re: [R] adding infrequent date labels to x-axis

2010-02-22 Thread emorway
Hello Hrishi, The command you suggested plotted the years on the x-axis only. It was the same exact plot as the one I included in the original post. Respectfully, Eric -- View this message in context: http://n4.nabble.com/adding-infrequent-date-labels-to-x-axis-tp1564804p1564875.html Sent fr

[R] Bootstrap Multivariate Times Series Forecast

2010-02-22 Thread Ricardo Gonçalves Silva
Dear Users, Consider a multivariate time series model: a_1*y(t)-...-a_k*y(t-k)=b+[c_1*z(t)-...-c_j*z(t-j)] i.e., a simple multivariate time series model with one exogenous variable. I would like to know what package can I use to do the following, using R: 1) Select k and j jointly; 2) Estimate t

Re: [R] Re : how to plot select points in preexisting persp plot

2010-02-22 Thread Hrishi Mittal
Preeti, Could you please post a sample of your data and all the code? It is hard to answer your question without knowing what "coords" and "act" are - Try http://prettygraph.com Pretty Graph , the easiest way to make R-powered graphs on the web. -- View this message in context: http://n4

Re: [R] [BioC] Empirical Pvalue

2010-02-22 Thread Hervé Pagès
Hi Mohamed, Please read the posting guides for both, the R-help and BioC mailing lists, so you learn the basics of how to ask questions on these lists (and software related mailing lists in general). Following those guides will greatly increase the chance that someone will actually be willing (an

Re: [R] counting repeating sequence lengths in a vector

2010-02-22 Thread Henrique Dallazuanna
Try this: rle(x)$length On Mon, Feb 22, 2010 at 2:27 PM, Larson, TR wrote: > Hello, > I have a very long (~50,000) sequence of repeating numbers.  The first 100 > are: > > [1]     0     0     0     0     0     0     0     0     0     0     0   429 >  [13]   429   429   429   429   429   429   4

[R] counting repeating sequence lengths in a vector

2010-02-22 Thread Larson, TR
Hello, I have a very long (~50,000) sequence of repeating numbers. The first 100 are: [1] 0 0 0 0 0 0 0 0 0 0 0 429 [13] 429 429 429 429 429 429 429 858 858 858 858 858 [25] 858 1287 1287 1287 1287 1287 1716 2145

Re: [R] Extract information from S4 object

2010-02-22 Thread Jay
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 wrote: > On Feb 22, 2010, at 8:05 AM, Jay wrote: > > > > > The function prediction() returns this: > > > Formal clas

Re: [R] Accessing values of a matrix

2010-02-22 Thread statquant
OKKK Thanks a lot for letting me know the subset function Cheers -- View this message in context: http://n4.nabble.com/Accessing-values-of-a-matrix-tp1561932p1564724.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.or

[R] Re : how to plot select points in preexisting persp plot

2010-02-22 Thread Preeti Iyer
Hello, I am trying to plot selected data points to a preexisting persp plot that satisfies a condition. I used the following statement - text(coords[,1], coords[,2], names(act[which(act > 8.75)]), cex=0.7) But I get all the points labeled, instead of the points that satisfy the condition speci

[R] adding infrequent date labels to x-axis

2010-02-22 Thread emorway
I'm sure there is a clever way to do the following, but I've been unable to find it on this forum or by writing my own functions. I have 8 years worth of weekly data but would like to restrict the labels on the x-axis to months only. I've included the first year's worth of data below. My line

Re: [R] SEM error

2010-02-22 Thread Jarrett Byrnes
I have often found this to happen if the scale of one variable is orders of magnitude different than the scale of other variables. Have you tried inspecting the covariance matrix and log transforming any such variables? On Feb 22, 2010, at 8:14 AM, Uwe Ligges wrote: On 20.02.2010 08:51

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

2010-02-22 Thread Marianne Promberger
> > gsub(paste(x, collapse = "|"), "something else", y) > [1] "something else blah" "something else blah" "something else blah" > [4] "something else blah" Many thanks! I didn't know about "collapse". Should have thought about reading up in ?paste ... Thanks Marianne -- Marianne Promberger P

Re: [R] Count between interval

2010-02-22 Thread Peter Ehlers
And as long as the interval is symmetric about zero: sum(abs(x) <= 2) [1] 6 -Peter Ehlers On 2010-02-22 9:18, Jorge Ivan Velez wrote: Hi, Here is a suggeston: x<- c(-1.3, 1, -1.5, -1, 1.5, -2.5, 3, -0.5) sum(x>=-2& x<=2) [1] 6 HTH, Jorge On Mon, Feb 22, 2010 at 9:49 AM, Bosken<> wro

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

2010-02-22 Thread Marc Schwartz
On Feb 22, 2010, at 10:25 AM, Marianne Promberger wrote: > Dear list, > > I have two vectors: > > x <- c("one","two") > y <- paste(rep(x,2),"blah") > > I want to replace all occurrences of each element of x in y with > something else, so that y looks like this: > > y > [1] "something else bla

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

2010-02-22 Thread Benilton Carvalho
gsub("one|two", "something else", y) ? On Mon, Feb 22, 2010 at 4:25 PM, Marianne Promberger wrote: > Dear list, > > I have two vectors: > > x <- c("one","two") > y <- paste(rep(x,2),"blah") > > I want to replace all occurrences of each element of x in y with > something else, so that y looks lik

Re: [R] replicate matrix

2010-02-22 Thread William Dunlap
> -Original Message- > From: r-help-boun...@r-project.org > [mailto:r-help-boun...@r-project.org] On Behalf Of wendy > Sent: Saturday, February 20, 2010 6:58 PM > To: r-help@r-project.org > Subject: [R] replicate matrix > > > Hi all, > > I have a matrix, for example > [,1] [,2] > [

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

2010-02-22 Thread Marianne Promberger
Dear list, I have two vectors: x <- c("one","two") y <- paste(rep(x,2),"blah") I want to replace all occurrences of each element of x in y with something else, so that y looks like this: y [1] "something else blah" "something else blah" "something else blah" [4] "something else blah" I can do

Re: [R] Error with KNN

2010-02-22 Thread Uwe Ligges
On 20.02.2010 19:10, Nancy Adam wrote: Hi all I receive this error message: “Error in mymodel$MSE : $ operator is invalid for atomic vectors” When I use KNN. This is the code: library(e1071) train<- rbind(iris3[,,1], iris3[,,2], iris3[,,3]) cl<- factor(c(rep("s",50), rep("c",50), rep("v"

Re: [R] change email subscription

2010-02-22 Thread Uwe Ligges
See the footer of your message. Uwe Ligges On 22.02.2010 13:27, Wendy Qiao wrote: Dear R team, As this is a university email address with very small inbox limit, can I please change my R email subscription to another email address instead of this one? My new email address is wendy2.q...@gma

Re: [R] Count between interval

2010-02-22 Thread Jorge Ivan Velez
Hi, Here is a suggeston: > x <- c(-1.3, 1, -1.5, -1, 1.5, -2.5, 3, -0.5) > sum(x >=-2 & x<=2) [1] 6 HTH, Jorge On Mon, Feb 22, 2010 at 9:49 AM, Bosken <> wrote: > > Hi all, > > I take a simple dataset like this: > > -1.3 > 1 > -1.5 > -1 > 1.5 > -2.5 > 3 > -0.5 > > Now I want to count how man

  1   2   >