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

[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

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

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] reading excel into R

2010-04-01 Thread milton ruser
Dear Cheba Please, install the package xlsReadWrite I suppose that read.csv is to read csv files, not xls ones. cheers milton On Thu, Apr 1, 2010 at 6:19 AM, cheba meier cheba.me...@googlemail.comwrote: Dear all, I am new R user and I am sure that this question has been asked quite often

[R] barplot with error bar in lattice

2010-04-01 Thread milton ruser
Dear all, I have a data.frame like below and I need to plot horizonal with error bar only for upper limit. On the code below I am able to plot the bars within groups, but I need (1) change from vertical to horizonal plot and (2) add the error bar. any hint are welcome. milton

[R] qplot(ggplo2): axis label sizes/colors

2010-03-23 Thread milton ruser
Dear all, How can I change the size and the color of axis? I would like that the xlab to be larger, xvalues to be larger, and xvalues in black instead of grey. x=runif(10) y=runif(10) require(ggplo2) qplot(x,y) cheers milton [[alternative HTML version deleted]]

Re: [R] labels gone

2010-03-15 Thread milton ruser
Hi Zhongyi, I must confess I not understood completely what you need, but... Tau-seq(0.05,0.95,0.05); Pi - seq(0.19,0.01,-0.01); par(cex.axis=0.8,ps=9,mar=c(1.5,1,0.5,1), oma=c(1,1,0.2,1) ,tck=-0.01); plot(Tau,Pi, type='l', xlab=Tau,ylab=Pi,col=4, xaxt=n, yaxt=n); axis(1,labels=F)

Re: [R] barplot with a probability density curve

2010-03-14 Thread milton ruser
Hi Olga, may be you can work around with the suggestions below, while other provide the best solution for your: barplot(c(1,2,3,5,2,3,1),names.arg=c(100,200,300,400,500,600,700), xlab=diameter,ylab=flow,main=some title,space=0.1) lines(lowess(c(1,2,3,5,2,3,1),f=1), type=b, lwd=2, lty=3,

Re: [R] modify the levels of a variable

2010-03-10 Thread milton ruser
Try data2$X-factor(data2$X) By the way, avoid call data as data, matrix as matrix, vector as vector. cheers milton On Wed, Mar 10, 2010 at 4:35 PM, clee cheel...@gmail.com wrote: If anyone can help with this, I would greatly appreciate it. When I subset a dataframe, and then create a

Re: [R] how to convert character variables into numeric variables directly

2010-03-08 Thread milton ruser
Try: as.numeric(as.character( typec)) milton On Mon, Mar 8, 2010 at 12:55 PM, Xumin Zeng xumin.z...@abbott.com wrote: Here is the example. age=18:29 height=c(76.1,77,78.1,78.2,78.8,79.7,79.9,81.1,81.2,81.8,82.8,83.5) type=c(A, B, C, D,A, B, C, D,A, B, C, D)

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

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

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] 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] 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] 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] 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] 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] 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] 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] (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] 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] 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] 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,

[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

[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

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),]

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

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

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

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

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

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] 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] 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] Count number of zeros in a collumn

2009-09-27 Thread milton ruser
Hi Marcio, How about if (min(dart[,1977])==0)... cheers milton On Sun, Sep 27, 2009 at 4:54 PM, Marcio Resende mresende...@yahoo.com.brwrote: I have a matrix 700x2000 which is sampled in each cycle from another matrix 788x2000 with the numbers 0,1 and 9 There is one specific collumn of

Re: [R] Teach me how to transpose in R

2009-09-26 Thread milton ruser
*Hum* bbb=t(as.matrix(data2)) ? good luck milton On Sun, Sep 27, 2009 at 12:39 AM, Hyo Lee totem...@gmail.com wrote: Hi guys, I need your help!! My goal is to make a csv file from ncdf file. This is the code i've used : hyo=open.ncdf(C:/CRUTEM3.nc) hyo [1] file C:/CRUTEM3.nc has

Re: [R] Create directory and copy files in R

2009-09-22 Thread milton ruser
Tammy Try avoid directories with spaces... may be this solve your problem. Case not, write us again. bests milton On Mon, Sep 21, 2009 at 10:05 AM, Tammy Ma metal_lical...@live.com wrote: HI, All R users, My problem is: fn [1] C:/Documents and

Re: [R] Superimposing (aligning) two graphs

2009-09-22 Thread milton ruser
Hi Maggie: 1) plot() combined with lines() may do the job. 2) google R graph gallery bests milton On Tue, Sep 22, 2009 at 5:23 PM, Maggie la.f...@gmail.com wrote: Hello, I need to superimpose two graphs (plotted representations of waver files) so I can see how much they align. Is there a

