Re: [R] Deleting a column in a dataframe by name

2009-09-29 Thread milton ruser
May be this: FRAME - FRAME[-c(NAME1, NAME2)] or FRAME- subset(FRAME, select=(-NAME1, -NAME2)) bests milton On Tue, Sep 29, 2009 at 4:12 PM, Dennis Fisher fis...@plessthan.com wrote: Colleagues, Hopefully a simple problem: I want to delete a column with a known name from a dataframe. I

Re: [R] Deleting a column in a dataframe by name

2009-09-29 Thread milton ruser
x=runif(12) y=runif(12) w=runif(12) mydf-data.frame(cbind(x,y,w)) head(mydf) mydf-subset(mydf, select=c(-x,-w)) head(mydf) On Tue, Sep 29, 2009 at 4:27 PM, Rolf Turner r.tur...@auckland.ac.nzwrote: On 30/09/2009, at 9:15 AM, milton ruser wrote: May be this: FRAME - FRAME[-c(NAME1

Re: [R] Deleting a column in a dataframe by name

2009-09-29 Thread milton ruser
Peace on your heart... it is still Tuesday, and Friday is away... :-) On Tue, Sep 29, 2009 at 4:54 PM, Rolf Turner r.tur...@auckland.ac.nzwrote: On 30/09/2009, at 9:32 AM, milton ruser wrote: x=runif(12) y=runif(12) w=runif(12) mydf-data.frame(cbind(x,y,w)) head(mydf) mydf-subset

Re: [R] Can You Recommend Books for Linear Mixed Models in R

2009-10-01 Thread milton ruser
Hi Pat, I am happy with the book Mixed Effects Models and Extensions in Ecology with R by Zuur, Ieno, Walker, Saveliev and Smith. In fact I am reading chapter by chapter, and the reading is very digestible, as well as the examples are quite easy to understand and to be reference for our real

Re: [R] is that possible to graph 4 dimention plot

2009-10-04 Thread milton ruser
How about : google r graph gallery 4d ? abs On Sun, Oct 4, 2009 at 3:31 PM, Duncan Murdoch murd...@stats.uwo.ca wrote: On 04/10/2009 3:14 PM, gcheer3 wrote: Suppose there are 4 variables d is a function of a , b and c I want to know how a, b and c change will make d change It will be

Re: [R] Filtering on a dataframe- newbie question

2009-10-18 Thread milton ruser
Hi there, Try ?subset ?str may also be useful.. bests milton On Sun, Oct 18, 2009 at 11:10 PM, ANJAN PURKAYASTHA anjan.purkayas...@gmail.com wrote: Hi, newbie question. I have a data-frame with 3 named columns: Name, Obs1, Obs2. The Name column members are made of alphanumeric

[R] updating columns using other column as reference

2009-10-19 Thread milton ruser
Dear R-gurus, Just supose I have a dara.frame that looks like myDF-read.table(stdin(),head=T,sep=,) codID,namesp,k1,k2,k3,k4 1,spA,2,5,6,3 2,spB,4,5,4,6 3,spC,2,1,5,6 4,spC,5,4,3,2 5,spD,1,2,3,4 6,spE,2,4,3,1 I need to update the columns k1-k4 with the namesp, but considering the math between

Re: [R] updating columns using other column as reference

2009-10-19 Thread milton ruser
Amazing. Thanks a lot milton On Mon, Oct 19, 2009 at 12:55 PM, Henrique Dallazuanna www...@gmail.comwrote: Milton, Try this: myDF[grep(k, names(myDF))] - sapply(myDF[grep(k, names(myDF))], function(idx)myDF$namesp[idx]) On Mon, Oct 19, 2009 at 2:03 PM, milton ruser milton.ru

Re: [R] rotate levelplot

2009-10-24 Thread milton ruser
Hi Kang, Could you send a reproducible sample-code? Bests miltinho On Sat, Oct 24, 2009 at 11:32 PM, Kang Min ngokang...@gmail.com wrote: Hi all, I have a matrix with 5 rows and 10 columns, which represent the grids on a rectangular map. I used the code below to plot, but it gives me

Re: [R] x-y plot as an histogram

2009-10-29 Thread milton ruser
Hi Federico, Have you checked google R graph library? bests milton On Thu, Oct 29, 2009 at 3:12 PM, Federico Abascal fedeabas...@yahoo.eswrote: Hi, I am investigating a problem for which I found no solution. I have a matrix with two columns. I have plotted it as an x-y plot

Re: [R] counting frequencies across two columns

2009-11-01 Thread milton ruser
Hi Jason, As your example is not reproducible, may be something like: myFreq-data.frame(table(articleID, author)) if you want to know only those articles with 1 author, you can try subset(myFreq, Freq==1) or something like. bests milton On Sun, Nov 1, 2009 at 2:20 AM, Jason Priem

[R] trend thresholds on time series

2009-11-07 Thread milton ruser
Dear all, I have several response variables estimated from some simulations, and I would like to identify the thresholds for trend changes. Fro example, below I forced two different response behaviours and on x is time unit. x-1:1500 y-x/exp(x^0.2) smaller15-y[y15]

Re: [R] LINEAR MIXED EFFECT

2009-11-11 Thread milton ruser
Dear Ana Golveia, It is completelly impossible someone realise what kind or help you need or what is happening. I suggest you give a look on the posting guide, mainly that part about a minimum reproducible code with self explaining information, etc. Cheers milton On Wed, Nov 11, 2009 at 7:22

Re: [R] LINEAR MIXED EFFECT

2009-11-13 Thread milton ruser
. cheers milton On Fri, Nov 13, 2009 at 8:30 AM, ANARPCG a.gouvei...@imperial.ac.uk wrote: Douglas Bates-2 wrote: On Thu, Nov 12, 2009 at 10:14 AM, milton ruser milton.ru...@gmail.com wrote: Hi Ana, I am not quite sure if it is the problem, but if you call your data.frame as exp

Re: [R] refactoring in R

2009-11-13 Thread milton ruser
Hi Peng, If that information is preliminary, so I guess you have a more clear problem and may be you are able to state a minimally reproducible code/example with what you really need. Bests milton On Fri, Nov 13, 2009 at 9:30 PM, Peng Yu pengyu...@gmail.com wrote: I'm wondering if there

Re: [R] Presentation of data in Graphical format

2009-11-15 Thread milton ruser
Google R graph grallery Google R ggplot2 Google R lattice and good luck milton On Sun, Nov 15, 2009 at 7:48 AM, Sunita22 sunita...@gmail.com wrote: Hello My data contains following columns: 1st column: Posts (GM, Secretary, AM, Office Boy) 2nd Column: Dept (Finance, HR, ...) 3rd column:

Re: [R] sensitivity analysis, input factors

2010-04-13 Thread milton ruser
Hi Szalai I had used only src function, and on that case you need to have a vector with your Y variable, and a data-frame with all your X (i.e. explanatory) variables. I have interest on stay in touch with others that have been using sensitivity package! bests milton On Tue, Apr 13, 2010 at

Re: [R] rcolorbrewer

2010-05-03 Thread milton ruser
?rev On Mon, May 3, 2010 at 4:06 PM, Steve Hempell sthemp...@gmail.com wrote: Can you reverse the color scheme order in rcolorbrewer. If so, how? [[alternative HTML version deleted]] __ R-help@r-project.org mailing list

[R] avoiding reinstall already installed library

2010-05-18 Thread milton ruser
Dear R-experts, I am installing new libraries using install.packages(ggplot2,dependencies=T). But I perceive that many dependencies are already installed. As I am using a low-band internet, how can avoid reinstall installed libraries? cheers milton [[alternative HTML version deleted]]

Re: [R] Re: Change class factor to numeric

2010-05-18 Thread milton ruser
Try arcilla-as.numeric(as.character(clay)) best milton On Tue, May 18, 2010 at 12:36 PM, Arantzazu Blanco Bernardeau aramu...@hotmail.com wrote: sorry I had a mistake sending my question without a subject. I do resend again. Please excuse me. Hello I have a data array with soil

Re: [R] avoiding reinstall already installed *package*

2010-05-18 Thread milton ruser
will not be reinstalled. cheers milton On Tue, May 18, 2010 at 12:49 PM, Martin Maechler maech...@stat.math.ethz.ch wrote: { I've modified the subject; I can't stand it hitting square into my face ... } mr == milton ruser milton.ru...@gmail.com on Tue, 18 May 2010 12:36:23 -0300 writes

Re: [R] avoiding reinstall already installed *package*

2010-05-18 Thread milton ruser
:25 PM To: milton ruser Cc: r-help@r-project.org Subject: Re: [R] avoiding reinstall already installed *package* On Tue, May 18, 2010 at 18:06, milton ruser milton.ru...@gmail.com wrote: Hi Martin, thanks for your reply, and very thanks for your kind tips about package

[R] (ecological) model sensitivity analysis book

2009-11-20 Thread milton ruser
Dear all, I know that my question may be so general, but I would like to know if is there a good book for model sensitivity analysis (applied to ecology?) that use R? Is there any package dedicated to model sensitivity analysis on R? Cheers milton [[alternative HTML version deleted]]

Re: [R] 11 distinguishable colors

2009-11-24 Thread milton ruser
Hi there, I picked it up from r.colors' grass manual: http://colorbrewer2.org/ The code bellow cold be usefull for you explore simbols and colors. x11(900,500) par(mfrow=c(1,2)) x-runif(20) y-runif(20) plot(y~x, type=n) Number.of.symbols-20 for (i in 1:Number.of.symbols) { points(y[i]~x[i],

Re: [R] Arrhenius plot with lattice

2009-11-25 Thread milton ruser
Hi Markus, PLEASE do read the posting guide http://www.R-project.org/posting-guide.htmlhttp://www.r-project.org/posting-guide.html Your example is not reproducible. Help your helper :-) bests miltinho On Wed, Nov 25, 2009 at 5:54 PM, Markus Häge markus_ha...@gmx.de wrote: hello there, I

Re: [R] Ggplot2: Elegant Graphics for Data Analysis (Use R) (2009 Paperback)

2009-11-30 Thread milton ruser
I agree with Ista. The book is very interesting. I put my 2cents too. milton On Mon, Nov 30, 2009 at 6:20 PM, Ista Zahn istaz...@gmail.com wrote: Depends on what you want. It's possible to find most of the information contained in the book in other places (the ggplto2 website, Hadley's

Re: [R] Normal tests disagree?

2009-12-01 Thread milton ruser
Hi there, Could you provide a minimum reproducible code, please. Bests milton On Tue, Dec 1, 2009 at 6:11 PM, rkevinbur...@charter.net wrote: If I have data that I feed into shapio.test and jarque.bera.test yet they seem to disagree. What do I use for a decision? For my data set I have

Re: [R] sort a data frame by a vector

2009-12-01 Thread milton ruser
not ellegant.. but... vecDF = data.frame(A1=c(C, A, B)) vecDF$A1.order=1:dim(vecDF) vecDF dataDF = data.frame(A1 = c(B, A, C), A2 = c(1,2,3)) dataDF2-merge(vecDF, dataDF, by=intersect(colnames(vecDF),colnames(dataDF))) dataDF2 dataDF2.ord-dataDF2[order(dataDF2$A1.order),] dataDF2.ord milton On

Re: [R] how to seperate a matrix

2009-12-03 Thread milton ruser
Hi there, mymat-matrix(c(6,8,1, 5,9,2, 20,10,3, 7,11,4, 8,12,5, 25,13,6, 14,14,7), byrow=T, ncol=3) mymat1-mymat[mymat[,1]mymat[,2],] mymat1 mymat2-mymat[mymat[,1]=mymat[,2],] mymat2 bests milton On Thu, Dec 3, 2009 at 10:12 PM, aegea gche...@gmail.com wrote: Hello, I am working on

Re: [R] length of a density curve (or any curve)

2009-12-04 Thread milton ruser
hi Sylvain, did you try ?density regards milton On Fri, Dec 4, 2009 at 7:19 AM, sylvain willart sylvain.will...@gmail.comwrote: Hello R users, When I type d - density(MyData$x) I obtain a density object I can plot, But I wonder if there is a way to easily compute the length of the

Re: [R] 4D plotting

2009-12-06 Thread milton ruser
Hi Jared, google R graph library 4d or http://addictedtor.free.fr/graphiques/RGraphGallery.php?graph=90 bests milton On Sun, Dec 6, 2009 at 7:31 AM, Jared Nance nanc...@phys.washington.eduwrote: Hello list, Thanks in advance for whatever help you can give me. I have a data set that I want

Re: [R] hi there!

2009-12-06 Thread milton ruser
Hi Dave, May be a good idea you (1) use subjects avoiding help, Hi there, problem, etc. The list is archived for future search, and nobody with problems with MCMC, Metropholis Hasting algorithm will search Hi there; (2) how about minimum reproducible problem? If you are needing use this kind of

Re: [R] Population Histogram

2009-12-09 Thread milton ruser
1. Read the posting guide: http://www.R-project.org/posting-guide.htmlhttp://www.r-project.org/posting-guide.html 2. Did you installed R? 2.1. Case yes, go to the help, click on Manual in PDF bests miltinho On Wed, Dec 9, 2009 at 1:47 PM, terry johnson terry.johnson@gmail.comwrote:

Re: [R] how to creat a matrix

2009-12-11 Thread milton ruser
Hi Him. Did you read the Help/Manual in PDF/A Instroduction to R.PDF? MyVect-scan() 0 0 0 0 1 1 0 0 0 0 0 1 0 0 0 0 0 1 0 0 0 0 0 1 0 MyVect MyMat-matrix(MyVect, ncol=5, byrow=T) MyMat bests milton On Fri, Dec 11, 2009 at 10:47 AM, Moohwan Kim kmhl...@gmail.com wrote: Dear R family I

Re: [R] help with graphing -- Points in my graph are not apparent, always displayed in steps

2009-12-12 Thread milton ruser
Hi Philip, I must confess that I not understood what is the problem. Could you clarify it a little bit more? Cheers miltinho brazil=toronto On Sat, Dec 12, 2009 at 12:39 AM, philip robinson robin...@students.wwu.edu wrote: I am trying to graphically represent a large set of data who's

Re: [R] how to import data from excel to R

2009-12-17 Thread milton ruser
Hi Sarath, try C:\\Documents and Settings\\admin\\Desktop\\data.txt or C:/Documents and Settings/admin/Desktop/data.txt or yet setwd(your path) mydata-read.table(data.txt) cheers miltinho On Thu, Dec 17, 2009 at 12:24 PM, sta_2...@spectrum.net.in wrote: Hi, I am using R and I want to

Re: [R] write.csv and col.names=F

2009-12-17 Thread milton ruser
?write.table On Thu, Dec 17, 2009 at 6:52 PM, kayj kjaj...@yahoo.com wrote: Hi All, I always have a problem with write.csv when I want the column names to be ignored, when I specify col.names=F, I get a header of V1 V2 V3 V4 etc. for example I tried write.csv(mydata, file=data.csv,

Re: [R] Help,Suggest me some methods to identify tr aining set and test set!!!

2009-12-21 Thread milton ruser
Not ellegant.. but... MyDF-data.frame(cbind(x=runif(10), y=rnorm(10))) TrainingSize=5 TrainingSize_list-sample(1:nrow(MyDF))[1:TrainingSize] TrainingSize_list MyDF.training-MyDF[(1:nrow(MyDF) %in% TrainingSize_list),] MyDF.training MyDF.test-MyDF[ ! (1:nrow(MyDF) %in% TrainingSize_list),]

[R] removing rows with NAs in anywere

2009-12-22 Thread milton ruser
Dear all, I have a data.frame with some NAs that can happens anywere, and I would like to keep only rows without NAs. Thanks in advance for your help, and Merry Xmas. myvect-runif(100) myvect[sample(1:100)[1:5]]-NA myDF-data.frame(matrix(myvect,ncol=5)) myDF miltinho [[alternative HTML

[R] prcomp : plotting only explanatory axis arrows

2009-12-23 Thread milton ruser
Dear all, I have a very large dataset (1712351 , 20) and would like to plot only the arrows that represent the contribution of each variables. On the sample below I woild like to plot only the explanatory variables (Murder, Assault..) and not the sites. prcomp(USArrests) # inappropriate

Re: [R] aggregate binary response data

2009-12-24 Thread milton ruser
?aggregate ? On Thu, Dec 24, 2009 at 7:56 AM, Graham Leask grahamle...@btinternet.comwrote: Dear list I have a response variable coded 0/1 i.e. a binary response. There are 20,000 individual responses that I would like to aggregate into numbers of each category (i.e. 0/1) by group called

Re: [R] aggregate binary response data

2009-12-24 Thread milton ruser
may be ?table ? On Thu, Dec 24, 2009 at 7:56 AM, Graham Leask grahamle...@btinternet.comwrote: Dear list I have a response variable coded 0/1 i.e. a binary response. There are 20,000 individual responses that I would like to aggregate into numbers of each category (i.e. 0/1) by group

Re: [R] help in merging

2009-12-24 Thread milton ruser
Hi there, You can add a order column on x or y and after use this field to order z. Like z-z[order(z$orderfield),] To generate a order on x or y you can do something like x$xorder-1:nrow(x) cheers milton On Thu, Dec 24, 2009 at 2:26 PM, utkarsh.sing...@global-analytics.comwrote: Hi All,

Re: [R] help: creating a unified histogram

2009-12-27 Thread milton ruser
Hi Yael, Have you tried to find it on google R graph gallery ? bests milton On Sun, Dec 27, 2009 at 2:45 PM, Yael Yohai yael.yo...@gmail.com wrote: Good evening, I would like to put the histograms of several matrices on the same graph, together on a unified histogram (instead of having

Re: [R] 2D array of strings

2009-12-28 Thread milton ruser
Hi Francesco, Be carefull with create a object named str, because you crash str() function. I don't know if it have implications on any package or functions. bests milton On Mon, Dec 28, 2009 at 4:32 PM, Francesco Napolitano franap...@gmail.comwrote: Sorry for the dumb question, but I

Re: [R] (no subject)

2009-12-29 Thread milton ruser
Hi Hao, I suggest you try again, starting by read posting guide at footnote of this email. How about a title for the message? How about identify yourself? bests milton On Tue, Dec 29, 2009 at 2:59 PM, Hao Cen h...@andrew.cmu.edu wrote: Hi, I wonder how to pass several functions and their

Re: [R] pass functions and arguments to function

2009-12-29 Thread milton ruser
Hi Hao, Ok. Sorry for my last post. bests milton On Tue, Dec 29, 2009 at 3:08 PM, Hao Cen h...@andrew.cmu.edu wrote: Hi, I wonder how to pass several functions and their arguments as arguments to a function. For example, the main function is f = function(X ) { process1(X) ...

Re: [R] Negbin Error Warnings

2009-12-30 Thread milton ruser
Dear Clara, could you, please, share the summary() of your data. bests milton 2009/12/30 Clara Brück clara_bru...@web.de Hi, I ran a negative binomial regression (NBR) using the Zelig-package and the negbin model. When I then try to use the simumlation approach using the setx () and

[R] Fwd: Negbin Error Warnings

2009-12-30 Thread milton ruser
Error Warnings To: milton ruser milton.ru...@gmail.com Dear Milton, Thanks for your email. I ran the NBR-model and then tried to compute the first differences of some expected values: m3.x1 -setx(m3, advisory=1) m3.x2 -setx(m3, advisory=0) m3.s1 -sim (m3,x=m3.x1, x1=m3.x2) However, I then get

Re: [R] joining two lists...

2009-12-31 Thread milton ruser
Dear Steve, If I undertood you have a list of data.frame?. Have you tried: NewData-rbind(data.frame(myData[1]),data.frame(myData[2])) Please, send us a str(object) so the experts can help you bests milton On Thu, Dec 31, 2009 at 1:05 AM, stevemew steve@alphaedge.org wrote: My R

Re: [R] How x[, 'colname1'] is implemented?

2009-12-31 Thread milton ruser
Hi Peng, If I undertood your point, try this: x-runif(10) y-runif(10) z-runif(10) w-runif(10) myDF-data.frame(cbind(x,y,z,w)) myDF myDF[,c(w,z)] Happy new year miltinho On Thu, Dec 31, 2009 at 5:15 PM, Peng Yu pengyu...@gmail.com wrote: I don't see where describes the implementation of

Re: [R] How to use read.table with Hebrew column names ?

2010-01-01 Thread milton ruser
Hi Tal, You want the colnames in hebrew? Case not, may be you can use colnames() to attrib new colnames. bests milton On Fri, Jan 1, 2010 at 4:47 PM, Tal Galili tal.gal...@gmail.com wrote: Hello dear R help group, I am trying to read a .txt file, with Hebrew column names, while keeping

Re: [R] How to not to terminate read.table if the input file is empty?

2010-01-01 Thread milton ruser
?try() On Fri, Jan 1, 2010 at 3:41 PM, Peng Yu pengyu...@gmail.com wrote: read.table terminates the program if the input file is empty. Is there way to let the program continue and return me a NULL instead of terminating the program? $ Rscript read_empty.R read.table(empty_data.txt)

Re: [R] How to calculate density function of Bivariate binomial distribution

2010-01-01 Thread milton ruser
Hi Nykee, I checked out ??bivariate on my R installed libraries, and found about a hundred of occurrences of bivariate. If nobody reply you with a streigth answer, give a look at: VGAM; splancs; sm; prada; geoR; fUtilities; fBasics; akima; adehabitat; vegan; ade4 and spatstat packages. bests

Re: [R] Boxplot- input the median point and the median value

2010-01-01 Thread milton ruser
Dear Elaine, The posting guide is http://www.R-project.org/posting-guide.htmlhttp://www.r-project.org/posting-guide.html bests milton On Fri, Jan 1, 2010 at 8:53 PM, elaine kuo elaine.kuo...@gmail.com wrote: Dear, I am a newbie to R. Now I am learning to draw boxplot using graphics().

Re: [R] Question on Reduce + rollmean

2010-01-02 Thread milton ruser
Dear M.Rahiz, Unfortunatelly I can't reproduce your example. PLEASE do read the posting guide http://www.R-project.org/posting-guide.htmlhttp://www.r-project.org/posting-guide.html bests milton On Sat, Jan 2, 2010 at 9:26 AM, Muhammad Rahiz muhammad.ra...@ouce.ox.ac.uk wrote: Hello useRs,

Re: [R] Artificial Neural Networks

2010-01-03 Thread milton ruser
Hi Alex, I am not a NN user, but I would suggestion you give more information about what are you planning to do. May be it helps the helper. Happy New Year. miltinho On Sun, Jan 3, 2010 at 3:54 PM, Alex Olafson alex.olaf...@yahoo.com wrote: Hi! I am studying to use some R libraries which are

Re: [R] function in aggregate applied to specific columns only

2010-01-03 Thread milton ruser
You want this? basicSubMean - aggregate(basicSub[c(score)], by=list(basicSub$student), FUN=mean, na.rm=TRUE) basicSubMean Group.1 score 1 1 55.0 2 2 60.0 3 3 67.5 bests milton On Sun, Jan 3, 2010 at 10:46 PM, david hilton shanabrook dhsha...@acad.umass.edu wrote: I

Re: [R] how to draw abline correctly?

2010-01-03 Thread milton ruser
Hi there, May be part of your problem is solved. size=c(80, 40, 20, 16, 10, 8, 4, 2) estimate=c(0.1353, 0.1337, 0.1237, 0.1164, 0.1058, 0.0983, 0.0759 , 0.0607) plot(size, estimate, type=o, lwd=1.5, col=blue, ylim=c(0, 0.2), xlim=c(0,80), pch=19, axes=F) abline(h=0.1361, col=red, lwd=1.5)

Re: [R] how to draw abline correctly?

2010-01-04 Thread milton ruser
You can also try: axis(1) axis(2, las=1) axis(2,at=c(0.1361), las=1) bests milton On Sun, Jan 3, 2010 at 9:26 PM, aegea gche...@gmail.com wrote: Hello, I am frastruated with this graph, just cannot get what I need. Thank you for any suggestions or help. I really appreciate it. I wrote

Re: [R] Fwd: string

2010-01-06 Thread milton ruser
Hi Robert, I use something like this with grass-R: myvar-paste(questionspan style=,color: #FF; font-weight: bold; font-size: 14pt,,sep=\) myvar May be it works for you. bests milton On Wed, Jan 6, 2010 at 6:22 AM, robert-mcfad...@o2.pl wrote: I would like to assign a variable y the

Re: [R] February 2010***New R Courses*** by XLSolutions Corp at 9 USA Cities: San Francisco, New York, Washington DC, Houston, Boston, Las Vegas, Seattle, etc

2010-01-07 Thread milton ruser
Dear Sue, Happy new your to you too. You will offer the course on R or S-Plus? If S-plus is the case, how deal with license? bests milton On Wed, Jan 6, 2010 at 3:24 PM, Sue Turner s...@xlsolutions-corp.com wrote: Happy New Year ! XLSolutions February 2010 R courses schedule is now

Re: [R] Linear regression - R^2

2010-01-09 Thread milton ruser
Dear Csanad, PLEASE do read the posting guide http://www.R-project.org/posting-guide.htmlhttp://www.r-project.org/posting-guide.html Reproducible code/problem are welcome. By the way, who you believe is right? bests milton On Sat, Jan 9, 2010 at 5:40 PM, Csanád Bertók

[R] xls.getshlib() error in xlsReadWrite package under W7

2010-07-06 Thread milton ruser
Dear all, I installed xlsReadWrite package under Windows 7, and got error when trying to run xls.getshlib(). Please, see the errors below: require(xlsReadWrite) Loading required package: xlsReadWrite xlsReadWrite version (cran shlib) Copyright (C) 2010 Hans-Peter Suter, Treetron, Switzerland.

Re: [R] (no subject)

2010-01-15 Thread milton ruser
Saleem, I have no idea about this, but may be you can use other data.frame name instead of data because data is a name of a function (see ?data). good luck milton On Fri, Jan 15, 2010 at 3:33 PM, saleem mukhtar saleem...@yahoo.com wrote: Hello, I am running R on cygwin for windows. File

Re: [R] How to detect and exclude outliers in R?

2010-01-19 Thread milton ruser
Hi V.S., Did you search first on r-repositories about this issue prior to ask? May be not. RSiteSearch(outliers) bests milton On Tue, Jan 19, 2010 at 1:08 AM, vikrant vikrant.shi...@tcs.com wrote: Suppose I am reading data from a file and the data contains some outliers. I want to know

Re: [R] how to open excel 2007 (.xlsx) file in R

2010-01-19 Thread milton ruser
Dear V.K., 1. http://www.R-project.org/posting-guide.htmlhttp://www.r-project.org/posting-guide.html 2. it is - in general - good someone that write to the list identify your self. It is a polite way of participate of the list. bests milton On Tue, Jan 19, 2010 at 12:15 AM, vikrant

Re: [R] concatenate 2 column vectors

2010-01-20 Thread milton ruser
Hi Magali, I confess that I still not figured out what you want, but try this: datatry2=transform(datatry,DIS=as.numeric(as.character(paste(datatry[,1],datatry[,2],sep= bests milton On Wed, Jan 20, 2010 at 11:14 PM, teurlai magali teurlai...@yahoo.frwrote: Sorry I forgot the subject in

Re: [R] 3d trend line

2010-01-21 Thread milton ruser
Hi mlcarte3 (it is you?) ?contour Try algo google R graph gallery bests milton On Thu, Jan 21, 2010 at 2:43 PM, mlcarte3 mlcar...@ncsu.edu wrote: I have a 3d scatter plot and I need to add a trend line. I can add a regression plane, but I need a line, not a plane. Thanks. -- View this

Re: [R] What is the difference between S and R?

2010-01-21 Thread milton ruser
Hi There, If you are starting with S or R, don't be worried. If you just started (I suppose with R), don't be worried too. Just enjoy and be happy :-) bests milton On Thu, Jan 21, 2010 at 10:44 AM, Peng Yu pengyu...@gmail.com wrote: On Thu, Jan 21, 2010 at 9:07 AM, Giovanni Petris

Re: [R] How to write '' to a csv with the default setting of write.csv?

2010-01-21 Thread milton ruser
Hi Peng, May be write.table(...,quote=F) bests milton On Thu, Jan 21, 2010 at 6:29 PM, Peng Yu pengyu...@gmail.com wrote: Please see the following example. I can not write '' to a csv file successfully. Could somebody let me if it is possible to write '' to a csv file with the default

Re: [R] How to write '' to a csv with the default setting of write.csv?

2010-01-21 Thread milton ruser
Peng, Are you having doble quotes? It is what you want? Did you tried the suggestion? bests milton On Thu, Jan 21, 2010 at 6:38 PM, Peng Yu pengyu...@gmail.com wrote: On Thu, Jan 21, 2010 at 5:32 PM, milton ruser milton.ru...@gmail.com wrote: Hi Peng, May be write.table(...,quote=F

Re: [R] Logistic regression

2010-01-21 Thread milton ruser
Dear someone, Take a look at drc and Hmisc/Design packages. Next time could be polite you identify yourself bests milton On Thu, Jan 21, 2010 at 4:17 PM, DispersionMap frenc...@btinternet.comwrote: can you do Logistic regression in R, if so how do you do it and how do you test the fit of a

Re: [R] R Output and ArcGIS

2010-01-25 Thread milton ruser
Hi Steve, As I have no problem with this because I run windows, I prefere to use .txt files. But may be on your case you install xlsReadWrite (or WriteRead), because .xls file works fine for ArcGis. Remember not use complex (like having dots) or large names because some times you get error on

Re: [R] How to choose seed in set.seed(seed, ...)?

2010-02-03 Thread milton ruser
Hi Michael, Experts will answer you better than me, but IMHO there is no best practices. In general you deal with set.seed() when you want that others (or yourself) reproduce the same sequence of random numbers generation as you run on a moment. Best regards milton On Wed, Feb 3, 2010 at

Re: [R] release memory

2010-02-08 Thread milton ruser
Try gc() best milton On Mon, Feb 8, 2010 at 5:16 PM, dys0022 dys0...@yahoo.com wrote: I am a new R user with the latest Ubuntu release. My executions consume large amount of memory (up to 1.6 GB). When I try to release the memory using rm(list=ls()), R still occupies 1.6GB. (I also tried

Re: [R] Separating columns, and sorting by rows

2010-02-14 Thread milton ruser
Hi Raging Jim may be this is a starting point. myDF-read.table(stdin(),head=T,sep=,) mm,Rainfall 1977-02,17.4 1977-03,34.0 1977-04,26.2 1977-05,42.6 1977-06,58.6 1977-07,23.2 1977-08,26.8 1977-09,48.4 1977-10,47.0 1977-11,37.2 1977-12,15.0 1978-01,2.6 1978-02,6.8 1978-03,9.0 1978-04,46.6

[R] parsing strings between [ ] in columns

2010-02-18 Thread milton ruser
Dear all, I have a data.frame with a column like the x shown below myDF-data.frame(cbind(x=c([[1, 0, 0], [0, 1]], [[1, 1, 0], [0, 1]],[[1, 0, 0], [1, 1]], [[0, 0, 1], [0, 1]]))) myDF x 1 [[1, 0, 0], [0, 1]] 2 [[1, 1, 0], [0, 1]] 3 [[1, 0, 0], [1, 1]] 4 [[0, 0, 1], [0,

Re: [R] [R-sig-Geo] Comparing spatial point patterns - Syrjala test

2008-02-09 Thread milton ruser
Dear Jean, I have no idea of how to solve this issue, but I suggest you write to the authors of spatstat package. I think they could help you very much. Another thing is that there was a threhead on this many times ago. May be that Alan Swanson cold also help you.

Re: [R] Length problem

2008-02-11 Thread milton ruser
[1] 0.9941125 1.7931375NA 1.0840625 1.1558000 length(coppie[1]) [1] 1 1 # Why ?? Thank you for your help Paolo Italia milton ruser wrote: Ciao Paolo, How about you show some row of your data? How many columns have your data.frame? One

Re: [R] Length problem

2008-02-11 Thread milton ruser
Ciao Paolo, How about you show some row of your data? How many columns have your data.frame? One? By the way data is not a so good name for your data frame. We will be very happy to help you Kindly, Miltinho Brasile On 2/11/08, Paolo Grillo [EMAIL PROTECTED] wrote: Hi all I have this

Re: [R] How to run one-way anova R?

2008-02-12 Thread milton ruser
Hi Kes, May be: my.df-data.frame(FT=rep(c(Fa,Fb),10),Y=runif(20)) aov(Y ~ FT, data = my.df ) # fixed effects Kind regards, miltinho On 2/12/08, Kes Knave [EMAIL PROTECTED] wrote: Dear all, How do I run a basic one-way anova in R? Regards Kes [[alternative HTML version deleted]]

[R] changing format of y-axis (barplot)

2008-02-25 Thread milton ruser
Dear all, I have some data.frame and I would like (1) change the format of y=axis to escape of scientific format (like 1e+05...) and (2) label the value or other variable on the top of the bar. df-data.frame(x=1:6,y=c(100,500,1000,25000,50,100)) barplot(df$y,names.arg=df$x) Thanks in

[R] plot y1 and y2 on one graph

2008-02-27 Thread milton ruser
Dear all I have a code like x-1:10 y1-x+runif(10)*2 y2-seq(0,50,length.out=10)+rnorm(10)*10 par(mfrow=c(1,2)) plot(y1~x) plot(y2~x) Now I would like to plot y1 and y2 on the same graph, with its two scales (y1 on left and y2 on rigth side). Any help are welcome. Kind regards Miltinho

Re: [R] plot y1 and y2 on one graph

2008-02-27 Thread milton ruser
Thanks to all, I reached up to my needs. miltinho Brazil. On 2/27/08, Gabor Grothendieck [EMAIL PROTECTED] wrote: There is an example in library(zoo) example(plot.zoo) On Wed, Feb 27, 2008 at 5:05 PM, milton ruser [EMAIL PROTECTED] wrote: Dear all I have a code like x-1:10

Re: [R] (performance) time in Windows vs Linux

2009-06-26 Thread milton ruser
Hi there, I have both systems on a DELL 64bit machine. I compiled R 2.9.0 on both systems, to get 64bits capability. Surpriselly, on Linux (Ubuntu with I installed 3 month ago) I spent 41s to run the same test you did, and less time (35s) under Vista. In fact I had noticed that I not have gained

Re: [R] applying a function to a pair of components for each row of a list

2009-06-27 Thread milton ruser
Hi Kavitha, I must confess you that I not understood well what you are looking for. But.. mylist-list(x=1:25, y=runif(n=25)) plot(mylist, type=n) points(mylist, type=p, col=mylist$x) Hth, miltinho On Sun, Jun 28, 2009 at 1:03 AM, Kavitha Venkatesan kavitha.venkate...@gmail.com wrote: Hi,

Re: [R] (performance) time in Windows vs Linux

2009-06-29 Thread milton ruser
Really? In fact I have a quadcore. But how can I know if Linux are really using only one core, and how can I setup it to use the 4cores? Thanks a lot, milton On Mon, Jun 29, 2009 at 1:46 AM, Patrick Connolly p_conno...@slingshot.co.nz wrote: On Fri, 26-Jun-2009 at 04:37PM -0400, milton

Re: [R] How to use subset in lm function

2009-06-29 Thread milton ruser
Like this? x-1:100 y-runif(100) df-data.frame(cbind(x,y)) mycoef.list-NULL for (i in 10:65) { mycoef-coef(lm(y~x, data=subset(df, x %in% 1:i)))[2] mycoef.list -c(mycoef.list, mycoef) } hist(mycoef.list) regards milton brazil=toronto On Mon, Jun 29, 2009 at 2:50 PM, Hongwei Dong

Re: [R] How do I change which R Graphics Device is active?

2009-06-30 Thread milton ruser
How about: dev.cur() dev.list() dev.next(which = dev.cur()) dev.prev(which = dev.cur()) dev.off(which = dev.cur()) dev.set(which = dev.next()) dev.new(...) graphics.off() dev.cur() dev.list() dev.next(which = dev.cur()) dev.prev(which = dev.cur()) dev.off(which = dev.cur()) dev.set(which =

Re: [R] (no subject)

2009-07-02 Thread milton ruser
Hi Andriy Fetsun, How about include a subject title on the subject of the email like length funcion not work. By the say, send us the result of str(your_object). Bests milton brazil=toronto On Thu, Jul 2, 2009 at 9:01 AM, Andriy Fetsun fet...@googlemail.com wrote: Hi Guys, It is very

Re: [R] output too large to display all

2009-07-06 Thread milton ruser
Hi Christine, How about try again? Cheers milton On Mon, Jul 6, 2009 at 11:34 AM, Christine Griffiths christine.griffi...@bristol.ac.uk wrote: Hi R Users, Hopefully a very simple solution, but I am stumped nevertheless. I am running glmer in which the output is too large so that not

Re: [R] How to separate the string?

2009-07-07 Thread milton ruser
hi hema may be strsplit can help on the job. bests. milton On Tue, Jul 7, 2009 at 3:54 AM, Hemavathi Ramulu hema.ram...@gmail.comwrote: Hi everyone, Hi want to separate the string(column1) for example column1 column2 column3 column4 column5 column6 bear b e

Re: [R] how to transpose a dataframe

2009-07-09 Thread milton ruser
Hi there, ?reshape do the job. bests milton On Thu, Jul 9, 2009 at 9:31 PM, roachyang roachy...@gmail.com wrote: I want to transpose a dataframe like level 2006 2007 2008 A B . C into levelyear A

Re: [R] how to transpose a dataframe

2009-07-10 Thread milton ruser
Hi roachyang, you need to be more precise on your email. But believe that you tryed one of the funtions like reshape or transform and the output has been printed on your R session. If your dataset are large, is better you assign your output on a object - like myobj-reshape(...) - and then save

Re: [R] drawing hmms

2009-07-10 Thread milton ruser
Hi Rajesh, As you was little clean on your asking, try give a look at one of those pages. http://finzi.psych.upenn.edu/R/library/msm/html/msm.html http://finzi.psych.upenn.edu/R/library/msm/html/hmm-dists.html http://finzi.psych.upenn.edu/R/library/mhsmm/html/hsmmspec.html

Re: [R] drawing hmms

2009-07-10 Thread milton ruser
, Jul 10, 2009 at 4:04 AM, rajesh j akshay.raj...@gmail.com wrote: Hi, I actually dont need something thats about hmm's.Something like the plotmat in diagram package would suffice.But I need to be able to include plots beside the graph nodes. On Fri, Jul 10, 2009 at 12:32 PM, milton ruser

Re: [R] sort data.frame by specific date column

2009-07-10 Thread milton ruser
Hi Mark, Like this? mydf-read.table(stdin(), head=T) Trade PosType EnDate EnTime ExDate ExTime PL_Pos 1 1 1040107915 1040107 1300164 2 1 1040108909 1040108 1300184 3 1 1040115921 1040115 1300 64 4 1 1040120 1134 1040120 1300124 5

Re: [R] maximum/minimum value of a function

2009-07-10 Thread milton ruser
Hi Derek, like this? myfunction-function(xmin, xmax, number_of_values) { xvalues-seq(from=xmin, to=xmax, by = ((xmax - xmin)/(number_of_values - 1))) outputvalues-xvalues^2+(xmax/xvalues)^3*runif(number_of_values) return(outputvalues) } my.values-myfunction(1,10,100) min(my.values)

Re: [R] mann-whitney U-test - U-value

2009-07-11 Thread milton ruser
Martin, If I am not wrong the difference is that if you has two levels you get U-test, and W-test for more than two levels. But the test is almost the same. cheers milton On Sat, Jul 11, 2009 at 2:13 PM, Martin Batholdy batho...@googlemail.comwrote: Hi, I know that I can perform a

  1   2   3   >