Re: [R] Problem reading mixed CSV file

2012-03-26 Thread Ashish Agarwal
Why does the output in the following say 2 and not 6? count.fields(textConnection(LL1532Ap,ABC# Depot-A+,,1971,8,2 + LL1532Ap,ABC# Depot-A+,Bhutan,1971,6,1 + LL1532Ap,ABC# Depot-A+,China,1971,17,1 + LL1532Ap,ABC# Depot-A+,China,1971,33,1 + LL1532Ap,ABC# Depot-A+,HongKong,1971,16,2 +

Re: [R] Problem reading mixed CSV file

2012-03-26 Thread Berend Hasselman
On 26-03-2012, at 08:16, Ashish Agarwal wrote: Why does the output in the following say 2 and not 6? count.fields(textConnection(LL1532Ap,ABC# Depot-A+,,1971,8,2 + LL1532Ap,ABC# Depot-A+,Bhutan,1971,6,1 + LL1532Ap,ABC# Depot-A+,China,1971,17,1 + LL1532Ap,ABC# Depot-A+,China,1971,33,1 +

Re: [R] Problem reading mixed CSV file

2012-03-26 Thread Ashish Agarwal
OMG. I think it uses comment character # as default in the argument. comment.char = # How do I turn it off? On Mon, Mar 26, 2012 at 11:56 AM, Berend Hasselman b...@xs4all.nl wrote: On 26-03-2012, at 08:16, Ashish Agarwal wrote: Why does the output in the following say 2 and not 6?

Re: [R] Problem reading mixed CSV file

2012-03-26 Thread Berend Hasselman
On 26-03-2012, at 08:26, Berend Hasselman wrote: On 26-03-2012, at 08:16, Ashish Agarwal wrote: Why does the output in the following say 2 and not 6? count.fields(textConnection(LL1532Ap,ABC# Depot-A+,,1971,8,2 + LL1532Ap,ABC# Depot-A+,Bhutan,1971,6,1 + LL1532Ap,ABC#

Re: [R] Problem reading mixed CSV file

2012-03-26 Thread Ashish Agarwal
comment.char = NULL does not work. Is there any way to make it NULL rather than having a specific character like '%'? On Mon, Mar 26, 2012 at 12:06 PM, Berend Hasselman b...@xs4all.nl wrote: comment.char = # How do I turn it off? ??? How about comment.car=% for example? Berend

Re: [R] Problem reading mixed CSV file

2012-03-26 Thread Berend Hasselman
On 26-03-2012, at 08:40, Ashish Agarwal wrote: comment.char = NULL does not work. Is there any way to make it NULL rather than having a specific character like '%'? Why don't you try something? comment.char= looks quite obvious and worth a try. Berend

Re: [R] Problem reading mixed CSV file

2012-03-26 Thread Berend Hasselman
On 26-03-2012, at 08:33, Ashish Agarwal wrote: OMG. I think it uses comment character # as default in the argument. comment.char = # How do I turn it off? ??? How about comment.car=% for example? Berend __ R-help@r-project.org mailing

Re: [R] Problem reading mixed CSV file

2012-03-26 Thread Ashish Agarwal
thanks. works very well. On Mon, Mar 26, 2012 at 12:12 PM, Berend Hasselman b...@xs4all.nl wrote: On 26-03-2012, at 08:40, Ashish Agarwal wrote: comment.char = NULL does not work. Is there any way to make it NULL rather than having a specific character like '%'? Why don't you try

Re: [R] reading header in txt file and making histogram

2012-03-26 Thread chuck.01
dbh is the column name. What is the name of your data? Lets assume your data is called YOUR_DATA then try: with(YOUR_DATA, hist(dbh)) OR hist(YOUR_DATA$dbh) OR hist(YOUR_DATA[, dbh]) etc... bamboohydraulics wrote Dear all I am a BEGINNER and have R on my Mac. I saved my excel

Re: [R] R Error : DATA to MATRIX

2012-03-26 Thread David Winsemius
On Mar 25, 2012, at 2:43 PM, ritwi...@isical.ac.in wrote: Thanks David, your suggestion works fine.btw I have another question..If I set (n,m) little bit large, say (n=20,m=10), R cannot handle the large data frame generated through expand.grid.Is there any way to increase

[R] What does package 'RDCOMClient' is not installed for 'arch=x64' exactly mean?

2012-03-26 Thread dthomas
Hi, I'm trying to use the excel.link package to write data to excel spreadsheets. I've installed the RDCOMClient package as required but get the error: package 'RDCOMClient' is not installed for 'arch=x64' I'm on Rx64 2.13.0. I assume it means the RDCOMClient package does not work on the x64

[R] Updating a Markov Chain

2012-03-26 Thread stivi
Hello, my question is if anyone has any good ideas how to create a Markov Chain from ordered data. So, I have some sort of time series, and if value1 happens as time1 and value2 happens at time2 I record this as an update to the probability transition matrix. The problem is that I cannot

Re: [R] How to test omitted level from a multiple level factor against overall mean in regression models?

2012-03-26 Thread Rolf Turner
The test you are requesting is ***MEANINGLESS***. The ``effect value'' of a single level is ill-defined (or in the more usual parlance, not estimable). The dummy.coef() procedure suggested by Gabor gives you point estimates *subject to the constraints* imposed by the contrasts used. The

Re: [R] Multivariate function from univariate functions

2012-03-26 Thread physicistintheory
Petr, thanks! I've now got a function I can work with. Regards! -- View this message in context: http://r.789695.n4.nabble.com/Multivariate-function-from-univariate-functions-tp4502670p4504599.html Sent from the R help mailing list archive at Nabble.com.

[R] sm.density kernel estimation for points

2012-03-26 Thread Matěj Plch
Hi! I have two dimensional dataset which has and I need to decide if a point lies in some confidence level. If a point has low confidence/density it can be anomaly which I need to find. For example: #load library library(sm) #get some data x.locs = c(74, 74.5, 75, 77,74.5) y.locs = c(64,

[R] Re : ROC Analysis

2012-03-26 Thread Pascal Oettli
Hi Camille, Probably you have to check wether there is any infinte value in x. Or calculate something like that for your x-axis: x[1:(ll-1)]+diff(x)/2 Regards, Pascal - Mail original - De : Camille Leclerc camille.lecl...@ymail.com À : r-help@r-project.org Cc : Envoyé le : Lundi 26

[R] row, col function but for a list (probably very easy question, cannot seem to find it though)

2012-03-26 Thread MBoersma
Hi guys, I'm quite new to R but quite enthousiastic. I'm trying to rewrite a bit of code in order to make it faster, so instead of nesting for loops I'm trying some apply functions. Since it's a combination of lists of matrices and other matrices, I would like to use lapply() on the list and

[R] The problem of using library(bigmemory)

2012-03-26 Thread yeheng...@gmail.com
To my understanding, library(bigmemory) should allow us to create a matrix larger than 2GB (say 17,000 by 17,000) by temporarily storing the data on harddisk. But I try the following codes to create such a matrix. library(bigmemory) A - big . matrix (17000 , 17000 , type =double, init = 0)

[R] Simple question regarding domain restrictions/piecewise functions in R

2012-03-26 Thread chad.mills
I am a novice R user. I would like to be able to graph some simple piecewise functions/functions with domain restrictions in R, but I'm having trouble defining such functions. For example, I would like to define the following function: f(x)={x^2 if -1xx; 1 if 2x3} Notably, the function is

[R] relative L1 bound

2012-03-26 Thread yx78
In the package lasso2, there is a Prostate Data. To find coefficients in the prostate cancer example we could impose L1 constraint on the parameters. code is: data(Prostate) p.mean - apply(Prostate, 5,mean) pros - sweep(Prostate, 5, p.mean, -) p.std - apply(pros, 5, var) pros - sweep(pros,

Re: [R] svycoxph and test statistics

2012-03-26 Thread Chirag Patel
Thank you, Dr. Therneau... I got a similar answer from Dr. Lumley: On Mar 24, 2012, at 4:05 PM, Thomas Lumley wrote: As far as I know there isn't any theoretical justification for the t-distribution but it empirically works better. You can get tests with a t or F reference distribution

Re: [R] relative L1 bound

2012-03-26 Thread yx78
by the way 0.44 is the optimal choice of lasso constrain. -- View this message in context: http://r.789695.n4.nabble.com/relative-L1-bound-tp4504620p4504678.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing

[R] reading header in txt file and making histogram

2012-03-26 Thread bamboohydraulics
Dear all I am a BEGINNER and have R on my Mac. I saved my excel file as .txt file, I have just one column with first row as the column name. My file when read by R looks like this. After reading the table I try to make a histogram by hist(dbh), it says object dbh not found. What am I doing

Re: [R] What does package 'RDCOMClient' is not installed for 'arch=x64' exactly mean?

2012-03-26 Thread Prof Brian Ripley
It 'exactly means' what it says. We have no idea who installed it from the sources, and you need to take that up with them (or yourself, if it were you). Please do follow the posting guide and give the 'at a minimum' information required. (Your R is seriously obsolete, and you were asked

Re: [R] Problem reading mixed CSV file

2012-03-26 Thread Ashish Agarwal
Need to fix up the file having 6 and 7 columns to be read as 6 columns only. Here is the working. Can somebody please let me know how do I maintain the order in which rows were read and append the two files into one: count.fields(textConnection(LL1532Ap,ABC# Depot-A+,,1971,8,2 LL1532Ap,ABC#

Re: [R] Handling 8GB .txt file in R?

2012-03-26 Thread Rainer M Krug
On 24/03/12 09:08, iliketurtles wrote: Hi, I am mediocre at R, maybe 1000 hours experience, but I received an 8GB dataset and I don't know what to do with it. I have to do extensive analysis over it for my Honours thesis. I can't even import it. I've tried; - Splitting it up using the

Re: [R] Simple question regarding domain restrictions/piecewise functions in R

2012-03-26 Thread Rolf Turner
Does f - function(x){ ifelse((-1 x x 1) | (2 x x 3),x^2,NA) } plot(f,xlim=c(-3,5)) give you what you want? cheers, Rolf Turner On 26/03/12 11:08, chad.mills wrote: I am a novice R user. I would like to be able to graph some simple

Re: [R] trellis plot

2012-03-26 Thread Peter Ehlers
On 2012-03-25 21:20, Sebastián Daza wrote: Hi everyone, I am just trying to figure out how to do a xyplot where in addition to dots and lines I can change dots' colors according to an individual variable (e.g., marital disruption across time, a dummy 0/1). When I use groups specification (see

[R] Enhance Pairs Scatterplot Matrix

2012-03-26 Thread Johannes Radinger
Hi, I am trying to do plot a scatterplot matrix using pairs() from the package graphics. There are now a few things I'd like to improve for a better understanding. 1) My scatterplot uses two log-scaled axis...as I don't know how to set them in pairs() I calculated the log before but now the

Re: [R] row, col function but for a list (probably very easy question, cannot seem to find it though)

2012-03-26 Thread Berend Hasselman
On 26-03-2012, at 00:17, MBoersma wrote: Hi guys, I'm quite new to R but quite enthousiastic. I'm trying to rewrite a bit of code in order to make it faster, so instead of nesting for loops I'm trying some apply functions. Since it's a combination of lists of matrices and other matrices,

Re: [R] R Error : DATA to MATRIX

2012-03-26 Thread Ashish Agarwal
On similar lines, how do I assign the position of the an occurence of a value say 7 in a data vector to another data vector? nFields [1] 6 6 6 6 6 6 7 7 6 6 6 7 6 6 6 6 I need the output data vector showing occurence of 7 as: [1] 7 8 12 I tried following but both has errors xFields7 - NULL for

[R] filter and merge 2 big Df's by=t_chr in sec.

2012-03-26 Thread Eugeniusz Kałuża
Please help to solve that problem: filter and merge 2 big Df's by=t_chr in sec. require(zoo) t1=2012-03-25 17:00:00.0 t2=2012-03-25 17:00:05.0 t1_POSIXlt-strptime(t1, %Y-%m-%d %H:%M:%OS)# t2_POSIXlt-strptime(t2, %Y-%m-%d %H:%M:%OS)# #

Re: [R] R Error : DATA to MATRIX

2012-03-26 Thread Berend Hasselman
On 26-03-2012, at 11:34, Ashish Agarwal wrote: On similar lines, how do I assign the position of the an occurence of a value say 7 in a data vector to another data vector? nFields [1] 6 6 6 6 6 6 7 7 6 6 6 7 6 6 6 6 I need the output data vector showing occurence of 7 as: [1] 7 8 12 I

[R] Error in par(pin = p) : invalid value

2012-03-26 Thread Aurelie Cosandey Godin
Dear list, I am getting an error message with pin and I am not quiet sure how to fix it. library(maps) op-par(mar=c(1,1,1,1)) map(xlim=c(-40,-45),ylim=c(60,65), boundary=T, fill=T, col='gray95') Error in par(pin = p) : invalid value specified for graphical parameter pin par(pin) par(pin)

Re: [R] Struggling with zoo and aggregate

2012-03-26 Thread Gabor Grothendieck
On Sun, Mar 25, 2012 at 10:20 PM, Thomas Adams thomas.ad...@noaa.gov wrote: Gabor, Thank you for your help -- it did help me a lot. However, with my data:    lead_time cycle    r_squared  fcst_date 1  6 0 5.405095e-02 07/31/2010 2 12 0 5.521620e-06 07/31/2010 3  

Re: [R] Error in par(pin = p) : invalid value

2012-03-26 Thread Duncan Murdoch
On 12-03-26 7:17 AM, Aurelie Cosandey Godin wrote: Dear list, I am getting an error message with pin and I am not quiet sure how to fix it. library(maps) op-par(mar=c(1,1,1,1)) map(xlim=c(-40,-45),ylim=c(60,65), boundary=T, fill=T, col='gray95') Error in par(pin = p) : invalid value

Re: [R] Error in par(pin = p) : invalid value

2012-03-26 Thread Aurelie Cosandey Godin
Should of seen that! It works now. Thank you Duncan, Aurelie On 2012-03-26, at 8:39 AM, Duncan Murdoch wrote: On 12-03-26 7:17 AM, Aurelie Cosandey Godin wrote: Dear list, I am getting an error message with pin and I am not quiet sure how to fix it. library(maps)

[R] circles()

2012-03-26 Thread John D. Muccigrosso
I cannot for the life of me figure this out: What's the parameter to fill in with color circles made with circles()? col changes the line color, but all I see in the help is a reference to additional graphic parameters, and no examples via google. Thanks! John Muccigrosso

[R] NA in R package randomForest

2012-03-26 Thread silje skår
I have a question regarding NA in randomForest (in R). I have a dataset which include both numerical and non-numerical variables, and the data includes some NA. I tried to use na.roughfix but then i get an error message na.roughfix only works for numeric or factor. I also tried rfImpute but this

Re: [R] NA in R package randomForest

2012-03-26 Thread Weidong Gu
Your data set probably has character variables. You can transform them into factors before running na.roughfix. Weidong Gu On Mon, Mar 26, 2012 at 8:17 AM, silje skår silje@gmail.com wrote: I have a question regarding NA in randomForest (in R). I have a dataset which include both numerical

Re: [R] Struggling with zoo and aggregate

2012-03-26 Thread Thomas Adams
Gabor, That does it! I can't thank you enough… Many thanks, Tom On Mon, Mar 26, 2012 at 7:22 AM, Gabor Grothendieck ggrothendi...@gmail.com wrote: On Sun, Mar 25, 2012 at 10:20 PM, Thomas Adams thomas.ad...@noaa.gov wrote: Gabor, Thank you for your help -- it did help me a lot.

[R] Dependence Hypothesis Test

2012-03-26 Thread Andrea Sica
Dear all, I have a data-set of 91 variables. It is related with a qualitative analysis I made to analyse a few phenomena. For this reason the questions are grouped in 10 clusters, since they express different phenomena of the analysis. I made the subset to collect the questions and represent the

Re: [R] circles()

2012-03-26 Thread Duncan Murdoch
On 26/03/2012 8:09 AM, John D. Muccigrosso wrote: I cannot for the life of me figure this out: What's the parameter to fill in with color circles made with circles()? col changes the line color, but all I see in the help is a reference to additional graphic parameters, and no examples via

Re: [R] how to cluster rows of words in a text file

2012-03-26 Thread Alekseiy Beloshitskiy
Hello, I didn't quite understand what you need, but maybe you can have a look here: www.slideshare.net/whitish/textmining-with-r R code fragments are in appendixes of the presentation. Hope this will help, -Alex From: r-help-boun...@r-project.org

[R] how to scale tokens

2012-03-26 Thread Alekseiy Beloshitskiy
Hi All, I need to scale variable of tokens (several words per observation) to be able to use it for svm(). For example, i have varible x4=(how,grow,tree) Any ideas how to scale to use in svm()? Thank you, -Alex [[alternative HTML version deleted]]

[R] SPSS R-Menu for Ordinal Factor Analysis

2012-03-26 Thread Till Below
Dear all, I am trying to conduct an enhanced version of factor analysis with a SPSS interface that allows to use R. This approach has been suggested in the recent article: Basto, M. and J.M. Pereira An SPSS R-Menu for Ordinal Factor Analysis. Journal of Statistical Software 46, pp. 1-29.

[R] Circular statistics- range.circular

2012-03-26 Thread Benedikt Niesterok
Hi everyone, does anyone know, which test of uniformity is run under the command range.circular() in the circular statistics package? It didn't say anything in the description of this command. Your help would be appreciated a lot. Benedikt [[alternative HTML version deleted]]

[R] copy the columns based on the code

2012-03-26 Thread MSousa
Hello, I have two different dataset, and wanted to join the two. For example I have a table of codes of cities, and other with with the codes of travels, [source for the destine]. what he wanted was to have a new data.frame with all the information city-data.frame(city=Barcelona,cod=1)

[R] SEM: Dependent binary: impact estimating wrong standard errors with hetcor()

2012-03-26 Thread Tryntsje Hoving-Wesselius
Hi, I'm using the SEM package to estimate a model with a binary variabele as dependent variable. In the literature I have to use then the correlation matrix, made by function hetcor(). Literature also says that the standard errors are not correct then. My question is if somebody knows what the

[R] R_unload_[PACKAGENAME] not called

2012-03-26 Thread Henrik Alsing Friberg
Hi R-help users.. I mistakenly posted a now rejected bug report (https://bugs.r-project.org/bugzilla3/show_bug.cgi?id=14855), but the subject still bugs me..! If I have a package called MyPackage with C code inside, then library or require(MyPackage) will call loadNamespace() which will look for

Re: [R] Seeming failure of options(width=60)

2012-03-26 Thread John C Nash
Yihui is correct that an error is generated. I should have mentioned that this is something I wanted to show in a vignette, and I got 103 characters as the length of the maximal line from a sink() file. Yihui's use of nchar is more elegant. Thus the question remains of how to avoid margin

Re: [R] copy the columns based on the code

2012-03-26 Thread Sarah Goslee
It isn't quite clear to me that cod in data frame city and pos in travel are are actually the same index, but if so then you can easily use merge() for this task. Sarah On Mon, Mar 26, 2012 at 5:41 AM, MSousa ricardosousa2...@clix.pt wrote:  Hello, I have two different dataset, and wanted to

Re: [R] reading header in txt file and making histogram

2012-03-26 Thread John Kane
Please provide your code as is requested in the posting guidelines. At a guess, you have read in your data into a data.frame called dhb with a variable called dhb. Try hist(dhb$dhb) and see what happens. To see a bit more of what you've got try str(dhb) to see what the sdtructure of the data

Re: [R] R Error : DATA to MATRIX

2012-03-26 Thread Petr PIKAL
On similar lines, how do I assign the position of the an occurence of a value say 7 in a data vector to another data vector? nFields [1] 6 6 6 6 6 6 7 7 6 6 6 7 6 6 6 6 I need the output data vector showing occurence of 7 as: [1] 7 8 12 which(nFields==7) Regards Petr I tried

Re: [R] R_unload_[PACKAGENAME] not called

2012-03-26 Thread Duncan Murdoch
On 26/03/2012 6:24 AM, Henrik Alsing Friberg wrote: Hi R-help users.. I mistakenly posted a now rejected bug report (https://bugs.r-project.org/bugzilla3/show_bug.cgi?id=14855), but the subject still bugs me..! If I have a package called MyPackage with C code inside, then library or

Re: [R] SPSS R-Menu for Ordinal Factor Analysis

2012-03-26 Thread Duncan Murdoch
On 26/03/2012 7:32 AM, Till Below wrote: Dear all, I am trying to conduct an enhanced version of factor analysis with a SPSS interface that allows to use R. This approach has been suggested in the recent article: Basto, M. and J.M. Pereira An SPSS R-Menu for Ordinal Factor Analysis. Journal

[R] many curves to the same plot

2012-03-26 Thread Alaios
Dear all, I would like to plot under the same plot many curves like mycdf-ecdf(runif(100)) plot(mycdf) curve(mycdf,from=-1,to=1) mycdf-ecdf(runif(200)+3) plot(mycdf) curve(mycdf,from=-1,to=1) as you can see each new call to curve redraws (erases) the window. How I can have an effect like

Re: [R] many curves to the same plot

2012-03-26 Thread R. Michael Weylandt
You can use lines() with ecdf objects or use the not-at-all-surprisingly-named add = TRUE parameter of curve(). Michael On Mon, Mar 26, 2012 at 10:44 AM, Alaios ala...@yahoo.com wrote: Dear all, I would like to plot under the same plot many curves like mycdf-ecdf(runif(100)) plot(mycdf)

Re: [R] SPSS R-Menu for Ordinal Factor Analysis

2012-03-26 Thread peter dalgaard
On Mar 26, 2012, at 16:31 , Duncan Murdoch wrote: I think you will need to ask this question in an SPSS help forum: - Version 2.10.x of R is quite old, and is no longer supported. - The error message is coming from SPSS, not R and actually, it doesn't look like an error message at all,

Re: [R] many curves to the same plot

2012-03-26 Thread Petr PIKAL
Hi If you went through help page you probably could find it yourself in shorter time then you spend writing email and waited for others to write back. see ?plot.stepfun Dear all, I would like to plot under the same plot many curves like mycdf-ecdf(runif(100)) plot(mycdf)

Re: [R] Seeming failure of options(width=60)

2012-03-26 Thread Duncan Murdoch
On 26/03/2012 9:30 AM, John C Nash wrote: Yihui is correct that an error is generated. I should have mentioned that this is something I wanted to show in a vignette, and I got 103 characters as the length of the maximal line from a sink() file. Yihui's use of nchar is more elegant. Thus the

Re: [R] SPSS R-Menu for Ordinal Factor Analysis

2012-03-26 Thread Eik Vettorazzi
Hi Till, you need a plug-in for SPSS to run R, which is bound to a fixed R version. SPSS 19 uses R2.10 and SPSS 20 uses R2.12 (both not up-to-date R versions). So you have to download this plugin (after registering on IBM Support) for your spss version and the according R version - or start with a

Re: [R] row, col function but for a list (probably very easy question, cannot seem to find it though)

2012-03-26 Thread David Winsemius
On Mar 25, 2012, at 6:17 PM, MBoersma wrote: Hi guys, I'm quite new to R but quite enthousiastic. I'm trying to rewrite a bit of code in order to make it faster, so instead of nesting for loops I'm trying some apply functions. It is an urban myth that using 'apply' functions will

Re: [R] R numerical integration

2012-03-26 Thread Hans W Borchers
casperyc casperyc at hotmail.co.uk writes: I don't know what is wrong with your Maple calculations, but I think you should check them carefully, because: (1) As Petr explained, the value of the integral will be 0.5 (2) The approach of Peter still works and returns : 0.4999777 (3) And the same

[R] Pareto frontier plots in three dimensions

2012-03-26 Thread Robbie Morrison
Hello all This is my first posting for some years. I am back using R again and must say I do like the language (regarding scripting, I also use matlab, perl, and bash). My question involves plotting a Pareto frontier in three dimensions. This is strictly a exercise in visualization, I make no

Re: [R] trellis plot

2012-03-26 Thread Sebastián Daza
Thank you Peter. The problem with you solution is that it doesn't represent the actual values of disruption. Look this example: person  - rep(1:2, each=4) income  - c(100, 120, 150, 200, 90, 100,120, 150) disruption  - c(0,0,0,1,0,1,1,0) time  - rep(c(1:4),2) dat  -

Re: [R] 'names' attribute must be the same length as the vector

2012-03-26 Thread Uwe Ligges
On 25.03.2012 22:10, jiefan wrote: I have met into this problem when I tried to run panel regression by plm. My code: library(plm) indus- read.csv(file=full.csv,header=TRUE) industry-as.data.frame(indus) reg-lm(LnTSO2 ~ LnPGDP + LnPGDP2 + LnSOES + LnCOES + LnLIMD + LnSHOLD + LnPRIV +

Re: [R] many curves to the same plot

2012-03-26 Thread Alaios
Thanks a lot :) That helped ! Alex From: R. Michael Weylandt michael.weyla...@gmail.com Cc: R help R-help@r-project.org Sent: Monday, March 26, 2012 4:48 PM Subject: Re: [R] many curves to the same plot You can use lines() with ecdf objects or use the

[R] Predefined set of gray scale colors

2012-03-26 Thread Alaios
Dear all, I am printing  7 different data sets to the same plots, if I use some colors like (red,blue,yelow...) the different datasets are easily recognized. I would like to use some gray scale colors for distinguising the datasets. Is there any ready set with grayscale colors (usually the

[R] normalization of multi-value string variable

2012-03-26 Thread Alekseiy Beloshitskiy
Hi All, I need to normalize/scale string variable which represents interests of customers (e.g., 'cycling, rollerblading, swimming' etc). Does anybody know how to do this, I want then use it along with other numeric variables for SVM classification. Appreciate for any advice. -Alex

Re: [R] Predefined set of gray scale colors

2012-03-26 Thread Sarah Goslee
You might start with ?gray Sarah On Mon, Mar 26, 2012 at 12:49 PM, Alaios ala...@yahoo.com wrote: Dear all, I am printing  7 different data sets to the same plots, if I use some colors like (red,blue,yelow...) the different datasets are easily recognized. I would like to use some gray

[R] writing output directly to file; sink?

2012-03-26 Thread Diann Prosser
Dear all, I am working with large matrices (19.6 million elements * 1000 simulations) and am trying to get around memory problems and vector length issues. I’ve split the inputs so that the output vector length will not exceed 2^31. Working on a 64bit machine with 80GB RAM, I still get close to

Re: [R] copy the columns based on the code

2012-03-26 Thread MSousa
Thanks for the reply. I am using the function you gave me. complete.travel-merge(travel, city, by.x = Source, by.y = cod, all = TRUE) complete.travel-merge(travel, city, by.x = Destine, by.y = cod, all = TRUE The problem is that it gives the result that I want The

Re: [R] reading header in txt file and making histogram

2012-03-26 Thread bamboohydraulics
this is what I do, IT READS THE TABLE, But then I am stuck, I removed the header read.table(nigeria slr misc/Nigeria India /DBHHISTOGRAM.txt, header=FALSE) hist(nigeria slr misc/Nigeria India /DBHHISTOGRAM.txt$V1) Error in nigeria slr misc/Nigeria India /DBHHISTOGRAM.txt$V1 : $ operator is

Re: [R] reading header in txt file and making histogram

2012-03-26 Thread jaclea
If there is a header in your .txt-file you should have header=TRUE, if there is no header you should have header=FALSE. You have to save your dataset in a variable, for example: myData-read.table(data.txt,header=TRUE) Then you can make a histogram by: hist(myData) -- View this message in

Re: [R] reading header in txt file and making histogram

2012-03-26 Thread chuck.01
You really need to read one of many beginners guides to R; but I will say this much: YOUR_DATA - read.table(nigeria slr misc/Nigeria India /DBHHISTOGRAM.txt, header=TRUE) then you should be able to plot it like I said previously. If that still doesn't work, you need to read the posting guide

Re: [R] Extracting numbers from a character variable of different types

2012-03-26 Thread irene
It worked perfectly! Thank you -- View this message in context: http://r.789695.n4.nabble.com/Extracting-numbers-from-a-character-variable-of-different-types-tp4482248p4505914.html Sent from the R help mailing list archive at Nabble.com. __

Re: [R] copy the columns based on the code

2012-03-26 Thread MSousa
thanks, is working -- View this message in context: http://r.789695.n4.nabble.com/copy-the-columns-based-on-the-code-tp4505253p4506128.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing list

[R] create waveform sawtooth

2012-03-26 Thread mail me
Hi: I am trying to create a sawtooth waveform. I used the following x - runif(500, min = -2, max = 2) y - (1 -abs(x3))* ((x3) = 1) combined - data.frame(x = x3, y = y3) plot(combined) and I get a triangular waveform, not sawtooth. Can someone give a solution to create a sawtooth waveform?

Re: [R] string substitution for argument in function

2012-03-26 Thread Prof. Dr. Pedro Martinez Arbizu
Me again, what I really dont understand is the behaivour of R here. notice, elem is not written in brackets, why it is interpreted as string and not as a variable? I would expect an error elem not found because it is not defined as variable, but in contrary it is accepted as column name, see

Re: [R] copy the columns based on the code

2012-03-26 Thread Rui Barradas
Hello, The problem is that it gives the result that I want Oh dear, it shouldn't? The idea is based on the column of source and intended Identify the cities and put a new data structure The idea is something like this. pos Source city Destine city_destine 1 1

[R] Building GLM Model with dependent variable as a proportion for 4 categories

2012-03-26 Thread prachi mehta
Hi All, I am trying to develop a GLM model in which the dependent data is a proportion data with 4 categories. I want to develop this model as an alternative method to Multinomial Logit Model for comparison of results. Here is the sample data: CategorySharesY x1 x2 1

Re: [R] copy the columns based on the code

2012-03-26 Thread R. Michael Weylandt
The problem is that it gives the result that I want I'm not sure that's what one would usually identify as a problemcan you say a little more specifically what you are looking to do? Michael On Mon, Mar 26, 2012 at 10:12 AM, MSousa ricardosousa2...@clix.pt wrote: Thanks for the reply.    

[R] How to find best parameter values using deSolve n optim() ?

2012-03-26 Thread mhimanshu
Hello, Can someone please help me out with the optim() function. # parameters pars- c(a1= 0.9, a2= 0.7, a3= 0.06, a4=0.02) y- c(Y=0.2, Z=0.1) Ymax- c(0.8) fucntion deriv derivs - function(time, y, pars) { with (as.list(c(y, pars)), { dy = a1*Y*(1-Y/Ymax) - a2*(Y+0.001) dz = a3*Y- a4*Z;

[R] confidence intervals for non-linear models

2012-03-26 Thread katem
Hi, Does anyone have any suggestions for plotting confidence intervals onto fitted non-linear models? Thank you for your help! Kate -- View this message in context: http://r.789695.n4.nabble.com/confidence-intervals-for-non-linear-models-tp4506095p4506095.html Sent from the R help mailing list

Re: [R] copy the columns based on the code

2012-03-26 Thread MSousa
thanks for the reply, But does not the results I need. What is confusing is that when making the first merge it gives only two cities when it should take three travel-data.frame(pos=1,Source=1,Destine=2) travel-rbind(travel,data.frame(pos=1,*Source=1*,Destine=3))

Re: [R] string substitution for argument in function

2012-03-26 Thread Prof. Dr. Pedro Martinez Arbizu
HI, thanks Weidong, Henrique, this works for the example (may be a bad example, but it should be as simple as possible), This is however only a workaround to name the columns of a dataframe, but is not addressing the problem itself, that is, how can I substitute a string (used as argument in

Re: [R] Simple question regarding domain restrictions/piecewise functions in R

2012-03-26 Thread chad.mills
Yes! Thanks. It was just the NA value instead of the as.null that does the trick. Correct code for the original piecewise I stated (for those who might be looking later) is: f - function(x){ ifelse((-1 x x 1),x^2,ifelse((2xx3),1,NA)) } plot(f,xlim=c(-1,3))

[R] Drape 3D scatterplot

2012-03-26 Thread John Benning
Hi, *Warning: R newb here. * I've got a 3D scatterplot ( http://www.fleetwoodswoodworks.com/scatterPlot.jpg) which I want to drape in color, with the draped plane corresponding to the mean z values on the plot (similar to: http://www.fleetwoodswoodworks.com/draped3Dscatterplot.jpg). Does anyone

[R] Year of data collection for 'diamonds' dataset in ggplot2

2012-03-26 Thread Marina Doucerain
Hello, I'm wondering what was the year (or year range) of collection for the data included in the 'diamonds' dataset in ggplot2. This information would be very helpful in interpreting the 'price' variable. Thank you! Marina Doucerain __

Re: [R] copy the columns based on the code

2012-03-26 Thread Sarah Goslee
The problem is that it gives the result that I want. That's a new sort of problem. You show two different merge() commands. What do you expect to happen, vs what does happen? Sarah On Mon, Mar 26, 2012 at 10:12 AM, MSousa ricardosousa2...@clix.pt wrote: Thanks for the reply.    I am using

Re: [R] create waveform sawtooth

2012-03-26 Thread R. Michael Weylandt
Perhaps something like this, x - seq(0, 10, length = 1000) y - x %%1 Michael On Mon, Mar 26, 2012 at 10:18 AM, mail me mailme...@googlemail.com wrote: Hi: I am trying to create a sawtooth waveform. I used the following x - runif(500, min = -2, max = 2) y - (1 -abs(x3))* ((x3) = 1)

Re: [R] string substitution for argument in function

2012-03-26 Thread R. Michael Weylandt
'elem' is an attribute (the name) of the 'data' variable (bad name for a variable), not something that is independently accessible. By very (exceptionally!) loose analogy, think of it as being something like a field/property of an object -- you can't get to it directly (and that's a good thing)

Re: [R] Updating a Markov Chain

2012-03-26 Thread ilai
On Sun, Mar 25, 2012 at 3:50 PM, stivi muhame...@wp.pl wrote: Hello, my question is if anyone has any good ideas how to create a Markov Chain from ordered data. So, I have some sort of time series, and if value1 happens as time1 and value2 happens at time2 I record this as an update to the

Re: [R] copy the columns based on the code

2012-03-26 Thread Bert Gunter
Fortunes candidate?! -- Bert On Mon, Mar 26, 2012 at 10:24 AM, Sarah Goslee sarah.gos...@gmail.com wrote: The OP wrote The problem is that it gives the result that I want. Sarah's reply: That's a new sort of problem. -- Bert Gunter Genentech Nonclinical Biostatistics Internal Contact

Re: [R] copy the columns based on the code

2012-03-26 Thread Sarah Goslee
Then there's something missing in what you tell us. # this is an awkward and inefficient way of constructing a data frame city-data.frame(city=Barcelona,cod=1) city-rbind(city,data.frame(city=Madrid,cod=2)) city-rbind(city,data.frame(city=Lisbon,cod=3))

[R] Completely Off Topic:Link to IOM report on use of -omics tests in clinical trials

2012-03-26 Thread Bert Gunter
Warning: This has little directly to do with R, although R and related tools (e.g. sweave and other reproducible research tools) have a natural role to play. The IOM report: http://www.iom.edu/Reports/2012/Evolution-of-Translational-Omics.aspx that arose out of the Duke Univ. genomics testing

Re: [R] Drape 3D scatterplot

2012-03-26 Thread Duncan Murdoch
On 26/03/2012 1:34 PM, John Benning wrote: Hi, *Warning: R newb here. * I've got a 3D scatterplot ( http://www.fleetwoodswoodworks.com/scatterPlot.jpg) which I want to drape in color, with the draped plane corresponding to the mean z values on the plot (similar to:

Re: [R] Predefined set of gray scale colors

2012-03-26 Thread Alaios
Great :) From: Sarah Goslee sarah.gos...@gmail.com Cc: R help R-help@r-project.org Sent: Monday, March 26, 2012 6:51 PM Subject: Re: [R] Predefined set of gray scale colors You might start with ?gray Sarah Dear all, I am printing  7 different data sets

[R] Different result with kruskal.test and post-hoc analysis with Nemenyi-Damico-Wolfe-Dunn test implemented in the help page for oneway_test in the coin package that uses multcomp

2012-03-26 Thread gianni lavaredo
Dear Researchers, Sorry for this email but I am not a statistician, and for this I have this problem to understand. Thanks in Advance for help and suggestions. Gianni I have 21 classes (00, 01, 02, 04, ,020) with different length. I did a kruskal wall test in R with the following code

  1   2   >