[R] About available datasets on PC

2010-11-23 Thread Stephen Liu
Hi folks, Win7 On running; data(package = .packages(all.available = TRUE)) it displays a list of datasets under:- Data sets in package ‘AER’: But I couldn't call/load all of them on the list. DJFranses Error: object 'DJFranses' not found CreditCard Error: object 'CreditCard' not found

[R] More detail in chart axis?

2010-11-23 Thread Noah Silverman
Hi, I have a series of data (about 80,000 pairs of x,y). Plotting it shows a great chart. However, R has randomly chosen about 6 labels for my x axis. Now, clearly I can't show them all, but would like some control over the granularity of what is displayed. I can't find anything in the

[R] Error: cannot allocate vector of size x Gb (64-bit ... yet again)

2010-11-23 Thread derek eder
Hello, I am facing the dreaded Error: cannot allocate vector of size x Gb and don't understand enough about R (or operating system) memory management to diagnose and solve the problem -- despite studying previous posts and relevant R help -- e.g.: Error messages beginning cannot allocate

Re: [R] About available datasets on PC

2010-11-23 Thread Jeff Newmiller
You need to load the package and then load the data: library(AER) data(DJFranses) then it will be available for you to work with. Stephen Liu wrote: Hi folks, Win7 On running; data(package = .packages(all.available = TRUE)) it displays a list of datasets under:- Data sets in package

Re: [R] More detail in chart axis?

2010-11-23 Thread Jim Lemon
On 11/23/2010 07:24 PM, Noah Silverman wrote: Hi, I have a series of data (about 80,000 pairs of x,y). Plotting it shows a great chart. However, R has randomly chosen about 6 labels for my x axis. Now, clearly I can't show them all, but would like some control over the granularity of what is

Re: [R] how to get rid of unused space on all 4 borders in plot() render

2010-11-23 Thread Petr PIKAL
Hi r-help-boun...@r-project.org napsal dne 21.11.2010 19:31:33: I have looked into par documentation, and only setting for size of the plot area was pin. But this setting sets the area as inflexible, that is no R intro manual Chapter 12: Graphical procedures page, 71 A typical figure

Re: [R] how to loop through variables in R?

2010-11-23 Thread Ivan Calandra
Hi! You haven't got any answer probably because you didn't provide a reproducible example. You can do it by copy/pasting the output of dput(d) or dput(df). Moreover, before this email, I couldn't really understand what you were trying to do. It's not crystal clear now, but I think I got it.

[R] Odp: save a regression model that can be used later

2010-11-23 Thread Petr PIKAL
Hi r-help-boun...@r-project.org napsal dne 22.11.2010 16:02:20: Hi everyone I have a question about how to save a regression model in R and how to retrieve it for making predictions in a new session. To be more specific, I fitted a multilevel logistic regression model using the lmer

Re: [R] Lost in POSIX

2010-11-23 Thread Jeff Newmiller
Dimitri Shvorob wrote: Nor would I call this much of an improvement in clarity... what about min? You want to know the minimum? LOL. (And apologies for the insensitivity). Thank you for help, Jeff. This works, but I am still curious to see a solution based on trunc, if anyone can find it.

Re: [R] About available datasets on PC

2010-11-23 Thread Stephen Liu
Hi Jeff, Tks for your advice. I got it B.R. satimis - Original Message From: Jeff Newmiller jdnew...@dcn.davis.ca.us To: Stephen Liu sati...@yahoo.com Cc: r-help@r-project.org Sent: Tue, November 23, 2010 4:41:57 PM Subject: Re: [R] About available datasets on PC You need to load

Re: [R] Calculating correlation

2010-11-23 Thread Tal Galili
Hi Gireesh, I seem to be missing something. Could you give me the formula you think should be used to calculate this? I don't get how you wish to get a correlation for only one pair of numbers. (or maybe I didn't understood what you explained - please try a simpler example)

[R] Reading parts of data files

2010-11-23 Thread fbielejec
Dear, I'm doing analysis where I need to work on relatively large (50-60 MB) text files, though I'm really interested only in parts with binary variables (named indicators1, indicators2, ... etc.) Every text file contains other numeric columns, but not always the same and not always in the same

Re: [R] Calculating correlation

2010-11-23 Thread gireesh bogu
in one output x - S1 i.e. 112 with all combinations ex: 112 vs 3 [ x-S1 vs a -S1] 112 vs 4 [ x-S1 vs a- S2 ] 112 vs 45 112 vs 34 112 vs 23 112 vs 112 …. ….. 112 vs 1 in second output x - S2 i.e. 0 with all . ex: 0 vs 3 0 vs 4 ……… ……… 0 vs 1 in next output x-S3

Re: [R] Error: cannot allocate vector of size x Gb (64-bit ... yet again)

