Re: [R] save conditions in a list

2012-07-06 Thread Christof Kluß
Hi Rui Barradas thank you very much, that's what I searched for result - lapply(conds, fun, DF) works, if day - DF$day val - DF$val Thanks Christof Am 02-07-2012 20:44, schrieb Rui Barradas: Hello, I'm not sure if this is what you want. #-- Make up a dataset

[R] function on strsplit output

2012-07-06 Thread Sarah Auburn
Hi, I am trying to format some data (example matrix m below) for which each data point has 2 associated values separated by a comma. I want to replace values 3 with 0 to give the example output below. I have got as far as: out-lapply(strsplit(m,,),as.numeric) Failed to identify anything along

Re: [R] Please help

2012-07-06 Thread Pascal Oettli
Hello, I did a mistake in my script: hour - day*24; hour - c(hour[1],cumsum(hour[1:length(hour)-1])) should be replaced by: hour - day*24; hour - c(0,cumsum(hour[1:length(hour)-1])) Regards Le 05/07/2012 10:49, Pascal Oettli a écrit : Hello, A short code to include real monthly

Re: [R] i need help in documentation

2012-07-06 Thread purushothaman
Hi, i tried with inlinedocs package function extract.docs.file. it have exactly what i need but it return as list. can u help me how to convert list to html documentation. Thanks B.Purushothaman -- View this message in context:

Re: [R] Graph showing fitted values obtained by binomial GLM

2012-07-06 Thread linda.kate
Thank you Simon! -- View this message in context: http://r.789695.n4.nabble.com/Graph-showing-fitted-values-obtained-by-binomial-GLM-tp4635573p4635581.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing list

[R] Definition of AIC (Akaike information criterion) for normal error models

2012-07-06 Thread Kunio takezawa
Dear R users (r-help@r-project.org), The definition of AIC (Akaike information criterion) for normal error models has just been changed. Please refer to the paper below on this matter. Eq.(22) is the new definition. The essential part is RSS(n+q+1)/(n-q-3); it is close to GCV. The paper

Re: [R] function on strsplit output

2012-07-06 Thread Jessica Streicher
sap-sapply(strsplit(m,,),as.numeric) sap[which(sap = 3)]-0 mNew-matrix(apply(sap,2,FUN=function(x){paste(x,collapse=,)}),ncol=4) works? On 06.07.2012, at 08:47, Sarah Auburn wrote: Hi, I am trying to format some data (example matrix m below) for which each data point has 2 associated values

Re: [R] function on strsplit output

2012-07-06 Thread Sarah Auburn
Perfect, thank you! From: Jessica Streicher j.streic...@micromata.de To: Sarah Auburn saub...@yahoo.com Cc: R help r-help@r-project.org Sent: Friday, 6 July 2012, 17:48 Subject: Re: [R] function on strsplit output sap-sapply(strsplit(m,,),as.numeric) sap[which(sap = 3)]-0

[R] Which is Best function to find Maximum Likelihood Estimates in R?

2012-07-06 Thread chamilka
I read there are several R functions to get the Maximum Likelihood Estimates of a probability distribution with the available data in R. For Example... http://stat.ethz.ch/R-manual/R-devel/library/stats4/html/mle.html mle {stats4} http://stat.ethz.ch/R-manual/R-devel/library/stats/html/optim.html

Re: [R] Problem using Lag on time series

2012-07-06 Thread Jeff Newmiller
You have made it extraordinarily difficult to help you. 1) Using a library function without indicating which library it is from. This is really a special case of failing to provide a reproducible example as requested in the posting guide. 2) You provide data that does not appear to be a result

Re: [R] function on strsplit output

2012-07-06 Thread Jessica Streicher
Could also do it in one apply i realized.. apply(m,1:2,FUN=function(x) {y-as.numeric(strsplit(x,split=,)[[1]]);y[which(y = 3)]-0;paste(y,collapse=,)}) On 06.07.2012, at 10:18, Jessica Streicher wrote: sap-sapply(strsplit(m,,),as.numeric) sap[which(sap = 3)]-0

[R] Different level set when predicting with e1071's Naive Bayes classifier