Re: [R] Correlate two time series

2009-09-22 Thread milton ruser
Have you looked at ts() ? On Tue, Sep 22, 2009 at 5:49 PM, Maggie la.f...@gmail.com wrote: I was wondering how to correlate two time series in R? I have to plotted waver files I need to correlate to one another to see how well they align.. Any guidance would be very much appreciated!

Re: [R] plotting least-squares regression against x-axis

2009-09-19 Thread milton ruser
By the way, how about: plot(y,mylm$residuals) So you can have an idea of your error distribution across your response variable domain. :-) milton toronto=brazil On Sat, Sep 19, 2009 at 9:02 PM, Sunil Suchindran sunilsuchind...@gmail.com wrote: x - seq(50) y - 10 + x * 2 + rnorm(50,0,10)

Re: [R] (no subject)

2009-09-13 Thread milton ruser
Hi Karen, As you are just starting (I guess), I suggest you take a shortcut: 1. copy and paste the data from Excel to the notepad, and save it as .TXT 2. If all is ok, you will have a tab-delimited file. 3. Be sure that your variables have its names on first line 4. read your data using

Re: [R] fitting nonlinear model

2009-09-09 Thread milton ruser
Hi Bill, I am not sure what you want, but... mydf-read.table(stdin(), head=T, sep=,) Np,year 96,2 91,5 89,7 85,10 plot(Np~year, data=mydf) mymodel-lm(Np~year, data=mydf) abline(mymodel, col=red) bests milton On Wed, Sep 9, 2009 at 2:27 PM, Bill Hyman billhym...@yahoo.com wrote: My data look

Re: [R] fitting nonlinear model

2009-09-09 Thread milton ruser
) cheers milton On Wed, Sep 9, 2009 at 3:01 PM, milton ruser milton.ru...@gmail.com wrote: Hi Bill, I am not sure what you want, but... mydf-read.table(stdin(), head=T, sep=,) Np,year 96,2 91,5 89,7 85,10 plot(Np~year, data=mydf) mymodel-lm(Np~year, data=mydf) abline(mymodel, col=red

Re: [R] Joining Characters in R {issue with paste}

2009-09-09 Thread milton ruser
You not tryed ?paste :-) paste(a,b,sep=) bests milton On Wed, Sep 9, 2009 at 5:08 PM, Abhishek Pratap abhishek@gmail.comwrote: Hi Guys I am want to join to strings in R. I am using paste but not getting desirable result. For the sake of clarity, a quick example: a=Bio b=iology

Re: [R] (no subject)

2009-09-06 Thread milton ruser
Hi Roslina, x-runif(1000,min=0, max=600) hist(x) x2-ifelse(x=300,x,300) summary(x2) hist(x2, ,breaks = c(0,50,100,150,200,250,300)) good luck milton On Mon, Sep 7, 2009 at 12:51 AM, Roslina Zakaria zrosl...@yahoo.com wrote: Hi r-users, I would like to know how to put all the data that is

Re: [R] (no subject)

2009-09-06 Thread milton ruser
May be you also want: x-runif(1000,min=0, max=600) x3-x[x=300] summary(x3) hist(x3, ,breaks = c(0,50,100,150,200,250,300)) bests milton On Mon, Sep 7, 2009 at 1:01 AM, milton ruser milton.ru...@gmail.com wrote: Hi Roslina, x-runif(1000,min=0, max=600) hist(x) x2-ifelse(x=300,x,300

Re: [R] to istall packages such as 'stats'

2009-09-03 Thread milton ruser
Hi Ajay, To install new packages try: install.packages(package_name, dependencies=T) But stats is part of base R. :-) bests miltinho On Thu, Sep 3, 2009 at 6:01 AM, Ajay Singh ajayvi...@yahoo.co.in wrote: Dear All, I have to install 'stats' package in my windows machine, could you let me

Re: [R] variable selection in logistic

2009-09-03 Thread milton ruser
Hi Annie, What kind of data (response and explanatory) you have? As a ecological modeller, I always think first what is my (our) mind models. After that I analyze a set of concorrent - *but with ecological meaning* - models. It helps me to test hypothesis as well as to discuss the results. If you

Re: [R] Easy way to get top 2 items from vector

2009-09-03 Thread milton ruser
Hi Noah, Next time try, please, send a short reproducible code/example. May be this can be (not so elegant but) helpfull myDF-data.frame(cbind(a=runif(5),b=runif(5))) myDF N=2 a.order-rev(order(myDF$a))[1:N] b.order-rev(order(myDF$b))[1:N] myDF.max2a-myDF[a.order,] myDF.max2a

  1   2   3   >