2010-11-23 Thread Uwe Ligges
On 23.11.2010 09:26, derek eder wrote: Hello, I am facing the dreaded Error: cannot allocate vector of size x Gb and don't understand enough about R (or operating system) memory management to diagnose and solve the problem -- despite studying previous posts and relevant R help -- e.g.: Error

Re: [R] How to call web service in R

2010-11-23 Thread 夏高
Actually RCurl is not for calling web service. But thank you all the same. Does anybody know how to call web service in R? Thanks! 2010/11/23 Steve Lianoglou mailinglist.honey...@gmail.com Hi, On Mon, Nov 22, 2010 at 4:51 AM, ÏÄ¸ß xiagao1...@gmail.com wrote: Hello everyone, I would like

Re: [R] Calculating correlation

2010-11-23 Thread Liviu Andronic
Hello On Tue, Nov 23, 2010 at 7:57 AM, gireesh bogu girishb...@gmail.com wrote: Hi guys I have an input file with multiple columns and and rows. Is it possible to calculate correlation of certain value of certain No (For example x of S1 = 112) with all other values (for example start with x

[R] extracted residuals as single column

2010-11-23 Thread Jennifer Dodd
Dear list, I am looking for some help with the format of extracted residual values output from: residuals (model1) [1] 0.74140534 0.72464895 0.74873373 0.78313611 0.78836470 0.76934676 0.77448847 0.78465077 0.79278282 [10] 0.80852056 0.81442047 0.80890838 0.78335232 0.83127343

Re: [R] How to call web service in R

2010-11-23 Thread Barry Rowlingson
On Tue, Nov 23, 2010 at 8:09 AM, 夏高 xiagao1...@gmail.com wrote: Actually RCurl is not for calling web service. But thank you all the same. Does anybody know how to call web service in R? Thanks! Use RCurl. It can be used for calling web services. A web service is just an HTTP endpoint that

[R] help

2010-11-23 Thread Francial Giscard LIBENGUE
S'il vous plait, j'ai un problème avec la fonction gamma je n'arrive pas à bienécrire l'exposant j'ai utilisé ^ et ça ne marche pas. cordialemont -- Francial Giscard LIBENGUE Doctorant en Mathématiques Appliquées ;Option : Statistique Université de Franche-Comté - UFR Sciences et Techniques

Re: [R] extracted residuals as single column

2010-11-23 Thread Barry Rowlingson
On Tue, Nov 23, 2010 at 10:36 AM, Jennifer Dodd j.d...@bio.gla.ac.uk wrote: Dear list, I am looking for some help with the format of extracted residual values output from: residuals (model1) [1]  0.74140534 0.72464895 0.74873373 0.78313611 0.78836470 0.76934676 0.77448847

Re: [R] How to call web service in R

2010-11-23 Thread Mike Marchywka
Date: Tue, 23 Nov 2010 11:28:15 + From: b.rowling...@lancaster.ac.uk To: xiagao1...@gmail.com CC: r-help@r-project.org Subject: Re: [R] How to call web service in R On Tue, Nov 23, 2010 at 8:09 AM, 夏高 wrote: Actually RCurl is not for

[R] format y-axis values in ggplot

2010-11-23 Thread Jack Johnson
Hi, I'd like to change scientific number formatting, 1,E+06, in y-axis values to 1 000 000 (number with 1000 separator). How to do this in ggplot..? I tried to use scale_discrete but couldn't get it to work. Below is a code I have this far: p-ggplot(Data,aes(x=Date,y=PRINCIPAL)) p+geom_point()

Re: [R] Lattice and Quartz

2010-11-23 Thread Dennis Murphy
Perhaps you're a victim of FAQ 7.22. Try putting a print () around your plot call - e.g., print(histogram(~ dollars)). And please read said FAQ to understand why... HTH, Dennis On Mon, Nov 22, 2010 at 10:11 PM, Sarah Berry escaloso...@gmail.com wrote: I ran this script in a source file on my

[R] cross validation using e1071:SVM

2010-11-23 Thread Neeti
Hi everyone I am trying to do cross validation (10 fold CV) by using e1071:svm method. I know that there is an option (“cross”) for cross validation but still I wanted to make a function to Generate cross-validation indices using pls: cvsegments method.

Re: [R] cross validation using e1071:SVM

2010-11-23 Thread Francial Giscard LIBENGUE
Hi everyone, Can you help me to plot Gamma(x/h+1) and Beta(x/h+1,(1-x)/h+1)?I want write x-seq(0,3,0.1) thank 2010/11/23 Neeti nikkiha...@gmail.com Hi everyone I am trying to do cross validation (10 fold CV) by using e1071:svm method. I know that there is an option (“cross”) for cross

Re: [R] format y-axis values in ggplot

2010-11-23 Thread Dennis Murphy
Hi: Borrowing an idea from a post by Baptiste Auguie a few weeks ago on the ggplot2 list, try this: http://groups.google.com/group/ggplot2/browse_thread/thread/46a78ab3315fdb7/f2ad305ae935a862?lnk=gstq=formatter#f2ad305ae935a862 mysep - function(x, ...) format(x, big.mark = ' ', trim = TRUE,

[R] read.table / type.convert: Unexpected column class

2010-11-23 Thread Andreas Borg
Hi all, I stumbled across the following, in my opinion surprising, behaviour of read.table: data - read.table(stdin(), sep=,) 0: 1, 2, 3 # Some comment 1: sapply(data, class) V1V2V3 integer integer numeric I was surprised to find that V3 has class numeric despite its

Re: [R] Reading parts of data files

2010-11-23 Thread jim holtman
?file - how to use connections ?read.table'skip' parameter, colClasses to only read columns you want That is not a large file. Read the whole thing in and then extract the data you need. On Tue, Nov 23, 2010 at 6:05 AM, fbielejec fbiele...@gmail.com wrote: Dear, I'm doing analysis

[R] Catching a RCurl error?

2010-11-23 Thread Tal Galili
Hi all, I'm running a complex script which accesses the internet, and sometimes it stops with the error: Error in curlPerform(url = url, headerfunction = header$update, curl = curl, : Failure when receiving data from the peer Is there a way to make the script wait longer, or not crash

[R] the first. from SAS in R

2010-11-23 Thread Joel
Is there any similar function in R to the first. in SAS? What it dose is: Lets say we have this table: a b c 1 1 5 1 0 2 2 0 2 2 0 NA 2 9 2 3 1 3 and then I want do to do one thing the first time the number 1 appers in a and something else the secund time 1 appers in a

Re: [R] Sweave Dynamic Graph Question

2010-11-23 Thread cameron
thx guys, this is what i came up with. it works fine, but too many lines of code. thanks again. Cameron \pagebreak \section{General Value Graph Daily} \setkeys{Gin}{width=1.1\textwidth} results=tex,echo=FALSE= rng=range(time(pf4[[1]]))

[R] Three-way ANOVA shows me two-way results

2010-11-23 Thread Marius 't Hart
Hi all, I'm doing a 3-way ANOVA like this: summary(aov(formula('FP ~ (lum * obj * man)^3 - Error(vp/(lum * obj * man)^3)'),data=dataf)) But in the output I only get 1- and 2-way effects, like this one: Error: vp:obj:man Df Sum Sq Mean Sq F value Pr(F) obj:man 1 1.5291e-34 1.5291e-34 5.7011

Re: [R] Lattice and Quartz

2010-11-23 Thread David Winsemius
On Nov 23, 2010, at 1:11 AM, Sarah Berry wrote: I ran this script in a source file on my Mac: library(lattice) year - 1900:2000 dollars - (year-1899)^2 plot(year,dollars) quartz(year) histogram(~dollars) The first plot appears in Quartz 2. The second quartz window, named year, opens but

Re: [R] Three-way ANOVA shows me two-way results

2010-11-23 Thread Joel
Marius think you wanted to do a new post and not reply to mine :P -- View this message in context: http://r.789695.n4.nabble.com/the-first-from-SAS-in-R-tp3055417p3055459.html Sent from the R help mailing list archive at Nabble.com. __

Re: [R] format y-axis values in ggplot

2010-11-23 Thread Jack Johnson
Dennis, Thanks for your answer. That is working nicely! Br, Jack 2010/11/23 Dennis Murphy djmu...@gmail.com Hi: Borrowing an idea from a post by Baptiste Auguie a few weeks ago on the ggplot2 list, try this:

Re: [R] the first. from SAS in R

2010-11-23 Thread Gustavo Carvalho
Perhaps something like this: a$d - ifelse(duplicated(a$a), 0, 1) On Tue, Nov 23, 2010 at 1:33 PM, Joel joda2...@student.uu.se wrote: Is there any similar function in R to the first. in SAS? What it dose is: Lets say we have this table:  a b  c  1 1  5  1 0  2  2 0  2  2 0 NA  2 9  2

Re: [R] the first. from SAS in R

2010-11-23 Thread David Winsemius
On Nov 23, 2010, at 8:33 AM, Joel wrote: Is there any similar function in R to the first. in SAS? What it dose is: Lets say we have this table: a b c 1 1 5 1 0 2 2 0 2 2 0 NA 2 9 2 3 1 3 and then I want do to do one thing the first time the number 1 appers in a and

Re: [R] compare GLM coefficients

2010-11-23 Thread Ben Bolker
Michael Bedward michael.bedward at gmail.com writes: Hello Kayce, My (very basic) understanding is that you can't directly compare the coefficients across models that have different response variables, nor could you use AIC and similar metrics of model goodness of fit. Instead, I think

[R] RDA Triplot

2010-11-23 Thread Danielwc
Hi Im using the VEGAN package to do an RDA ordination. In my plot I get my environmental scores as arrows/vectors, but my species scores as points. I would like to get the species scores as arrows as well. I tried the envfit function, and this creates arrows, but the lenght is not proportionate

Re: [R] question about constraint minimization

2010-11-23 Thread Ravi Varadhan
You need to spend some time and learn the basics of R. There are several excellent tutorials available on the CRAN R-project website. This will be a very good investment of your time. After you have spent a couple of days learning the basics, try using my package to solve your problem. Best,

Re: [R] Reading parts of data files

2010-11-23 Thread Gabor Grothendieck
On Tue, Nov 23, 2010 at 6:05 AM, fbielejec fbiele...@gmail.com wrote: Dear, I'm doing analysis where I need to work on relatively large (50-60 MB) text files, though I'm really interested only in parts with binary variables (named indicators1, indicators2, ... etc.) Every text file contains

[R] Three-way ANOVA shows me two-way results

2010-11-23 Thread Marius 't Hart
Hi all, I'm doing a 3-way ANOVA like this: summary(aov(formula('FP ~ (lum * obj * man)3 - Error(vp/(lum * obj * man)3)'),data=dataf)) But in the output I only get 1- and 2-way effects, like this one: Error: vp:obj:man Df Sum Sq Mean Sq F value Pr(F) obj:man 1 1.5291e-34 1.5291e-34 5.7011

Re: [R] help

2010-11-23 Thread Spencer Graves
La langue de cette liste ci est desolais anglais. S'il vous plait, nous devon un example. J'ai trouve: 2^3 [1] 8 PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. Hope this helps. Spencer On

Re: [R] Sporadic errors when training models using CARET

2010-11-23 Thread Max Kuhn
Kendric, I've seen these too and traceback() usually goes back to ksvm(). This doesn't mean that the error is there, but the results fo traceback() from you would be helpful. thanks, Max On Mon, Nov 22, 2010 at 6:18 PM, Kendric Wang kendr...@interchange.ubc.ca wrote: Hi. I am trying to

[R] Executing multiple .R files

2010-11-23 Thread Santosh Srinivas
Hello R-Helpers, I have a directory with some .R files that I execute every day. I want to write a script that executes each one of time sequentially. Is there a statement for this? Thank you. __ R-help@r-project.org mailing list

Re: [R] how to add frequencies to barplot

2010-11-23 Thread Juliet Hannah
Also check out the following post: http://permalink.gmane.org/gmane.comp.lang.r.general/210897 On Sat, Nov 20, 2010 at 4:32 PM, casperyc caspe...@hotmail.co.uk wrote: Hi, I have count data x2=rep(c(0:3),c(13,80,60,27)) x2  0  1  2  3 13 80 60 27 I want to graph to be ploted as

Re: [R] quot;helpquot;

2010-11-23 Thread Ben Bolker
Francial Giscard LIBENGUE libengue at gmail.com writes: S'il vous plait, j'ai un problème avec la fonction gamma je n'arrive pas à bienécrire l'exposant j'ai utilisé ^ et ça ne marche pas. cordialemont ici on parle anglais (je ne sais pas s'il y a une autre liste ou on pourrait poser

[R] Searching dataset

2010-11-23 Thread Stephen Liu
Hi folks, help.search(item_name) displays a list of dataset containing the item_name. If I need to search containing more than one item_name. What will be the command TIA B.R. Stephen L __ R-help@r-project.org mailing list

Re: [R] the first. from SAS in R

2010-11-23 Thread Dennis Murphy
Interesting. Check this out: u - sample(c(TRUE, FALSE), 10, replace = TRUE) u [1] FALSE FALSE TRUE FALSE FALSE TRUE FALSE FALSE FALSE FALSE class(u) [1] logical u + 0 [1] 0 0 1 0 0 1 0 0 0 0 0 + u [1] 0 0 1 0 0 1 0 0 0 0 v - rpois(10, 3) !duplicated(v) [1] TRUE FALSE TRUE TRUE

Re: [R] the first. from SAS in R

2010-11-23 Thread Charles C. Berry
On Tue, 23 Nov 2010, Joel wrote: Is there any similar function in R to the first. in SAS? What it dose is: Lets say we have this table: a b c 1 1 5 1 0 2 2 0 2 2 0 NA 2 9 2 3 1 3 and then I want do to do one thing the first time the number 1 appers in a and something else the

Re: [R] Executing multiple .R files

2010-11-23 Thread jim holtman
?source On Tue, Nov 23, 2010 at 10:04 AM, Santosh Srinivas santosh.srini...@gmail.com wrote: Hello R-Helpers, I have a directory with some .R files that I execute every day. I want to write a script that executes each one of time sequentially. Is there a statement for this? Thank you.

[R] Slow update(insert) on Data.frame

2010-11-23 Thread Joel
Hi When I try to update an number in a large data.frame by its pos It's really slow it take almost a sec to do this and I wonder why and if where is any faster way to update a number in a data.frame ive tried DF$col[POS]-number DF[xPOS,yPOS]-number Thx //Joel -- View this message in

[R] Finding the rows with the duplicated index

2010-11-23 Thread Santosh Srinivas
Hello Group, I have a huge time series dataset with sample below. I am basically trying to read it into a zoo object with columns 1:6 to index. Zoo issues a warning that some of the rows have duplicated index. dput(z) structure(list(TrdTimestamp = structure(list(sec = c(19, 19, 18, 10, 12, 43,

Re: [R] the first. from SAS in R

2010-11-23 Thread Charles C. Berry
On Tue, 23 Nov 2010, Dennis Murphy wrote: Interesting. Check this out: u - sample(c(TRUE, FALSE), 10, replace = TRUE) u [1] FALSE FALSE TRUE FALSE FALSE TRUE FALSE FALSE FALSE FALSE class(u) [1] logical u + 0 [1] 0 0 1 0 0 1 0 0 0 0 0 + u [1] 0 0 1 0 0 1 0 0 0 0 v - rpois(10, 3)

Re: [R] Three-way ANOVA shows me two-way results

2010-11-23 Thread Tal Galili
It seems your model defines more parameters then your data permits (e.g: singular). What you can do is try something like this: summary(aov(formula('FP ~ (lum + obj + man)^2 - Error(vp/(lum + obj + man)^2)),data=dataf)) (notice the use of + and ^) This is something your data might permit

Re: [R] the first. from SAS in R

2010-11-23 Thread David L Lorenz
It all has to do with the precedence of the ! operator. Compare !duplicated(v) + 0 with (!duplicated(v)) + 0 Dave From: Charles C. Berry cbe...@tajo.ucsd.edu To: Dennis Murphy djmu...@gmail.com Cc: r-help@r-project.org Date: 11/23/2010 10:08 AM Subject: Re: [R] the first. from SAS in R Sent

Re: [R] Slow update(insert) on Data.frame

2010-11-23 Thread Charles C. Berry
On Tue, 23 Nov 2010, Joel wrote: Hi When I try to update an number in a large data.frame by its pos It's really slow it take almost a sec to do this and I wonder why and if where is any faster way to update a number in a data.frame ive tried DF$col[POS]-number DF[xPOS,yPOS]-number See

Re: [R] Finding the rows with the duplicated index

2010-11-23 Thread Gabor Grothendieck
On Tue, Nov 23, 2010 at 11:04 AM, Santosh Srinivas santosh.srini...@gmail.com wrote: Hello Group, I have a huge time series dataset with sample below. I am basically trying to read it into a zoo object with columns 1:6 to index. Zoo issues a warning that some of the rows have duplicated

Re: [R] the first. from SAS in R

2010-11-23 Thread David Winsemius
On Nov 23, 2010, at 11:04 AM, Charles C. Berry wrote: On Tue, 23 Nov 2010, Dennis Murphy wrote: Interesting. Check this out: u - sample(c(TRUE, FALSE), 10, replace = TRUE) u [1] FALSE FALSE TRUE FALSE FALSE TRUE FALSE FALSE FALSE FALSE class(u) [1] logical u + 0 [1] 0 0 1 0 0 1 0 0 0

[R] How to bold items on a data frame?

2010-11-23 Thread Stephen Liu
Hi folks, I ran following commands to create a data file: Year - c(1983, 1984, 1985, 1986, 1987, 1988, 1989, 1990) Yield - c(60, 50, 70, 70, 80, 50, 60, 40) Concentration% - c(30, 40, 50, 60, 70, 80, 90, 100) Temperature - c(56, 47, 53, 53, 56, 47, 44, 44) datafile -

Re: [R] question about constraint minimization

2010-11-23 Thread Hans W Borchers
dhacademic at gmail.com dhacademic at gmail.com writes: Hi, I have struggled on this bound optimization with equality constraint by using optim function for two days, but still fail to prepare a good input. Can anyone help to prepare the input for my specific case? Many thanks. Best,

Re: [R] How to bold items on a data frame?

2010-11-23 Thread Joshua Wiley
Hi Stephen, bold is presentation type idea, you can do it in reports (e.g., via Sweave in LaTeX with \textbf{} or \mathbf{}), or graphs from R (e.g., make your title bold: expression(bold(Your Title))), but it is not something you do in R itself. Some front ends to R or text editors may provide

Re: [R] Finding the rows with the duplicated index

2010-11-23 Thread Santosh Srinivas
Thanks. Sorry about the naming (just cobbled the example). The rows are duplicate in the actual data. I am looking for a way to find the number of duplicate rows for each unique rowindex combination which is based on columns 1:6. -Original Message- From: Gabor Grothendieck

Re: [R] Finding the rows with the duplicated index

2010-11-23 Thread Gabor Grothendieck
On Tue, Nov 23, 2010 at 11:49 AM, Santosh Srinivas santosh.srini...@gmail.com wrote: Thanks. Sorry about the naming (just cobbled the example). The rows are duplicate in the actual data. I am looking for a way to find the number of duplicate rows for each unique rowindex combination which is

[R] Question about list function

2010-11-23 Thread Guido Leoni
Dear List I'm a newbie R user. I'm utilizing the list function in order to make a var like this: clusters-list(a=var1,b=var2) My problem is that the total numer of variables that I need to include in my list is up to 200. I've the text string with the complete list of my variables but is too

[R] problem with parameters

2010-11-23 Thread alkalinetrio
Hi, I just downloaded R, but I can't create a matrix with a parameter p. Could someone help me please? Andrea -- View this message in context: http://r.789695.n4.nabble.com/problem-with-parameters-tp3055682p3055682.html Sent from the R help mailing list archive at Nabble.com.

Re: [R] Three-way ANOVA shows me two-way results

2010-11-23 Thread Mike Marchywka
From: tal.gal...@gmail.com Date: Tue, 23 Nov 2010 18:14:39 +0200 To: m.t.h...@ai.rug.nl CC: r-help@r-project.org Subject: Re: [R] Three-way ANOVA shows me two-way results It seems your model defines more parameters then your data permits (e.g: singular). if you do google search

Re: [R] problem with parameters

2010-11-23 Thread Joshua Wiley
Hi Andrea, Help us help you. Please read this: http://www.R-project.org/posting-guide.html and then reply to the list with more details on what you have, what you have tried, and what you expect the results to be. Then we can get you started (and maybe even finished). To borrow from Sherlock

Re: [R] cross validation using e1071:SVM

2010-11-23 Thread Neeti
@Francial Giscard LIBENGUE please post your query again so that with different subject -- View this message in context: http://r.789695.n4.nabble.com/cross-validation-using-e1071-SVM-tp3055335p3055831.html Sent from the R help mailing list archive at Nabble.com.

Re: [R] cross validation using e1071:SVM

2010-11-23 Thread Neeti
could anyone help me with my last problem. if the question is not clear please let me know thank you Hi everyone I am trying to do cross validation (10 fold CV) by using e1071:svm method. I know that there is an option (“cross”) for cross validation but still I wanted to make a

[R] PROC COMPARE

2010-11-23 Thread Jim Edd Jones
Is it possible to use R to compare two datasets to look for discrepancies, as one would with the SAS procedure PROC COMPARE? Any help on this would be greatly appreciated. James E. Jones 135 Salina St. Lafayette, CO 80026 [[alternative HTML version deleted]]

Re: [R] PROC COMPARE

2010-11-23 Thread Henrique Dallazuanna
See all.equal, a toy example: all.equal(iris, Orange) On Tue, Nov 23, 2010 at 3:17 PM, Jim Edd Jones jejones...@comcast.netwrote: Is it possible to use R to compare two datasets to look for discrepancies, as one would with the SAS procedure PROC COMPARE? Any help on this would be greatly

Re: [R] Question about list function

2010-11-23 Thread Henrique Dallazuanna
Try this: as.list(scan(your_txt_file)) On Tue, Nov 23, 2010 at 3:05 PM, Guido Leoni guido.le...@gmail.com wrote: Dear List I'm a newbie R user. I'm utilizing the list function in order to make a var like this: clusters-list(a=var1,b=var2) My problem is that the total numer of variables

Re: [R] Question about list function

2010-11-23 Thread jim holtman
Try this: input - textConnection(var1 a + var2 b + var3 c + var4 d + var5 e) x - read.table(input, as.is = TRUE) close(input) # create a list xList - as.list(x$V2) names(xList) - x$V1 xList $var1 [1] a $var2 [1] b $var3 [1] c $var4 [1] d $var5 [1] e On Tue, Nov 23, 2010 at 12:05

Re: [R] the first. from SAS in R

2010-11-23 Thread Charles C. Berry
On Tue, 23 Nov 2010, David Winsemius wrote: On Nov 23, 2010, at 11:04 AM, Charles C. Berry wrote: On Tue, 23 Nov 2010, Dennis Murphy wrote: Interesting. Check this out: u - sample(c(TRUE, FALSE), 10, replace = TRUE) u [1] FALSE FALSE TRUE FALSE FALSE TRUE FALSE FALSE FALSE FALSE

[R] effect that cluster sampling has on count data (extra-multinomial or extra-binomial variation?)

2010-11-23 Thread Powers, Randall - BLS
I'm wondering if R has a package for generating count data subject to extra-multinomial variation (or extra-binomial variation?)?? In other words, the effect that cluster sampling has on count data. For example: Have a multinomial (p,n) distribution with p=(p11, p12,p21,p22) and n=number of

Re: [R] PROC COMPARE

2010-11-23 Thread Joshua Wiley
Hi James, Here is another option: all.equal() is more general, the sapply() option I wrote assumes both data frames are basically the same, but that maybe a few elements are different (e.g., you have two data sets and are not sure which is more recent). y - x - iris # x y copies of iris ##

[R] Plot two zoo object with different indexes

2010-11-23 Thread Manta
Dear R community, I have the following two zoo objects: MONTHLY CPI plot(z) par(usr) [1] 1977.76333 2011.15333 70.39856 227.03744 z=zooreg(cpius$Value,as.yearmon(1979-11),frequency=12) str(z) ‘zooreg’ series from Nov 1979 to Oct 2010 Data: num [1:372] 76.2 77 77.8 78.5 79.5 80.3 81.1

Re: [R] Plot two zoo object with different indexes

2010-11-23 Thread Gabor Grothendieck
On Tue, Nov 23, 2010 at 12:53 PM, Manta mantin...@libero.it wrote: Dear R community, I have the following two zoo objects: MONTHLY CPI plot(z) par(usr) [1] 1977.76333 2011.15333   70.39856  227.03744 z=zooreg(cpius$Value,as.yearmon(1979-11),frequency=12) str(z) ‘zooreg’ series from Nov

Re: [R] Plot two zoo object with different indexes

2010-11-23 Thread Manta
Thanks as always Gabor. What I was looking for was to plot the daily series on the same graph of the monthly using the command 'lines'. That is why I wanted to change the index of one of the two, as it seems the reason why at the moment it is not working! -- View this message in context:

Re: [R] Plot two zoo object with different indexes

2010-11-23 Thread Gabor Grothendieck
On Tue, Nov 23, 2010 at 1:11 PM, Manta mantin...@libero.it wrote: Thanks as always Gabor. What I was looking for was to plot the daily series on the same graph of the monthly using the command 'lines'. That is why I wanted to change the index of one of the two, as it seems the reason why at

Re: [R] the first. from SAS in R

2010-11-23 Thread Seeliger . Curt
Is there any similar function in R to the first. in SAS? ?duplicated a$d - ifelse( duplicated( a$a ), 0 , 1 ) a$d.2 - as.numeric( !duplicated( a$a ) ) Actually, duplicated does not duplicate SAS' first. operator, though it may suffice for the OP's needs. To illustrate,

[R] specifying colours in a ggplot2 piechart

2010-11-23 Thread John Kane
Someone was asking how to do a 16 category piechart in OpenOffice Calc and it appears that it can not be done (which we, probably, should be happy about) but I thought that I'd try it in ggplot2. It works but I then thought I'd like to make the colours more distinctive but fro some reason I

[R] Barplot and plot means with confidence intervals in the same plot

2010-11-23 Thread Ignacio Alborés
I want to obtain a barplot with error bars and a mean plot with error bars with other scale on the left in the same plot. I need help to get the same x axis (centered) when overlay two plots (barplot2 and plotCI (or lineplot.CI) with errorbars), with diferent y axis. I use par (new=T) but the

Re: [R] specifying colours in a ggplot2 piechart

2010-11-23 Thread baptiste auguie
fixing the various typos in your code, this works, lat_plot() + scale_fill_manual(value=mycolours) HTH, baptiste On 23 November 2010 19:47, John Kane jrkrid...@yahoo.ca wrote: Someone was asking how to do a 16 category piechart in OpenOffice Calc and it appears that it can not be done

Re: [R] the first. from SAS in R

2010-11-23 Thread Charles C. Berry
On Tue, 23 Nov 2010, seeliger.c...@epamail.epa.gov wrote: Is there any similar function in R to the first. in SAS? ?duplicated a$d - ifelse( duplicated( a$a ), 0 , 1 ) a$d.2 - as.numeric( !duplicated( a$a ) ) Actually, duplicated does not duplicate SAS' first. operator, though

Re: [R] specifying colours in a ggplot2 piechart

2010-11-23 Thread John Kane
Obviously I cannot type or spell. Your solution was great but do you have any idea why scale_colour_manual(value=mycolours) does not work. More typos on my part or am I misunderstanding what scale_colour_manual(value=mycolours) is supposed to do? Thanks --- On Tue, 11/23/10, baptiste

Re: [R] specifying colours in a ggplot2 piechart

2010-11-23 Thread baptiste auguie
Colour affects the lines and points in ggplot2, whilst fill is for areas (polygons, etc.) By the way, I can't type either -- I meant last_plot(). baptiste On 23 November 2010 20:18, John Kane jrkrid...@yahoo.ca wrote: Obviously I cannot type or spell. Your solution was great but do you

[R] Possible to pretty-printing using str()?

2010-11-23 Thread Janko Thyson
Dear list, I'm looking for a suitable way to sort of one-line-pretty-print an arbitrary R object in some of my log outputs. Consider this: cat(paste(The object/value is: , x, ., sep=), sep=\n) No problem if x is of class: - character (length=1) - numeric (length=1) - logical

Re: [R] the first. from SAS in R

2010-11-23 Thread Seeliger . Curt
... # I've found I need to lag a column to mimic SAS' first. # operator, thusly, though perhaps someone else knows # differently. Note this does not work on unordered # dataframes! lag.k1 - c(NA, tt$k1[1:(nrow(tt) - 1)]) tt$r.first.k1 - ifelse(is.na(lag.k1), 1, tt$k1 != lag.k1)

Re: [R] Possible to pretty-printing using str()?

2010-11-23 Thread Janko Thyson
A short annotation: I already considered deparsing the object via deparse(x), but that's not exactly the info I'd like. I'm only interested in a short summary like str() would give me, not in the deparsed object. Thx, Janko ## SYSTEM INFO ## Windows XP SP3 R 2.12.0 (patched as

Re: [R] Executing multiple .R files

2010-11-23 Thread Henrik Bengtsson
See R.utils package and sourceDirectory(), e.g. sourceDirectory(R/); You can also specify that you only want to reload files that have been modified, e.g. sourceDirectory(R/, modifiedOnly=TRUE); /Henrik On Tue, Nov 23, 2010 at 7:04 AM, Santosh Srinivas santosh.srini...@gmail.com wrote: Hello

Re: [R] Possible to pretty-printing using str()?

2010-11-23 Thread Duncan Murdoch
On 23/11/2010 2:35 PM, Janko Thyson wrote: A short annotation: I already considered deparsing the object via deparse(x), but that's not exactly the info I'd like. I'm only interested in a short summary like str() would give me, not in the deparsed object. I think you need to give more details.

Re: [R] Possible to pretty-printing using str()?

2010-11-23 Thread Janko Thyson
Aha: capture.output(str(x)) Sorry for the two previous eMails, but sometimes posting to the list helps me thinking somehow ;-) Greetz, Janko -Ursprüngliche Nachricht- Von: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] Im Auftrag von Janko Thyson Gesendet:

[R] Help with grouped barplot

2010-11-23 Thread Art Burke
Given the data structure below, how can I create a bar plot for the values of disc for each area grouped by year? bar -structure(list(year = c(2003, 2003, 2003, 2003, 2003, 2003, 2003, 2007, 2007, 2007, 2007, 2007, 2007, 2007), area = structure(c(6L, 4L, 1L, 2L, 3L, 5L, 7L, 6L, 4L, 1L, 2L, 3L,

Re: [R] Help with grouped barplot

2010-11-23 Thread Gabor Grothendieck
On Tue, Nov 23, 2010 at 2:49 PM, Art Burke art.bu...@educationnorthwest.org wrote: Given the data structure below, how can I create a bar plot for the values of disc for each area grouped by year? bar -structure(list(year = c(2003, 2003, 2003, 2003, 2003, 2003, 2003, 2007, 2007, 2007, 2007,

Re: [R] overlay histograms on map at map coordinates

2010-11-23 Thread Greg Snow
Look at the subplot function in the TeachingDemos package. -- Gregory (Greg) L. Snow Ph.D. Statistical Data Center Intermountain Healthcare greg.s...@imail.org 801.408.8111 -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r- project.org] On Behalf Of

Re: [R] how to add frequencies to barplot

2010-11-23 Thread Greg Snow
I would suggest reading the discussion (all the posts) starting at: http://tolstoy.newcastle.edu.au/R/e2/help/07/08/22858.html This gives reasons for not doing what you are asking along with some alternatives. -- Gregory (Greg) L. Snow Ph.D. Statistical Data Center Intermountain Healthcare

Re: [R] Barplot and plot means with confidence intervals in the same plot

2010-11-23 Thread Yuliya Matveyeva
2010/11/23 Ignacio Alborés nachoalbo...@gmail.com I want to obtain a barplot with error bars and a mean plot with error bars with other scale on the left in the same plot. I need help to get the same x axis (centered) when overlay two plots (barplot2 and plotCI (or lineplot.CI) with

Re: [R] save a regression model that can be used later

2010-11-23 Thread Greg Snow
One possibility for what is happening is that when you load the object it does not automatically load the package as well, so you get an error when working with the object. Try loading lme4 package in a new session, then load the saved object and see if things work for you. -- Gregory (Greg)

  1   2   >