2012-07-06 Thread Joao Azevedo
Hi! I think I had some issues with the charset on the previous message, so I'm sending this again. Sorry for the double post. I'm using the Naive Bayes classifier provided by the e1071 package ( http://cran.r-project.org/web/packages/e1071) and I've noticed that the predict function has a

[R] Problem with Apriori

2012-07-06 Thread Jolinda Bartlett
Good morning When I try to run an Apriori, I keep often the following message: checking subsets of size 1Error in apriori(tr, parameter = list(supp = 3e-05, conf = 0.05, minlen = 2, : internal error in trio library Other time it works fine on the same data. What can be the reason for this?

Re: [R] function curve() swap axes

2012-07-06 Thread Boudewijn Verkooijen
Thank you guys! The solution turned out to be very simple, but I couldn't find it anywhere (not even in the help of the function curve()). I think there's not much difference in both methods, but with curve you can easily adjust the amount of a values with 'n = ...' option, to make a smooth line.

Re: [R] Problem with Apriori

2012-07-06 Thread Pascal Oettli
Good morning, 1. What is your OS? 2. Which package do you use to perform the Apriori? 3. Without providing any reproducible code, it will probably be complicated to answer you. Best Regards Le 06/07/2012 18:15, Jolinda Bartlett a écrit : Good morning When I try to run an Apriori, I keep

Re: [R] i need help in documentation

2012-07-06 Thread Jim Lemon
On 07/05/2012 10:37 PM, purushothaman wrote: Hi, i am new in R development , so i need help in documentation from comments in R file for example sample-function(filepath) { ##title read csv file ##author Purushoth ##description this function is used to read a csv

[R] How to compute hazard function using coxph.object

2012-07-06 Thread Yashwanth M.R
My question is, how to compute hazard function(H(t)) after building the coxph model. I even aware of the terminology that differs from hazard function(H(t)) and the hazard rate(h(t)). Here onward I wish to calculate both. Here what I have done in two different methods;

Re: [R] save conditions in a list

2012-07-06 Thread arun
Hi, You could also try: set.seed(1) n - 1e2 DF - data.frame(year=2010 + sample(3, n, TRUE),     day=sample(365, n, TRUE),     val=sample(100, n, TRUE)) a = day 100; b = val 50; c = year == 2012 conds - list(a=a, b=b, c=c)#To check individual conditions eval(parse(text=conds[[1]]))

[R] How to compute hazard function using coxph.object

2012-07-06 Thread Yashwanth M.R
My question is, how to compute hazard function(H(t)) after building the coxph model. I even aware of the terminology that differs from hazard function(H(t)) and the hazard rate(h(t)). Here onward I wish to calculate both. Here what I have done in two different methods;

[R] coefficient of lasso and lasso plot

2012-07-06 Thread guguspelangi
i have try lasso use lars package with real data.. when i running the data and plot it (plot(lasso)).. in plot showed|beta|/max|beta, when i calculated beta from coefficient lasso (coef(lasso)).. there are some differentiation between the calculation result and the plot.. what cause it? must i

Re: [R] Error in hclust?

2012-07-06 Thread Jean V Adams
Mateus, I see what you mean. I can't figure it out. I found a message from a few years ago that mentions the same problem. http://tolstoy.newcastle.edu.au/R/e8/help/09/12/8165.html This may be a bug. Since hclust() is in the stats package, I am ccing the maintainer, r-c...@r-project.org, on

[R] S3 / S4 Method for print and frames

2012-07-06 Thread Julien-Laferriere Alice
Hello, i am trying to understand the difference between S3 function print and S4 method. My main problem is that i need to count the environment (frames) (i use /eval/() afterwards with /parent.frame/() ). If I use a S4 method/print/, a 'hidden' environment seem to appears. It is not the

Re: [R] Error in hclust?

2012-07-06 Thread Mateus Teixeira
Thank you, Jean. Mateus 2012/7/6 Jean V Adams jvad...@usgs.gov Mateus, I see what you mean. I can't figure it out. I found a message from a few years ago that mentions the same problem. http://tolstoy.newcastle.edu.au/R/e8/help/09/12/8165.html This may be a bug. Since hclust() is in

[R] How to do goodness-of-fit diagnosis and model checking for rlm in R?

2012-07-06 Thread Michael
Hi all, I am reading the MASS book but it doesn't give examples about the diagnosis and model checking for rlm... My data is highly non-Gaussian so I am using rlm instead of lm. My questions are: 0. Are goodness-of-fit and model-checking using rlm completely the same as usual regression?

Re: [R] automatic completion of object names

2012-07-06 Thread peter dalgaard
On Jul 6, 2012, at 15:36 , knallg...@gmx.com wrote: Hello there, I just upgraded to R 2.15 (from R 2.12) on a Windows XP machine and noticed some puzzling behaviour (that in my opinion did not exist in R 2.12). It is possible now to call objects without spelling out the full object

Re: [R] How to do goodness-of-fit diagnosis and model checking for rlm in R?

2012-07-06 Thread Bert Gunter
On Fri, Jul 6, 2012 at 7:19 AM, Michael comtech@gmail.com wrote: Hi all, I am reading the MASS book but it doesn't give examples about the diagnosis and model checking for rlm... My data is highly non-Gaussian so I am using rlm instead of lm. Well, rlm() is not necessarily going

Re: [R] automatic completion of object names

2012-07-06 Thread Jeff Newmiller
This has always been how the $ operator has worked. See ?[[, and use that operator directly if you want exact matching. --- Jeff NewmillerThe . . Go Live...

Re: [R] automatic completion of object names

2012-07-06 Thread David Winsemius
On Jul 6, 2012, at 9:36 AM, knallg...@gmx.com wrote: Hello there, I just upgraded to R 2.15 (from R 2.12) on a Windows XP machine and noticed some puzzling behaviour (that in my opinion did not exist in R 2.12). It is possible now to call objects without spelling out the full object

Re: [R] Error in hclust?

2012-07-06 Thread Jean V Adams
Mateus, Apparently I cc-ed the wrong group (see snippet below). I think you found a bug. I suggest you write something up and send it to r-b...@r-project.org Follow the directions on How to report a bug http://www.stat.ucl.ac.be/ISdidactique/Rhelp/library/base/html/bug.report.html Jean

[R] Anova Type II and Contrasts

2012-07-06 Thread mails
the study design of the data I have to analyse is simple. There is 1 control group (CTRL) and 2 different treatment groups (TREAT_1 and TREAT_2). The data also includes 2 covariates COV1 and COV2. I have been asked to check if there is a linear or quadratic treatment effect in the data. I

[R] Mann-Whitney by group

2012-07-06 Thread David Chertudi
Hello, I'm having some issues grouping cases for some Mann-Whitney U tests I'm attempting to run. I'm willing to use wilcox.test if it'll work; I've also tried wilcox_test() from the coin package. Here's the deal: for each column (A through H), I would like to run the two-sample independent

[R] convert a table

2012-07-06 Thread Amanduh320
I have my data in a table table - table(test2$Filename, test2$PREDICT) I need to convert this table so it keeps the same structure, but is a different format. The current output is count data by Filename and I want to get the max for each Filename. Columns are: Filename, 1, 2, 3, 4, 5, 6, 7

[R] Maximum number of patterns and speed in grep

2012-07-06 Thread mdvaan
Hi, I am using R's grep function to find patterns in vectors of strings. The number of patterns I would like to match is 7,700 (of different sizes). I noticed that I get an error message when I do the following: data - array() for (j in 1:length(x)) { array[j] -

[R] help in R programming

2012-07-06 Thread umesh khatri
Dear everyone I'm dealing with a problem related to Poisson Ridge Regression. If anyone can help me in this regard by telling if any changes in the source code of glm.fit may help -- Regards Umesh Khatri [[alternative HTML version deleted]]

Re: [R] Maximum number of patterns and speed in grep

2012-07-06 Thread Sarah Goslee
Hi, Given that you can't provide a full example, please at least provide str() on your data, more complete information on the problem, and ideally a small toy example that demonstrates precisely what you are doing. For instance, you tell us that you get an error message but you never tell us

[R] Based-Density cluster algorithms, Help where find some algorithms

2012-07-06 Thread gianni lavaredo
Dear Researches, Did anyone know where to find the following Based-Density cluster algorithms in R or function implemented? OPTICS (Ordering Points To Identify the Clustering Structure) DECLUN (DENsity CLUstering) Thanks in Advance Gianni [[alternative HTML version deleted]]

Re: [R] convert a table

2012-07-06 Thread andrija djurovic
Hi. Here are some examples that could be useful here: set.seed(1) x - sample(1:2, 100, rep=T) y - sample(1:3, 100, rep=T) TAB - table(x,y) #check the object class class(TAB) apply(TAB, 2, max) apply(TAB[,2:3], 2, max) apply(TAB, 1, max) DF - as.data.frame.matrix(TAB) class(DF) sapply(DF, max)

[R] Compare date Oracle with Sys.time

2012-07-06 Thread cindy.dol
I would like to import only datas of my table where DATEtoday-7days. But my datas in Oracle are 'dates' and in R are 'characters'. now_7-format(Sys.time()-(7*60*60*24), %Y-%m-%d 00:00:00) How to do? -- View this message in context:

Re: [R] convert a table

2012-07-06 Thread David L Carlson
A table is a matrix Filename - sample(LETTERS[1:5], 50, replace=TRUE) PREDICT - sample(1:7, 50, replace=TRUE) To get a data frame that preserves the row/column structure of the table use dfm - as.data.frame.matrix(tbl) But you can do what you want on the table directly: rowmx - apply(tbl, 1,

Re: [R] convert a table

2012-07-06 Thread David L Carlson
I managed to leave out one key line: Filename - sample(LETTERS[1:5], 50, replace=TRUE) PREDICT - sample(1:7, 50, replace=TRUE) tbl - table(Filename, PREDICT) # this one. Then the rest follows. -- David -Original Message- From:

Re: [R] Mann-Whitney by group

2012-07-06 Thread David L Carlson
Can you describe how your data is organized. It is clear there are eight columns, but it is not clear how the groups are represented, a Group column or do the groups have to be assembled from information in another column (a column with CD8.14, etc)? Create a small version of the data and use

Re: [R] help in R programming

2012-07-06 Thread David L Carlson
Please hold. The r-help hive mind is attempting to probe your brain for further details of your problem related to Poisson Ridge Regression. -- David L Carlson Associate Professor of Anthropology Texas AM University College Station, TX 77843-4352

Re: [R] Maximum number of patterns and speed in grep

2012-07-06 Thread Gabor Grothendieck
On Fri, Jul 6, 2012 at 10:45 AM, mdvaan mathijsdev...@gmail.com wrote: Hi, I am using R's grep function to find patterns in vectors of strings. The number of patterns I would like to match is 7,700 (of different sizes). I noticed that I get an error message when I do the following: data -

Re: [R] convert a table

2012-07-06 Thread Jeff Newmiller
If you don't know what format your data is in, how are we supposed to know? Have you tried looking at your data file with a text editor? Perhaps reading your data using read.table or read.csv will help? Use the str function to learn its structure? use dput as the posting guide recommends so we

Re: [R] help in R programming

2012-07-06 Thread Jeff Newmiller
Ding! Semantic error : problem not defined. Return to beginning of Posting Guide mentioned at bottom of every email on this list, process contents, and reformulate question. --- Jeff NewmillerThe

Re: [R] Anova Type II and Contrasts

2012-07-06 Thread David Winsemius
Dear Peter; This is an exact duplicate of a question posted on SO. Cross-posting is deprecated on Rhelp. -- David. On Jul 6, 2012, at 11:06 AM, mails wrote: the study design of the data I have to analyse is simple. There is 1 control group (CTRL) and 2 different treatment groups (TREAT_1

Re: [R] How to do goodness-of-fit diagnosis and model checking for rlm in R?

2012-07-06 Thread Michael
Thanks a lot Bert. Yet still, I would like to know does the diagnosis and model-checking for rlm follow the same procedure as lm ... Thank you! On Fri, Jul 6, 2012 at 9:27 AM, Bert Gunter gunter.ber...@gene.com wrote: On Fri, Jul 6, 2012 at 7:19 AM, Michael comtech@gmail.com

Re: [R] Anova Type II and Contrasts

2012-07-06 Thread John Fox
Dear Peter, Because your model is additive, type-II and type-III tests are identical, and the t-tests for the linear and quadratic coefficients are interpretable. I hope this helps, John John Fox Sen. William McMaster Prof. of Social Statistics

[R] How to import SAS data in R?

2012-07-06 Thread C W
Hi all I have a large SAS data set, how do I get it read in R? The data is too big (about 400,000 rows by 100 columns) to be saved as an Excel file. How should I get it read in R? Any packages? I don't seem to find any. Thanks, Mike [[alternative HTML version deleted]]

Re: [R] How to import SAS data in R?

2012-07-06 Thread Jeff Newmiller
You need to learn how to search. RSiteSearch(SAS) or Google : R read SAS --- Jeff NewmillerThe . . Go Live... DCN:jdnew...@dcn.davis.ca.usBasics: ##.#. ##.#. Live

Re: [R] How to import SAS data in R?

2012-07-06 Thread Duncan Murdoch
On 06/07/2012 1:11 PM, C W wrote: Hi all I have a large SAS data set, how do I get it read in R? The data is too big (about 400,000 rows by 100 columns) to be saved as an Excel file. How should I get it read in R? Any packages? I don't seem to find any. You could write it out in some plain

Re: [R] Compare date Oracle with Sys.time

2012-07-06 Thread Marc Schwartz
On Jul 6, 2012, at 10:50 AM, cindy.dol wrote: I would like to import only datas of my table where DATEtoday-7days. But my datas in Oracle are 'dates' and in R are 'characters'. now_7-format(Sys.time()-(7*60*60*24), %Y-%m-%d 00:00:00) How to do? Huh? How are you getting the data from Oracle

Re: [R] How to import SAS data in R?

2012-07-06 Thread Marc Schwartz
On Jul 6, 2012, at 12:11 PM, C W wrote: Hi all I have a large SAS data set, how do I get it read in R? The data is too big (about 400,000 rows by 100 columns) to be saved as an Excel file. How should I get it read in R? Any packages? I don't seem to find any. Thanks, Mike Strangely

[R] Is R BOOT package available for 2.15.1 version for Mac OS?

2012-07-06 Thread SergeyHC
Hi everyone! I am trying to install boot package in order to do bootstrap analysis with R. I am trying to use the following command: install.packages(/var/folders/++/++53BE++6+0++4RjPqRgNE+-lHI/-Tmp-//Rtmpvmc3aE/downloaded_packages) But get the following message: package

[R] automatic completion of object names

2012-07-06 Thread knallgrau
Whoops! Thanks to everyone and sorry for asking something so trivial and obvious for everyone. I really hadn't noticed it before. Cheers! i. - Original Message - From: David Winsemius Sent: 07/06/12 05:32 PM To: knallg...@gmx.com Subject: Re: [R] automatic completion of

[R] Poisson Ridge Regression

2012-07-06 Thread umesh khatri
Dear everyone I'm dealing with a problem related to Poisson Ridge Regression. If anyone can help me in this regard by telling if any changes in the source code of glm.fit may help -- Regards Umesh Khatri __ R-help@r-project.org mailing list

Re: [R] Maximum number of patterns and speed in grep

2012-07-06 Thread mdvaan
Thanks for the quick response. I should phrase my question differently because everything is working fine, I am just trying to find a more efficient approach: 1. What's the maximum size of the patterns argument in grep? Can't find it online. 2. I am trying to match 7,700 character strings to

Re: [R] convert a table

2012-07-06 Thread Amanduh320
Fantastic, thank you so much! -- View this message in context: http://r.789695.n4.nabble.com/convert-a-table-tp4635615p4635627.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing list

[R] Describing Multistage PPS Survey Designs

2012-07-06 Thread david rae
Hi all, I am wondering if the survey package has the capability to describe a multistage survey with PPS clusters. Speficially I have a national survey that is intially stratified by region, then by PPS cluster by locality. I'm not sure if this needs to be done with an FPC statement or if I need

[R] Tables extraction in R ?

2012-07-06 Thread Greeknovice
Hi, I 'm a novice user of R statistics and my hands-on experience with it is minimal. I want to create a table for my MBA course assignment that looks like the ones that SPSS and MS Excel produces ,the data that the table has to include are the following : table(agec) agec 1 2 3 749 160

Re: [R] Compare date Oracle with Sys.time

2012-07-06 Thread arun
Hello, Not quite understand about ur problem.  Could you dput your data? The below should be general format to convert to class date  now_7-format(Sys.time()-(7*60*60*24), %Y-%m-%d 00:00:00) now_7 [1] 2012-06-29 00:00:00 class(now_7) [1] character now_7_new-as.Date(Sys.time()-(7*60*60*24),

Re: [R] how to do a graph with tree different colors??

2012-07-06 Thread John Kane
Thanks Dennise, Now I see what you have. It was not obvious from the original data set that you posted that is is a time series. I am afraid my knowledge of time series is about zero but I am sure others here can help now that we have a good idea of what the data really looks

Re: [R] Is R BOOT package available for 2.15.1 version for Mac OS?

2012-07-06 Thread Marc Schwartz
On Jul 6, 2012, at 10:46 AM, SergeyHC wrote: Hi everyone! I am trying to install boot package in order to do bootstrap analysis with R. I am trying to use the following command: install.packages(/var/folders/++/++53BE++6+0++4RjPqRgNE+-lHI/-Tmp-//Rtmpvmc3aE/downloaded_packages) But

Re: [R] Compare date Oracle with Sys.time

2012-07-06 Thread arun
Hello, On further thinking about the issue, I guess this may be close to what you are looking for: date1-c(2012-06-12,2012-06-13,2012-06-14,2012-06-27,2012-07-01,2012-07-02,2012-07-03,2012-07-04)  class(date1) [1] character #converting to date class date2-as.Date(date1,format=%Y-%m-%d)

Re: [R] How to import SAS data in R?

2012-07-06 Thread C W
Quick question, what the memory size in R? I converted to CSV, but only 53300 of the 1,000,000 rows were read in. Did R run out of memory? If so, is there a work around? Thanks, Mike On Fri, Jul 6, 2012 at 1:24 PM, Duncan Murdoch murdoch.dun...@gmail.comwrote: On 06/07/2012 1:11 PM, C W

[R] Error in plot.new() : figure margins too large

2012-07-06 Thread Karly Harrod
Hello All, I am running the following code in RStudio, and I keep on getting an error message that says: Error in plot.new() : figure margins too large Is there something that I am doing wrong? # Import Data nba - read.csv(http://datasets.flowingdata.com/ppg2008.csv;, sep=,) nba #Sort Data

Re: [R] Tables extraction in R ?

2012-07-06 Thread John Kane
Have a look at the xtables package. I have not used it in some time but I think it may do what you want. A google search R statistics xtables should bring up some useful information on this. John Kane Kingston ON Canada -Original Message- From: sggk...@gmail.com Sent: Fri, 6 Jul

Re: [R] Poisson Ridge Regression

2012-07-06 Thread John Kane
Probabaly not since you do not tell us what the problem is. Please read the posting notes at the bottom of this email John Kane Kingston ON Canada -Original Message- From: khatriumes...@gmail.com Sent: Fri, 6 Jul 2012 21:13:32 +0530 To: r-help@r-project.org Subject: [R] Poisson

Re: [R] Compare date Oracle with Sys.time

2012-07-06 Thread Jeff Newmiller
Yes, this request is woefully incomplete, but it looked to me like the question was how to select the data so that a limited amount of data could be imported, rather than loading the whole table and subsetting in R. If so, the discussion should be about constructing appropriate SQL (with or

Re: [R] Is R BOOT package available for 2.15.1 version for Mac OS?

2012-07-06 Thread David Winsemius
On Jul 6, 2012, at 11:46 AM, SergeyHC wrote: Hi everyone! I am trying to install boot package in order to do bootstrap analysis with R. I am trying to use the following command: install.packages(/var/folders/++/++53BE++6+0++4RjPqRgNE+-lHI/-Tmp-// Rtmpvmc3aE/downloaded_packages) But

Re: [R] Error in plot.new() : figure margins too large

2012-07-06 Thread John Kane
Your code works okay in a plan R terminal. It may be an Rstudio problem or perhaps you already have a graphics device open and are trying to draw into it? John Kane Kingston ON Canada -Original Message- From: khar...@lbl.gov Sent: Fri, 6 Jul 2012 10:39:41 -0700 To:

[R] Help with multinomial regression to survey

2012-07-06 Thread Diana Marcela Martinez Ruiz
I would like to know which function to use for modeling a multinomial logistic regression model in the package survey. thanks [[alternative HTML version deleted]] __ R-help@r-project.org mailing list

Re: [R] How to import SAS data in R?

2012-07-06 Thread David Winsemius
On Jul 6, 2012, at 1:39 PM, C W wrote: Quick question, what the memory size in R? I converted to CSV, but only 53300 of the 1,000,000 rows were read in. Did R run out of memory? If so, is there a work around? You probably have mismatched quotes. Consider using quote=. Also consider

Re: [R] Error in plot.new() : figure margins too large

2012-07-06 Thread mlell08
On 06.07.2012 20:02, John Kane wrote: Your code works okay in a plan R terminal. It may be an Rstudio problem or perhaps you already have a graphics device open and are trying to draw into it? John Kane Kingston ON Canada -Original Message- From: khar...@lbl.gov Sent: Fri,

Re: [R] How to import SAS data in R?

2012-07-06 Thread David Winsemius
On Jul 6, 2012, at 2:49 PM, C W wrote: Hey, David table(count.fields()) is telling me have 11 columns, but I have way more, more like 30 columns. table(count.fields(persistency.csv)) 1 2 3 4 5 6 7 8 9 11 439384 39617 16130 21993 12556

Re: [R] How to import SAS data in R?

2012-07-06 Thread C W
Has anyone read SAS file .sas7bdat into R. The above suggestions don't work. -M On Fri, Jul 6, 2012 at 3:06 PM, David Winsemius dwinsem...@comcast.netwrote: On Jul 6, 2012, at 2:49 PM, C W wrote: Hey, David table(count.fields()) is telling me have 11 columns, but I have way more, more

Re: [R] How to import SAS data in R?

2012-07-06 Thread Duncan Murdoch
On 06/07/2012 3:30 PM, C W wrote: Has anyone read SAS file .sas7bdat into R. The above suggestions don't work. They work if you follow them. Duncan Murdoch -M On Fri, Jul 6, 2012 at 3:06 PM, David Winsemius dwinsem...@comcast.netwrote: On Jul 6, 2012, at 2:49 PM, C W wrote: Hey,

Re: [R] how to do a graph with tree different colors??

2012-07-06 Thread William Dunlap
David W. wrote Ø It also looses it year labels and you may need to suppress the x-axis Ø and substitute your own. You can use plot(yourTimeSeries, type=n) to set up the usual time series axis and then use lines() or points() to add the data. E.g., plot(serieONI, type=n)

Re: [R] How to import SAS data in R?

2012-07-06 Thread Daniel Nordlund
-Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of C W Sent: Friday, July 06, 2012 12:30 PM To: David Winsemius Cc: r-help help Subject: Re: [R] How to import SAS data in R? Has anyone read SAS file .sas7bdat into R. The

Re: [R] Help with multinomial regression to survey

2012-07-06 Thread Thomas Lumley
On Sat, Jul 7, 2012 at 5:50 AM, Diana Marcela Martinez Ruiz dianamm...@hotmail.com wrote: I would like to know which function to use for modeling a multinomial logistic regression model in the package survey. thanks There isn't an implementation of multinomial regression in the survey package.

Re: [R] Plotting the probability curve from a logit model with 10 predictors

2012-07-06 Thread Greg Snow
Look at the Predict.Plot and TkPredict functions in the TeachingDemos package. These will not plot all 11 dimensions at once, but will plot 2 of the dimensions conditioned on the others. You can then change the conditioning to see relationships. These use base rather than ggplot graphics. On

[R] estimating NA values against selected slots

2012-07-06 Thread eliza botto
Dear R Users, Could you please help me on the following issue? I have a real large yearly data set. For each year I have 365 flow values. Some of the flow values are not known and that’s why you will see NA written in those slots. I wanted to know, is there a way that I can estimate those

Re: [R] Compare date Oracle with Sys.time

2012-07-06 Thread Cindy Dolomieu
Yes, it's more like that, to limit the amout of data imported in R from Oracle. - Mail original - De: Jeff Newmiller jdnew...@dcn.davis.ca.us À: Marc Schwartz marc_schwa...@me.com, cindy.dol cindy.dolom...@insa-lyon.fr Cc: r-help@r-project.org Envoyé: Vendredi 6 Juillet 2012 20:00:37

Re: [R] NADA Data Frame Format: Wide or Long?

2012-07-06 Thread MacQueen, Don
Hi Rich, So what you're faced with is that the cenros() function has no built-in methods for grouping or subsetting -- unlike some other R methods, especially those that work with the lattice package, or the many modeling functions like lm() that have a subset argument or employ a conditioning

[R] differences between survival models between STATA and R

2012-07-06 Thread JPF
Dear Community, I have been using two types of survival programs to analyse a data set. The first one is an R function called aftreg. The second one an STATA function called streg. Both of them include the same analyisis with a weibull distribution. Yet, results are very different. Shouldn't

Re: [R] Plotting the probability curve from a logit model with 10 predictors

2012-07-06 Thread Abraham Mathew
Ok, so let's say I have a logit equation outlined as Y= 2.5 + 3X1 + 2.3X2 + 4X3 + 3.6X4 + 2.2X5 So a one unit increase in X2 is associated with a 2.3 increase in Y, regardless of what the other predictor values are. So I guess instead of trying to plot of curve with all the predictors accounted

Re: [R] How to import SAS data in R?

2012-07-06 Thread David Winsemius
On Jul 6, 2012, at 3:30 PM, C W wrote: Has anyone read SAS file .sas7bdat into R. The above suggestions don't work. help(package=foreign) And read the Import/Export Manual -- -M On Fri, Jul 6, 2012 at 3:06 PM, David Winsemius dwinsem...@comcast.net wrote: On Jul 6, 2012, at

Re: [R] Mann-Whitney by group

2012-07-06 Thread Oxenstierna
Hi David, Thank you for the insight: I could have sworn I added a picture of the data, but providing the actual data is worlds easier to deal with, I'm sure. I've never used dput(), so I entered it using the dataframe in question as the object, and I've pasted the results below. Essentially, I

[R] How to compare stacked histograms/datasets

2012-07-06 Thread Atulkakrana
Hello All, I have a couple of stacked histograms which I need to compare/evaluate for similarity or difference. http://r.789695.n4.nabble.com/file/n4635668/Selection_011.png I believe rather than evaluating histograms is will be east to work with dataset used to plot these stacked histograms,

Re: [R] Is R BOOT package available for 2.15.1 version for Mac OS?

2012-07-06 Thread SergeyHC
David, thanks for your helpful response. Now i figured out that I do have boot command and do not need to install any additional packages. With regards t R GUI, I guess I would like to use it, but I could not figure out how to install it. For now I downloaded R-2.15.1.pkg file, which I could

[R] Plotting rpart trees with long list of class members

2012-07-06 Thread MarkBeauchene
I have a class with 732 members, so using rpart.plot is giving me a tiny plot in the middle of the window. Is there a good way to modify the plot, or replace the long list with something like group1? -- View this message in context:

[R] How to intall car package on linux

2012-07-06 Thread Chenwei Gao
I got this error when I tried to install the package in order to use the logit(logit transformation) function on linux. install.packages(car) Error in install.packages(car) : object 'car' not found I was using the Terminal in my macbook accessing the server in the lab. I could successfully

Re: [R] differences between survival models between STATA and R

2012-07-06 Thread Joshua Wiley
Hi J, You have not provided nearly enough information for us to evaluate whether the results should be similar. You are talking about two completely different packages, with no information on your data, only a small amount of information about your model (the same analysis) but clearly the

Re: [R] How to intall car package on linux

2012-07-06 Thread R. Michael Weylandt
Install.packages needs a quoted argument -- the error message suggests that its looking for a variable named `car` and not finding it -- in fact, most functions taking package names do, but library() allows some trickery, which is convenient, but occasionally leads to mishaps like this.

Re: [R] How to intall car package on linux

2012-07-06 Thread Michael Weylandt
Please do keep messages on list unless there's a particular reason to take them private: it increases your visibility to the wider community and, by extension, your chance of getting a good reply. On Jul 6, 2012, at 10:43 PM, Chenwei Gao ga...@umich.edu wrote: Thank you so much Michael !!